Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,090,126 C→T 5.5% intergenic (+80/‑66) hisL → / → hisG his operon leader peptide/ATP phosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,090,1260CT 5.5% 343.9 / 5.9 108intergenic (+80/‑66)hisL/hisGhis operon leader peptide/ATP phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base C (39/63);  new base T (4/2);  total (43/65)
Fisher's exact test for biased strand distribution p-value = 2.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.77e-01

Only 100 of 109 total aligned reads displayed.

CATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGACTCACGCGAATTGCTGGCGCGCTGT  >  NC_000913/2089979‑2090275
                                                                                                                                                   |                                                                                                                                                     
cATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCt                                                                                                                                                     >  2:537138/1‑150 (MQ=255)
cATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCt                                                                                                                                                     >  2:220326/1‑150 (MQ=255)
 aTTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCtt                                                                                                                                                    >  1:864373/1‑150 (MQ=255)
    cacaGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCg                                                                                                                                                 <  2:151424/150‑1 (MQ=255)
    cacaGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCg                                                                                                                                                 >  1:1347375/1‑150 (MQ=255)
       agagACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCgggg                                                                                                                                              <  1:724102/150‑1 (MQ=255)
        gagaCTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCggggg                                                                                                                                             <  1:1279667/150‑1 (MQ=255)
        gagaCTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGACGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGGTGATCTTCCggggg                                                                                                                                             <  1:48793/150‑1 (MQ=255)
         agaCTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGc                                                                                                                                            <  2:979085/150‑1 (MQ=255)
            cTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCttt                                                                                                                                         >  1:773643/1‑150 (MQ=255)
              tttATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCttttt                                                                                                                                       >  1:1157292/1‑150 (MQ=255)
               ttATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCtttttt                                                                                                                                      <  1:363233/150‑1 (MQ=255)
                      acGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTgcgc                                                                                                                               <  1:458761/150‑1 (MQ=255)
                        gcgTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCgg                                                                                                                             >  2:710850/1‑150 (MQ=255)
                        gcgTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCgg                                                                                                                             <  1:300458/150‑1 (MQ=255)
                          gTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGtt                                                                                                                           >  2:67683/1‑150 (MQ=255)
                          gTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGtt                                                                                                                           >  1:118758/1‑150 (MQ=255)
                          gTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCACCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGtt                                                                                                                           <  1:1427170/150‑1 (MQ=255)
                           ttCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTg                                                                                                                          <  2:56457/150‑1 (MQ=255)
                             cAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGAt                                                                                                                        >  1:840288/1‑150 (MQ=255)
                             cAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGAt                                                                                                                        <  2:261127/150‑1 (MQ=255)
                                tttAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAAc                                                                                                                     <  2:1510631/150‑1 (MQ=255)
                                 ttAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACg                                                                                                                    >  1:1193638/1‑150 (MQ=255)
                                 ttAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACg                                                                                                                    <  1:1516451/150‑1 (MQ=255)
                                  tAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACgg                                                                                                                   <  2:1528556/150‑1 (MQ=255)
                                  tAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACgg                                                                                                                   <  1:1387282/150‑1 (MQ=255)
                                  tAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACgg                                                                                                                   <  1:1200887/150‑1 (MQ=255)
                                   aaaCACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCATTTTTATTTCGCGGTTGATAAAGGt                                                                                                                  >  2:1025038/1‑150 (MQ=255)
                                       accaccATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAg                                                                                                              <  1:649608/150‑1 (MQ=255)
                                       accaccATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAg                                                                                                              >  2:1177307/1‑150 (MQ=255)
                                        ccaccaTCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGa                                                                                                             >  2:386961/1‑150 (MQ=255)
                                         caccatcatcaccatcatCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGAc                                                                                                            <  2:257563/150‑1 (MQ=255)
                                            catcatcaccatcatcCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAgg                                                                                                         <  2:111711/150‑1 (MQ=255)
                                              tcatcaccatcatcCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCTGTTGATAACGGTTCAGACAGttt                                                                                                       <  2:1251878/150‑1 (MQ=255)
                                              tcatcaccatcatcCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGtt                                                                                                       <  1:171138/150‑1 (MQ=255)
                                              tcatcaccatcatcCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGtt                                                                                                       <  1:710850/150‑1 (MQ=255)
                                              tcatcaccatcatcCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGtt                                                                                                       <  1:386961/150‑1 (MQ=255)
                                                  caccatcatcCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTaaa                                                                                                   >  2:593015/1‑150 (MQ=255)
                                                    ccatcatcCTGACTAGTCTTTCAGGCTATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAaga                                                                                                 <  1:1212746/150‑1 (MQ=255)
                                                    ccatcatcCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAaga                                                                                                 >  2:1106111/1‑150 (MQ=255)
                                                       tcatcCTGAGTAGTCTTTCTGGCGATGTGTGCTGGAAGACTTTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGATGATCTCCTTCCGGGGGCTTTTTTATTGAGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGa                                                                                              <  2:821357/150‑1 (MQ=255)
                                                       tcatcCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGa                                                                                              <  1:673053/150‑1 (MQ=255)
                                                          tcCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATa                                                                                           <  1:492035/150‑1 (MQ=255)
                                                          tcCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATa                                                                                           <  2:1326235/150‑1 (MQ=255)
                                                          tcCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATa                                                                                           <  2:657199/150‑1 (MQ=255)
                                                           ccTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATaa                                                                                          <  1:48312/150‑1 (MQ=255)
                                                            cTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAAc                                                                                         >  1:445125/1‑150 (MQ=255)
                                                                cTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACaaaa                                                                                     <  2:892647/150‑1 (MQ=255)
                                                                cTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACaaaa                                                                                     <  2:110409/150‑1 (MQ=255)
                                                                cTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACaaaa                                                                                     <  1:537138/150‑1 (MQ=255)
                                                                cTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACaaaa                                                                                     >  1:375846/1‑150 (MQ=255)
                                                                  aGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATg                                                                                   <  1:342899/150‑1 (MQ=255)
                                                                         cAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGGTGATCTTCCggggg                                                                                                                                             <  2:444614/85‑1 (MQ=38)
                                                                         cAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGGTGATCTTCCggggg                                                                                                                                             >  1:444614/1‑85 (MQ=38)
                                                                           ggCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAAt                                                                                            >  1:1423486/1‑132 (MQ=255)
                                                                           ggCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAAt                                                                                            <  2:1423486/132‑1 (MQ=255)
                                                                                 gtgtgCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGt                                                                    <  1:1353222/150‑1 (MQ=255)
                                                                                 gtgtgCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGt                                                                    >  2:1426213/1‑150 (MQ=255)
                                                                                  tgtgCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGtt                                                                   >  1:604344/1‑150 (MQ=255)
                                                                                  tgtgCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGtt                                                                   >  1:1272241/1‑150 (MQ=255)
                                                                                     gCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTAc                                                                >  2:441089/1‑150 (MQ=255)
                                                                                     gCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTAc                                                                >  1:1188990/1‑150 (MQ=255)
                                                                                     gCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTAc                                                                >  1:1505300/1‑150 (MQ=255)
                                                                                      cTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACg                                                               <  1:441089/150‑1 (MQ=255)
                                                                                        ggAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCa                                                             >  2:42747/1‑150 (MQ=255)
                                                                                         gAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCAt                                                            <  2:807366/150‑1 (MQ=255)
                                                                                           aGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAgg                                                                                                         >  2:74853/1‑103 (MQ=255)
                                                                                           aGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAgg                                                                                                         <  1:74853/103‑1 (MQ=255)
                                                                                            gACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGc                                                         <  2:1157292/150‑1 (MQ=255)
                                                                                            gACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGc                                                         >  1:924338/1‑150 (MQ=255)
                                                                                              cATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTa                                                       <  1:220326/150‑1 (MQ=255)
                                                                                                ttCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATg                                                     <  2:1193638/150‑1 (MQ=255)
                                                                                                 tCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGc                                                    <  2:79190/150‑1 (MQ=255)
                                                                                                  cAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCa                                                   >  1:1505327/1‑150 (MQ=255)
                                                                                                   aGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGacaac                                                                          >  2:58045/1‑126 (MQ=255)
                                                                                                   aGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGacaac                                                                          <  1:58045/126‑1 (MQ=255)
                                                                                                   aGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAg                                                  <  1:42747/150‑1 (MQ=255)
                                                                                                      tCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCATTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGaaa                                               >  1:700027/1‑150 (MQ=255)
                                                                                                      tCTTCCAGTGCTGCATGAACGCATGAGAAAGCCCAAGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGATGATAACGGTTAAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGATATGAAGaaa                                               >  2:1143767/1‑150 (MQ=255)
                                                                                                       cTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAAt                                              >  2:302765/1‑150 (MQ=255)
                                                                                                          ccAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCg                                           >  2:177654/1‑150 (MQ=255)
                                                                                                           cAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCgg                                          <  1:1426213/150‑1 (MQ=255)
                                                                                                           aaGTGGTGAATTAACGCATGATAAAGCCCCCGTAATATAACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCgg                                          <  1:1025038/149‑1 (MQ=255)
                                                                                                                   cATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTaa                                  >  2:1273174/1‑150 (MQ=255)
                                                                                                                         cGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGtgatg                            <  1:314508/150‑1 (MQ=255)
                                                                                                                         cGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGtgatg                            <  2:840288/150‑1 (MQ=255)
                                                                                                                          gCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGtgatga                           >  1:1449544/1‑150 (MQ=255)
                                                                                                                           cATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGAc                          >  1:333430/1‑150 (MQ=255)
                                                                                                                            ataaaaaaGCCCCCGGAAGGTGATCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGACt                         >  2:1525084/6‑150 (MQ=255)
                                                                                                                             tGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGACTc                        <  2:1188990/150‑1 (MQ=255)
                                                                                                                              gagaAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGACTCa                       <  1:207974/150‑1 (MQ=255)
                                                                                                                                 aaaGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGACTCAcgc                    <  2:864373/150‑1 (MQ=255)
                                                                                                                                  aaGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGACTCAcgcg                   <  2:1505300/150‑1 (MQ=255)
                                                                                                                                  aaGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGACTCAcgcg                   <  2:439611/150‑1 (MQ=255)
                                                                                                                                      ccccGGAAGGTGATCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGACTCACGCGAAtt               >  1:157371/1‑150 (MQ=255)
                                                                                                                                       cccGGAAGGTGATCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGACTCACGCGAATTGc             >  1:1220880/1‑150 (MQ=255)
                                                                                                                                         cGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGACTCACGCGAATTGCt            <  1:177654/150‑1 (MQ=255)
                                                                                                                                          ggAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGACTCACGCGAATTGCTg           <  2:1272241/150‑1 (MQ=255)
                                                                                                                                           gAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGACTCACGCGAATTGCTgg          <  1:302765/150‑1 (MQ=255)
                                                                                                                                                   ccTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGACTCACGCGAATTGCTGGCGCGCtgt  >  2:520780/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGACTCACGCGAATTGCTGGCGCGCTGT  >  NC_000913/2089979‑2090275

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