Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,534,406 A→G 5.3% Q128R (CAG→CGG)  nhoA → N‑hydroxyarylamine O‑acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,534,4060AG 5.3% 238.7 / 2.8 75Q128R (CAG→CGG) nhoAN‑hydroxyarylamine O‑acetyltransferase
Reads supporting (aligned to +/- strand):  ref base A (36/35);  new base G (2/2);  total (38/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.29e-01

TTACGCGAGCTGGGGTTTAACGTTCGCAGCTTGTTAGGGCGCGTAGTGTTATCAAATCCGCCAGCATTACCGCCGCGCACCCATCGTTTGCTGTTGGTGGAACTGGAAGAGGAAAAATGGATTGCTGATGTCGGTTTCGGTGGGCAGACGCTAACCGCGCCGATTCGTTTAGTTTCCGATCTCGTGCAGACCACGCCACACGGAGAGTATCGGTTGTTGCAGGAGGGTGATGATTGGGTGTTGCAGTTTAATCATCATCAGCATTGGCAGTCGATGTACCGTTTTGATCTCTGC  >  NC_000913/1534261‑1534554
                                                                                                                                                 |                                                                                                                                                    
ttACGCGAGCTGGGGTTTAACGTTCGCAGCTTGTTAGGGCGCGTAGTGTTATCAAATCCGCCAGCATTACCGCCGCGCACCCATCGTTTGCTGTTGGTGGAACTGGAAGAGGAAAAATGGATTGCTGATGTCGGTTTCGGTGGGCAGa                                                                                                                                                    <  2:403026/148‑1 (MQ=255)
           ggggTTTAACGTTCGCAGCTTGTTAGGGCGCGTAGTGTTATCAAATCCGCCAGCATTACCGCCGCGCACCCATCGTTTGCTGTTGGTGGAACTGGAAGAGGAAAAATGGATTGCTGATGTCGGTTTCGGTGGGCAGACGCTAACCgcgc                                                                                                                                        >  2:274422/1‑149 (MQ=255)
              gTTTAACGTTCGCAGCTTGTTAGGGCGCGTAGTGTTATCAAATCCGCCAGCATTACCGCCGCGCACCCATCGTTTGCTGTTGGTGGAACTGGAAGAGGAAAAATGGATTGCTGATGTCGGTTTCGGTGGGCAGACGCTAACCGCGCCGa                                                                                                                                     >  2:235950/1‑149 (MQ=255)
              gTTTAACGTTCGCAGCTTGTTAGGGCGCGTAGTGTTATCAAATCCGCCAGCATTACCGCCGCGCACCCATCGTTTGCTGTTGGGGGAACTGGAAGAGGAAAAATGGATTGCTGAGGTCGGTTTCGGTGGGCAGACGCTAACCGCGCCGa                                                                                                                                     >  8:368029/1‑149 (MQ=255)
               tttAACGTTCGCAGCTTGTTAGGGCGCGTAGTGTTATCAAATCCGCCAGAATTACCGCCGCGCACCCATCGTTTGCTGTTGGTGGAACTGGAAGAGGAAAAATGGATTGCTGATGTCGGTTTCGGTGGGCAGACGCTAACCGCGCCGAt                                                                                                                                    >  7:211141/1‑149 (MQ=255)
                 tAACGTTCGCAGCTTGTTAGGGCGCGTAGTGTTATCAAATCCGCCAGCATTACCGCCGCGCACCCATCGTTTGCTGTTGGTGGAACTGGAAGAGGAAAAATGGATTGCTGATGTCGGTTTCGGTGGGCAGACGCTAACCGCGCCGATTc                                                                                                                                  >  2:114085/1‑149 (MQ=255)
                   aCGTTCGCAGCTTGTTAGGGCGCGTAGTGTTATCAAATCCGCCAGCATTACCGCCGCGCACCCATCGTTTGCTGTTGGTGGAACTGGAAGTGGAAAAATGGATTGCTGATGTCGGTTTCGGTGGGCAGACGCTAACCGCGCCGATTCGt                                                                                                                                >  2:400820/1‑149 (MQ=255)
                   aCGTTCGCAGCTTGTTAGGGCGCGTAGTGTTATCAAATCCGCCAGCATTACCGCCGCGCACCCATCGTTTGCTGTTGGTGGAACTGGAAGAGGAAAAATGGATTGCTGATGTCGGTTTCGGTGGGCAGACGCTAACCGCGCCGATTCGt                                                                                                                                <  1:354158/149‑1 (MQ=255)
                      ttCGCAGCTTGTTAGGGCGCGTAGTGTTATCAAATCCGCCAGCATTACCGCCGCGCACCCATCGTTTGCTGTTGGTGGAACTGGAAGAGGAAAAATGGATTGCTGATGTCGGTTTCGGTGGGCAGACGCTAACCGCGCCGATTCGTTTa                                                                                                                             <  1:274422/149‑1 (MQ=255)
                         gCAGCTTGTTAGGGCGCGTAGTGTTATCAAATCCGCCAGCATTACCGCCGCGCACCCATCGTTTGCTGTTGGTGGAACTGGAAGAGGAAAAATGGATTGCTGCTGTCGGTTTCGGTGGGCAGACGCTAACCGCGCCGATTCGTTTAGt                                                                                                                           <  6:2165/148‑1 (MQ=255)
                         gCAGCTTGTTAGGGCGCGTAGTGTTATCAAATCCGCCAGCATTACCGCCGCGCACCCATCGTTTGCTGTTGGTGGAACTGGAAGAGGAAAAATGGATTGCTGATGTCGGGTTCGGTGGGCAGACGCTAACCGCGCCGATTCGTTTAGt                                                                                                                           <  3:292395/148‑1 (MQ=255)
                         gCAGCTAGTTAGGGCGCGTAGTGTTATCAAATCCGCCAGCAATACCGCCGCGCACCCATCGTTTGCTGTTGGTGGAACTGGAAGAGGAAAAATGGATTGCTGATGTCGTTTTCGTTAGGCAGACGCTAACCGCGCCGATTCGTATAGt                                                                                                                           >  6:219607/1‑148 (MQ=255)
                                 ttAGGGCGCGTAGTGTTATCAAATCCGCCAGCATTACCGCCGCGCACCCATCGTTTGCTGTTGGTGGAACTGGAAGAGGAAAAATGGATTGCTGATGTCGGTTTCGGTGGGCAGACGCTAACCGCGCCGATTCGTTTAGTTTCCGAtc                                                                                                                   <  7:128359/148‑1 (MQ=255)
                                   aGGGCGCGTAGTGTTATCAAATCCGCCAGCATTACCGCCGCGCACCCATCGTTTGCTGTTGGTGGAACTGGAAGAGGAAAAATGGATTGCTGATGTCGGTTTCGGTGGGCAGACGCTAACCGCGCCGATTCGTTTAGTTTCCGATCTCg                                                                                                                >  7:75289/1‑149 (MQ=255)
                                      gcgcGTAGTGTTATCAAATCCGCCAGCATTACCGCCGCGCACCCATCGTTTGCTGTTGGTGGAACTGGAAGTGGAAAAATGGATTGCTGATGTCGGTTTCGGTGGGCAGACGCTAACCGCGCCGATTCGTTTAGTTTCCGATCTCGTGc                                                                                                             <  1:400820/149‑1 (MQ=255)
                                      gcgcGTAGTGTTATCAAATCCGCCAGCATTACCGCCGCGCACCCATCGTTTGCTGTTGGTGGAACTGGAAGAGGAAAAATGGATTGCTGATGTCGGTTTCGGTGGGCAGACGCTAACCGCGCCGATTCg                                                                                                                                 <  4:230684/129‑1 (MQ=255)
                                      gcgcGTAGTGTTATCAAATCCGCCAGCATTACCGCCGCGCACCCATCGTTTGCTGTTGGTGGAACTGGAAGAGGAAAAATGGATTGCTGATGTCGGTTTCGGTGGGCAGACGCTAACCGCGCCGATTCg                                                                                                                                 >  3:230684/1‑129 (MQ=255)
                                       cgcgTAGTGTTATCAAATCCGCCAGCATTACCGCCGCGCACCCATCGTTTGCTGTTGGTGGAACTGGAAGAGGAAAAATGGATTGCTGATGTCGGTTTCGGTGGGCAGACGCTAACCGCGCCGATTCGTTTAGTTTCCGATCTCGTGCa                                                                                                            <  1:4394/149‑1 (MQ=255)
                                       cgcgTAGTGTTATCAAATCCGCCAGCATTACCGCCGCGCACCCATCGTTTGCTGTTGGTGGAACTGGAAGAGGAAAAATGGATTGCTGATGTCGGTTTCGGTGGGCAGACGCTAACCGCGCCGATTCGTTTAGTTTCCGATCTCGTGCa                                                                                                            <  1:114085/149‑1 (MQ=255)
                                            aGTGTTATCAAATCCGCCAGCATTACCGCCGCGCACCCATCGTTTGCTGTTGGTGGAACTGGAAGAGGAAAAATGGATTGCTGATGTCGGTTTCGGTGGGCAGACg                                                                                                                                                  >  6:112164/1‑106 (MQ=255)
                                            aGTGTTATCAAATCCGCCAGCATTACCGCCGCGCACCCATCGTTTGCTGTTGGTGGAACTGGAAGAGGAAAAATGGATTGCTGATGTCGGTTTCGGTGGGCAGACg                                                                                                                                                  <  5:112164/106‑1 (MQ=255)
                                              tgtTATCAAATCCGCCAGCATTACAGCCGCGCAACCATCGTTTGCTGTTGGTGGATCTGGAAGAGGAAAAATGGATTGCTGATGTCGGTTTCGGTGGGCAGACGCTAACCGCGCCGATTCGTTTAGTTTCCGATCTCGTGCAGACCAc                                                                                                      <  5:219607/148‑1 (MQ=255)
                                                ttATCAAATCCGCCAGCATTACCGCCGCGCACCCATCGTTTGCTGTTGGTGGAACTGGAAGAGGAAAAATGGATTGCTGATGTCGGTTTCGGTGGGCGGACGCTAACCGCGCCGATTCGTTTAGTTTCCGATCTCGTGCAGCCCACGcc                                                                                                   >  5:6070/1‑149 (MQ=255)
                                                  aTCAAATCCGCCAGCATTACCGCCGCGCACCCATCGTTTGCTGTTGGTGGAACTGGAAGAGGAAAAATGGATTGCTGATGTCGGTTTCGGTGGGCAGACGc                                                                                                                                                 <  1:204869/101‑1 (MQ=255)
                                                  aTCAAATCCGCCAGCATTACCGCCGCGCACCCATCGTTTGCTGTTGGTGGAACTGGAAGAGGAAAAATGGATTGCTGATGTCGGTTTCGGTGGGCAGACGc                                                                                                                                                 >  2:204869/1‑101 (MQ=255)
                                                  aTCAAATCCGCCAGCATTACCGCCGCGCACCCATCGTTTGCTGTTGGTGGAACTGGAAGAGGAAAAATGGATTGCTGATGTCGGTTTCGGTGGGCAGACGCTAACCGCGCCGATTCGTTTAGTTTCCGATCTCGTGCAGACCACGCca                                                                                                  >  4:187540/1‑148 (MQ=255)
                                                     aaaTCCGCCAGCATTACCGCCGCGCACCCATCGTTTGCTGTTGGTGGAACTGGAAGAGGAAAAATGGATTGCTGATGTCGGTTTCGGTGGGCAGACGCTAACCGCGCCGATTCGTTTa                                                                                                                             >  7:10225/1‑118 (MQ=255)
                                                     aaaTCCGCCAGCATTACCGCCGCGCACCCATCGTTTGCTGTTGGTGGAACTGGAAGAGGAAAAATGGATTGCTGATGTCGGTTTCGGTGGGCAGACGCTAACCGCGCCGATTCGTTTAGTTTCCGATCCCGTGCAGACCACGCCACACg                                                                                              >  1:102856/1‑149 (MQ=255)
                                                      aaTCCGCCAGCATTACCGCCGCGCACCCATCGTTTGCTGTTGGTGGAACTGGAAGAGGAAAAATGGATTGCTGATGTCGGTTTCGGTGGGCAGACGCTAACCGCGCCGATTCGTTTAGTTTCCGAtctc                                                                                                                 >  5:56208/1‑129 (MQ=255)
                                                      aaTCCGCCAGCATTACCGCCGCGCACCCATCGTTTGCTGTTGGTGGAACTGGAAGAGGAAAAATGGATTGCTGATGTCGGTTTCGGTGGGCAGACGCTAACCGCGCCGATTCGTTTAGTTTCCGAtctc                                                                                                                 <  6:56208/129‑1 (MQ=255)
                                                      aaTCCGCCAGCATTACCGCCGCGCACCCATCGTTTGCTGTTGGTGGAACTGGAAGAGGAAAAATGGATTGCTGATGTCGGTTTCGGTGGGCAGACGCTAACCGCGCCGATTCGTTTAGTTTCCGATCTCGTGCAGACCACg                                                                                                     >  5:328025/1‑141 (MQ=255)
                                                      aaTCCGCCAGCATTACCGCCGCGCACCCATCGTTTGCTGTTGGTGGAACTGGAAGAGGAAAAATGGATTGCTGATGTCGGTTTCGGTGGGCAGACGCTAACCGCGCCGATTCGTTTAGTTTCCGATCTCGTGCAGACCACg                                                                                                     <  6:328025/141‑1 (MQ=255)
                                                      aaTCCGCCAGCATTACCGCCGCGCACCCATCGTTTGCTGTTGGTGGAACTGGAAGAGGAAAAATGGATTGCTGATGTCGGTTTCGGTGGGCAGACGCTAACCGCGCCGATTCGTTTAGTTTCCGATCTCGTGCAGACCACGCCACACg                                                                                              >  4:179942/1‑148 (MQ=255)
                                                        tCCGCCAGCATTACCGCCGCGCACCCATCGTTTGCTGTTGGTGGAACTGGAAGAGGAAAAATGGATTGCTGATGTCGGTTTCGGTGGGCAGACGCTAACCGCGCCGATTCGTTTAGTTTCCGATCTCGTGCAGACCACGCCacac                                                                                               >  5:408039/1‑145 (MQ=255)
                                                        tCCGCCAGCATTACCGCCGCGCACCCATCGTTTGCTGTTGGTGGAACTGGAAGAGGAAAAATGGATTGCTGATGTCGGTTTCGGTGGGCAGACGCTAACCGCGCCGATTCGTTTAGTTTCCGATCTCGTGCAGACCACGCCacac                                                                                               <  6:408039/145‑1 (MQ=255)
                                                           gCCAGCATTACCGCCGCGCACCCATCGTTTGCTGTTGGTGGAACTGGAAGAGGAAAAATGGATTGCTGATGTCGGTTTCGGTGGGCAGACGCTAACCGCGCCGATTCGTTTAGTTTCCGAtc                                                                                                                   >  4:280385/1‑122 (MQ=255)
                                                           gCCAGCATTACCGCCGCGCACCCATCGTTTGCTGTTGGTGGAACTGGAAGAGGAAAAATGGATTGCTGATGTCGGTTTCGGTGGGCAGACGCTAACCGCGCCGATTCGTTTAGTTTCCGAtc                                                                                                                   <  3:280385/122‑1 (MQ=255)
                                                               gCATTACCGCCGCGCACCCATCGTTTGCTGTTGGTGGAACTGGAAGAGGAAAAATGGATTGCTGATGTCGGTTTCGGTGGGCAGACGCTAACCGCGCCGATTCGTTTAGTTTCCGATCTCGTGCAGACCACGCCACACGGAGAGTATCg                                                                                    >  4:94411/1‑149 (MQ=255)
                                                                cATTACCGCCGCGCACCCATCGTTTGCTGTTGGTGGAACTGGAAGAGGAAAAATGGATTGCTGATGTCGGTTTCGGTGGGCAGACGCTAACCGCGCCGATTCGTTTAGTTTCCGATCTCGTGCAGACCACGCCACACGGAGAGTATCgg                                                                                   >  7:104201/1‑149 (MQ=255)
                                                                  ttACCGCCGCGCACCCATCGTTTGCTGTTGGTGGAACTGGAAGAGGAAAAATGGATTGCTGATGTCGGTTTCGGTGGGCAGACGCTAACCGCGCAGATTCGTTTAGTTTCCGATCTCGTGCAGACCACGACACACGGAGAGTATCGgtt                                                                                 <  8:75289/149‑1 (MQ=255)
                                                                       gccgcGCACCCATCGTTTGCTGTTGGTGGAACTGGAAGAGGAAAAATGGATTGCTGATGTCGGTTTCGGTGGGCAGACGCTAACCGCGCCGAtt                                                                                                                                   <  7:240570/94‑1 (MQ=255)
                                                                       gccgcGCACCCATCGTTTGCTGTTGGTGGAACTGGAAGAGGAAAAATGGATTGCTGATGTCGGTTTCGGTGGGCAGACGCTAACCGCGCCGAtt                                                                                                                                   >  8:240570/1‑94 (MQ=255)
                                                                              aCCCATCGTTTGCTGTTGGTGGAACTGGAAGAGGAAAAATGGATTGCTGATGTCGGTTTCGGTGGGCGGACGCTAACCGCGCCGATTCGTTTAGTTTCCGa                                                                                                                     >  6:221609/1‑101 (MQ=255)
                                                                              aCCCATCGTTTGCTGTTGGTGGAACTGGAAGAGGAAAAATGGATTGCTGATGTCGGTTTCGGTGGGCGGACGCTAACCGCGCCGATTCGTTTAGTTTCCGa                                                                                                                     <  5:221609/101‑1 (MQ=255)
                                                                                ccATCGTTTGCTGTTGGTGGAACTGGAAGAGGAAAAATGGATTGCTGATGTCGGTTTCGGTGGGCAGACGCTAACCGCGCCGATTCGTTTAGTTTCCGATCTCGTGCAGACCACGCCACACGGAGAGTATCGGTTGTTGCAGGAGGGt                                                                    >  2:362251/1‑148 (MQ=255)
                                                                                ccATCGTTTGCTGTTGGTGGAACTGGAAGAGGAAAAATGGATTGCTGATGTAGGTTTCGGTGGGCAGACGCTAACCGCGCCGATTCGTTTAGTTTCCGATCTCGTGCAGAAAACGCCACACGGAGAGTATAGGTTGTTGCAGGAGGGt                                                                    >  5:205105/1‑148 (MQ=255)
                                                                                   tCGTTTGCTGTTGGTGGAACTGGAAGAGGAAAAATGGATTGCTGATGTCGGTTTCGGTGGGCAGACGCTAACCGCGCCGATTCGTTTAGTTTCCGATCTCGTGCAGACCACGCCACACGGAGAGTATCGGTTGTTGCAGGAGGGtgat                                                                 >  3:312708/1‑148 (MQ=255)
                                                                                      tttGCTGTTGGTGGAACTGGAAGAGGAAAAATGGATTGCTGATGTCGGTTTCGGTGGGCAGACGCTAACCGCGCAGATTAGTTCAGTTTCCGATCTCGTGCAGACCACGCCACACGGAGAGTATCGGTTGTTGCAGGAGGGTGAGGAtt                                                             >  6:352059/1‑149 (MQ=255)
                                                                                          cTGTTGGTGGAACTGGATGAGGATAAATTGATTGCTGATTTCGGTTTCGGTGGGAAGACGATAACCGCGCCGATTCGTTTAGTTTCCGATCTCGTGCAGACCACTCCAAACTGAGAGTATCTGTTGTTGCTGTAGGGTGATGATTGGgt                                                         <  8:211141/149‑1 (MQ=255)
                                                                                          cTGTTGGTGGAACTGGAAGAGGAAAAATGGATTGCTGATGTCGGTTTCGGTGGGCAGACGCTAACCGCGCCGATTCGTTTAGTTTCCGATCTCGTGCAGACCACGCCACACGGAGAGTATCGGTTGTTGCAGGAGGGTGATGATTGGgt                                                         <  1:235950/149‑1 (MQ=255)
                                                                                           tGTTGGTGGAACTGGAAGAGGAAAAATGGATTGCTGATGTCGGTTTCGGTGGGCAGACGCTAACCGCGCCGATTCGTTTAGTTTCCGATCTCGTGCAGACCACGCCCCACGGAGAGTATCGGTTGTTGCAGGAGGGTGATGATTGGgtg                                                        >  8:355557/1‑149 (MQ=255)
                                                                                            gTTGGTGGAACTGGAAGAGGAAAAATGGATTGCTGATGTCGGTTTCGGTGGGCAGACGCTAACCGCGCCGATTCGTTTAGTTTCCGATCTCGTGCAGACCACGCCACACGGAGAGTATCGGTTGTTGCAGGAGGGTGATGATTGGgtgt                                                       >  3:186352/1‑149 (MQ=255)
                                                                                              tggtggAACTGGAAGAGGAAAAATGGATTGCTGATGTCGGTTTCGGTGGGCGGACGCTAACCGCGCCGATTCGTTTAGTTTCCGATCTCGTGCAGACCACGCCACACTGAGAGTATCGGTTGTTGCAGGAGGGTGATGATTGGGTGtt                                                      <  6:6070/148‑1 (MQ=255)
                                                                                               ggtggAACTGGAAGAGGAAAAATGGATTGCTGATGTCGGTTTCGGTGGGCAGACGCTAACCGCGCCGATTCGTTTAGTTTCCGATCTCGTGCAGACCACGCCACACGGAGAGTATCGGTTGTTGCAGGAGGGTGATGATTGGGTGTTg                                                     <  6:250497/148‑1 (MQ=255)
                                                                                                gtggAACTGGAAGAGGAAAAATGGATTGCTGATGTCGGTTTCGGTGGGCAGACGCTAACCGCGCCGATTCGTTTAGTTTCCGATCTCGTGCAGACCACGCCACACGGAGAGTATCGGTTGTTGc                                                                            >  5:225220/1‑124 (MQ=255)
                                                                                                gtggAACTGGAAGAGGAAAAATGGATTGCTGATGTCGGTTTCGGTGGGCAGACGCTAACCGCGCCGATTCGTTTAGTTTCCGATCTCGTGCAGACCACGCCACACGGAGAGTATCGGTTGTTGc                                                                            <  6:225220/124‑1 (MQ=255)
                                                                                                gtggAACTGGAAGAGGAAAAATGGATTGCTGATGTCGGTTTCGGTGGGCAGACGCTAACCGCGCCGATTCGTTTAGTTTCCGATCTCGTGCAGACCACGCCACACGGAGAGTATCGGTTGTTGCAGGAGGGtga                                                                  <  3:124322/134‑1 (MQ=255)
                                                                                                gtggAACTGGAAGAGGAAAAATGGATTGCTGATGTCGGTTTCGGTGGGCAGACGCTAACCGCGCCGATTCGTTTAGTTTCCGATCTCGTGCAGACCACGCCACACGGAGAGTATCGGTTGTTGCAGGAGGGtga                                                                  >  4:124322/1‑134 (MQ=255)
                                                                                                   gAACTGGACGAGGAAAAATGGATTGCTGATGTCGGTTTCGGTGGGCAGACGCTAACCGCGCCGATTCGTTTAGTTTCCGATCTCGTGCAGACCACGCCACACGGAGAGTATCGGTTGTTGCAGGAGGGTGATGATTGGGTGTTGCAGtt                                                <  1:270930/149‑1 (MQ=255)
                                                                                                         gAAGAGGAAAAATGGATTGCTGATGTCGGTTTCGGTGGGCAGACGCTAACCGCGCCGATTCGTTTAGTTTCCGATCTCGTGCAGACCACGCCACACGGAGAGTATCGGTTGTTGCAGGAGGGTGATGATTGGGTGTTGCAGTTTAatca                                          >  6:104099/1‑149 (MQ=255)
                                                                                                         gAAGAGGAAAAATGGATTGCTGATGTCGGTTTCGGTGGGCAGACGCTAACCGCGCCGATTCGTTTAGTTTCCGATCTCGTGCAGACCACGCCACACGGAGAGTATCGGTTGTTGCAGGAGGGTGATGATTGGGTGTTGCAGTTTAat                                            <  4:312708/147‑1 (MQ=255)
                                                                                                          aaGAGGAAAAATGGATTGCTGATGTCGGTTTCGGTGGGCAGACGCTAACCGCGCCGATTCGTTTAGTTTCCGATCTCGTGCAGACCACGCCACACGGAGAGTATCGGTTGTTGCAGGAGGGTGATGATTGGGTGTTGCAGTTTAatcat                                         <  7:355557/149‑1 (MQ=255)
                                                                                                           agagGAAAAATGGATTGCTGATGTCGGTTTCGGTGGGCAGACGCTAACCGCGCCGATTCGTTTAGTTTCCGATCTCGTGCAGACCACGCCACACGGAGAGTATCGGTTGTTGCAGGAGGGTGATGATTGGGTGTTGCAGTTTAatcat                                         >  3:324570/1‑148 (MQ=255)
                                                                                                                aaaaaTGGATTGCTGCTGTCGGGTGCGGTGGGCAGACGCTAACCGCGCCGATGCGTTTAGTTTCCGATCTCGTGCAGACCACGCCACACGGAGAGTATCGGTTGTTGCAGGAGGGTGATGATTGGGTGTTGCAGTTTAATCATCATCAg                                   <  4:324570/149‑1 (MQ=255)
                                                                                                                aaaaaTGGATTGCTGATGTCGGTTTCGGTGGGCAGACGCTAACCGCGCCGATTCGTTTAGTTTCCGATCTCGTGCAGACCACGCCACACGGAGAGTATCGGGTGTTGCAGGAGGGTGATGATTGGGTGTTGCAGTTTAATCATCATCAg                                   <  3:130968/149‑1 (MQ=255)
                                                                                                                      ggATTGCTGATGTCGGTTTCGGTGGGCAGACGCTAACCGCGCCGATTCGTTTAGTTTCCGATCTCGTGCAGACCACGCCACACGGAGAGTATCGGGTGTTGCAGGAGGGTGATGATTGGGTGTTGCAGTTTAATCATCATCAGCAGTgg                             >  6:254079/1‑149 (MQ=255)
                                                                                                                       gATTGCTGATGTCGGTTTCGGTGGGCAGACGCTAACCGCGCCGATTCGTTTAGTTTCCGATCTCGTGCAGACCACGCCACACGGCGAGTATCGGTTGTTGCAGGAGGGTGATGATTGGGTGTTGCAGTTTAATCATCATCAGCATTgg                             <  3:179942/148‑1 (MQ=255)
                                                                                                                       gATTGCTGATGTCGGTTTCGGTGGGCAGACGCTAACCGCGCCGATTCGTTTAGTTTCCGATCTCGTGCAGACCACGCCACACGGAGAGTATCGGTTGTTGCAGGAGGGTGATGATTGGGTGTTGCAGTTTAATCATCATCAGCATTgg                             <  5:256136/148‑1 (MQ=255)
                                                                                                                           gCTGATGTCGGTTTCGGTGGGCAGACGCTAACCGCGCCGATTCGTTTAGTTTCCGATCTCGTGCAGACCACGCCACACGGAGAGTATCGGTTGTTGCAGGAGGGTGATGATTGGGTGTTGCAGTTTAATCATCATCAGCATTGGCAGTc                        >  8:417486/1‑149 (MQ=255)
                                                                                                                               aTGTCGGTTTCGGTGGGCAGACGCTAACCGCGCCGATTCGTTTAGTTTCCGATCTCGTGCAGACCACGCCACACGGAGAGTATCGGTTGTTGCAGGAGGGTGATGATTGGGTGTTGCAGTTTAATCATCATCAGCATTGGCAGTCGATg                    <  8:104201/149‑1 (MQ=255)
                                                                                                                               aTGTCGGTTTCGGTGGGCAGACGCTAACCGCGCCGATTCGTTTAGTTTCCGATCCCGTGCAGACCACGCCACACGGAGAGTATCGGTTGTTGCAGGAGGGTGATGATTGGGTGTTGCAGTTTAATCATCATCAGCATTGGCAGTCGa                      <  2:102856/147‑1 (MQ=255)
                                                                                                                                 gTCGGTTTCGGTGGGCAGACGCTAACCGCGCCGATTCGTTTAGTTTCCGATCTCGTGGAGACCACGCCACACGGTGAGTATCGGTTGTTGCAGGAGGGTGATGATTGGGGGTTGCAGTTTAATCATCATCAGCATTGGCAGTCGATGt                   <  3:338000/148‑1 (MQ=255)
                                                                                                                                  tCGGTTTCGGTGGGCAGACGCTAACCGCGCCGATTCGTTTAGTTTCCGATCTCGTGCAGACCACGCCACACGGAGAGTATCGGTTGTTGCAGGAGGGTGATGATTGGGTGTTGAAGTTTAATAATCATCAGCAGTGGCAGTCGATGTAc                 >  1:334170/1‑149 (MQ=255)
                                                                                                                                           gTGGGCAGACGCTAACCGCGCCGATTCGTTTAGTTTCCGATCTCGTGCAGACCACGCCACACGGAGAGTATCGGTTGTTGCAGGAGGGTGATGATTGGGTGTTGCAGTTTAATCATCATCAGCATTTGCAGTCGATGTACCGTTTTGAt        <  5:104099/149‑1 (MQ=255)
                                                                                                                                               gCAGACGCTAACCGCGCCGATTCGTTTAGTTTCCGATCTCGTGCAGACCACGCCACACGGAGAGTATCGGTTGTTGCAGGAGGGTGATGATTGGGTGTTGCAGTTTAATCATCATCAGCATTGGCAGTCGATGTACCGTTTTGATctct    >  4:167584/1‑149 (MQ=255)
                                                                                                                                                 aGACGCTAACCTCGCCGATTCGTTTTGTTTCCGATCTCGTGCAGACCACGCCACACGGAGTGTATCGGTTGTTGCAGGAGGGTGATGATTGGGTGTTGCAGTTTAATCATCATAAGCATTGGCAGTCGATGTACCGTTTTGATCTCTga  <  6:205105/149‑2 (MQ=255)
                                                                                                                                                 aGACGCTAACCGCGCCGATTCGTTTAGTTTCCGATCTCGTGCAGACCACGCCACACGGAGAGTATCGGTTGTTGCAGGAGGGTGATGATTGGGTGTTGCAGTTTAATCATCATCAGCa                                 >  6:246925/1‑118 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
TTACGCGAGCTGGGGTTTAACGTTCGCAGCTTGTTAGGGCGCGTAGTGTTATCAAATCCGCCAGCATTACCGCCGCGCACCCATCGTTTGCTGTTGGTGGAACTGGAAGAGGAAAAATGGATTGCTGATGTCGGTTTCGGTGGGCAGACGCTAACCGCGCCGATTCGTTTAGTTTCCGATCTCGTGCAGACCACGCCACACGGAGAGTATCGGTTGTTGCAGGAGGGTGATGATTGGGTGTTGCAGTTTAATCATCATCAGCATTGGCAGTCGATGTACCGTTTTGATCTCTGC  >  NC_000913/1534261‑1534554

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑