Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,588,479 T→A 5.6% N121Y (AAT→TAT)  ydeS ← putative fimbrial protein YdeS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,588,4790TA 5.6% 384.4 / 7.7 120N121Y (AAT→TAT) ydeSputative fimbrial protein YdeS
Reads supporting (aligned to +/- strand):  ref base T (54/59);  new base A (5/2);  total (59/61)
Fisher's exact test for biased strand distribution p-value = 2.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 120 total aligned reads displayed.

TTGCCGAAGCCCTCACCAGTCCCGGATTGACGGACTTCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCACTTAACGTAATCTG  >  NC_000913/1588332‑1588626
                                                                                                                                                   |                                                                                                                                                   
ttGCCGAAGCCCTCACCAGTCCCGGATTGACGGACTTCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTg                                                                                                                                                   <  1:998344/150‑1 (MQ=255)
 tGCCGAAGCCCTCACCAGTCCCGGATTGACGGACTTCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGc                                                                                                                                                  <  2:1432343/150‑1 (MQ=255)
  gCCGAAGCCCTCACCAGTCCCGGATTGACGGACTTCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCg                                                                                                                                                 >  2:90745/1‑150 (MQ=255)
  gCCGAAGCCCTCACCAGTCCCGGATTGACGGACTTCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCg                                                                                                                                                 >  1:1487803/1‑150 (MQ=255)
     gAAGCCCTCACCAGTCCCGGATTGACGGACTTCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTcc                                                                                                                                              <  1:1230180/150‑1 (MQ=255)
         cccTCACCAGTCCCGGATTGACGGACTTCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAata                                                                                                                                          <  1:1385186/150‑1 (MQ=255)
           cTCACCAGTCCCGGATTGACGGACTTCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCAAATATGCGTccggtgaa                                                                                                                                        <  1:1561818/150‑7 (MQ=255)
            tCACCAGTCCCGGATTGACGGACTTCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATAttt                                                                                                                                       <  1:225371/150‑1 (MQ=255)
             cACCAGTCCCGGATTGACGGACTTCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCATTTTTATTCACTTCATGGCGTCATGATTATTATTCAGTTTCAAAGTATGTATATTTTCGTGAATTAttta                                                                                                                                      <  2:607401/150‑2 (MQ=255)
             cACCAGTCCCGGATTGACGGACTTCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTc                                                                                                                                      <  2:1389684/150‑1 (MQ=255)
             cACCAGTCCCGGATTGACGGACTTCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTc                                                                                                                                      >  1:98987/1‑150 (MQ=255)
             cACCAGTCCCGGATTGACGGACTTCTGAGTTGACTTCAGACGACCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGCCGCATATTTGCGTCAAATAtttg                                                                                                                                      >  2:833982/1‑149 (MQ=255)
              aCCAGTCCCGGATTGACGGACTTCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCt                                                                                                                                     <  2:1310721/150‑1 (MQ=255)
               ccAGTCCCGGATTGACGGACTTCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTc                                                                                                                                      >  1:1475272/1‑148 (MQ=255)
                   tCCCGGATTGACGGACTTCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATAcc                                                                                                                                <  1:681807/150‑1 (MQ=255)
                         aTTGACGGACTTCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAAcc                                                                                                                          <  2:29076/150‑1 (MQ=255)
                          ttGACGGACTTCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCg                                                                                                                         <  1:814315/150‑1 (MQ=255)
                             aCGGACTTCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAg                                                                                                                      <  2:1243868/150‑1 (MQ=255)
                              cGGACTTCTGAGTTGACTTCAGTCGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCTTCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAgg                                                                                                                     <  2:534444/150‑1 (MQ=255)
                                 aCTTCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCg                                                                                                                  <  1:1284794/150‑1 (MQ=255)
                                  cTTCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGt                                                                                                                 >  1:1035699/1‑150 (MQ=255)
                                    tCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTAt                                                                                                               >  1:1496279/1‑150 (MQ=255)
                                     cTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTAtt                                                                                                              <  1:849330/150‑1 (MQ=255)
                                     cTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTAtt                                                                                                              <  1:1457840/150‑1 (MQ=255)
                                     cTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTAtt                                                                                                              <  2:1371085/150‑1 (MQ=255)
                                      tGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTg                                                                                                             <  2:1370003/150‑1 (MQ=255)
                                        aGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGct                                                                                                           >  1:1578839/1‑150 (MQ=255)
                                              cTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCCCTGCATCCCGGCATGACGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGcttccttc                                                                                                     <  2:519993/150‑1 (MQ=255)
                                              cTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGcttccttc                                                                                                     >  1:1550522/1‑150 (MQ=255)
                                                tCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCAt                                                                                                   >  2:1292272/1‑150 (MQ=255)
                                                   gACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAgg                                                                                                                     >  1:900473/1‑129 (MQ=255)
                                                   gACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAgg                                                                                                                     <  2:900473/129‑1 (MQ=255)
                                                    aCGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAg                                                                                               >  2:165035/1‑150 (MQ=255)
                                                     cGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGt                                                                                              >  2:711859/1‑150 (MQ=255)
                                                        gCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTc                                                                                           <  2:1487803/150‑1 (MQ=255)
                                                        gCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTc                                                                                           >  2:1287727/1‑150 (MQ=255)
                                                          ggAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTcaa                                                                                         >  1:952803/1‑150 (MQ=255)
                                                          ggAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTcaa                                                                                         >  2:241209/1‑150 (MQ=255)
                                                          ggAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTAATTTTCACCTGACGCATATTTGAGTCCAATATTTCTATACCCAAAACGGATGCAGTATTGCTTCCTTCATCCAGTTTcaa                                                                                         >  1:451235/1‑150 (MQ=255)
                                                              tAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACaaa                                                                                     <  1:471271/150‑1 (MQ=255)
                                                                aGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGt                                                                                   >  2:949798/1‑150 (MQ=255)
                                                                 ggCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAgg                                                                                                                     >  2:874191/1‑115 (MQ=255)
                                                                 ggCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGta                                                                                  >  2:83598/1‑150 (MQ=255)
                                                                   cAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGtatt                                                                                >  2:1297543/1‑150 (MQ=255)
                                                                   cAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCAACGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGtatt                                                                                >  1:1202679/1‑150 (MQ=255)
                                                                      aaTATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTattatt                                                                             <  1:1280723/150‑1 (MQ=255)
                                                                        tatTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTAtttt                                                                           >  2:924336/1‑150 (MQ=255)
                                                                           tgttctgttcTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTg                                                                        <  2:438344/150‑1 (MQ=255)
                                                                               ctgttctgGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATc                                                                    >  1:1152277/1‑150 (MQ=255)
                                                                                 gttctgGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCt                                                                  >  2:148397/1‑150 (MQ=255)
                                                                                  ttctgGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTc                                                                 <  2:34525/150‑1 (MQ=255)
                                                                                    ctgGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCaa                                                               <  1:1287727/150‑1 (MQ=255)
                                                                                     tgGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAAt                                                              >  2:152022/1‑150 (MQ=255)
                                                                                       gTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATAc                                                            >  2:1608868/1‑150 (MQ=255)
                                                                                       gTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATAc                                                            <  1:165035/150‑1 (MQ=255)
                                                                                          ccAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGt                                                         >  1:908367/1‑150 (MQ=255)
                                                                                          ccAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGt                                                         <  2:1152277/150‑1 (MQ=255)
                                                                                           tgccgggATGCAGTGGATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGtt                                                        >  1:1083688/5‑150 (MQ=255)
                                                                                            aGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTa                                                       >  1:1094879/1‑150 (MQ=255)
                                                                                               gggATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAAt                                                    <  2:1270132/150‑1 (MQ=255)
                                                                                                ggATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATg                                                   >  1:1603066/1‑150 (MQ=255)
                                                                                                 gATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGc                                                  <  1:233365/150‑1 (MQ=255)
                                                                                                 gATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGc                                                  >  1:483553/1‑150 (MQ=255)
                                                                                                 gATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGc                                                  >  1:478005/1‑150 (MQ=255)
                                                                                                 gATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGc                                                  <  1:949798/150‑1 (MQ=255)
                                                                                                  aTCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCa                                                 <  2:1603066/150‑1 (MQ=255)
                                                                                                    ccACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAAc                                               >  1:181332/1‑150 (MQ=255)
                                                                                                      aCTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAAccc                                             <  2:1616843/150‑1 (MQ=255)
                                                                                                      aCTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAAccc                                             <  2:276370/150‑1 (MQ=255)
                                                                                                       cTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCg                                            <  1:1297543/150‑1 (MQ=255)
                                                                                                         gCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGaa                                          >  1:87354/1‑150 (MQ=255)
                                                                                                           aTCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAAcc                                        <  1:1004392/150‑1 (MQ=255)
                                                                                                           aTCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGATTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAAcc                                        <  2:1469887/150‑1 (MQ=255)
                                                                                                             cccGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAAAACCGTTAAATGCAACCCGAAcccc                                      <  1:273712/150‑1 (MQ=255)
                                                                                                                   aTGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGt                                <  1:90745/150‑1 (MQ=255)
                                                                                                                          tCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCGTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGCACCCCCGTTGTCCCTTtg                         >  2:564307/1‑150 (MQ=255)
                                                                                                                          tCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTtg                         >  1:643437/1‑150 (MQ=255)
                                                                                                                           cATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTtgc                        <  2:87354/150‑1 (MQ=255)
                                                                                                                            aTTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTtgct                       >  1:540343/1‑150 (MQ=255)
                                                                                                                            aTTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTtg                         >  1:1547519/1‑148 (MQ=255)
                                                                                                                            aTTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTtg                         <  2:1547519/148‑1 (MQ=255)
                                                                                                                             ttCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTtgctg                      <  1:711859/150‑1 (MQ=255)
                                                                                                                               cAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCa                    >  2:501812/1‑150 (MQ=255)
                                                                                                                               cAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCa                    >  2:1296196/1‑150 (MQ=255)
                                                                                                                                aGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCAt                   <  2:1525165/150‑1 (MQ=255)
                                                                                                                                  tttCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTc                 <  2:98987/150‑1 (MQ=255)
                                                                                                                                  tttCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTc                 <  1:83598/150‑1 (MQ=255)
                                                                                                                                    atattggACGCAAATATGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCAc               >  1:1247592/6‑150 (MQ=255)
                                                                                                                                     cACCGGACGCATATTTGCGTGCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGAAACCCGAACCCCCGTTGTCGCTTTGCTGCATTCACt              >  2:944061/1‑150 (MQ=255)
                                                                                                                                      aCCGGACGCATATATGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTACATCCAGTTTCAACAAAGTATTATTTACTGCATCCTCAAAACCGTTAAATCCAACCCAAACCCCCTTTCTCCCTTTACTGCATTCACtt             >  1:681186/1‑150 (MQ=255)
                                                                                                                                       ttggACGCAAATATGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCa                    >  2:64884/3‑142 (MQ=255)
                                                                                                                                       ttggACGCAAATATGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCa                    <  1:64884/140‑1 (MQ=255)
                                                                                                                                       ccGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCACTTa            >  1:1263899/1‑150 (MQ=255)
                                                                                                                                         ggACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCACTTAAc          <  2:1629146/150‑1 (MQ=255)
                                                                                                                                          gACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGACCCCCCGTTGTCCCTTTGCTGCATTCACTTAACg         >  2:960423/1‑150 (MQ=255)
                                                                                                                                          gACGCAAATATGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCACTTAACg         >  2:349281/1‑150 (MQ=255)
                                                                                                                                           aCGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCACTTAACGt        <  1:864903/150‑1 (MQ=255)
                                                                                                                                              cATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCACTTAACGTAAt     >  1:1279072/1‑150 (MQ=255)
                                                                                                                                              cATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCACTTAACGTAAt     >  1:1215771/1‑150 (MQ=255)
                                                                                                                                                 atTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCACTTAACGTAATCTg  <  2:952803/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TTGCCGAAGCCCTCACCAGTCCCGGATTGACGGACTTCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCACTTAACGTAATCTG  >  NC_000913/1588332‑1588626

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: ‑