Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,557,197 G→T 8.2% P281Q (CCG→CAG)  ddpF ← putative D,D‑dipeptide ABC transporter ATP‑binding subunit DdpF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,557,1970GT 8.2% 272.8 / 10.9 85P281Q (CCG→CAG) ddpFputative D,D‑dipeptide ABC transporter ATP‑binding subunit DdpF
Reads supporting (aligned to +/- strand):  ref base G (43/35);  new base T (2/5);  total (45/40)
Fisher's exact test for biased strand distribution p-value = 2.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.05e-01

AATCGTAAAGCATTGCCGGATGACGCGTCAGGCGCGTGAATGCCTGATGCGTTGTTAGCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCC  >  NC_000913/1557049‑1557343
                                                                                                                                                    |                                                                                                                                                  
aaTCGTAAAGCATTGCCGTATGACGCGTCAGGCGCGTGAATGCCTGATGCGTTGTTAGCATCTCTACAGCGCCCGCCAGCACTGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCgg                                                                                                                                                   <  2:286571/150‑1 (MQ=255)
   cGTAAGGCATTGCCGGATGACGCGTCAGGCGCGTGAATGCCTGATGCGTTGTTAGCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCa                                                                                                                                                <  1:56247/150‑1 (MQ=255)
     tAAAGCATTGCCGGATGACGCGTCAGGCGCGTGAATGCCTGATGCGTTGTTAGCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCCGCGGGCAAc                                                                                                                                              >  1:23895/1‑150 (MQ=255)
        aGCATTGCCGGAGGACGCGGCAGGCGCGTGAATGCCTGATGCGTTGTTAGCATCTCTACAGGGCCCGCCAGCACCGGAGCTCTCGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCCGCTGGCAAccca                                                                                                                                           <  2:1725332/150‑4 (MQ=255)
         gCATTGCCGGATGACGCGTCAGGCGCGTGAATGCCTGATGCGTTGTTAGCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCCGCTGGCAAcccac                                                                                                                                          <  1:483744/150‑5 (MQ=255)
         gCATTGCCGGATGACGCGTCAGGCGCGTGAATGCCTGATGCGTTGTTAGCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCCGCTGGCAAcccac                                                                                                                                          <  1:308632/150‑5 (MQ=255)
          cATTGCCGGATGACGCGTCAGGCGCGTGAATGCCTGATGCGTTGTTAGCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCCGCTGGCAAcccacg                                                                                                                                         <  2:1346019/150‑6 (MQ=255)
          cATTGCCGGATGACGCGTCAGGCGCGTGAATGCCTGATGCGTTGTTAGCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCCGCTGGCAAcccacg                                                                                                                                         <  1:1624922/150‑6 (MQ=255)
          cATTGCCGGATGACGCGTCAGGCGCGTGAATGCCTGATGCGTTGTTAGCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAAcccacg                                                                                                                                         <  1:29717/150‑6 (MQ=255)
          cATTGCCGGATGACGCGTCAGGCGCGTGAATGCCTGATGCGTTGTTAGCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCg                                                                                                                                         <  2:1416120/150‑1 (MQ=255)
          cATTGCCGGATGACGCGTCAGGCGCGTGAATGCCTGATGCGTTGTTAGCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCg                                                                                                                                         <  1:1199014/150‑1 (MQ=255)
              gCCGGATGACGCGTCAGGCGCGTGAATGCCTGATGCGTTGTTAGCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGa                                                                                                                                     >  1:1119770/1‑150 (MQ=255)
              gCCGGATGACGCGTCAGGCGCGTGAATGCCTGATGCGTTGTTAGCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGa                                                                                                                                     >  2:1739841/1‑150 (MQ=255)
                cGGATGACGCGTCAGGCGCGTGAATGCCTGATGCGTTGTTAGCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGaaa                                                                                                                                   <  1:944947/150‑1 (MQ=255)
                 ggATGACGCGTCAGGCGCGTGAATGCCTGATGCGTTGTTAGCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGaaaa                                                                                                                                  >  1:70784/1‑150 (MQ=255)
                      aCGCGTCAGGCGCGTGAATGCCTGATGCGTTGTTAGCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAg                                                                                                                             >  1:189778/1‑150 (MQ=255)
                      aCGCGTCAGGCGCGTGAATGCCTGATGCGTTGTTAGCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAg                                                                                                                             >  1:1212147/1‑150 (MQ=255)
                      aCGCGTCAGGCGCGTGAATGCCTGATGCGTTGTTAGCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAg                                                                                                                             >  1:254782/1‑150 (MQ=255)
                         cgTCAGGCGCGTGAATGCCTGATGCGTTGTTAGCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCt                                                                                                                          >  2:728204/1‑150 (MQ=255)
                             aGGCGCGTGAATGCCTGATGCGTTGTTAGCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCg                                                                                                                      <  2:1489916/150‑1 (MQ=255)
                               gcgcGTGAATGCCTGATGCGTTGTTAGCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGc                                                                                                                    <  2:647032/150‑1 (MQ=255)
                                     gAATGCCTGATGCGTTGTTAGCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTg                                                                                                              >  2:1999090/1‑150 (MQ=255)
                                           cTGATGCGTTGTTAGCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGttt                                                                                                        >  1:852728/1‑150 (MQ=255)
                                           cTGATGCGTTGTTAGCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGttt                                                                                                        >  2:171715/1‑150 (MQ=255)
                                             gATGCGTTGTTAGCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTcc                                                                                                      >  1:1174124/1‑150 (MQ=255)
                                              aTGCGTTGTTAGCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTccc                                                                                                     <  2:1050237/150‑1 (MQ=255)
                                               tGCGTTGTTAGCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCg                                                                                                    <  2:189778/150‑1 (MQ=255)
                                               tGCGTTGTTAGCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCg                                                                                                    >  1:574409/1‑150 (MQ=255)
                                                  gttgttAGCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCa                                                                                                 <  1:54716/150‑1 (MQ=255)
                                                    tgttAGCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGa                                                                                               >  1:20978/1‑150 (MQ=255)
                                                     gttAGCATCTCTACAGCGCCCGACAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGAt                                                                                              >  2:560359/1‑150 (MQ=255)
                                                          cATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGtt                                                                                         <  2:165486/150‑1 (MQ=255)
                                                           aTCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGttt                                                                                        >  1:1117611/1‑150 (MQ=255)
                                                           aTCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGttt                                                                                        >  2:1492555/1‑150 (MQ=255)
                                                              tctACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTAc                                                                                     >  1:1510252/1‑150 (MQ=255)
                                                                 aCAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTa                                                                                  >  1:396599/1‑150 (MQ=255)
                                                                      gcccgccAGCACCGGAGCTAACGTCCGTCCTCCCTTATTGCTAATGATTGACGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTAGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCATATCGCGTTTAAGTAATGcc                                                                             >  2:322933/1‑150 (MQ=255)
                                                                       cccgccAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGccc                                                                            >  2:1681971/1‑150 (MQ=255)
                                                                       cccgccAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACAGCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGccc                                                                            >  1:1747576/1‑150 (MQ=255)
                                                                           ccAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCAttc                                                                        <  2:574409/150‑1 (MQ=255)
                                                                            cAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCAttct                                                                       <  1:1739841/150‑1 (MQ=255)
                                                                              gCACCGGAGCTCACGTCCGTCCTCCCTTATTGGTAATGATTGCCGGAGTTCAGATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAAGGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCAttcttc                                                                     <  1:798193/150‑1 (MQ=255)
                                                                               cACCGGAGCTCACGTCTGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTtcc                                                                    <  2:845032/150‑1 (MQ=255)
                                                                               cACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTtcc                                                                    <  1:1742549/150‑1 (MQ=255)
                                                                                         cACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGtt                                                          >  2:172931/1‑150 (MQ=255)
                                                                                           cgtccgtcCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTa                                                        <  1:1306609/150‑1 (MQ=255)
                                                                                           cgtccgtcCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTa                                                        >  2:1818402/1‑150 (MQ=255)
                                                                                           cgtacgtcATCCCTAATTGCAAATTATTGAAGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTa                                                        >  2:792977/1‑150 (MQ=255)
                                                                                                gtccTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAAt                                                   <  2:70784/150‑1 (MQ=255)
                                                                                                  cctccCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATgg                                                 <  1:1966413/150‑1 (MQ=255)
                                                                                                    tccCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCg                                               <  2:552897/150‑1 (MQ=255)
                                                                                                    tccCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCg                                               >  2:1361243/1‑150 (MQ=255)
                                                                                                     cccTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCgg                                              >  2:146556/1‑150 (MQ=255)
                                                                                                          aTTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGgggagg                                         >  1:815210/1‑150 (MQ=255)
                                                                                                             gCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGaa                                      >  1:308587/1‑150 (MQ=255)
                                                                                                             gCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGaa                                      >  2:150955/1‑150 (MQ=255)
                                                                                                             gCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATACGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGaa                                      >  1:1748153/1‑150 (MQ=255)
                                                                                                               tAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATc                                    <  2:1119770/150‑1 (MQ=255)
                                                                                                               tAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATc                                    <  2:1117611/150‑1 (MQ=255)
                                                                                                                 aTGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCa                                  <  1:1917904/150‑1 (MQ=255)
                                                                                                                   gATTGCCGGACTTCACATCCTTGGGTTGCCAGCGGGCAACGTTCCTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGACCATTCTTCATCCAGCGTTTTATAAATGGCGGGGAGGGAATCCAGc                                >  2:862013/1‑150 (MQ=255)
                                                                                                                   gATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGccc                                                                            <  1:1276267/106‑1 (MQ=255)
                                                                                                                   gATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGccc                                                                            >  2:1276267/1‑106 (MQ=255)
                                                                                                                       gCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCaata                            >  2:647007/1‑150 (MQ=255)
                                                                                                                        ccGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCgggg                                            >  2:624528/1‑133 (MQ=255)
                                                                                                                        ccGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCgggg                                            <  1:624528/133‑1 (MQ=255)
                                                                                                                             cTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCggg                      <  2:1174124/150‑1 (MQ=255)
                                                                                                                             cTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCAGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCggg                      <  2:1702372/150‑1 (MQ=255)
                                                                                                                              ttCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGgggaggg                                        <  1:708978/131‑1 (MQ=255)
                                                                                                                              ttCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGgggaggg                                        >  2:708978/1‑131 (MQ=255)
                                                                                                                               tCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGgtg                    <  2:396599/150‑1 (MQ=255)
                                                                                                                               tCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGgtg                    <  1:1492555/150‑1 (MQ=255)
                                                                                                                               tCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGgtg                    >  2:1320285/1‑150 (MQ=255)
                                                                                                                                  caTCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTAt                 <  1:65460/150‑1 (MQ=255)
                                                                                                                                  caTCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTAt                 >  1:342681/1‑150 (MQ=255)
                                                                                                                                    tCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGgtgt                   >  1:1278264/1‑146 (MQ=255)
                                                                                                                                    tCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGgtgt                   <  2:1278264/146‑1 (MQ=255)
                                                                                                                                    tCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATgg               >  2:119429/1‑150 (MQ=255)
                                                                                                                                     ccGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGa              >  2:1052127/1‑150 (MQ=255)
                                                                                                                                      cgaacgTTGCCCGCTGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGAt             >  2:1643468/6‑150 (MQ=255)
                                                                                                                                      cgaacgTTGCCCGCTGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGAt             >  1:1590098/6‑150 (MQ=255)
                                                                                                                                           gTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCa        >  2:1633714/1‑150 (MQ=255)
                                                                                                                                           gTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCa        <  1:792977/150‑1 (MQ=255)
                                                                                                                                             tGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAgg      <  2:1791028/150‑1 (MQ=255)
                                                                                                                                               ccAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTg    >  1:455337/1‑150 (MQ=255)
                                                                                                                                                 aGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGcc  <  1:647007/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
AATCGTAAAGCATTGCCGGATGACGCGTCAGGCGCGTGAATGCCTGATGCGTTGTTAGCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCC  >  NC_000913/1557049‑1557343

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: