Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,286,142 T→A 5.0% intergenic (+6/‑69) 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,1420TA 5.0% 312.8 / 3.8 110intergenic (+6/‑69)yejM/proLessential inner membrane DUF3413 domain‑containing protein; lipid A production and membrane permeability factor/tRNA‑Pro
Reads supporting (aligned to +/- strand):  ref base T (49/55);  new base A (4/2);  total (53/57)
Fisher's exact test for biased strand distribution p-value = 4.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 119 total aligned reads displayed.

ACCCCGAAAAAGACGCTGGTACTGAACAATAACGGTAAATACCGCACTTACAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCG‑‑GTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAA  >  NC_000913/2285993‑2286291
                                                                                                                                                     |                                                                                                                                                       
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)
                                                                                                                                                     |                                                                                                                                                       
ACCCCGAAAAAGACGCTGGTACTGAACAATAACGGTAAATACCGCACTTACAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCG‑‑GTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAA  >  NC_000913/2285993‑2286291

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: