Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,090,122 A→G 5.3% intergenic (+76/‑70) 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,1220AG 5.3% 398.2 / 8.1 131intergenic (+76/‑70)hisL/hisGhis operon leader peptide/ATP phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base A (53/71);  new base G (2/5);  total (55/76)
Fisher's exact test for biased strand distribution p-value = 6.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 132 total aligned reads displayed.

AATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGACTCACGCGAATTGCTG  >  NC_000913/2089973‑2090266
                                                                                                                                                     |                                                                                                                                                
aaTAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCAACCTGACTATTCTTTAAGGGGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGa                                                                                                                                                  >  1:749875/1‑150 (MQ=255)
   aaaCATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCa                                                                                                                                               <  1:825842/150‑1 (MQ=255)
       aTTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACAATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCtt                                                                                                                                           >  1:554132/1‑150 (MQ=255)
        ttCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGACCGCATGAGAAAGCCCCCGGAAGATCACCTTc                                                                                                                                          >  1:604229/1‑150 (MQ=255)
           acaGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCgg                                                                                                                                       >  2:882807/1‑150 (MQ=255)
            caGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGGTGATCTTCCggg                                                                                                                                      <  1:223917/150‑1 (MQ=255)
            caGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCggg                                                                                                                                      <  2:888800/150‑1 (MQ=255)
            caGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCggg                                                                                                                                      <  2:782374/150‑1 (MQ=255)
             agagACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCgggg                                                                                                                                     <  1:410146/150‑1 (MQ=255)
             agagACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAAAAAGCCCCCGGAAGGTGATCTTCCgggg                                                                                                                                     <  1:1552180/150‑1 (MQ=255)
              gagaCTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCCTCCTTACTAGTCTTTCAGGCGATGTTTGCTGGAAGAACTTCACATCTCCCATTGGTGAATTAACGCAAGAGAAAGCCCACGGAAGATCACCTTCAggggg                                                                                                                                    >  1:142478/1‑150 (MQ=255)
              gagaCTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGGTGATCTTCCggggg                                                                                                                                    <  2:1227377/150‑1 (MQ=255)
               agaCTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGc                                                                                                                                   <  2:1728283/150‑1 (MQ=255)
                  cTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCttt                                                                                                                                <  2:1138342/150‑1 (MQ=255)
                  cTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCttt                                                                                                                                >  2:1391154/1‑150 (MQ=255)
                   ttttATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCtttt                                                                                                                               >  2:958503/1‑150 (MQ=255)
                      tATGACATGAGTTCAATTTTAACACCACTATGATCACCATCATCCTGAATAGTCTTTCAGGCGATGTGTGCTGTAAGACAGTCAGAGCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTa                                                                                                                            <  2:911403/150‑1 (MQ=255)
                         gACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTg                                                                                                                         <  2:1812429/150‑1 (MQ=255)
                         gACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTg                                                                                                                         <  2:503241/150‑1 (MQ=255)
                           cacGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTgcg                                                                                                                       >  2:1111221/1‑150 (MQ=255)
                             cgcgTTCAATTTAAACACCACCATCATCACCATCCTCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTGTTTTATTGcgcg                                                                                                                     >  2:901855/1‑150 (MQ=255)
                             cgcgTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGATGACGTTCAGATCTTCCAGTGGTGCATGAAAGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTGTTTATTGcgcg                                                                                                                     >  2:1701314/1‑150 (MQ=255)
                                 ttCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTg                                                                                                                 <  1:1436015/150‑1 (MQ=255)
                                   cAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGAt                                                                                                               <  1:560311/150‑1 (MQ=255)
                                     aTTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCCGTTGATaa                                                                                                             >  1:1152576/1‑150 (MQ=255)
                                        aaaaCACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACgg                                                                                                          <  1:1160613/149‑1 (MQ=255)
                                         aaaCACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGt                                                                                                         >  2:389104/1‑150 (MQ=255)
                                         aaaCACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGt                                                                                                         <  1:1630614/150‑1 (MQ=255)
                                         aaaCACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGt                                                                                                         <  2:1909700/150‑1 (MQ=255)
                                           acacCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGtta                                                                                                       >  1:823141/1‑149 (MQ=255)
                                            caccacCATCATCACCATCATCCTGACTAGTCTTTCAGGGGATGTGGGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCa                                                                                                      <  2:217310/150‑1 (MQ=255)
                                            caccacCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCa                                                                                                      <  1:1042609/150‑1 (MQ=255)
                                            caccacCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTACCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCa                                                                                                      >  2:1911954/1‑150 (MQ=255)
                                              ccaccaTCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGa                                                                                                    >  2:1253263/1‑150 (MQ=255)
                                               caccatcatcaccatcatCCTGACTAGTCTTACAGGCGATGTGTGCTGGAAGACATACAGATCTTCCAGTGGTCAATGAACGCATGAGAAAGCCCCCGGAAGAACACCTTCCGGGTGCTTTTTTATTGCGCGGTTCATAACGGTTCAGaa                                                                                                   >  1:1659206/1‑149 (MQ=255)
                                                 ccatcatcaccatcatcCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAg                                                                                                 >  1:230967/1‑150 (MQ=255)
                                                  catcatcaccatcatcCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAgg                                                                                                <  1:14229/150‑1 (MQ=255)
                                                  catcatcaccatcatcCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAgg                                                                                                <  1:1893920/150‑1 (MQ=255)
                                                  catcatcaccatcatcCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGAGCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAgg                                                                                                <  1:1678095/150‑1 (MQ=255)
                                                      atcaccatcatcCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTa                                                                                            >  2:1070290/1‑150 (MQ=255)
                                                        caccatcatcCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTaaa                                                                                          >  2:152300/1‑150 (MQ=255)
                                                          ccatcatcCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAaga                                                                                        >  1:986002/1‑150 (MQ=255)
                                                            atcatcCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGTTCTTCCAGTGTTGCATGAACGCATGAGAAAGCCCGCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAgg                                                                                      <  1:1279853/150‑1 (MQ=255)
                                                            atcatcCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAgg                                                                                      <  2:1111919/150‑1 (MQ=255)
                                                             tcatcCTGACTAGTCTTTCAGGCGATGTCTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGa                                                                                     <  2:1278244/150‑1 (MQ=255)
                                                              catcCTGACTAGTCTTTCAGGCGATGTGTGCTGGAACAAATTCAGATCTTCCTGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGTGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGaa                                                                                    <  1:1214876/150‑1 (MQ=255)
                                                                   tGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACa                                                                               <  2:1557844/150‑1 (MQ=255)
                                                                       tAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAAt                                                                           <  1:1635323/150‑1 (MQ=255)
                                                                       tAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAAt                                                                           <  2:338748/150‑1 (MQ=255)
                                                                        aGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATg                                                                          <  2:1829011/150‑1 (MQ=255)
                                                                        aGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATg                                                                          >  1:755316/1‑150 (MQ=255)
                                                                           ctgtCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGACCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATgaca                                                                       <  1:76795/147‑1 (MQ=255)
                                                                           cTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATgaca                                                                       <  1:1095290/150‑1 (MQ=255)
                                                                           cTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATgaca                                                                       >  1:737846/1‑150 (MQ=255)
                                                                            tttCAGGCGATGTGTGCTGGAAGACAATCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCACCCGGAAGATCACCTTCCGGGGGCTTATTTATTGCGCGTATGATAACAGTTCACACAGGGTTAAAGAGGAATTACAAAATgacag                                                                      >  1:248696/1‑150 (MQ=255)
                                                                               cAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGAGAGGTTTAAAGAGGAATAACAAAATgacagaca                                                                   <  1:469560/150‑1 (MQ=255)
                                                                               cAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATgacagaca                                                                   <  1:683016/150‑1 (MQ=255)
                                                                                 ggCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGacaac                                                                 >  2:1792892/1‑150 (MQ=255)
                                                                                        tgtgCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGtt                                                          >  2:1181008/1‑150 (MQ=255)
                                                                                           gCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTaa                                                       >  1:1671650/1‑149 (MQ=255)
                                                                                             tGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGc                                                     <  2:1236995/150‑1 (MQ=255)
                                                                                             tGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGc                                                     >  2:1732737/1‑150 (MQ=255)
                                                                                             tGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTGCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGc                                                     <  2:957111/150‑1 (MQ=255)
                                                                                             tGGAAGACATTCAGATCTTCCAGTGGTGCATCAACGCATGAGAAAGCCCCCGGAAGATCACCTGCCGGGGGCATTTTTATTGCGCGGATGATAACAGTTAATCAACGTTTAAAGAGGAATAACAAAATGCCAGACAAACCTCGTATACga                                                     >  1:1572358/1‑149 (MQ=255)
                                                                                              ggAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCa                                                    >  1:257017/1‑150 (MQ=255)
                                                                                                aaGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATa                                                  <  2:1765079/150‑1 (MQ=255)
                                                                                                    cATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTa                                              >  1:1583926/1‑150 (MQ=255)
                                                                                                     aTTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTTTTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACAATCGTTTACGCATAGCTAt                                             >  2:346530/1‑150 (MQ=255)
                                                                                                      ttCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATTACAAAATGTCAGACAACCCTCGTTTACGCATAGCTATg                                            >  1:259854/1‑150 (MQ=255)
                                                                                                      ttCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATg                                            <  2:503546/150‑1 (MQ=255)
                                                                                                        cAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCa                                          <  1:1760402/150‑1 (MQ=255)
                                                                                                        cAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTATATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAAAACTCGTTTACGCATAGCTATGCa                                          >  2:111996/1‑150 (MQ=255)
                                                                                                         aGATCTTCCAGTGGGGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAg                                         <  2:386127/150‑1 (MQ=255)
                                                                                                          gATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGa                                        >  2:402880/1‑150 (MQ=255)
                                                                                                             cTTCCAGTGGTGCATGAACGCATGAGAAAGCCACCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGATCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAAAACTCGTTTACGCATAGCTATGCAGAAAt                                     >  1:1492832/1‑150 (MQ=255)
                                                                                                              ttCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATc                                    <  2:628948/150‑1 (MQ=255)
                                                                                                               tCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATcc                                   <  1:1040748/150‑1 (MQ=255)
                                                                                                               tCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCAACCGGGGGCTTTTTTATTTCGCGGTTGATAACGGTTCACACAGCTTTAAAGAGGAAAAACAAAATGACAGACAACACACGGTTAAGCATAGCTATTCAGAAAtac                                   >  1:477724/1‑148 (MQ=255)
                                                                                                                     ggtgCATGAACGCATGAGAAAGCCCCCGGAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGtt                            >  2:1290926/1‑150 (MQ=255)
                                                                                                                      gtgCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTAAGAAAAGCTATGCAGAAATCCGGCCGtt                            >  1:1392236/1‑150 (MQ=255)
                                                                                                                        gCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGCCAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTa                          >  2:1233319/1‑150 (MQ=255)
                                                                                                                         cATGAACTCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATATCAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTaa                         <  2:491155/150‑1 (MQ=255)
                                                                                                                            gAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGtg                      <  1:422131/150‑1 (MQ=255)
                                                                                                                              aCGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCTCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGtgat                    <  1:1391154/150‑1 (MQ=255)
                                                                                                                               cGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGtgatg                   <  2:342724/150‑1 (MQ=255)
                                                                                                                               cGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGtgatg                   >  1:253350/1‑150 (MQ=255)
                                                                                                                                 cATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGAc                 >  2:1340894/1‑150 (MQ=255)
                                                                                                                                   tGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGACTc               >  2:820050/1‑150 (MQ=255)
                                                                                                                                      taaaGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGACTCAcg            <  2:612079/149‑1 (MQ=255)
                                                                                                                                      gaAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGACTCAcg            >  2:785295/1‑150 (MQ=255)
                                                                                                                                        aaGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGACTCAcgcg          <  1:958503/150‑1 (MQ=255)
                                                                                                                                           cccccGGAAGGTGATCTTCCGGGGGCTTTTTTATAGCGCGGATGAAAACGGTTCAGACAGGTGTAAAGAGGAATGACAAAATGACAGACAACACTCGTTTACGCATAGATATGCAGAAATCCGGACGTTTAAGTGATGACTCACGCGAAt       >  1:98107/1‑150 (MQ=255)
                                                                                                                                           cccccGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGACTCACGCGAAt       <  1:936678/150‑1 (MQ=255)
                                                                                                                                           cccccGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGACTCACGCGAAt       <  2:1440902/150‑1 (MQ=255)
                                                                                                                                            ccccGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGACTCACGCGAAtt      >  2:1500311/1‑150 (MQ=255)
                                                                                                                                             cccGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGACTCACGCGAATTg     <  1:1472332/150‑1 (MQ=255)
                                                                                                                                             cccGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGACTCACGCGAATTg     <  1:1041683/150‑1 (MQ=255)
                                                                                                                                             cacGGAAGGTGATCTTCCGGGGGCTTTTTTATTGCTCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCCTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGACTCACGCGAAtag     >  2:1115290/3‑148 (MQ=255)
                                                                                                                                                ggAAGATCACGTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGGTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTTATGACTCACGCGATTTGCTg  <  1:148071/150‑1 (MQ=255)
                                                                                                                                                ggAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGACTCACGCGAATTGCTg  <  2:677287/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                
AATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGACTCACGCGAATTGCTG  >  NC_000913/2089973‑2090266

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