Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,286,145 A→T 6.3% intergenic (+9/‑66) yejM → / → proL essential inner membrane DUF3413 domain‑containing protein; lipid A production and membrane permeability factor/tRNA‑Pro

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,286,1450AT 6.3% 295.8 / 6.4 109intergenic (+9/‑66)yejM/proLessential inner membrane DUF3413 domain‑containing protein; lipid A production and membrane permeability factor/tRNA‑Pro
Reads supporting (aligned to +/- strand):  ref base A (49/53);  new base T (5/2);  total (54/55)
Fisher's exact test for biased strand distribution p-value = 2.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.22e-01

Only 100 of 117 total aligned reads displayed.

CGAAAAAGACGCTGGTACTGAACAATAACGGTAAATACCGCACTTACAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCG‑‑GTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCC  >  NC_000913/2285997‑2286294
                                                                                                                                                    |                                                                                                                                                       
cGAAAAAGACGCTGGTACTGAACAATAACGGTAAATACCGCACTTACAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATaa                                                                                                                                                        <  2:1177223/150‑1 (MQ=255)
cGAAAAAGACGCTGGTACTGAACAATAACGGTAAATACCGCACTTACAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATaa                                                                                                                                                        >  1:1709243/1‑150 (MQ=255)
ataaaaaCAAGCTGGTACAGAACAATAACGTTAAATACCGCACTTACAACTTACGTTGTGAAAGAGTGAAAGATGAAAAAACACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGATGCGTTTTATCGCTAACTGATTAATTATaa                                                                                                                                                        <  1:309579/148‑1 (MQ=255)
 gAAAAAGACGCTGGTACTGAACAATAACGGTAAATACCGCACTTACAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATaaa                                                                                                                                                       <  1:1163387/150‑1 (MQ=255)
  aaaaaGACGCTGGTACTGAACAATAACGGTAAATACCGCACTTACAACTTACGTTGTGAAATAGTGAAATATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTTACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAAt                                                                                                                                                      <  2:1242487/150‑1 (MQ=255)
      aGACGCTGGTACTGAACAATAACGGTAAATACCGCACTTACAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGt                                                                                                                                                  <  2:767599/150‑1 (MQ=255)
         cGCTGGTACTGAACAATAACGGTAAATACCGCACTTACAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAg                                                                                                                                               <  1:2272/150‑1 (MQ=255)
          ccttctACAGAACAATAACGGAAAATACCGCAATTACAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTTATAATTAATCAGTTAGc                                                                                                                                              <  1:1584428/145‑1 (MQ=255)
           cTGGTACTGAACAATAACGGTAAATACCGCACTTACAACTTACGTTGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCg                                                                                                                                             <  2:866116/150‑1 (MQ=255)
              gTACTGAACAATAACGGTAAATACCGCACTTACAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGaaa                                                                                                                                          >  2:979368/1‑150 (MQ=255)
                 cTGAACAATAACGGTAAATACCGCACTTACAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAatat                                                                                                                                       >  1:511081/1‑150 (MQ=255)
                   gAACAATAACGGTAAATACCGCACTTACAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCt                                                                                                                                     >  2:1484020/1‑150 (MQ=255)
                     aCAATAACGGTAAATACCGCACTTACAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTa                                                                                                                                   <  2:754427/150‑1 (MQ=255)
                     aCAATAACGGTAAATACCGCACTTACAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTa                                                                                                                                   <  2:379461/150‑1 (MQ=255)
                     aCAATAACGGTAAATACCGCACTTACAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTa                                                                                                                                   >  1:374302/1‑150 (MQ=255)
                      cAATAACGGTAAATACCGCACTTACAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTAc                                                                                                                                  <  1:1077478/150‑1 (MQ=255)
                          aaCGGTAAATACCGCACTTACAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGc                                                                                                                              <  1:534062/150‑1 (MQ=255)
                              gTAAATACCGCACTTACAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATc                                                                                                                          <  1:1653431/150‑1 (MQ=255)
                              gTAAATACCGCACTTACAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATc                                                                                                                          >  2:448075/1‑150 (MQ=255)
                              gTAAATACCGCACTTACAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGATTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATc                                                                                                                          >  1:615747/1‑150 (MQ=255)
                               tAAATACCGCACTTACAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCg                                                                                                                         <  2:1407911/150‑1 (MQ=255)
                               tAAATACCGCACTTACAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCGGGCAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATCTCTTACTTGCAATCg                                                                                                                         >  2:88598/1‑150 (MQ=255)
                                aaaTACCGCACTTACAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCg‑‑g                                                                                                                      <  1:974618/150‑1 (MQ=255)
                                 aaTACCGCACTTACAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCG‑‑gt                                                                                                                     <  2:166089/150‑1 (MQ=255)
                                   tACCGCACTTACAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCG‑‑gtgt                                                                                                                   >  2:806180/1‑150 (MQ=255)
                                    aCCGCACTTACAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCG‑‑Gtgtg                                                                                                                  <  1:1118270/150‑1 (MQ=255)
                                    aCCGCACTTACAACTTACGTGGTGAAAGAGTCAAAGAATAAATACCACAGTTAATTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCG‑‑Gtgtg                                                                                                                  <  1:56515/150‑1 (MQ=255)
                                      ctcACTTAAAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTTACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCG‑‑GTGTGGa                                                                                                                <  1:1538267/148‑1 (MQ=255)
                                      cGCACTTACAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCG‑‑GTGTGGa                                                                                                                >  2:1683304/1‑150 (MQ=255)
                                       gCACTTACAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCG‑‑GTGTGGaa                                                                                                               <  2:1284887/150‑1 (MQ=255)
                                          cTTACAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCG‑‑GTGTGGAAAAc                                                                                                            <  1:1618471/150‑1 (MQ=255)
                                            tACAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCG‑‑GTGTGGAAAACgg                                                                                                          >  1:474308/1‑150 (MQ=255)
                                             aCAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCG‑‑GTGTGGAAAACGgt                                                                                                         <  1:919225/150‑1 (MQ=255)
                                               aaCTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCG‑‑GTGTGGAAAACGgtag                                                                                                       <  1:1660497/150‑1 (MQ=255)
                                                 cTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCG‑‑GTGTGGAAAACGgtagta                                                                                                     >  1:1221590/1‑150 (MQ=255)
                                                 cTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCG‑‑GTGTGGAAAACGgtagta                                                                                                     >  2:1435029/1‑150 (MQ=255)
                                                   tACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCG‑‑GTGTGGAAAACGGTAGTAtt                                                                                                   <  1:104136/150‑1 (MQ=255)
                                                       tggtgAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCG‑‑GTGTGGAAAACGGTAGTATTagca                                                                                               >  2:850606/1‑150 (MQ=255)
                                                                   gaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCG‑‑GTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCgg                                                                                   >  2:761909/1‑150 (MQ=255)
                                                                   gaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCG‑‑GTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCgg                                                                                   >  2:541298/1‑150 (MQ=255)
                                                                    aaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCG‑‑GTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGc                                                                                  <  1:806180/150‑1 (MQ=255)
                                                                      agaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCG‑‑GTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCAc                                                                                >  2:1224911/1‑150 (MQ=255)
                                                                      agaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCG‑‑GTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCAc                                                                                <  1:492415/150‑1 (MQ=255)
                                                                       gaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCG‑‑GTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACg                                                                               <  1:522326/150‑1 (MQ=255)
                                                                       gaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCG‑‑GTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACg                                                                               >  1:596625/1‑150 (MQ=255)
                                                                        aTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCG‑‑GTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGt                                                                              >  1:1665311/1‑150 (MQ=255)
                                                                         tGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCG‑‑GTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTa                                                                             >  1:126990/1‑150 (MQ=255)
                                                                         tGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGAGAAATATCTTACTTGCAATCG‑‑GTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTa                                                                             >  2:1162921/1‑150 (MQ=255)
                                                                          gAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCG‑‑GTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAg                                                                            <  2:1665311/150‑1 (MQ=255)
                                                                            aaaaCCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCG‑‑GTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAgcg                                                                          <  1:1449546/150‑1 (MQ=255)
                                                                             aaaCCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCG‑‑GTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAgcgc                                                                         <  2:143966/150‑1 (MQ=255)
                                                                             aaaCCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCG‑‑GTGTGGAAAACGGTAGTATTAGCAGCCACGACTCGGCACGTAgcgc                                                                         >  2:434657/1‑150 (MQ=255)
                                                                              aaCCACAGTTATTTTTGTTATTGCAAGTGCTGACAGACGTGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCG‑‑GTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCa                                                                        <  2:786524/150‑1 (MQ=255)
                                                                               aCCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCG‑‑GTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAg                                                                       >  2:808296/1‑150 (MQ=255)
                                                                                ccACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCG‑‑GTGTGGAAAACGGTAGTATTAGCAGCCACGAg                                                                                       <  2:671596/133‑1 (MQ=255)
                                                                                ccACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCG‑‑GTGTGGAAAACGGTAGTATTAGCAGCCACGAg                                                                                       >  1:671596/1‑133 (MQ=255)
                                                                                ccACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCG‑‑GTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGc                                                                      <  2:374302/150‑1 (MQ=255)
                                                                                   caGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCG‑‑GTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTg                                                                   >  1:461392/1‑150 (MQ=255)
                                                                                    aGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCG‑‑GTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTgg                                                                  <  2:1263136/150‑1 (MQ=255)
                                                                                     gTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCG‑‑GGGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGt                                                                 <  2:266254/150‑1 (MQ=255)
                                                                                             tGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCG‑‑GTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCAcc                                                         >  2:959754/1‑150 (MQ=255)
                                                                                               ttATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCG‑‑GTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGt                                                       <  1:1484020/150‑1 (MQ=255)
                                                                                                 aTTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGt                                                       <  2:182232/150‑1 (MQ=255)
                                                                                                    gCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCG‑‑GTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATgg                                                  >  2:80765/1‑150 (MQ=255)
                                                                                                    gCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCG‑‑GTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATgg                                                  >  2:1188485/1‑150 (MQ=255)
                                                                                                     cAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCG‑‑GTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATggg                                                 <  1:979368/150‑1 (MQ=255)
                                                                                                     cAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCG‑‑GTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATggg                                                 >  1:1575175/1‑150 (MQ=255)
                                                                                                        gTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCG‑‑GTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGgtg                                              >  2:599614/1‑150 (MQ=255)
                                                                                                          gCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCG‑‑GTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTc                                            >  2:1129983/1‑150 (MQ=255)
                                                                                                              aCAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCG‑‑GTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCgggg                                        <  2:1221590/150‑1 (MQ=255)
                                                                                                               cAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCG‑‑GTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCggggg                                       <  2:588422/150‑1 (MQ=255)
                                                                                                                 gACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCG‑‑GTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTc                                     >  2:1678501/1‑150 (MQ=255)
                                                                                                                   cGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCG‑‑GTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCgg                                   >  2:788474/1‑150 (MQ=255)
                                                                                                                      gaAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCG‑‑GTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAgg                                <  2:474308/150‑1 (MQ=255)
                                                                                                                      gaAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCG‑‑GTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAgg                                <  2:1575175/150‑1 (MQ=255)
                                                                                                                      gaAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCG‑‑GTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAgg                                <  2:1019247/150‑1 (MQ=255)
                                                                                                                         gCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCG‑‑GTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTc                             >  2:317163/1‑150 (MQ=255)
                                                                                                                           gTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCG‑‑GTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCaa                           >  1:1176500/1‑150 (MQ=255)
                                                                                                                              ttATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCG‑‑GTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATc                        <  2:461392/150‑1 (MQ=255)
                                                                                                                                ttCGCTAACTGATTTATAATTAATCAGTTAGCGAAATATCTTACTTGCAATCG‑‑GTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCct                      >  1:665639/2‑150 (MQ=255)
                                                                                                                                aTCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCG‑‑GTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCct                      <  1:1129983/150‑1 (MQ=255)
                                                                                                                                 tCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCG‑‑GTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCctc                     >  1:1119548/1‑150 (MQ=255)
                                                                                                                                  cGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCG‑‑GTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCctct                    <  2:1119548/150‑1 (MQ=255)
                                                                                                                                   gCTAACTGATTTATAATTAATCAGTTAGCGAAATATCTTACTTGCAATCG‑‑GTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCtctc                   >  2:802626/1‑150 (MQ=255)
                                                                                                                                   gCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCG‑‑GTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCtctc                   >  1:593342/1‑150 (MQ=255)
                                                                                                                                    cTAACTGATTTATAATTAATCAGTTAGCGAAATATCTTACTTGCAATCG‑‑GTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCg                  >  2:1202268/1‑150 (MQ=255)
                                                                                                                                    cTAACTGATTTATAATTAATCAGTTAGCGAAATATCTTACTTGCAATCG‑‑GTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCg                  >  2:1418943/1‑150 (MQ=255)
                                                                                                                                     tAACTGATTAATTATTAATCAGTTAGCGAAATATCTTACTTGCAATCG‑‑GTGTGGAAAAAGGTAGTATTAGCAGCCACGCGTCGGCACTTACAGCAGACTGGTAGAGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAAACCTCTCGt                 >  1:1599456/1‑150 (MQ=255)
                                                                                                                                     tAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCG‑‑GTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGt                 >  2:1578578/1‑150 (MQ=255)
                                                                                                                                      aaCTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCG‑‑GTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTg                <  1:99204/150‑1 (MQ=255)
                                                                                                                                      aaCTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCG‑‑GTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTg                >  1:984716/1‑150 (MQ=255)
                                                                                                                                       aCTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCG‑‑GGGGGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGc               <  2:1487114/150‑1 (MQ=255)
                                                                                                                                        cTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCG‑‑GTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGcc              <  2:332184/150‑1 (MQ=255)
                                                                                                                                         tGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCG‑‑GTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCg             <  2:504497/150‑1 (MQ=255)
                                                                                                                                          gATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCG‑‑GTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGAAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCgc            >  1:1358604/1‑149 (MQ=255)
                                                                                                                                           attaattaTAAATCAGTTAGCGAAATATCTTACTTGCAATCG‑‑GTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGAc           <  1:468644/150‑1 (MQ=255)
                                                                                                                                             taattaTAAATCAGTTAGCGAAATATCTTACTTGCAATCG‑‑GTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCa         >  1:770216/1‑150 (MQ=255)
                                                                                                                                              aattaTAAATCAGTTAGCGAAATATCTTTCTTGCAATCG‑‑GTGTGGAAAACGGTAGTATTAGCAGCAACGAGTCGGCACGTAGCGCCGCCAGGAAGCTCACCGACATGGGGTGACCGGGGTCGGAGGTTCATATCCTCTAGTGCCGCCCaa        >  1:1130738/1‑150 (MQ=255)
                                                                                                                                                   taAATCAGTTAGCGAAATATCTTACTTGCAATCG‑‑GTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATc   <  1:317163/150‑1 (MQ=255)
                                                                                                                                                    aaaTCAGTTAGCGAAATATCTTACTTGCAATCG‑‑GTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATcc  >  1:677678/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                       
CGAAAAAGACGCTGGTACTGAACAATAACGGTAAATACCGCACTTACAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCG‑‑GTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCC  >  NC_000913/2285997‑2286294

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