Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,420,667:1 +C 5.1% coding (534/1182 nt) yhdX → putative ABC transporter membrane subunit YhdX

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,420,6671.C 5.1% 410.5 / 16.5 98coding (534/1182 nt)yhdXputative ABC transporter membrane subunit YhdX
Reads supporting (aligned to +/- strand):  ref base . (47/46);  new base C (2/3);  total (49/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.45e-01

Only 100 of 105 total aligned reads displayed.

GAGATCTTCCGTAATATTCCCCCGCTGCTGCAAATCTTCTTCTGGTACTTTGCCGTGTTGCGCAATTTGCCCGGACCACGCCAGGCTGTGAGCGCGTTTGATCTGGCCTTTTTGAGCAATCGTGGGCTTTATATTCCGTCACCGCAGCTGGGAGACGGATTTATTGCGTTTATCCTGGCTGTTGTTATGGCTATAGTCCTTTCTGTTGGGCTATTCCGCTTTAATAAAACATACCAGATAAAGACCGGACAACTGCGCCGCACCTGGCCGATCGCCGCAGTGTTGATCATTGGTTTGCC  >  NC_000913/3420518‑3420816
                                                                                                                                                      |                                                                                                                                                     
gagaTCTTCCGTAATATTCCCCCGCTGCTGCAAATCTTCTTCTGGTACTTTGCCGTGTTGCGCAATTTGCCCGGACCACGCCAGGCTGTGAGCGCGTTTGATCTGGCCTTTTTGAGCAATCGTGGGCTTTATATTCCGTCACCGCAGCTg                                                                                                                                                        <  2:1435084/150‑1 (MQ=255)
 agaTCTTCCGTAATATTCCCCCGCTGCTGCAAATCTTCTTCTGGTACTTTGCCGTGTTGCGCAATTTGCCCGGACCACGCCAGGCTGTGAGCGCGTTTGATCTGGCCTCTTTGAGCAATCGTGGGCTTTATATTCCGTCACCGCAGCTgg                                                                                                                                                      <  2:804115/150‑1 (MQ=255)
  taTCTTCCGTAATATTCCCCCGCTGCTGCAAATCTTCTTCTGGTACTTTGCCGTGTTGCGCAATTTGCCCGGACCACGCCAGGCTGTGAGCGCGTTTGATCTGGCCTTTTTGAGCAATCGTGGGCTTTATATTCCGTCACCGCAGCTggg                                                                                                                                                     <  2:204912/149‑1 (MQ=255)
  gaTCTTCCGTAATATTCCCCCGCTGCTGCAAATCTTCTTCTGGTACTTTGCCGTGTTGCGCAATTTGCCCGGACCACGCCAGGCTGTGAGCGCGTTTGATCTGGCCTTTTTGAGCAATCGTGGGCTTTATATTCCGTCACCGCAGCTggg                                                                                                                                                     <  1:1043703/150‑1 (MQ=255)
  gaTCTTCCGTAATATTCCCCCGCTGCTGCAAATCTTCTTCTGGTACTTTGCCGTGTTGCGCAATTTGCCCGGACCACGCCAGGCTGTGAGCGCGTTTGATCTGGCCTTTTTGAGCAATCGTGGGCTTTATATTCCGTCACCGCAGCTggg                                                                                                                                                     <  1:156427/150‑1 (MQ=255)
    tCTTCCGTAATATTCCCCCGCTGCTGCAAATCTTCTTCTGGTACTTTGCCGTGTTGCGCAATTTGCCCGGACCACGCCAGGCTGTGAGCGCGTTTGATCTGGCCTTTTTGAGCAATCGTGGGCTTTATATTCCGTCACCGCAGCTGGgag                                                                                                                                                   >  1:627094/1‑150 (MQ=255)
        gcGTAATTTTCCCCCGCTGCTGCAAATCTTCTTCTGGTACTTTGCCGTGTTGCGCAATTTGCCCTGACCACGCCAGGCTGTGAGCGCGTTTGATCTGGCCTTTTTGAGCAATCGTGGGCTTTATATTCCGTCACCGCAGCTGGGAGACgg                                                                                                                                               <  2:1084444/149‑1 (MQ=255)
          gTAATATTCCCCCGCTGCTGCAAATCTTCTTCTGGTACTTTGCCGTGTTGCGCAATTTGCCCGGACCACGCCAGGCTGTGAGCGCGTTTGATCTGGCCTTTTTGAGCAATCGTGGGCTTTATATTCCGTCTCCCAGCTGCGGTGACGGaa                                                                                                                                             <  1:1344152/150‑2 (MQ=255)
          gTAATATTCCCCCCCTGCTGCAAATCTTCTTCTGGTACTTTGCCGTGTTGCGCAATTTGCCCGGACAACGACAGGCTGTGAGCGCGTTTGATCTGGACTTTTTGAGCAATCGTTGGCTTTATATTCCGTCACCGCAGCTGGGAGACGGAt                                                                                                                                             >  2:1231132/1‑150 (MQ=255)
           tAATATTCCCCCGCTGCTGCAAATCTTCTTCTGGTACTTTGCCGTGTTGCGCAATTTGCCCGGACCACGCCAGGCTGTGAGCGCGTTTGATCTGGCCTTTTTGAGCAATCGTGGGCTTTATATTCCGTCTCCCAGCTGCGGTGACGGaat                                                                                                                                            <  1:745105/150‑3 (MQ=255)
             atatTCCCCCGCTGCTGCAAATCTTCTTCTGGTACTTTGCCGTGTTGCGCAATTTGCCCGGACCACGCCAGGCTGTGAGCGCGTTTGATCTGGCCTTTTTGAGCAATCGTGGGCTTTATATTCCGTCACCGCAGCTGGGAGACGGATTTa                                                                                                                                          <  1:541737/150‑1 (MQ=255)
              tatTCCCCCGCTGCTGCAAATCTTCTTCTGGTACTTTGCCGTGTTGCGCAATTTGCCCGGACCACGCCAGGCTGTGAGCGCGTTTGATCTGGCCTTTTTGAGCAATCGTGGGCTTTATATTCCGTCTCCCAGCTGCGGTGACGGaatata                                                                                                                                         <  1:1266776/150‑6 (MQ=255)
              tatTCCCCCGCTGCTGCAAATCTTCTTCTGGTACTTTGCCGTGTTGCGCAATTTGCCCGGACCACGCCAGGCTGTGAGCGCGTTTGATCTGGCCTTTTTGAGCAATCGTGGGCTTTATATTCCGTCACCGCAGCTGGGAGACGGATTTAt                                                                                                                                         >  1:1526692/1‑150 (MQ=255)
                   ccccGCTGCTGCAAATCTTCTTCTGGTACTTTGCCGTGTTGCGCAATTTGCCCGGACCACGCCAGGCTGTGAGCGCGTTTGATCTGGCCTTTTTGAGCAATCGTGGGCTTTATATTCCGTCACCGCAGCTGGGAGACGGATTTATTGCGt                                                                                                                                    >  2:767100/1‑150 (MQ=255)
                   ccccGCTGCTGCAAATCTTCTTCTGGTACTTTCCCGTGTTGCGCAATTTGCCCGGACCACGCCAGGCTGTGAGCGCGTTTGATCTGGCCTTTTTGAGCAATCGTGGGCTTTATATTCCGTCACCGCAGCTGGGACACGGATTTATTGCGt                                                                                                                                    >  2:1122524/1‑150 (MQ=255)
                    cccGCTGCTGCAAATCTTCTTCTGGTACTTTGCCGTGTTGCGCAATTTGCCCGGACCACGCCAGGCTGTGAGCGCGTTTGATCTGGCCTTTTTGAGCAATCGTGGGCTTTATATTCCGTCACCGCAGCTGGGAGACGGATTTATTGCGtt                                                                                                                                   <  2:349545/150‑1 (MQ=255)
                           ctgcAAATCTTCTTCTGGTACTTTGCCGTGTTGCGCAATTTGCCCGGACCACGCCAGGCTGTGAGCGCGTTTGATCTGGCCTTTTTGAGCAATCGTGGGCTTTATATTCCGTCACCGCAGCTGGGAGACGGATTTATTGCGTTTATCCTg                                                                                                                            <  1:1422303/150‑1 (MQ=255)
                             gcAAATCTTCTTCTGGTACTTTGCCGTGTTGCGCAATTTGCCCGGACCACGCCAGGCTGTGAGCGCGTTTGATCTGGCCTTTTTGAGCAATCGTGGGCTTTATATTCCGTCACCGCAGCTGGGAGACGGATTTATTGCGTTTATCCTGGc                                                                                                                          <  1:105516/150‑1 (MQ=255)
                                 aTCTTCTTCTGGTACTTTGCCGTGTTGCGCAATTTGCCCGGACCACGCCAGGCTGTGAGCGCGTTTGATCTGGCCTTTTTGAGCAATCGTGGGCTTTATATTCCGTCACCGCAGCTGGGAGACGGATTTATTGCGTTTATCCTGGCtgtt                                                                                                                      >  1:53323/1‑150 (MQ=255)
                                 aTCTTCTTCTGGTACTTTGCCGTGTTGCGCAATTTGCCCGGACCACGCCAGGCTGTGAGCGCGTTTGATCTGGCCTTTTTGAGCAATCGTGGGCTTTATATTCCGTCACCGCAGCTGGGAGACGGATTTATTGCGTTTATCCTGGCtgtt                                                                                                                      >  1:672128/1‑150 (MQ=255)
                                  tcttcttctGGTACTTTGCCGTGTTGCGCAATTTGCCCGGACCACGCCAGGCTGTGAGCGCGTTTGATCTGGCCTTTTTGAGCAATCGTGGGCTTTATATTCCGTCACCGCAGCTGGGAGACGGATTTATTGCGTTTATCCTGGCtgttg                                                                                                                     <  2:607360/150‑1 (MQ=255)
                                    ttcttctGGTACTTTGCCGTGTTGCGCAATTTGCCCGGACCACGCCAGGCTGTGAGCGCGTTTGATCTGGCTTTTTGAGCAATCGTGGGCTTTATATTCCGTCACCGCAGCTGGGAGACGGATTTATTGCGTTTATCCTGGCTGTTGTTa                                                                                                                  <  1:463678/150‑1 (MQ=255)
                                     tcttctGGTACTTTGCCGTGTTGCGCAATTTGCCCGGACCACGCCAGGCTGTGAGCGCGTTTGATCTGGCCTTTTTGAGCAATCGTGGGCTTTATATTCCGTCACCGCAGCTGGGAGACGGATTTATTGCGTTTATCCTGGCTGTTGTTa                                                                                                                  >  2:265585/1‑150 (MQ=255)
                                     tcttctGGTACTTTGCCGTGTTGCGCAATTTGCCCGGACCACGCCAGGCTGTGAGCGCGTTTGATCTGGCCTTTTTGAGCAATCGTGGGCTTTATATTCCGTCACCGCAGCTGGGAGACGGATTTATTGCGTTTATCCTGGCTGTTGTTa                                                                                                                  <  1:442547/150‑1 (MQ=255)
                                        tctGGTACTTTGCCGTGTTGCGCAATTTGCCCGGACCACGCCAGGCTGTGAGCGCGTTTGATCTGGCCTTTTTGAGCAATCGTGGGCTTTATATTCCGTCACCGCAGCTGGGAGACGGATTTATTGCGTTTATCCTGGCTGTTGTTATgg                                                                                                               >  1:928941/1‑150 (MQ=255)
                                         ctGGTACTTTGCCGTGTTGCGCAATTTGCCCGGACCACGCCAGGCTGTGAGCGCGTTTGATCTGGCCTTTTTGAGCAATCGTGGGCTTTATATTCCGTCACCGCAGCTGGGAGACGGATTTATTGCGTTTATCCTGGCTGTTGTTATGGc                                                                                                              <  1:1331108/150‑1 (MQ=255)
                                          tGGTACTTTGCCGTGTTGCGCAATTTGCCCGGACCACGCCAGGCTGTGAGCGCGTTTGATCTGGCCTTTTTGAGCAATCGTGGGCTTTATATTCCGTCACCGCAGCTGGGAGACGGATTTATTGCGTTTATCCTGGCTGTTGTTATGGCt                                                                                                             >  1:936652/1‑150 (MQ=255)
                                          tGGTACTTTGCCGTGTTGCGCAATTTGCCCGGACCACGCCAGGCTGTGAGCGCGTTTGATCTGGCCTTTTTGAGCAATCGTGGGCTTTATATTCCGTCACCGCAGCTGGGAGACGGATTTATTGCGTTTATCCTGGCTGTTGTTATGGCt                                                                                                             <  1:918551/150‑1 (MQ=255)
                                           ggTACTTTGCCGTGTTGCGCAATTTGCCCGGACCACGCCAGGCTGTGAGCGCGTTTGATCTGGCCTTTTTGAGCAATCGTGGGCTTTATATTCCGTCACCGCAGCTGGGAGACGGATTTATTGCGTTTATCCTGGCTGTTGTTATGGCta                                                                                                            >  1:1006711/1‑150 (MQ=255)
                                            gTACTTTGCCGTGTTGCGCAATTTGCCCGGACCACGCCAGGCTGTGAGCGCGTTTGATCTGGCCTTTTTGAGCAATCGTGGGCTTTATATTCCGTCACCGCAGCTGGGAGACGGATTTATTGCGTTTATCCTGGCTGTTGTTATGGCtat                                                                                                           <  1:672534/150‑1 (MQ=255)
                                             tACTTTGCCGTGTTGCGCAATTTGCCCGGACCACGCCAGGCTGTGAGCGCGTTTGATCTGGCCTTTTTGAGCAATCGTGGGCTTTATATTCCGTCACCGCAGCTGGGAGACGGATTTATTGCGTTTATCCTGGCTGTTGTTATGGCtata                                                                                                          >  2:242357/1‑150 (MQ=255)
                                                tttGCCGTGTTGCGCAATTTGCCCGGACCACGCCAGGCTGTGAGCGCGTTTGATCTGGCCTTTTTGAGCAATCGTGGGCTTTATATTCCGTCACCGCAGCTGGGAGACGGATTTATTGCGTTTATCCTGGCTGTTGTTATGGCTATAGTc                                                                                                       >  2:473028/1‑150 (MQ=255)
                                                   gCCGTGTTGCGCAATTTGCCCGGACCACGCCAGGCTGTGAGCGCGTTTGATCTGGCCTTTTTGAGCAATCGTGGGCTTTATATTCCGTCACCGCAGCTGGGAGACGGATTTATTGCGTTTATCCTGGCTGTTGTTATGGCTATAGTCCtt                                                                                                    >  1:808354/1‑150 (MQ=255)
                                                   gCCGTGTTGCGCAATTTGCCCGGACCACGCCAGGCTGTGAGCGCGTTTGATCTGGCCTTTTTGAGCAATCGTGGGCTTTATATTCCGTCACCGCAGCTGGGAGACGGATTTATTGCGTTTATCCTGGCTGTTGTTATGGCTATAGTCCtt                                                                                                    <  1:1071307/150‑1 (MQ=255)
                                                      gggtTGCGCAATTTGCCCGGACCACGCCAGGCTGTGAGCGCGTTTGATCTGGCCTTTTTGAGCGATCGTGGGCTTTATATTCCGTCCCCGCAGCTGGGAGACGGATTTATTGCGTTTATCCTGGCTGTTGTTATGGCTATAGTCGTTTCt                                                                                                 >  2:559177/3‑150 (MQ=255)
                                                         ttGCGCAATTTGCCCGGACCACGCCAGGCTGTGAGCGCGTTTGATCTGGCCTTTTTGAGCAATCGTGGGCTTTATATTCCGTCACCGCAGCTGGGAGACGGATTTATTGCGTTTATCCTGGCTGTTGTTATGGCTATAGTCCTTTCTGtt                                                                                              <  2:683749/150‑1 (MQ=255)
                                                         atGCGCAATTTGCCCGGACCACGCCAGGCTGTGAGCGCGTTTGATCTGGCCTTTTTGAGCAATCGTGGGCTTTATATTCCGTCACCGCAGCTGGGAGACGGATTTATTGCGTTTATCCTGGCTGTTGTTATGGCTATAGTCCTTTCTGtt                                                                                              >  1:583226/2‑150 (MQ=255)
                                                          tGCGCAATTTGCCCGGACCACGCCAGGCTGTGAGCGCGTTTGATCTGGCCTTTTTGAGCAATCGTGGGCTTTATATTCCGTCACCGCAGCTGGGAGACGGATTTATTGCGTTTATCCTGGCTGTTGTTATGGCTATAGTCCTTTCTGTTg                                                                                             >  1:556807/1‑150 (MQ=255)
                                                          tGCGCAATTTGCCCGGACCACGCCAGGCTGTGAGCGCGTTTGATCTGGCCTTTTTGAGCAATCGTGGGCTTTATATTCCGTCACCGCAGCTGGGAGACGGATTTATTGCGTTTATCCTGGCTGTTGTTATGGCTATAGTCCTTTCTGTTg                                                                                             <  1:272769/150‑1 (MQ=255)
                                                           gcgcAATTTGCCCGGACCACGCCAGGCTGTGAGCGCGTTTGATCTGGCCTTTTTGAGCAATCGTGGGCTTTATATTCCGTCACCGCAGCTGGGAGACGGATTTATTGCGTTTATCCTGGCTGTTGTTATGGCTATAGTCCTTTCt                                                                                                 >  2:732282/1‑145 (MQ=255)
                                                           gcgcAATTTGCCCGGACCACGCCAGGCTGTGAGCGCGTTTGATCTGGCCTTTTTGAGCAATCGTGGGCTTTATATTCCGTCACCGCAGCTGGGAGACGGATTTATTGCGTTTATCCTGGCTGTTGTTATGGCTATAGTCCTTTCt                                                                                                 <  1:732282/145‑1 (MQ=255)
                                                           gcgcAATTTGCCCGGACCACGCCAGGCTGTGAGCGCGTTTGATCTGGCCTTTTTGAGCAATCGTGGGCTTTATATTCCGTCACCGCAGCTGGGAGACGGATTTATTGCGTTTATCCTGGCTGTTGTTATGGCTATAGTCCTTTCTGTTgg                                                                                            >  2:949143/1‑150 (MQ=255)
                                                           gcgcAATTTGCCCGGACCACGCCAGGCTGTGAGCGCGTTTGATCTGGCCTTTTTGAGCAATCGTGGGCTTTATATTCCGTCACCGCAGCTGGGAGACGGATTTATTGCGTTTATCCTGGCTGTTGTTATGGCTATAGTCCTTTCTGTTgg                                                                                            <  1:1113516/150‑1 (MQ=255)
                                                           gcgcAATTTGCCCGGACCACGCCAGGCTGTGAGCGCGTTTGATCTGGCCTTTTTGAGCAATCGTGGGCTTTATATTCCGTCACCGCAGCTGGGAGACGGATTTATTGCGTTTATCCTGGCTGTTGTTATGGCTATAGTCCTTTCTGTTgg                                                                                            >  1:532987/1‑150 (MQ=255)
                                                           gcgcAATTTGCCCGGACCACGCCAGGCTGTGAGCGCGTTTGATCTGGCCTTTTTGAGCAATCATGGGCTTTATATTCCGTCACCGCAGCTGGGAGACGGATTTATTGCGTTTATCCTGGCTGTTGTTATGGCTATAGTCCTTTCTGTTgg                                                                                            >  1:620792/1‑150 (MQ=255)
                                                            cgcAATTTGCCCGGACCACGCCAGGCTGTGAGCGCGTTTGATCTGGCCTTTTTGAGCAATCGTGGGCTTTATATTCCGTCACCGCAGCTGGGAGACGGATTTATTGCGTTTATCCTGGCTGTTGTTATGGCTATAGTCCTTTCTGTTggg                                                                                           <  2:906578/150‑1 (MQ=255)
                                                             gcAATTTGCCCTGACCACGCCAGGGTGTGAGAGCGTTTGATCTGGCCTTTTTGAGCAATCGTGGGTTTTATATTCCGTCACCTCTGCTGTGTGAAGGATTTATTGCGTTTATCCTGGCTGTTGTTATGGCTATAGTCCTTTCTGTTGGGc                                                                                          <  2:790043/150‑1 (MQ=255)
                                                             gcAATTTGCCCGGACCACGCCAGGCTGTGAGCGCGTTTGATCTGGCCTTTTTGAGTAATCGTGGGCTTTATATTCCGTCACCGCAGCTGGGAGACGGATTTATTGCGTTTATCCTGGCTGTTGTTATGGCTATAGTCCTTTCTGTTGGGc                                                                                          <  1:1317360/150‑1 (MQ=255)
                                                                  ttGCCCGGACCACGCCAGGCTGTGAGCGCGTTTGATCTGGCCTTTTTGAGCAATCGTGGGCTTTATATTCCGTCACCGCAGCTGGGAGACGGATTTATTGCGTTTATCCTGGCTGTTGTTATGGCTATAGTCCTTTCTGTTGGGCTATTc                                                                                     <  1:242357/150‑1 (MQ=255)
                                                                    gCCCGGACCACGCCAGGCTGTGAGCGCGTTTGATCTGGCCTTTTTGAGCAATCGTGGGCTTTATATTCCGTCACCGCAGCTGGGAGACGGATTTATTGCGTTTATCCTGGCTGTTGTTATGGCTATAGTCCTTTCTGTTGGGCTATTCCg                                                                                   >  1:125655/1‑150 (MQ=255)
                                                                      ccGGACCACGCCAGGCTGTGAGCGCGTTTGATCTGGCCTTTTTGAGCAATCGTGGGCTTTATATTCCGTCACCGCAGCTGGGAGACGGATTTATTGCGTTTATCCTGGCTGTTGTTATGGCTATAGTCCTTTCTGTTGGGCTATTCCGCt                                                                                 >  1:247002/1‑150 (MQ=255)
                                                                       cGGACCACGCCAGGCTGTGAGCGCGTTTGATCTGGCCTTTTTGAGCAATCGTGGGCTTTATATTCCGTCACCGCAGCTGGGAGACGGATTTATTGCGTTTATCCTGGCTGTTGTTATGGCTATAGTCCTTTCTGTTGGGCTATTCCGCtt                                                                                <  2:139731/150‑1 (MQ=255)
                                                                        ggACCACGCCAGGCTGTGAGCGCGTTTGATCTGGCCTTTTTGAGCAATCGTGGGCTTTATATTCCGTCACCGCAGCTGGGAGACGGATTTATTGCGTTTATCCTGGCTGTTGTTATGGCTATAGTCCTTTCTGTTGGGCTATTCCGCttt                                                                               <  2:909777/150‑1 (MQ=255)
                                                                        ggACCACGCCAGGCTGTGAGCGCGTTTGATCTGGCCTTTTTGAGCAATCGTGGGCTTTATATTCCGTCACCGCAGCTGGGAGACGGATTTATTGCGTTTATCCTGGCTGTTGTTATGGCTATAGTCCTTTCTGTTGGGCTATTCCGCttt                                                                               <  1:381115/150‑1 (MQ=255)
                                                                        ggACCACGCCAGGCTGTGAGCGCGTTTGATCTGGCCTTTTTGAGCAATCGTGGGCTTTATATTCCGTCACCGCAGCTGGGAGACGGATTTATTGCGTTTATCCTGGCTGTTGTTATGGCTATAGTCCTTTCTGTTGGGCTATTCCGCttt                                                                               <  1:753989/150‑1 (MQ=255)
                                                                        ggACCACGCCAGGCTGTGAGCGCGTTTGATCTGGCCTTTTTGAGCAATCGTGGGCTTTATATTCCGTCACCGCAGCTGGGAGACGGATTTATTGCGTTTATCCTGGCTGTTGTTATGGCTATAGTCCTTTCTGTTGGGCTATTCCGCttt                                                                               <  1:473028/150‑1 (MQ=255)
                                                                         gACCACGCCAGGCTGTGAGCGCGTTTCATCTGGCCTTTTTGAGCAATCGTGGGCTTTATATTCCGTCACCGCAGCTGGGAGACGGATTTATTGCGTTTATCCTGGCTGTTGTTATGGCTATAGTCCTTTCTGTTGGGCTATTCCGCTTta                                                                              >  2:1010985/1‑150 (MQ=255)
                                                                             aCGCCAGGCTGTGAGCGCGTTTGATCTGGCCTTTTTGAGCAATCGTGGGCTTTATATTCCGTCACCGCAGCTGGGAGACGGATTTATTGCGTTTATCCTGGCTGTTGTTATGGCTATAGTCCTTTc                                                                                                  >  2:1204652/1‑126 (MQ=255)
                                                                             aCGCCAGGCTGTGAGCGCGTTTGATCTGGCCTTTTTGAGCAATCGTGGGCTTTATATTCCGTCACCGCAGCTGGGAGACGGATTTATTGCGTTTATCCTGGCTGTTGTTATGGCTATAGTCCTTTc                                                                                                  <  1:1204652/126‑1 (MQ=255)
                                                                             aCGCCAGGCTGTGAGCGCGTTTGATCTGGCCTTTTTGAGCAATCGTGGGCTTTATATTCCGTCACCGCAGCTGGGAGACGGATTTATTGCGTTTATCCTGGCTGTTGTTATGGCTATAGTCCTTTCTGTTGGGCTATTCCGCTTtaataa                                                                          >  2:626859/1‑150 (MQ=255)
                                                                                        tgAGCGCGTTTGATCTGGCCTTTTTGAGCAATCGTGGGCTTTATATTCCGTCACCGCAGCTGGGAGACGGATTTATTGCGTTTATCCTGGCTGTTGTTATGGCTATAGTCCTTTCTGTTGGGCTATTCCGCTTTAATAAAACATACCAGa                                                               >  2:1192272/1‑150 (MQ=255)
                                                                                         gAGCGCGTTTGATCTGGCCTTTTTGAGCAATCGTGGGCTTTATATTCCGTCACCGCAGCTGGGAGACGGATTTATTGCGTTTATCCTGGCTGTTGTTATGGCTATAGTCCTTTCTGTTGGGCTATTCCGCTTTAATAAAACATACCAGAt                                                              <  1:767100/150‑1 (MQ=255)
                                                                                          aGCGCGTTTGATCTGGCCTTTTTGAGCAATCGTGGGCTTTATATTCCGTCACCGCAGCTGGGAGACGGATTTATTGCGTTTATCCTGGCTGTTGTTATGGCTATAGTCCTTTCTGTTGGGCTATTCCGCTTTAATAAAACATACCAGATa                                                             >  1:569544/1‑150 (MQ=255)
                                                                                           gcgcGTTTGATCTGGCCTTTTTGAGCAATCGTGGGCTTTATATTCCGTCACCGCAGCTGGGAGACGGATTTATTGCGTTTATCCTGGCTGTTGTTATGGCTATAGTCCTTTCTGTTGGGCTATTCCGCTTTAATAAAACATACCAGATaa                                                            >  2:1220886/1‑150 (MQ=255)
                                                                                             gcgTTTGATCTGGCCTTTTTGAGCAATCGTGGGCTTTATATTCCGTCACCGCAGCTGGGAGACGGATTTATTGCGTTTATCCTGGCTGTTGTTATGGCTATAGTCCTTTCTGTTGGGCTATTCCGCTTTAATAAAACATACCAGATAAAg                                                          <  1:393514/150‑1 (MQ=255)
                                                                                              cgTTTGATCTGGCCTTTTTGAGCAATCGTGGGCTTTATATTCCGTCACCGCAGCTGGGAGACGGATTTATTGCGTTTATCCTGGCTGTTGTTATGGCTATAGTCCTTTCTGTTGGGCTATTCCGCTTTAATAAAACATACCAGATAAAGa                                                         >  2:189457/1‑150 (MQ=255)
                                                                                                      cTGGCCTTTTTGAGCAATCGTGGGCTTTATATTCCGTCACCGCAGCTGGGAGACGGATTTATTGCGTTTATCCTGGCTGTTGTTATGGCTATAGTCCTTTCTGTTGGGCTATTCCGCTTTAATAAAACATACCAGATAAAGACCGGACaa                                                 <  2:1006711/150‑1 (MQ=255)
                                                                                                          ccTTTTTGAGCAATCGTGGGCTTTATATTCCGTCACCGCAGCTGGGAGACGGATTTATTGCGTTTATCCTGGCTGTTGTTATGGCTATAGTCCTTTCTGTTGGGCTATTCCGCTTTAATAAAACATACCAGATAAAGACCGGACAACTgc                                             >  1:95425/1‑150 (MQ=255)
                                                                                                          ccTTTTTGAGCAATCGTGGGCTTTATATTCCGTCACCGCAGCTGGGAGACGGATTTATTGCGTTTATCCTGGCTGTTGTTATGGCTATAGTCCTTTCTGTTGGGCTATTCCGCTTTAATAAAACATACCAGATAAAGACCGGACAACTgc                                             >  1:1064575/1‑150 (MQ=255)
                                                                                                           cTTTTTGAGCAATCGTGGGCTTTATATTCCGTCACCGCAGCTGGGAGACGGATTTATTGCGTTTATCCTGGCTGTTGTTATGGCTATAGTCCTTTCTGTTGGGCTATTCCGCTTTAATAAAACATACCAGATAAAGACCGGACAAc                                                >  2:410512/1‑146 (MQ=255)
                                                                                                           cTTTTTGAGCAATCGTGGGCTTTATATTCCGTCACCGCAGCTGGGAGACGGATTTATTGCGTTTATCCTGGCTGTTGTTATGGCTATAGTCCTTTCTGTTGGGCTATTCCGCTTTAATAAAACATACCAGATAAAGACCGGACAAc                                                <  1:410512/146‑1 (MQ=255)
                                                                                                            tttttGAGCAATCGTGGGCTTTATATTCCGTCACCGCAGCTGGGAGACGGATTTATTGCGTTTATCCTGGCTGTTGTTATGGCTATAGTCCTTTCTGTTGGGCTATTCCGCTTTAATAAAACATACCAGATAAAGACCGGACAACTgcgc                                           >  2:33279/1‑150 (MQ=255)
                                                                                                             ttttGAGCAATCGTGGGCTTTATATTCCGTCACCGCAGCTGGGAGACGGATTTATTGCGTTTATCCTGGCTGTTGTTATGGCTATAGTCCTTTCTGTTGGGCTATTCCGCTTTAATAAAACATACCAGATAAAGACCGGACAACTGcgcc                                          >  2:1441370/1‑150 (MQ=255)
                                                                                                               ttGAGCAATCGTGGGCTTTATATTCCGTCACCGCAGCTGGGAGACGGATTTATTGCGTTTATCCTGGCTGTTGTTATGGCTATAGTCCTTTCTGTTGGGCTATTCCGCTTTAATAAAACATACCAGATAAAGACCGGACAACTGcgccgc                                        >  1:767793/1‑150 (MQ=255)
                                                                                                               ttGAGCAATCGTGGGCTTTATATTCCGTCACCGCAGCTGGGAGACGGATTTATTGCGTTTATCCTGGCTGTTGTTATGGCTATAGTCCTTTCTGTTGGGCTATTCCGCTTTAATAAAACATACCAGATAAAGACCGGACAACTGcgccgc                                        <  1:33279/150‑1 (MQ=255)
                                                                                                                tGAGCAATCGTGGGCTTTATATTCCGTCACCGCAGCTGGGAGACGGATTTATTGCGTTTATCCTGGCTGTTGTTATGGCTATAGTCCTTTCTGTTGGGCTATTCCGCTTTAATAAAACATACCAGATAAAGACCGGACAACTGCGCCGCa                                       <  1:1024402/150‑1 (MQ=255)
                                                                                                                tGAGCAATCGTGGGCTTTATATTCCGTCACCGCAGCTGGGAGACGGATTTATTGCGTTTATCCTGGCTGTTGTTATGGCTATAGTCCTTTCTGTTGGGCTATTCCGCTTTAATAAAACATACCAGATAAAGACCGGACAACTGCGCCGCa                                       <  1:484147/150‑1 (MQ=255)
                                                                                                                      aTCGTGGGCTTTATATTCCGTCACCGCAGCTGGGAGACGGATTTATTGCGTTTATCCTGGCTGTTGTTATGGCTATAGTCCTTTCTGTTGGGCTATTCCGCTTTAATAAAACATACCAGATAAAGACCGGACAACTGCGCCGCACCTGGc                                 >  1:1073736/1‑150 (MQ=255)
                                                                                                                        cGTGGGCTTTATATTCCGTCACCGCAGCTGGGAGACGGATTTATTGCGTTTATCCTGGCTGTTGTTATGGCTATAGTCCTTTCTGTTGGGCTATTCCGCTTTAATAAAACATACCAGATAAAGACCGGACAACTGCGCCGCACCTGGCCg                               <  2:672128/150‑1 (MQ=255)
                                                                                                                           tggCTTTTTATTCCGTCACCGCAGCTGGGAGACGGATTTATTGCGTTTATCCTGGCTGTTGTTATGGCTATAGTCCTTTCTGTTGGGCTATTCCGCTTTAATAAAACATACCAGATAAAGACCGGACAACTGCGCCGCACCTGGCCGATc                            <  1:1122524/149‑1 (MQ=255)
                                                                                                                            ggCTTTATATTCCGTCACCGCAGCTGGGAGACGGATTTATTGCGTTTATCCTGGCTGTTGTTATGGCTATAGTCCTTTCTGTTGGGCTATTCCGCTTTAATAAAACATACCAGATAAAGACCGGACAACTGCGCCGCACCTGGCCGATcg                           <  2:323205/150‑1 (MQ=255)
                                                                                                                            ggCTTTATATTCCGTCACCGCAGCTGGGAGACGGATTTATTGCGTTTATCCTGGCTGTTGTTATGGCTATAGTCCTTTCTGTTGGGCTATTCCGCTTTAATAAAACATACCAGATAAAGACCGGACAACTGCGCCGCACCTGGCCGATcg                           <  2:1329899/150‑1 (MQ=255)
                                                                                                                            ggCTTTATATTCCGTCACCGCAGCTGGGAGACGGATTTATTGCGTTTATCCTGGCTGTTGTTATGGCTATAGTCCTTTCTGTTGGGCTATTCCGCTTTAATAAAACATACCAGATAAAGACCGGACAACTGCGCCGCACCTGGCCGATcg                           <  1:1366162/150‑1 (MQ=255)
                                                                                                                                 tataTTCCGTCACCGCAGCTGGGAGACGGATTTATTGCGTTTATCCTGGCTGTTGTTATGGCTATAGTCCTTTCTGTTGGGCTAt                                                                                       >  2:1489645/1‑85 (MQ=255)
                                                                                                                                 tataTTCCGTCACCGCAGCTGGGAGACGGATTTATTGCGTTTATCCTGGCTGTTGTTATGGCTATAGTCCTTTCTGTTGGGCTAt                                                                                       <  1:1489645/85‑1 (MQ=255)
                                                                                                                                  atatTCCGTCACCGCAGCTGGGAGACGGATTTATTGCGTTTATCCTGGCTGTTGTTATGGCTATAGTCCTTTCTGTTGGGCTATTCCGCTTTAATAAAACATACCAGATAAAGACCGGACAACTGCGCCGCACATGGCCGATCGCCGCAg                     >  2:849287/1‑150 (MQ=255)
                                                                                                                                    atTCCGTCACCGCAGCTGGGAGACGGATTTATTGCGTTTATCCTGGCTGTTGTTATGGCTATAGTCCTTTCTGTTGGGCTATTCCGCTTTAATAAAACATACCAGATAAAGACCGGACAACTGCGCCGCACCTGGCCGATCGCCGCAgtg                   >  2:1025811/1‑150 (MQ=255)
                                                                                                                                    aatCCGTCTCCCAGCTGCGGTGACGGATTTATTGCGTTTATCCTGGCTGTTGTTATGGCTATAGTCCTTTCTGTTGGGCTATTCCGCTTTAATAAAACATACCAGATAAAGACCGGACAACTGCGCCGCACCTGGCCGATCGCCGCAgtg                   >  2:1360048/3‑150 (MQ=255)
                                                                                                                                     ttCCGTCACCGCAGCTGGGAGACGGATTTATTGCGTTTATCCTGGCTGTTGTTATGGCTATAGTCCTTTCTGTTGGGCTATTCCGCTTTAATAAAACATACCAGATAAAGACCGGACAACTGCGCCGCACCTGGCCGATCGCCGCAgtgt                  >  1:856560/1‑150 (MQ=255)
                                                                                                                                     atCCGTCTCCCAGCTGCGGTGACGGATTTATTGCGTTTATCCTGGCTGTTGTTATGGCTATAGTCCTTTCTGTTGGGCTATTCCGCTTTAATAAAACATACCAGATAAAGACCGGACAACTGCGCCGCACCTGGCCGATCGCCGCAgtgt                  >  2:45030/2‑150 (MQ=255)
                                                                                                                                       ccGTCACCGCAGCTGGGAGACGGATTTATTGCGTTTATCCTGGCTGTTGTTATGGCTATAGTCCTTTCTGTTGGGCTATTCCGCTTTAATAAAACATACCAGATAAAGACCGGACAACTGCGCCGCACCTGGCCGATCGCCGCAGTGTTg                <  2:620792/150‑1 (MQ=255)
                                                                                                                                         gTCACCGCAGCTGGGAGACGGATTTATTGCGTTTATCCTGGCTGTTGTTATGGCTATAGTCCTTTCTGTTGGGCTATTCCGCTTTAATAAAACATACCAGATAAAGACCGGACAACTGCGCCGCACCTGGCCGATCGCCGCAGTGTTGAt              <  1:629813/150‑1 (MQ=255)
                                                                                                                                         gTCACCGCAGCTGGGAGACGGATTTATTGCGTTTATCCCGGCTGTTGTTATGGCTATAGTCCTTTCTGTTGGGCTATTCCGCTTTAATAAAACATACCAGATAAAGACCGGACAACTGCGCCGCACCTGGCCGATCGCCGCAGTGTTGAt              >  1:417890/1‑150 (MQ=255)
                                                                                                                                             ccGCAGCTGGGAGACGGATTTATTGCGTTTATCCTGGCTGTTGTTATGGCTATAGTCCTTTCTGTTGGGCTATTCCGCTTTAATAAAACATACCAGATAAAGACCGGACAACTGCGCCGCACCTGGCCGATCGCCGCAGTGTTGATCAtt          >  1:965044/1‑150 (MQ=255)
                                                                                                                                               gCAGCTGGGAGACGGATTTATTGCGTTTATCCTGGCTGTTGTTATGGCTATAGTCCTTTCTGTTGGGCTATTCCGCTTTAATAAAACATACCAGATAAAGACCGGACAACTGCGCCGCACCTGGCCGATCACCGCAGTGTTGATCATTgg        <  1:265585/150‑1 (MQ=255)
                                                                                                                                                 aGCGGGCAGACGGATTTATAGCGTGTATCCTGGCTGTGGTTATGGCTCTAGTCCTTTCTGTGGGGCGATTCCGCTTTAATAAAACCTACCACATAAAGACCGGCCAACCGCGCCGCACCTGGCCGAGCGCCGCAGTGTTGATCATTGGtt      <  2:329303/150‑1 (MQ=255)
                                                                                                                                                    tGGGAGACGGATTTATTGCGTTTATCCTGGCTGTTGTTATGGCTATAGTCCTTTCTGTTGGGCTATTCCGCTTTAATAAAACATACCAGATAAAGACCGGACAACTGCGCCGCACCTGGCCGATCGCCGCAGTGTTGATCATTGGTTTGc   >  2:466749/1‑150 (MQ=255)
                                                                                                                                                    tGGGAGACGGATTTATTGCGTTTATCCTGGCTGTTGTTATGGCTATAGTCATTTCTGTTGGGCTATTCCGCTTTAATAAAAAATACCAGATAAAGACCGGACAACTGCGCCGCACCTGGCCCATCGCCGCAGTGTTGATCATTGGTTTGc   >  2:78858/1‑150 (MQ=255)
                                                                                                                                                     gggAGACGGATTTATTGCGTTTATCCTGGCTGTTGTTATGGCTATAGTCCTTTCTGTTGGGCTATTCCGCTTTAATAAAACATACCAGATAAAGACCGGACAACTGCGCCGCACCTGGCCGATCGCCGCAGTGTTGATCATTGGTTTGcc  <  1:294805/150‑1 (MQ=255)
                                                                                                                                                     gggAGACGGATTTATTGCGTTTATCCTGGCTGTTGTTATGGCTATAGTCCTTTCTGTTGGGCTATTCCGCTTTAATAAAACATACCAGATAAAGACCGGACAACTGCGCCGCACCTGGCCGATCGCCGCAGTGTTGATCATTGGTTTGcc  >  2:533305/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
GAGATCTTCCGTAATATTCCCCCGCTGCTGCAAATCTTCTTCTGGTACTTTGCCGTGTTGCGCAATTTGCCCGGACCACGCCAGGCTGTGAGCGCGTTTGATCTGGCCTTTTTGAGCAATCGTGGGCTTTATATTCCGTCACCGCAGCTGGGAGACGGATTTATTGCGTTTATCCTGGCTGTTGTTATGGCTATAGTCCTTTCTGTTGGGCTATTCCGCTTTAATAAAACATACCAGATAAAGACCGGACAACTGCGCCGCACCTGGCCGATCGCCGCAGTGTTGATCATTGGTTTGCC  >  NC_000913/3420518‑3420816

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: