Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,286,139 A→T 5.3% intergenic (+3/‑72) 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,1390AT 5.3% 297.2 / 5.3 113intergenic (+3/‑72)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 (51/56);  new base T (4/2);  total (55/58)
Fisher's exact test for biased strand distribution p-value = 4.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.09e-01

Only 100 of 119 total aligned reads displayed.

CACCCCGAAAAAGACGCTGGTACTGAACAATAACGGTAAATACCGCACTTACAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCG‑‑GTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAA  >  NC_000913/2285992‑2286288
                                                                                                                                                   |                                                                                                                                                       
 accCCGAAAAAGACGCTGGTACTGAACAATAACGGTAAATACCGCACTTACAACTTACGTCGTGAAAGAGGGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGattaatt                                                                                                                                                      <  1:699789/150‑1 (MQ=255)
  ccccGAAAAAGACGCTGGTACTGAACAATAACGGTAAATACCGCACTTACAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGattaatta                                                                                                                                                     <  2:560110/150‑1 (MQ=255)
   cccTAAAAAGACGATTGTACTGAACAATAACGGTAAATACCTCACTTACAATTTACGTGGTGAATTATTTAAAGATGAAAAACCACTGTTAAGTTTGTTATTTCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATtat                                                                                                                                                    <  1:1335260/150‑1 (MQ=255)
   cccGAAAAAGACGCTGGTACTGAACAATAACGGTAAATACCGCACTTACAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTTACTGATTAATtat                                                                                                                                                    >  2:50377/1‑150 (MQ=255)
   cccGAAAAAGACGCTGGTACTGAACAATAACGGTAAATACCGCACTTACAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTCTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATtat                                                                                                                                                    >  1:266254/1‑150 (MQ=255)
     cGAAAAAGACGCTGGTACTGAACAATAACGGTAAATACCGCACTTACAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATaa                                                                                                                                                  <  2:1177223/150‑1 (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)
              ggCTGGTAATGAACAATAACGGTAAATACCGCACTTACAACTTACGTTGTGAAAGAGTGAAAGATGAAAAACCTCAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTTATAATTAATCAGTTAg                                                                                                                                         <  2:1527881/149‑1 (MQ=255)
              cGCTGGTACTGAACAATAACGGTAAATACCGCACTTACAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAg                                                                                                                                         <  1:2272/150‑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:379461/150‑1 (MQ=255)
                          aCAATAACGGTAAATACCGCACTTACAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTa                                                                                                                             <  2:754427/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)
                                 cGGAAAATACCGCACTTACAACTTACGTGGTCAAAGAGTGAAAAAAGAAAAACCACAGTTAAGATTGTTATTGCAAGTGCTGACATACCAGAAGCGATTAAACGCTAACTGATAAATTATAAATCAGTTAGCGAAATATCTTACTTGCaa                                                                                                                      >  1:310717/1‑150 (MQ=255)
                                   gTAAATACCGCACTTACAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATc                                                                                                                    <  1:1653431/150‑1 (MQ=255)
                                   gTAAATACCGCACTTACAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATc                                                                                                                    >  2:448075/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)
                                          ccGCACTTACAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCG‑‑GTGTgg                                                                                                           <  1:1478636/150‑1 (MQ=255)
                                           ctcACTTAAAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTTACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCG‑‑GTGTGGa                                                                                                          <  1:1538267/148‑1 (MQ=255)
                                           cGCACTTACAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCG‑‑GTGTGGa                                                                                                          >  2:1683304/1‑150 (MQ=255)
                                               cTTACAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCG‑‑GTGTGGAAAAc                                                                                                      <  1:1618471/150‑1 (MQ=255)
                                                 tACAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCG‑‑GTGTGGAAAACgg                                                                                                    >  1:504497/1‑150 (MQ=255)
                                                 tACAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCG‑‑GTGTGGAAAACgg                                                                                                    >  1:474308/1‑150 (MQ=255)
                                                    aaCTTACGTGGTGAAATAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGGGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCG‑‑GTGTGGAAAACGgtag                                                                                                 <  2:554313/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:541298/1‑150 (MQ=255)
                                                                        gaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCG‑‑GTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCgg                                                                             >  2:761909/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:596625/1‑150 (MQ=255)
                                                                             aTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCG‑‑GTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGt                                                                        >  1:1665311/1‑150 (MQ=255)
                                                                              tGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCG‑‑GTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACgaa                                                                       >  2:769578/1‑148 (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)
                                                                                    aCCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCG‑‑GTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAg                                                                 >  2:808296/1‑150 (MQ=255)
                                                                                     ccACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCG‑‑GTGTGGAAAACGGTAGTATTAGCAGCCACGAg                                                                                 <  2:671596/133‑1 (MQ=255)
                                                                                     ccACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCG‑‑GTGTGGAAAACGGTAGTATTAGCAGCCACGAg                                                                                 >  1:671596/1‑133 (MQ=255)
                                                                                        caGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCG‑‑GTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTg                                                             >  2:42593/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)
                                                                                                    ttATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCG‑‑GTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGt                                                 <  1:850606/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:1575175/1‑150 (MQ=255)
                                                                                                          cAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCG‑‑GTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATggg                                           <  1:979368/150‑1 (MQ=255)
                                                                                                               gCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCG‑‑GTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTc                                      >  2:1129983/1‑150 (MQ=255)
                                                                                                                   aCAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCG‑‑GTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCgggg                                  <  2:1221590/150‑1 (MQ=255)
                                                                                                                    cAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCG‑‑GTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCggggg                                 >  1:1620401/1‑150 (MQ=255)
                                                                                                                    cAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCG‑‑GTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCggggg                                 <  2:588422/150‑1 (MQ=255)
                                                                                                                     aGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCG‑‑GTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGt                                >  2:1343609/1‑150 (MQ=255)
                                                                                                                        cGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCG‑‑GTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCgg                             >  2:788474/1‑150 (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)
                                                                                                                                 ttttATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCG‑‑GTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCaaa                    <  1:599614/150‑1 (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:984716/1‑150 (MQ=255)
                                                                                                                                           aaCTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCG‑‑GTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTg          <  1:99204/150‑1 (MQ=255)
                                                                                                                                             cTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCG‑‑GTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGcc        <  2:332184/150‑1 (MQ=255)
                                                                                                                                              tGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCG‑‑GTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCg       <  2:504497/150‑1 (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)
                                                                                                                                                   |                                                                                                                                                       
CACCCCGAAAAAGACGCTGGTACTGAACAATAACGGTAAATACCGCACTTACAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCG‑‑GTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAA  >  NC_000913/2285992‑2286288

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: