Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,202,959 T→A 5.9% intergenic (‑26/‑65) ymfK ← / → ymfT e14 prophage; putative repressor protein YmfK/e14 prophage; putative DNA‑binding transcriptional regulator YmfT

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,202,9590TA 5.9% 510.1 / 14.0 162intergenic (‑26/‑65)ymfK/ymfTe14 prophage; putative repressor protein YmfK/e14 prophage; putative DNA‑binding transcriptional regulator YmfT
Reads supporting (aligned to +/- strand):  ref base T (76/76);  new base A (8/2);  total (84/78)
Fisher's exact test for biased strand distribution p-value = 1.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 165 total aligned reads displayed.

CTCCCATAAACCTGCTGGTCTGGGTTGGCTCGCGATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGCTCAACACCTGGCTACTTACGGC  >  NC_000913/1202810‑1203108
                                                                                                                                                     |                                                                                                                                                     
ctcCCATAAACCTGCTGGTCTGGGTTGGCTCGCGATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGtt                                                                                                                                                       >  2:1277747/1‑150 (MQ=255)
  cccATAAACCTGCTGGTCTGGGTTGGCTCGCGATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTaa                                                                                                                                                     >  1:967244/1‑150 (MQ=255)
  cccATAAACCTGCTGGTCTGGGTTGGCTCGCGATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTaa                                                                                                                                                     >  2:1539218/1‑150 (MQ=255)
     aTAAACCTGCTGGTCTGGGTTGGCTCGCGATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGt                                                                                                                                                  <  1:273370/150‑1 (MQ=255)
      tAAACCTGCTGGTCTGGGTTGGCTCGCGATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTa                                                                                                                                                 >  2:73623/1‑150 (MQ=255)
        aaCCTGCTGGTCTGGGTTGGCTCGCGATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACa                                                                                                                                               <  1:2150817/150‑1 (MQ=255)
          ccTGCTGGTCTGGGTTGGCTCGCGATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAg                                                                                                                                             <  2:989275/150‑1 (MQ=255)
           ctgctgGTCTGGGTTGGCTCGCGATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCTTAACGGTACAAgg                                                                                                                                            <  2:1427225/150‑2 (MQ=255)
           ctgctgGTCTGGGTTGGCTCGCGATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGt                                                                                                                                            >  2:1778329/1‑150 (MQ=255)
                 gTCTGGGTTGGCTCGCGATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTGAAACCCTTGTACCGTTAAGGTACAAGTATCTTg                                                                                                                                      <  1:1470513/150‑1 (MQ=255)
                   cTGGGTTGGCTCGCGATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGaa                                                                                                                                    <  1:2471152/150‑1 (MQ=255)
                     gggTTGGCTCGCGATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAg                                                                                                                                             >  1:1316441/1‑139 (MQ=255)
                      ggTTGGCTCGCGATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGt                                                                                                                                 >  2:1820859/1‑150 (MQ=255)
                         tGGCTCGCGATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTGAGGTACAAGTATCTTGAAGGTTCa                                                                                                                              >  2:2165038/1‑150 (MQ=255)
                           gCTCGCGATCAATCATAGTGGCAAAGGAATAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCAtt                                                                                                                            <  1:1850907/150‑1 (MQ=255)
                            cTCGCGATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCAttt                                                                                                                           >  1:2243945/1‑150 (MQ=255)
                            cTCGCGATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCAttt                                                                                                                           >  2:1226587/1‑150 (MQ=255)
                                    caatcaTAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATg                                                                                                                   >  2:370549/1‑150 (MQ=255)
                                     aatcaTAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGt                                                                                                                  >  1:223742/1‑150 (MQ=255)
                                     aatcaTAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGt                                                                                                                  >  1:740049/1‑150 (MQ=255)
                                     aatcaTAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGt                                                                                                                  <  1:1380957/150‑1 (MQ=255)
                                      atcaTAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTa                                                                                                                 >  1:1376613/1‑150 (MQ=255)
                                      atcaTAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTATTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTa                                                                                                                 <  2:145801/150‑1 (MQ=255)
                                       tcaTAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTaa                                                                                                                >  1:1144224/1‑150 (MQ=255)
                                              ggCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTac                                                                                                         >  1:1390707/1‑150 (MQ=255)
                                               gCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTaca                                                                                                        >  1:13537/1‑150 (MQ=255)
                                                cAAAGGAAGAATTCCCGCAAACACCATCTCTCAGTTTTCTGGCGTTAGCCCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATTTGTacac                                                                                                       >  2:2392011/1‑150 (MQ=255)
                                                 aaaGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCGTTTCAATCATGTAATATGTACAcc                                                                                                      <  1:2452048/150‑1 (MQ=255)
                                                   aGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCgg                                                                                                    <  1:2156243/150‑1 (MQ=255)
                                                   aGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCgg                                                                                                    >  2:549928/1‑150 (MQ=255)
                                                        gaaTTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTa                                                                                               <  2:607905/150‑1 (MQ=255)
                                                            tCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACata                                                                                           <  1:1922704/150‑1 (MQ=255)
                                                               cGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATAttt                                                                                        >  1:399278/1‑149 (MQ=255)
                                                                 ccAACACCATCTCTCAGTTTTCTTGCTTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTa                                                                                      <  1:1224147/150‑1 (MQ=255)
                                                                 ccAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACCTCCGGTGTACATAttaca                                                                                      <  2:509213/150‑4 (MQ=255)
                                                                 ccAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTa                                                                                      <  1:2283513/150‑1 (MQ=255)
                                                                 ccAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTa                                                                                      <  1:447533/150‑1 (MQ=255)
                                                                 ccAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTa                                                                                      <  2:168227/150‑1 (MQ=255)
                                                                  cAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACCTCCGGTGTACATAttacatg                                                                                    <  2:2374169/150‑6 (MQ=255)
                                                                    acacCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACCTCCGGTGTACATATTACATGa                                                                                   <  2:662801/150‑1 (MQ=255)
                                                                       ccATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAg                                                                                <  2:423175/150‑1 (MQ=255)
                                                                               cAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTgg                                                                        <  1:92404/150‑1 (MQ=255)
                                                                               cAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTgg                                                                        >  2:1418114/1‑150 (MQ=255)
                                                                                aGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTggg                                                                       <  2:770358/150‑1 (MQ=255)
                                                                                 gTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGa                                                                      >  2:964193/1‑150 (MQ=255)
                                                                                  ttttCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGAc                                                                     <  1:1778329/150‑1 (MQ=255)
                                                                                   tttCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGAct                                                                    >  1:212089/1‑150 (MQ=255)
                                                                                    ttCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCCATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGActc                                                                   <  2:1193122/150‑1 (MQ=255)
                                                                                    ttCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGActc                                                                   <  2:216134/150‑1 (MQ=255)
                                                                                      cTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCtt                                                                 >  1:1602682/1‑150 (MQ=255)
                                                                                        ggCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTa                                                               >  1:2123914/1‑150 (MQ=255)
                                                                                           ggtacaCCGCCGGATGTCATGGATTGTTTTCATCACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGTTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAAcc                                                            <  2:2315222/145‑1 (MQ=255)
                                                                                           gTTAGCCCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTATTTAAcc                                                            >  2:478417/1‑150 (MQ=255)
                                                                                           gTTAGACCGCCGTATTTCATGGATTGTTTTCATAACGAAATTAAAAGACTTGTACCTTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAAcc                                                            <  1:315738/150‑1 (MQ=255)
                                                                                           gTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAAcc                                                            <  2:42187/150‑1 (MQ=255)
                                                                                           gTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAAcc                                                            <  2:13537/150‑1 (MQ=255)
                                                                                               gACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAg                                                        >  2:2291458/1‑150 (MQ=255)
                                                                                                 ccgccgGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGaa                                                      >  1:1166116/1‑150 (MQ=255)
                                                                                                   gccgGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAAca                                                    <  1:726221/150‑1 (MQ=255)
                                                                                                    ccgGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAAcag                                                   <  2:763674/150‑1 (MQ=255)
                                                                                                     cgGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAAcagc                                                  <  1:2063872/150‑1 (MQ=255)
                                                                                                       gATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAAcagcag                                                >  2:562314/1‑150 (MQ=255)
                                                                                                          gTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAg                                                        <  2:2470509/139‑1 (MQ=255)
                                                                                                          gTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGc                                             >  2:2036034/1‑150 (MQ=255)
                                                                                                          gTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGc                                             <  1:1277747/150‑1 (MQ=255)
                                                                                                             aTGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAg                                          >  1:536450/1‑150 (MQ=255)
                                                                                                                   tGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGcc                                    <  1:2036034/150‑1 (MQ=255)
                                                                                                                    gTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCg                                   <  1:243450/150‑1 (MQ=255)
                                                                                                                    gTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCg                                   <  2:1941331/150‑1 (MQ=255)
                                                                                                                    gTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCg                                   >  2:1281859/1‑150 (MQ=255)
                                                                                                                     ttttCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCgg                                  <  1:1571436/150‑1 (MQ=255)
                                                                                                                     ttttCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTAGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCgg                                  <  2:213156/150‑1 (MQ=255)
                                                                                                                      tttCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGa                                 >  2:1571359/1‑150 (MQ=255)
                                                                                                                          aTAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAg                             <  1:1820859/150‑1 (MQ=255)
                                                                                                                           tAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGt                            >  2:728475/1‑150 (MQ=255)
                                                                                                                           tAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGt                            >  1:2447163/1‑150 (MQ=255)
                                                                                                                           tAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCGTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGt                            >  1:373804/1‑150 (MQ=255)
                                                                                                                            aaCGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGtt                           >  2:153488/1‑150 (MQ=255)
                                                                                                                             aCGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTg                          >  2:2191432/1‑150 (MQ=255)
                                                                                                                             aCGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTg                          >  1:2323597/1‑150 (MQ=255)
                                                                                                                             aCGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTAAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTg                          >  1:232502/1‑150 (MQ=255)
                                                                                                                               gAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGc                        >  1:1066538/1‑150 (MQ=255)
                                                                                                                               gAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGc                        <  2:967244/150‑1 (MQ=255)
                                                                                                                                 aaTTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGCTc                      <  1:2205398/150‑1 (MQ=255)
                                                                                                                                  aTTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGCTCa                     <  2:877684/150‑1 (MQ=255)
                                                                                                                                   ttAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGCTCaa                    >  1:2408089/1‑150 (MQ=255)
                                                                                                                                     aaaaCCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGCTCAaca                  <  2:1637001/150‑1 (MQ=255)
                                                                                                                                       gatacTTGTACCTTAACGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGCTCAACAcc                >  1:621736/5‑150 (MQ=255)
                                                                                                                                       aaCCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGCTCAACAcc                >  2:1207519/1‑150 (MQ=255)
                                                                                                                                         tacTTGTACCTTAACGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGGTCAACACCTg              >  2:753991/3‑150 (MQ=255)
                                                                                                                                         cccTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGCTCAACACCTg              <  2:1108131/150‑1 (MQ=255)
                                                                                                                                          acTTGTACCTTAACGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGCTCAACACCTgg             >  1:2029756/2‑150 (MQ=255)
                                                                                                                                          acTTGTACCTTAACGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGCTCAACACCTgg             >  2:1030511/2‑150 (MQ=255)
                                                                                                                                          acTTGTACCTTAACGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGCTCAACACCTgg             >  2:2010468/2‑150 (MQ=255)
                                                                                                                                           cTTGTACCTTAACGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGCTCAACACCTGGc            >  2:2271399/1‑150 (MQ=255)
                                                                                                                                           cTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGCTCAACACCTGGc            <  1:378051/150‑1 (MQ=255)
                                                                                                                                                  cGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGCTCAACACCTGGCTACTTAc     <  1:948535/150‑1 (MQ=255)
                                                                                                                                                   gTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGCTCAACACCTGGCTACTTACg    >  2:466267/1‑150 (MQ=255)
                                                                                                                                                     tAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGCTCAACACCTGGCTACTTACGGc  <  2:1066538/150‑1 (MQ=255)
                                                                                                                                                     tAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGCTCAACACCTGGCTACTTACGGc  <  2:292724/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CTCCCATAAACCTGCTGGTCTGGGTTGGCTCGCGATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGCTCAACACCTGGCTACTTACGGC  >  NC_000913/1202810‑1203108

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