Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,383,959 A→T 10.4% intergenic (+35/+36) elaD → / ← yfbK protease ElaD/IPR002035/DUF3520 domain‑containing protein YfbK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,383,9590AT10.4% 303.2 / 20.8 106intergenic (+35/+36)elaD/yfbKprotease ElaD/IPR002035/DUF3520 domain‑containing protein YfbK
Reads supporting (aligned to +/- strand):  ref base A (38/57);  new base T (6/5);  total (44/62)
Fisher's exact test for biased strand distribution p-value = 5.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.15e-01

Only 100 of 111 total aligned reads displayed.

CAGAAATTGTTATTTAGATAAATATGGCGATGCAAATATCAATCATTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGCCTGCTGAGCCCACTGTTTGATCTGCTG  >  NC_000913/2383810‑2384108
                                                                                                                                                     |                                                                                                                                                     
cAGAAATTGTTATTTAGATAAATATGGCGATGCAAATATCAATCATTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAata                                                                                                                                                       <  2:300236/150‑1 (MQ=255)
  gAAATTGTTATTTAGATAAATATGGCGATGCAAATATCAATCATTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCATATTTg                                                                                                                                                     <  1:1199388/150‑1 (MQ=255)
    aaTTGTTATTTAGATAAATATGGCGATGCAAATATCAATCATTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCt                                                                                                                                                   <  1:67245/150‑1 (MQ=255)
     aTTGTTATTTAGATAAATATGGCGATGCAAATATCAATCATTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTg                                                                                                                                                  <  2:954384/150‑1 (MQ=255)
     aTTGTTATTTAGATAAATATGGCGATGCAAATATCAATCATTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTg                                                                                                                                                  <  2:948735/150‑1 (MQ=255)
         ttATTTAGATAAATATGGCGATGCAAATATCAATCATTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCATATTTGCTGACTg                                                                                                                                              <  2:722985/150‑1 (MQ=255)
         ttATTTAGATAAATATGGCGATGCAAATATCAATCATTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTg                                                                                                                                              <  1:1174839/150‑1 (MQ=255)
           aTTTAGATAAATATGGCGATGCAAATATCAATCATTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCATATTTGCTGACTGGt                                                                                                                                            <  2:64556/150‑1 (MQ=255)
           aTTTAGATAAATATGGCGATGCAAATATCAATCATTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGt                                                                                                                                            >  2:420154/1‑150 (MQ=255)
           aTTTAGATAAATATGGCGATGCAAATATCAATCATTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGt                                                                                                                                            >  1:515062/1‑150 (MQ=255)
            tttAGATAAATATGGCGATGCAAATATCAATCATTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTc                                                                                                                                           <  1:331793/150‑1 (MQ=255)
              tAGATAAATATGGCGATGCAAATATCAATCATTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCATATTTGCTGACTGGTCGa                                                                                                                                         <  2:1049071/150‑1 (MQ=255)
              tAGATAAATATGGCGATGCAAATATCAATCATTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGa                                                                                                                                         >  2:901257/1‑150 (MQ=255)
              tAGATAAATATGGCGATGCAAATATCAATCATTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGa                                                                                                                                         >  2:48967/1‑150 (MQ=255)
              tAGATAAATATGGCGATGCAAATATCAATCATTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGa                                                                                                                                         >  1:263263/1‑150 (MQ=255)
               aGATAAATATGGCGATGCAAATATCAATCATTATTATAGAACCTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGaa                                                                                                                                        >  2:1167788/1‑150 (MQ=255)
                  tAAATATGGCGATGCAAATATCAATCATTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACAt                                                                                                                                     >  1:614468/1‑150 (MQ=255)
                    aaTATGGCGATGCAAATATCAATCATTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAg                                                                                                                                   <  1:725082/150‑1 (MQ=255)
                    aaTATGGCCATGCAAATATCAATCATTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAg                                                                                                                                   <  2:311453/150‑1 (MQ=255)
                     atatGGCGATGCAAATATCAATCATTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGt                                                                                                                                  <  2:935991/150‑1 (MQ=255)
                      tatGGCGATGCAAATATCAATCATTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCATATTTGCTGACTGGTCGaattaatta                                                                                                                                 <  2:1041754/150‑8 (MQ=255)
                         ggCGATGCAAATATCAATCATTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGc                                                                                                                              <  2:282023/150‑1 (MQ=255)
                         ggCGATGCAAATATCAATCATTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGc                                                                                                                              <  1:1324514/150‑1 (MQ=255)
                          gCGATGCAAATATCAATCATTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCt                                                                                                                             <  2:1308676/150‑1 (MQ=255)
                             aTGCAAATATCAATCATTATTATAGAAACTTAGAAATAAAAAACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTaaa                                                                                                                          >  2:1375694/1‑150 (MQ=255)
                              tGCAAATATCAATCATTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAAc                                                                                                                         <  1:968740/150‑1 (MQ=255)
                               gCAAATATCAATCATTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACa                                                                                                                        <  2:1376880/150‑1 (MQ=255)
                                cAAATATCAATCATTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAg                                                                                                                       <  2:537632/150‑1 (MQ=255)
                                 aaaTATCAATCATTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGt                                                                                                                      <  2:62357/150‑1 (MQ=255)
                                   atatCAATCATTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGtca                                                                                                                    <  2:1290471/150‑1 (MQ=255)
                                      tcaatcaTTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGCGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGtcatca                                                                                                                 <  2:1380256/150‑1 (MQ=255)
                                       caatcaTTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCac                                                                                                                >  1:1514084/1‑150 (MQ=255)
                                           caTTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCactgac                                                                                                            <  1:1488619/150‑1 (MQ=255)
                                            attattATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGTCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCactgact                                                                                                           <  1:1511256/150‑1 (MQ=255)
                                             ttattaTAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCactgactg                                                                                                          <  2:515062/150‑1 (MQ=255)
                                             ttattaTAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCactgactg                                                                                                          <  1:638030/150‑1 (MQ=255)
                                             ttattaTAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCactgactg                                                                                                          <  1:258728/150‑1 (MQ=255)
                                              tattatAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGAAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCActgactga                                                                                                         >  2:1174993/1‑150 (MQ=255)
                                               attatAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGAt                                                                                                        >  2:1343416/1‑150 (MQ=255)
                                                ttatAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATg                                                                                                       <  2:1367280/150‑1 (MQ=255)
                                                ttatAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATg                                                                                                       <  2:118139/150‑1 (MQ=255)
                                                  ataGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATgtc                                                                                                     <  2:522041/150‑1 (MQ=255)
                                                  ataGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATgtc                                                                                                     >  1:496554/1‑150 (MQ=255)
                                                  ataGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATgtc                                                                                                     <  2:1530156/150‑1 (MQ=255)
                                                  ataGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATgtc                                                                                                     <  1:119971/150‑1 (MQ=255)
                                                    aGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGATTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATgtcag                                                                                                   <  1:865423/150‑1 (MQ=255)
                                                     gAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATgtcagt                                                                                                  >  1:358967/1‑150 (MQ=255)
                                                     gAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATgtcagt                                                                                                  >  1:379450/1‑150 (MQ=255)
                                                     gAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATgtcagt                                                                                                  >  1:19348/1‑150 (MQ=255)
                                                      aaaCTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATgtcagtc                                                                                                 <  2:488096/150‑1 (MQ=255)
                                                        aCTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTcac                                                                                               <  2:808660/150‑1 (MQ=255)
                                                          ttAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCacac                                                                                             >  2:774729/1‑150 (MQ=255)
                                                           tAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACAcc                                                                                            <  1:947002/150‑1 (MQ=255)
                                                             gAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATgtc                                                                                                     >  1:763798/1‑139 (MQ=255)
                                                             gAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATgtc                                                                                                     <  2:763798/139‑1 (MQ=255)
                                                               aaTAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATcc                                                                                        >  2:46312/1‑150 (MQ=255)
                                                                  aaaaCACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGcc                                                                                     <  2:379450/150‑1 (MQ=255)
                                                                   aaaCACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAg                                                                                                                       <  2:1242506/115‑1 (MQ=255)
                                                                   aaaCACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAg                                                                                                                       >  1:1242506/1‑115 (MQ=255)
                                                                          cacaacacaaAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAAc                                                                                                                                       <  2:1019449/92‑1 (MQ=255)
                                                                          cacaACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAAc                                                                                                                                       >  1:1019449/1‑92 (MQ=255)
                                                                          cacaACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAAt                                                                             >  1:473923/1‑150 (MQ=255)
                                                                          cacaACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAAt                                                                             >  2:627481/1‑150 (MQ=255)
                                                                          cacaACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGATGACTGCTCGAACATAGTTTGCTAAAAAGTCATAACTGACTGAAGTCAGTCACACAATACGCCAGTTCAAt                                                                             >  1:737295/1‑150 (MQ=255)
                                                                              aCCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAgg                                                                         <  1:848064/150‑1 (MQ=255)
                                                                              aCCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAgg                                                                         >  2:1491869/1‑150 (MQ=255)
                                                                               cccAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGc                                                                        <  2:414199/150‑1 (MQ=255)
                                                                                    aaaTAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATa                                                                   <  2:1064832/150‑1 (MQ=255)
                                                                                       tAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAAt                                                                >  1:431730/1‑150 (MQ=255)
                                                                                              tCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGccc                                                         <  1:1478523/150‑1 (MQ=255)
                                                                                              tCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGccc                                                         <  1:901257/150‑1 (MQ=255)
                                                                                                     aaaGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAAcc                                                  >  1:906364/1‑150 (MQ=255)
                                                                                                      aaGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAAccc                                                 >  1:514854/1‑150 (MQ=255)
                                                                                                      aaGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAAccc                                                 <  1:1343416/150‑1 (MQ=255)
                                                                                                          gtgagtTAACATCATTTTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTAATCACTGACTGATGTCAGTCACAACATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCtgtg                                             <  2:432337/150‑1 (MQ=255)
                                                                                                                 aaCATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCAGCCAGTTCAATCAGGCGAATAAATTACGCCCTGTAACCCTGTGGATCTTc                                      >  1:476773/1‑150 (MQ=255)
                                                                                                                   catcatTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTccc                                    >  1:717931/1‑150 (MQ=255)
                                                                                                                   catcatTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTccc                                    >  1:979928/1‑150 (MQ=255)
                                                                                                                         tGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTg                              <  1:228491/150‑1 (MQ=255)
                                                                                                                         tGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTg                              <  1:774729/150‑1 (MQ=255)
                                                                                                                           ttaattaaTTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGcc                            >  2:1063324/1‑150 (MQ=255)
                                                                                                                           ttaattaaTTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGcc                            <  2:19348/150‑1 (MQ=255)
                                                                                                                             aattaattCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGCctg                          >  2:976316/1‑150 (MQ=255)
                                                                                                                                taattCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGCctgctg                       >  1:1333785/1‑150 (MQ=255)
                                                                                                                                       aCCAGTCAGCATATTTGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGCCTGCTGCGCCCAc                >  2:526034/1‑150 (MQ=255)
                                                                                                                                        ccAGTCAGCATATTTGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGCCTGCTGAGCCCACt               >  2:428362/1‑150 (MQ=255)
                                                                                                                                        ccAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGCCTGCTGAGCCCACt               >  2:264343/1‑150 (MQ=255)
                                                                                                                                         cAGTCAGCATATTTGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGCCTGCTGAGCCCACTg              >  1:423675/1‑150 (MQ=255)
                                                                                                                                         cAGTCAGCATATTTGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGCCTGCTGAGCCCACTg              >  1:275077/1‑150 (MQ=255)
                                                                                                                                          aGTCAGCATATTTGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGCCTGCTGAGCCCACTGt             >  2:262628/1‑150 (MQ=255)
                                                                                                                                             cAGCATATTTGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGCCTGCTGAGCCCACTGTTTg          >  2:237175/1‑150 (MQ=255)
                                                                                                                                               gCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGCCTGCTGAGCCCACTGTTTGAt        <  1:420154/150‑1 (MQ=255)
                                                                                                                                                cAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGCCTGCTGAGCCCACTGTTTGATc       >  2:60599/1‑150 (MQ=255)
                                                                                                                                                cAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGCCTGCTGAGCCCACTGTTTGATc       <  1:216115/150‑1 (MQ=255)
                                                                                                                                                 aaaTATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGCCTGCTGAGCCCACTGTTTGATct      <  2:473923/150‑1 (MQ=255)
                                                                                                                                                 aaaTATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGCCTGCTGAGCCCACTGTTTGATct      <  1:1174993/150‑1 (MQ=255)
                                                                                                                                                   atatGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGCCTGCTGAGCCCACTGTTTGATctgc    <  2:263263/150‑1 (MQ=255)
                                                                                                                                                   atatGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGCCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGCCTGCTGAGCCCACTGTTTGATctgc    <  1:1187928/150‑1 (MQ=255)
                                                                                                                                                    tatGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCCAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGCCTGCTGAGCCCACTGTTTGATctgct   <  2:540589/150‑1 (MQ=255)
                                                                                                                                                     atGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGCCTGCTGAGCCCACTGTTTGATctgctg  <  2:1128741/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CAGAAATTGTTATTTAGATAAATATGGCGATGCAAATATCAATCATTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGCCTGCTGAGCCCACTGTTTGATCTGCTG  >  NC_000913/2383810‑2384108

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: