Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,557,194 A→C 9.0% L282R (CTG→CGG)  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,1940AC 9.0% 248.5 / 12.0 86L282R (CTG→CGG) ddpFputative D,D‑dipeptide ABC transporter ATP‑binding subunit DdpF
Reads supporting (aligned to +/- strand):  ref base A (43/35);  new base C (3/5);  total (46/40)
Fisher's exact test for biased strand distribution p-value = 4.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTGAAATCGTAAAGCATTGCCGGATGACGCGTCAGGCGCGTGAATGCCTGATGCGTTGTTAGCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCC  >  NC_000913/1557045‑1557343
                                                                                                                                                     |                                                                                                                                                     
ttGAAATCGTAAAGCATTGCCGGATGACGCGTCAGGCGCGTGAATGCCTGATGCGTTGTTAGCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCa                                                                                                                                                       >  2:29717/1‑150 (MQ=255)
 tGAAATCGTAAAGCATTGCCGGATGACGCGTCAGGCGCGTGAATGCCTGATGCGTTGTTAGCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAg                                                                                                                                                      >  1:1702372/1‑150 (MQ=255)
    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                                                                                                        >  2:171715/1‑150 (MQ=255)
                                               cTGATGCGTTGTTAGCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGttt                                                                                                        >  1:852728/1‑150 (MQ=255)
                                                 gATGCGTTGTTAGCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTcc                                                                                                      >  1:1174124/1‑150 (MQ=255)
                                                  aTGCGTTGTTAGCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTccc                                                                                                     <  2:1050237/150‑1 (MQ=255)
                                                   tGCGTTGTTAGCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCg                                                                                                    >  1:574409/1‑150 (MQ=255)
                                                   tGCGTTGTTAGCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCg                                                                                                    <  2:189778/150‑1 (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                                                                                        >  2:1492555/1‑150 (MQ=255)
                                                               aTCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGttt                                                                                        >  1:1117611/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                                        >  2:708978/1‑131 (MQ=255)
                                                                                                                                  ttCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGgggaggg                                        <  1:708978/131‑1 (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                   <  2:1278264/146‑1 (MQ=255)
                                                                                                                                        tCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGgtgt                   >  1:1278264/1‑146 (MQ=255)
                                                                                                                                        tCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATgg               >  2:119429/1‑150 (MQ=255)
                                                                                                                                         ccGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGa              >  2:1052127/1‑150 (MQ=255)
                                                                                                                                          cgaacgTTGCCCGCTGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGAt             >  1:1590098/6‑150 (MQ=255)
                                                                                                                                          cgaacgTTGCCCGCTGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGAt             >  2:1643468/6‑150 (MQ=255)
                                                                                                                                               gTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCa        <  1:792977/150‑1 (MQ=255)
                                                                                                                                               gTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCa        >  2:1633714/1‑150 (MQ=255)
                                                                                                                                                 tGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAgg      <  2:1791028/150‑1 (MQ=255)
                                                                                                                                                   ccAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTg    >  1:455337/1‑150 (MQ=255)
                                                                                                                                                     aGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGcc  <  1:647007/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TTGAAATCGTAAAGCATTGCCGGATGACGCGTCAGGCGCGTGAATGCCTGATGCGTTGTTAGCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCC  >  NC_000913/1557045‑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: