Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,270,338 T→A 32.3% L187Q (CTG→CAG)  yejF → microcin C transporter, ATP‑binding subunit; ABC family

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,270,3380TA32.3% 401.5 / 290.3 272L187Q (CTG→CAG) yejFmicrocin C transporter, ATP‑binding subunit; ABC family
Reads supporting (aligned to +/- strand):  ref base T (83/101);  new base A (42/46);  total (125/147)
Fisher's exact test for biased strand distribution p-value = 6.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 275 total aligned reads displayed.

TTGATCGCGTTGGTATCCGCCAGGCGGCAAAACGGCTGACAGATTATCCGCATCAGCTCTCCGGCGGCGAACGGCAGCGGGTGATGATTGCGATGGCGCTGTTAACGCGACCGGAATTATTAATTGCCGATGAACCGACCACCGCACTGGACGTCTCTGTCCAGGCGCAGATTTTACAGCTGTTGCGCGAACTGCAAGGCGAGCTGAATATGGGCATGCTGTTTATTACTCATAACCTCAGCATTGTCAGAAAACTGGCCCACCGCGTGGCGGTAATGCAAAACGGTCGCTGTGTCG  >  CP009273/2270191‑2270487
                                                                                                                                                   |                                                                                                                                                     
ttGATCGCGTTGGTATCCGCCAGGCGGCAAAACGGCTGACAGATTATCCGCATCAGCTCTCCGGCGGCGAACGGCAGCGGGTGATGATTGCGATGGCGCTGTTAACGCGACCGGAATTATTAATTGCCGATGAACCGACCACCGCACTgg                                                                                                                                                     >  2:3263017/1‑150 (MQ=255)
  gATCGCGTTGGTATCCGCCAGGCGGCAAAACGGCTGACAGATTATCCGCATCAGCTCTCCGGCGGCGAACGGCAGCGGGTGATGATTGCGATGGCGCTGTTAACGCGACCGGAATTATTAATTGCCGATGAACCGACCACCGCACTGGAc                                                                                                                                                   >  2:2325830/1‑150 (MQ=255)
   aTCGCGTTGGTATCCGCCAGGCGGCAAAACGGCTGACAGATTATCCGCATCAGCTCTCCGGCGGCGAACGGCAGCGGGTGATGATTGCGATGGCGCTGTTAACGCGACCGGAATTATTAATTGCCGATGAACCGACCACCGCACTGGACg                                                                                                                                                  <  2:3917263/150‑1 (MQ=255)
        gTTGGTATCCGCCAGGCGGCAAAACGGCTGACAGATTATCCGCATCAGCTCTCCGGCGGCGAACGGCAGCGGGTGATGATTGCGATGGCGCTGTTAACGCGACCGGAATTATTAATTGCCGATGAACCGACCACCGCACTGGACGTctct                                                                                                                                             >  1:1118343/1‑150 (MQ=255)
        gTTGGTATCCGCCAGGCGGCAAAACGGCTGACAGATTATCCGCATCAGCTCTCCGGCGGCGAACGGCAGCGGGTGATGATTGCGATGGCGCTGTTAACGCGACCGGAATTATTAATTGCCGATGAACCGACCACCGCACTGGACGTctct                                                                                                                                             >  2:6432/1‑150 (MQ=255)
         ttGGTATCCGCCAGGCGGCAAAACGGCTGACAGATTATCCGCATCAGCTCTCCGGCGGCGAACGGCAGCGGGTGATGATTGCGATGGCGCTGTTAATGCGACCGGAATTATTAATTGCCGATGAACCGACCACCGCACTGGACGTCTCTg                                                                                                                                            <  1:231891/150‑1 (MQ=255)
         ttGGTATCCGCCAGGCGGCAAAACGGCTGACAGATTATCCGCATCAGCTCTCCGGCGGCGAACGGCAGCGGGTGATGATTGCGATGGCGCTGTTAACGCGACCGGAATTATTAATTGCCGATGAACCGACCACCGCACTGGACGTCTCTg                                                                                                                                            >  1:706497/1‑150 (MQ=255)
            gTATCCGCCAGGCGGCAAAACGGCTGACAGATTATCCGCATCAGCTCTCCGGCGGCGAACGGCAGCGGGTGATGATTGCGATGGCGCTGTTAACGCGACCGGAATTATTAATTGCCGATGAACCGACCACCGCACTGGACGTCTCTGTcc                                                                                                                                         <  2:1514010/150‑1 (MQ=255)
            gTATCCGCCAGGCGGCAAAACGGCTGACAGATTATCCGCATCAGCTCTCCGGCGGCGAACGGCAGCGGGTGATGATTGCGATGGCGCTGTTAACGCGACCGGAATTATTAATTGCCGATGAACCGACCACCGCACAGGACGTCTCTGTcc                                                                                                                                         <  1:3232882/150‑1 (MQ=255)
             tATCCGCCAGGCGGCAAAACGGCTGACAGATTATCCGCATCAGCTCTCCGGCGGCGAACGGCAGCGGGTGATGATTGCGATGGCGCTGTTAACGCGACCGGAATTATTAATTGCCGATGAACCGACCACCGCACAGGACGTCTCTGTCCa                                                                                                                                        <  1:3154210/150‑1 (MQ=255)
                ccGCCAGGCGGCAAAACGGCTGACAGATTATCCGCATCAGCTCTCCGGCGGCGAACGGCAGCGGGTGATGATTGCGATGGCGCTGTTAACGCGACCGGAATTATTAATTGCCGATGAACCGACCACCGCACTGGACGTCTCTGTCCAGgc                                                                                                                                     >  2:3562303/1‑150 (MQ=255)
                ccGCCAGGCGGCAAAACGGCTGACAGATTATCCGCATCAGCTCTCCGGCGGCGAACGGCAGCGGGTGATGATTGCGATGGCGCTGTTAACGCGACCGGAATTATTAATTGCCGATGAACCGACCACCGCACAGGACGTCTCTGTCCAGgc                                                                                                                                     <  1:2016274/150‑1 (MQ=255)
                ccGCCAGGCGGCAAAACGGCTGACAGATTATCCGCATCAGCTCTCCGGCGGCGAACGGCAGCGGGTGATGATTGCGATGGCGCTGTTAACGCGACCGGAATTATTAATTGCCGATGAACCGACCACCGCACAGGACGTCTCTGTCCAGgc                                                                                                                                     <  1:3113981/150‑1 (MQ=255)
                   ccAGGCGGCAAAACGGCTGACAGATTATCCGCATCAGCTCTCCGGCGGCGAACGGCAGCGGGTGATGATTGCGATGGCGCTGTTAACGCGACCGGAATTATTAATTGCCGATGAACCGACCACCGCACTGGACGTCTCTGTCCAGGCGCa                                                                                                                                  <  2:4482217/150‑1 (MQ=255)
                    cAGGCGGCAAAACGGCTGACAGATTATCCGCATCAGCTCTCCGGCGGCGAACGGCAGCGGGTGATGATTGCGATGGCGCTGTTAACGCGACCGGAATTATTAATTGCCGATGAACCGACCACCGCACAGGACGTCTCTGTCCAGGCGCAg                                                                                                                                 >  2:1430402/1‑150 (MQ=255)
                     aGGCGGCAAAACGGCTGACAGATTATCCGCATCAGCTCTCCGGCGGCGAACGGCAGCGGGTGATGATTGCGATGGCGCTGTTAACGCGACCGGAATTATTAATTGCCGATGAACCGACCACCGCACAGGACGTCTCTGTCCAGGCGCAGa                                                                                                                                >  2:4080905/1‑150 (MQ=255)
                         ggcAAAACGGCTGACAGATTATCCGCATCAGCTCTCCGGCGGCGAACGGCAGCGGGTGATGATTGCGATGGCGCTGTTAACGCGACCGGAATTATTAATTGCCGATGAACCGACCACCGCACTGGACGTCTCTGTCCAGGCGCAGAtttt                                                                                                                            <  2:1569645/150‑1 (MQ=255)
                         ggcAAAACGGCTGACAGATTATCCGCATCAGCTCTCCGGCGGCGAACGGCAGCGGGTGATGATTGCGATGGCGCTGTTAACGCGACCGGAATTATTAATTGCCGATGAACCGACCACCGCACTGGACGTCTCTGTCCAGGCGCAGAtttt                                                                                                                            >  2:3588526/1‑150 (MQ=255)
                         ggcAAAACGGCTGACAGATTATCCGCATCAGCTCTCCGGCGGCGAACGGCAGCGGGTGATGATTGCGATGGCGCTGTTAACGCGACCGGAATTATTAATTGCCGATGAACCGACCACCGCACTGGACGTCTCTGTCCAGGCGCAGAtttt                                                                                                                            >  1:2734265/1‑150 (MQ=255)
                          gcAAAACGGCTGACAGATTATCCGCATCAGCTCTCCGGCGGCGAACGGCAGCGGGTGATGATTGCGATGGCGCTGTTAACGCGACCGGAATTATTAATTGCCGATGAACCGACCACCGCACAGGACGTCTCTGTCCAGGCGCAGATTTTa                                                                                                                           >  2:3014830/1‑150 (MQ=255)
                            aaaaCGGCTGACAGATTATCCGCATCAGCTCTCCGGCGGCGAACGGCAGCGGGTGATGATTGCGATGGCGCTGTTAACGCGACCGGAATTATTAATTGCCGATGAACCGACCACCGCACTGGACGTCTCTGTCCAGGCGCAGATTTTACa                                                                                                                         >  2:785235/1‑150 (MQ=255)
                             aaaCGGCTGACAGATTATCCGCATCAGCTCTCCGGCGGCGAACGGCAGCGGGTGATGATTGCGATGGCGCTGTTAACGCGACCGGAATTATTAATTGCCGATGAACCGACCACCGCACTGGACGTCTCTGTCCAGGCGCAGATTTTACAg                                                                                                                        <  2:1096833/150‑1 (MQ=255)
                             aaaCGGCTGACAGATTATCCGCATCAGCTCTCCGGCGGCGAACGGCAGCGGGTGATGATTGCGATGGCGCTGTTAACGCGACCGGAATTATTAATTGCCGATGAACCGACCACCGCACTGGACGTCTCTGTCCAGGCGCAGATTTTACAg                                                                                                                        <  2:1262135/150‑1 (MQ=255)
                                cGGCTGACAGATTATCCGCATCAGCTCTCCGGCGGCGAACGGCAGCGGGTGATGATTGCGATGGCGCTGTTAACGCGACCGGAATTATTAATTGCCGATGAACCGACCACCGCACAGGACGTCTCTGTCCAGGCGCAGATTTTACAGCTg                                                                                                                     <  1:1353556/150‑1 (MQ=255)
                                 ggCTGACAGATTATCCGCATCAGCTCTCCGGCGGCGAACGGCAGCGGGTGATGATTGCGATGGCGCTGTTAACGCGACCGGAATTATTAATTGCCGATGAACCGACCACCGCACTGGACGTCTCTGTCCAGGCGCAGATTTTACAGCTGt                                                                                                                    <  1:4734277/150‑1 (MQ=255)
                                   cTGACAGATTATCCGCATCAGCTCTCCGGCGGCGAACGGCAGCGGGTGATGATTGCGATGGCGCTGTTAACGCGACCGGAATTATTAATTGCCGATGAACCGACCACCGCACTGGACGTCTCTGTCCAGGCGCAGATTTTACAGCTGTTg                                                                                                                  <  1:723274/150‑1 (MQ=255)
                                    tGACAGATTATCCGCATCAGCTCTCCGGCGGCGAACGGCAGCGGGTGATGATTGCGATGGCGCTGTTAACGCGACCGGAATTATTAATTGCCGATGAACCGACCACCGCACTGGACGTCTCTGTCCAGGCGCAGATTTTACAGCTGTTgc                                                                                                                 <  1:3764220/150‑1 (MQ=255)
                                    tGACAGATTATCCGCATCAGCTCTCCGGCGGCGAACGGCAGCGGGTGATGATTGCGATGGCGCTGTTAACGCGACCGGAATTATTAATTGCCGATGAACCGACCACCGCACTGGACGTCTCTGTCCAGGCGCAGATTTTACAGCTGTTgc                                                                                                                 <  2:2542530/150‑1 (MQ=255)
                                     gACAGATTATCCGCATCAGCTCTCCGGCGGCGAACGGCAGCGGGTGATGATTGCGATGGCGCTGTTAACGCGACCGGAATTATTAATTGCCGATGAACCGACCACCGCACAGGACGTCTCTGTCCAGGCGCAGATTTTACAGCTGTTgcg                                                                                                                >  2:3944213/1‑150 (MQ=255)
                                         gATTATCCGATTCTGCTCTCCGGCGGCGAACGGCAGCGGGTGTTGTTTGCGATGGCGCTGTTATCGCGACCTGAATTATTAATTGCCGTTGAACCGACCACCGCTCTGGACGTCTCTGTCCAGGCGCAGATTTTACAGCTGTTGCGCGaa                                                                                                            <  1:943395/150‑1 (MQ=255)
                                          aTTATCCGCATCAGCTCTCCGGCGGCGAACGGCAGCGGGTGATGATTGCGATGGCGCTGTTAACGCGACCGGAATTATTAATTGCCGATGAACCGACCACCGCACTGGACGTCTCTGTCCAGGCGCAGATTTTACAGCTGTTGCGCGAAc                                                                                                           <  1:3285451/150‑1 (MQ=255)
                                            tATCCGCATCAGCTCTCCGGCGGCGAACGGCAGCGGGTGATGATTGCGATGGCGCTGTTAACGCGACCGGAATTATTAATTGCCGATGAACCGACCACCGCACAGGACGTCTCTGTCCAGGCGCAGATTTTACAGCTGTTGCGCGAACTg                                                                                                         >  1:3350002/1‑150 (MQ=255)
                                             aTCCGCATCAGCTCTCCGGCGGCGAACGGCAGCGGGTGATGATTGCGATGGCGCTGTTAACGCGACCGGAATTATTAATTGCCGATGAACCGACCACCGCACTGGACGTCTCTGTCCAGGCGCAGATTTTACAGCTGTTGCGCGAACTGc                                                                                                        <  2:921715/150‑1 (MQ=255)
                                                    tCAGCTCTCCGGCGGCGAACGGCAGCGGGTGATGATTGCGATGGCGCTGTTAACGCGACCGGAATTATTAATTGCCGATGAACCGACCACCGCACTGGACGTCTCTGTCCAGGCGCAGATTTTACAGCTGTTGCGCGAACTGCAAGGCGa                                                                                                 >  1:4377595/1‑150 (MQ=255)
                                                       gCTCTCCGGCGGCGAACGGCAGCGGGTGATGATTGCGATGGCGCTGTTAACGCGACCGGAATTATTAATTGCCGATGAACCGACCACCGCACAGGACGTCTCTGTCCAGGCGCAGATTTTACAGCTGTTGCGCGAACTGCAAGGCGAGCt                                                                                              <  1:441247/150‑1 (MQ=255)
                                                        ctctCCGGCGGCGAACGGCAGCGGGTGATGATTGCGATGGCGCTGTTAACGCGACCGGAATTATTAATTGCCGATGAACCGACCACCGCACTGGACGTCTCTGTCCAGGCGCAGATTTTACAGCTGTTGCGCGAACTGCAAGGCGAGCTg                                                                                             >  1:4105865/1‑150 (MQ=255)
                                                        ctctCCGGCGGCGAACGGCAGCGGGTGATGATTGCGATGGCGCTGTTAACGCGACCGGAATTATTAATTGCCGATGAACCGACCACCGCACAGGACGTCTCTGTCCAGGCGCAGATTTTACAGCTGTTGCGCGAACTGCAAGGCGAGCTg                                                                                             <  1:341406/150‑1 (MQ=255)
                                                         tctcCGGCGGCGAACGGCAGCGGGTGATGATTGCGATGGCGCTGTTAACGCGACCGGAATTATTAATTGCCGATGAACCGACCACCGCACAGGACGTCTCTGTCCAGGCGCAGATTTTACAGCTGTTGCGCGAACTGCAAGGCGAGCTGa                                                                                            >  2:3138196/1‑150 (MQ=255)
                                                         tctcCGGCGGCGAACGGCAGCGGGTGATGATTGCGATGGCGCTGTTAACGCGACCGGAATTATCAATTGCCGATGAACCGACCACCGCACTGGACGTCTCTGTCCAGGCGCAGATTTTACAGCTGTTGCGCGAACTGCAAGGCGAGCTGa                                                                                            <  2:969752/150‑1 (MQ=255)
                                                           tcCGGCGGCGAACGGCAGCGGGTGATGATTGCGATGGCGCTGTTAACGCGACCGGAATTATTAATTGCCGATGAACCGACCACCGCACAGGACGTCTCTGTCCAGGCGCAGATTTTACAGCTGTTGCGCGAACTGCAAGGCGAGCTGAat                                                                                          <  2:2637139/150‑1 (MQ=255)
                                                                cggcgAACGGCAGCGGGTGATGATTGCGATGGCGCTGTTAACGCGACCGGAATTATTAATTGCCGATGAACCGACCACCGCACTGGACGTCTCTGTCCAGGCGCAGATATTACAGCTGTTGCGCGAACTGCAAGGCGAGCTGAATATggg                                                                                     >  1:3685592/1‑150 (MQ=255)
                                                                 ggcgAACGGCAGCGGGTGATGATTGCGATGGCGCTGTTAACGCGACCGGAATTATTAATTGCCGATGAACCGACCACCGCACAGGACGTCTCTGTCCAGGCGCAGATTTTACAGCTGTTGCGCGAACTGCAAGGCGAGCTGAATATGGGc                                                                                    >  1:897453/1‑150 (MQ=255)
                                                                    gAACGGCAGCGGGTGATGATTGCGATGGCGCTGTTAACGCGACCGGAATTATTAATTGCCGATGAACCGACCACCGCACTGGACGTCTCTGTCCAGGCGCAGATTTTACAGCTGTTGCGCGAACTGCAAGGCGAGCTGAATATGGGCATg                                                                                 >  2:2255820/1‑150 (MQ=255)
                                                                       cGGCAGCGGGTGATGATTGCGATGGCGCTGTTAACGCGACCGGAATTATTAATTGCCGATGAACCGACCACCGCACTGGACGTCTCTGTCCAGGCGCAGATTTTACAGCTGTTGCGCGAACTGCAAGGCGAGCTGAATATGGGCATGCTg                                                                              <  1:3465697/150‑1 (MQ=255)
                                                                         gCAGCGGGTGATGATTGCGATGGCGCTGTTAACGCGACCGGAATTATTAATTGCCGATGAACCGACCACCGCACTGGACGTCTCTGTCCAGGCGCAGATTTTACAGCTGTTGCGCGAACTGCAAGGCGAGCTGAATATGGGCATGCTGtt                                                                            <  1:4143360/150‑1 (MQ=255)
                                                                         gCAGCGGGTGATGATTGCGATGGCGCTGTTAACGCGACCGGAATTATTAATTGCCGATGAACCGACCACCGCACTGGACGTCTCTGTCCAGGCGCAGATTTTACAGCTGTTGCGCGAACTGCAAGGCGAGCTGAATATGGGCATGCTGtt                                                                            <  2:1833010/150‑1 (MQ=255)
                                                                          cAGCGGGTGATGATTGCGATGGCGCTGTTAACGCGACCGGAATTATTAATTGCCGATGAACCGACCACCGCACTGGACGTCTCTGTCCAGGCGCAGATTTTACAGCTGTTGCGCGAACTGCAAGGCGAGCTGAATATGGGCATGCTGttt                                                                           <  2:1321364/150‑1 (MQ=255)
                                                                               ggTGATGATTGCGATGGCGCTGTTAACGCGACCGGAATTATTAATTGCCGATGAACCGACCACCGCACTGGACGTCTCTGTCCAGGCGCAGATTTTACAGCTGTTGCGCGAACTGCAAg                                                                                                     >  1:3027618/1‑119 (MQ=255)
                                                                                gTGATGATTGCGATGGCGCTGTTAACGCGACCGGAATTATTAATTGCCGATGAACCGACCACCGCACAGGACGTCTCTGTCCAGGCGCAGATTTTACAGCTGTTGCGCGAACTGCAAGGCGAGCTGAATATGGGCATGCTGTTTATTACt                                                                     <  1:1019226/150‑1 (MQ=255)
                                                                                 tgatgaTTGCGATGGCGCTGTTAACGCGACCGGAATTATTAATTGCCGATGAACCGACCACCGCACTGGACGTCTCTGTCCAGGCGCAGATTTTACAGCTGTTGCGCGAACTGCAAGGCGAGCTGAATATGGGCATGCTGTTTATTACTc                                                                    <  2:4288221/150‑1 (MQ=255)
                                                                                   atgatTGCGATGGCGCTGTTAACGCGACCGGAATTATTAATTGCCGATGAACCGACCACCGCACAGGACGTCTCTGTCCAGGCGCAGATTTTACAGCTGTTGCGCGAACTGCAAGGCGAGCTGAATATGGGCATGCTGTTTATTACTCAt                                                                  <  1:2945041/150‑1 (MQ=255)
                                                                                     gatTGCGATGGCGCTGTTAACGCGACCGGAATTATTAATTGCCGATGAACCGACCACCGCACAGGACGTCTCTGTCCAGGCGCAGATTTTACAGCTGTTGCGCGAACTGCAAGGCGAGCTGAATATGGGCATGCTGTTTATTACTCATaa                                                                >  1:4345555/1‑150 (MQ=255)
                                                                                       ttGCGATGGCGCTGTTAACGCGACCGGAATTATTAATTGCCGATGAACCGACCACCGCACAGGACGTCTCTGTCCAGGCGCAGATTTTACAGCTGTTGCGCGAACTGCAAGGCGAGCTGAATATGGGCATGCTGTTTATTACTCATAAcc                                                              <  1:1263194/150‑1 (MQ=255)
                                                                                          cGATGGCGCTGTTAACGCGACCGGAATTATTAATTGCCGATGAACCGACCACCGCACTGGACGTCTCTGTCCAGGCGCAGATTTTACAGCTGTTGCGCGAACTGCAAGGCGAGCTGAATATGGGCATGCTGTTTATTACTCATAACCTCa                                                           <  2:1552519/150‑1 (MQ=255)
                                                                                           gATGGCGCTGTTAACGCGACCGGAATTATTAATTGCCGATGAACCGACCACCGCACTGGACGTCTCTGTCCAGGCGCAGATTTTACAGCTGTTGCGCGAACTGCAAGGCGAGCTGAATATGGGCATGCTGTTTATTACTCATAACCTCAg                                                          <  2:602880/150‑1 (MQ=255)
                                                                                            aTGGCGCTGTTAACGCGACCGGAATTATTAATTGCCGATGAACCGACCACCGCACTGGACGTCTCTGTCCAGGCGCAGATTTTACAGCTGTTGCGCGAACTGCAAGGCGAGCTGAATATGGGCATGCTGTTTATTACTCATAACCTCAGc                                                         >  2:9838/1‑150 (MQ=255)
                                                                                              ggCGCTGTTAACGCGACCGGAATTATTAATTGCCGATGAACCGACCACCGCACTGGACGTCTCTGTCCAGGCGCAGATTTTACAGCTGTTGCGCGAACTGCAAGGCGAGCTGAATATGGGCATGCTGTTTATTACTCATAACCTCAGCAt                                                       >  2:2180256/1‑150 (MQ=255)
                                                                                               gcgcTGTTAACGCGACCGGAATTATTAATTGCCGATGAACCGACCACCGCACTGGACGTCTCTGTCCAGGCGCAGATTTTACAGCTGTTGCGCGAACTGCAAGGCGAGCTGAATATGGGCATGCTGTTTATTACTCATAACCTCAGCAtt                                                      <  2:426023/150‑1 (MQ=255)
                                                                                                  cTGTTAACGCGACCGGAATTATTAATTGCCGATGAACCGACCACCGCACTGGACGTCTCTGTCCAGGCGCAGATTTTACAGCTGTTGCGCGAACTGCAAGGCGAGCTGAATATGGGCATGCTGTTTATTACTCATAACCTCAGCATTGTc                                                   >  1:1571485/1‑150 (MQ=255)
                                                                                                   tGTTAACGCGACCGGAATTATTAATTGCCGATGAACCGACCACCGCACTGGACGTCTCTGTCCAGGCGCAGATTTTACAGCTGTTGCGCGAACTGCAAGGCGAGCTGAATATGGGCATGCTGTTTATTACTCATAACCTCAGCATTGTCa                                                  <  2:1281051/150‑1 (MQ=255)
                                                                                                   tGTTAACGCGACCGGAATTATTAATTGCCGATGAACCGACCACCGCACTGGACGTCTCTGTCCAGGCGCAGATTTTACAGCTGTTGCGCGAACTGCAAGGCGAGCTGAATATGGGCATGCTGTTTATTACTCATAACCTCAGCATTGTCa                                                  <  2:3708529/150‑1 (MQ=255)
                                                                                                    gTTAACGCGACCGGAATTATTAATTGCCGATGAACCGACCACCGCACTGGACGTCTCTGTCCAGGCGCAGATTTTACAGCTGTTGCGCGAACTGCAAGGCGAGCTGAATATGGGCATGCTGTTTATTACTCATAACCTCAGCATTGTCAg                                                 >  1:2527592/1‑150 (MQ=255)
                                                                                                    gTTAACGCGACCGGAATTATTAATTGCCGATGAACCGACCACCGCACTGGACGTCTCTGTCCAGGCGCAGATTTTACAGCTGTTGCGCGAACTGCAAGGCGAGCTGAATATGGGCATGCTGTTTATTACTCATAACCTCAGCATTGTCAg                                                 <  2:476920/150‑1 (MQ=255)
                                                                                                      tAACGCGACCGGAATTATTAATTGCCGATGAACCGACCACCGCACTGGACGTCTCTGTCCAGGCGCAGATTTTACAGCTGTTGCGCGAACTGCAAGGCGAGCTGAATATGGGCATGCTGTTTATTACTCATAACCTCAGCATTGTCAGaa                                               >  1:548105/1‑150 (MQ=255)
                                                                                                      tAACGCGACCGGAATTATTAATTGCCGATGAACCGACCACCGCACAGGACGTCTCTGTCCAGGCGCAGATTTTACAGCTGTTGCGCGAACTGCAAGGCGAGCTGAATATGGGCATGCTGTTTATTACTCATAACCTCAGCATTGTCAGaa                                               <  2:4204863/150‑1 (MQ=255)
                                                                                                      tAACGCGACCGGAATTATTAATTGCCGATGAACCGACCACCGCACAGGACGTCTCTGTCCAGGCGCAGATTTTACAGCTGTTGCGCGAACTGCAAGGCGAGCTGAATATGGGCATGCTGTTTATTACTCATAACCTCAGCATTGTCAGaa                                               <  2:2866151/150‑1 (MQ=255)
                                                                                                       aaCGCGACCGGAATTATTAATTGCCGATGAACCGACCACCGCACAGGACGTCTCTGTCCAGGCGCAGATTTTACAGCTGTTGCGCGAACTGCAAGGCGAGCTGAATATGGGCATGCTGTTTATTACTCATAACCTCAGCATTGTCAGaaa                                              >  1:3099535/1‑150 (MQ=255)
                                                                                                        aCGCGACCGGAATTATTAATTGCCGATGAACCGACCACCGCACTGGACGTCTCTGTCCAGGCGCAGATTTTACAGCTGTTGCGCGAACTGCAAGGCGAGCTGAATATGGGCATGCTGTTTATTACTCATAACCTCAGCATTGTCAGaaaa                                             >  2:3462550/1‑150 (MQ=255)
                                                                                                        aCGCGACCGGAATTATTAATTGCCGATGAACCGACCACCGCACTGGACGTCTCTGTCCAGGCGCAGATTTTACAGCTGTTGCGCGAACTGCAAGGCGAGCTGAATATGGGCATGCTGTTTATTACTCATAACCTCAGCATTGTCAGaaaa                                             <  1:1761052/150‑1 (MQ=255)
                                                                                                          gcgACCGGAATTATTAATTGCCGATGAACCGACCACCGCACTGGACGTCTCTGTCCAGGCGCAGATTTTACAGCTGTTGCGCGAACTGCAAGGCGAGCTGAATATGGGCATGCTGTTTATTACTCATAACCTCAGCATTGTCAGAAAACt                                           <  2:1730500/150‑1 (MQ=255)
                                                                                                          gcgACCGGAATTATTAATTGCCGATGAACCGACCACCGCACTGGACGTCTCTGTCCAGGCGCAGATTTTACAGCTGTTGCGCGAACTGCAAGGCGAGCTGAATATGGGCATGCTGTTTATTACTCATAACCTCAGCATTGTCAGAAAACt                                           <  1:4074538/150‑1 (MQ=255)
                                                                                                           cgACCGGAATTATTAATTGCCGATGAACCGACCACCGCACAGGACGTCTCTGTCCAGGCGCAGATTTTACAGCTGTTGCGCGAACTGCAAGGCGAGCTGAATATGGGCATGCTGTTTATTACTCATAACCTCAGCATTGTCAGAAAACTg                                          <  1:4495613/150‑1 (MQ=255)
                                                                                                             aCCGGAATTATTAATTGCCGATGAACCGACCACCGCACAGGACGTCTCTGTCCAGGCGCAGATTTTACAGCTGTTGCGCGAACTGCAAGGCGAGCTGAATATGGGCATGCTGTTTATTACTCATAACCTCAGCATTTTCAGAAAACTGGc                                        >  1:3135262/1‑150 (MQ=255)
                                                                                                                ggAATTATTAATTGCCGATGAACCGACCACCGCACTGGACGTCTCTGTCCAGGCGCAGATTTTACAGCTGTTGCGCGAACTGCAAGGCGAGCTGAATATGGGCATGCTGTTTATTACTCATAACCTCAGCATTGTCAGAAAACTGGCCCa                                     >  2:860619/1‑150 (MQ=255)
                                                                                                                ggAATTATTAATTGCCGATGAACCGACCACCGCACTGGACGTCTCTGTCCAGGCGCAGATTTTACAGCTGTTGCGCGAACTGCAAGGCGAGCTGAATATGGGCATGCTGTTTATTACTCATAACCTCAGCATTGTCAGAAAACTGGCCCa                                     >  2:4066220/1‑150 (MQ=255)
                                                                                                                 gAATTATTAATTGCCGATGAACCGACCACCGCACAGGACGTCTCTGTCCAGGCGCAGATTTTACAGCTGTTGCGCGAACTGCAAGGCGAGCTGAATATGGGCATGCTGTTTATTACTCATAACCTCAGCATTGTCAGAAAACTGGCCCAc                                    <  1:3953931/150‑1 (MQ=255)
                                                                                                                  aaTTATTAATTGCCGATGAACCGACCACCGCACAGGACGTCTCTGTCCAGGCGCAGATTTTACAGCTGTTGCGCGAACTGCAAGGCGAGCTGAATATGGGCATGCTGTTTATTACTCATAACCTCAGCATTGTCAGAAAACTGGCCCAcc                                   >  2:512487/1‑150 (MQ=255)
                                                                                                                     tattaATTGCCGATGAACCGACCACCGCACTGGACGTCTCTGTCCAGGCGCAGATTTTACAGCTGTTGCGCGAACTGCAAGGCGAGCTGAATATGGGCATGCTGTTTATTACTCATAACCTCAGCATTGTCAGAAAACTGGCCCACcgcg                                <  1:3843500/150‑1 (MQ=255)
                                                                                                                       ttaATTGCCGATGAACCGACCACCGCACTGGACGTCTCTGTCCAGGCGCAGATTTTACAGCTGTTGCGCGAACTGCAAGGCGAGCTGAATATGGGCATGCTGTTTATTACTCATAACCTCAGCATTGTCAGAAAACTGGCCCACCGCGTg                              >  2:1349447/1‑150 (MQ=255)
                                                                                                                         aaTTGCCGATGAACCGACCACCGCACAGGACGTCTCTGTCCAGGCGCAGATTTTACAGCTGTTGCGCGAACTGCAAGGCGAGCTGAATATGGGCATGCTGTTTATTACTCATAACCTCAGCATTGTCAGAAAACTGGCCCACCGCGTGGc                            >  2:2146857/1‑150 (MQ=255)
                                                                                                                            tGCCGATGAACCGACCACCGCACTGGACGTCTCTGTCCAGGCGCAGATTTTACAGCTGTTGCGCGAACTGCAAGGCGAGCTGAATATGGGCATGCTGTTTATTACTCATAACCTCAGCATTGTCAGAAAACTGGCCCACCGCGTGGCGGt                         >  1:3652032/1‑150 (MQ=255)
                                                                                                                            tGCCGATGAACCGACCACCGCACAGGACGTCTCTGTCCAGGCGCAGATTTTACAGCTGTTGCGCGAACTGCAAGGCGAGCTGAATATGGGCATGCTGTTTATTACTCATAACCTCAGCATTGTCAGAAAACTGGCCCACCGCGTGGCGGt                         <  1:2871381/150‑1 (MQ=255)
                                                                                                                             gCCGATGAACCGACCACCGCACTGGACGTCTCTGTCCAGGCGCAGATTTTACAGCTGTTGCGCGAACTGCAAGGCGAGCTGAATATGGGCATGCGGTTTATTACTCATAACCTCAGCATTGTCAGAAAACTGGCCCACCGCGTGGCGGTa                        <  2:3655311/150‑1 (MQ=255)
                                                                                                                              ccGATGAACCGACCACCGCACTGGACGTCTCTGTCCAGGCGCAGATTTTACAGCTGTTGCGCGAACTGCAAGGCGAGCTGAATATGGGCATGCTGTTTATTACTCATAACCTCAGCATTGTCAGAAAACTGGCCCACCGCGTGGCGGTaa                       >  2:1127048/1‑150 (MQ=255)
                                                                                                                               cGATGAACCGACCACCGCACTGGACGTCTCTGTCCAGGCGCAGATTTTACAGCTGTTGCGCGAACTGCAAGGCGAGCTGAATATGGGCATGCTGTTTATTACTCATAACCTCAGCATTGTCAGAAAACTGGCCCACCGCGTGGCGGTAAt                      <  1:3768213/150‑1 (MQ=255)
                                                                                                                                 aTGAACCGACCACCGCACTGGACGTCTCTGTCCAGGCGCAGATTTTACAGCTGTTGCGCGAACTGCAAGGCGAGCTGAATATGGGCATGCTGTTTATTACTCATAACCTCAGCATTGTCAGAAAACTGGCCCACCGCGTGGCGGTAATGc                    <  1:3081436/150‑1 (MQ=255)
                                                                                                                                  tGAACCGACCACCGCACAGGACGTCTCTGTCCAGGCGCAGAGTTTACAGCTGTTGCGCGAACTGCAAGGCGAGCTGAATATGGGCATGCTGTTTATTACTCATAACCTCAGCATTGTCAGAAAACTGGCCCACCGCGTGGCGGTAATGCa                   >  1:2608499/1‑150 (MQ=255)
                                                                                                                                   gAACCGACCACCGCACTGGACGTCTCTGTCCAGGCGCAGATTTTACAGCTGTTGCGCGAACTGCAAGGCGAGCTGAATATGGGCATGCTGTTTATTACTCATAACCTCAGCATTGTCAGAAAACTGGCCCACCGCGTGGCGGTAATGCaa                  >  1:4597886/1‑150 (MQ=255)
                                                                                                                                   gAACCGACCACCGCACTGGACGTCTCTGTCCAGGCGCAGATTTTACAGCTGTTGCGCGAACTGCAAGGCGAGCTGAATATGGGCATGCTGTTTATTACTCATAACCTCAGCATTGTCAGAAAACTGGCCCACCGCGTGGCGGTAATGCaa                  >  1:4595481/1‑150 (MQ=255)
                                                                                                                                    aaCCGACCACCGCACAGGACGTCTCTGTCCAGGCGCAGATTTTACAGCTGTTGCGCGAACTGCAAGGCGAGCTGAATATGGGCATGCTGTTTATTACTCATAACCTCAGCATTGTCAGAAAACTGGCCCACCGCGTGGCGGTAATGCaaa                 >  2:3388431/1‑150 (MQ=255)
                                                                                                                                     aCCGACCACCGCACTGGACGTCTCTGTCCAGGCGCAGATTTTACAGCTGTTGCGCGAACTGCAAGGCGAGCTGAATATGGGCATGCTGTTTATTACTCATAACCTCAGCATTGTCAGAAAACTGGCCCACCGCGTGGCGGTAATGCaaaa                <  2:1959204/150‑1 (MQ=255)
                                                                                                                                        gACCACCGCACTGGACGTCTCTGTCCAGGCGCAGATTTTACAGCTGTTGCGCGAACTGCAAGGCGAGCTGAATATGGGCATGCTGTTTATTACTCATAACCTCAGCATTGTCAGAAAACTGGCCCACCGCGTGGCGGTAATGCAAAACgg             <  1:4206589/150‑1 (MQ=255)
                                                                                                                                        gACCACCGCACTGGACGTCTCTGTCCAGGCGCAGATTTTACAGCTGTTGCGCGAACTGCAAGGCGAGCTGAATATGGGCATGCTGTTTATTACTCATAACCTCAGCATTGTCAGAAAACTGGCCCACCGCGTGGCGGTAATGCAAAACgg             <  1:855305/150‑1 (MQ=255)
                                                                                                                                        gACCACCGCACTGGACGTCTCTGTCCAGGCGCAGATTTTACAGCTGTTGCGCGAACTGCAAGGCGAGCTGAATATGGGCATGCTGTTTATTACTCATAACCTCAGCATTGTCAGAAAACTGGCCCACCGCGTGGCGGTAATGCAAAACgg             >  1:1591266/1‑150 (MQ=255)
                                                                                                                                        gACCACCGCACAGGACGTCTCTGTCCAGGCGCAGATTTTACAGCTGTTGCGCGAACTGCAAGGCGAGCTGAATATGGGCATGCTGTTTATTACTCATAACCTCAGCATTGTCAGAAAACTGGCCCACCGCGTGGCGGTAATGCAAAACgg             <  2:304150/150‑1 (MQ=255)
                                                                                                                                            accGCACTGGACGTCTCTGTCCAGGCGCAGATTTTACAGCTGTTGCGCGAACTGCAAGGCGAGCTGAATATGGGCATGCTGTTTATTACTCATAACCTCAGCATTGTCAGAAAACTGGCCCACCGCGTGGCGGTAATGCAAAACGGTCGc         >  2:3829512/1‑150 (MQ=255)
                                                                                                                                             ccGCACTGGACGTCTCTGTCCAGGCGCAGATTTTACAGCTGTTGCGCGAACTGCAAGGCGAGCTGAATATGGGCATGCTGTTTATTACTCATAACCTCAGCATTGTCAGAAAACTGGCCCACCGCGTGGCGGTAATGCAAAACGGTCGCt        <  2:4558959/150‑1 (MQ=255)
                                                                                                                                             ccGCACTGGACGTCTCTGTCCAGGCGCAGATTTTACAGCTGTTGCGCGAACTGCAAGGCGAGCTGAATATGGGCATGCTGTTTATTACTCATAACCTCAGCATTGTCAGAAAACTGGCCCACCGCGTGGCGGTAATGCAAAACGGTCGCt        <  2:4105865/150‑1 (MQ=255)
                                                                                                                                              cGCACTGGACGTCTCTGTCCAGGCGCAGATTTTACAGCTGTTGCGCGAACTGCAAGGCGAGCTGAATATGGGCATGCTGTTTATTACTCATAACCTCAGCATTGTCAGAAAACTGGCCCACCGCGTGGCGGTAATGCAAAACGGTCGCtg       >  1:2280243/1‑150 (MQ=255)
                                                                                                                                                cACTGGACGTCTCTGTCCAGGCGCAGATTTTACAGCTGTTGCGCGAACTGCAAGGCGAGCTGAATATGGGCATGCTGTTTATTACTCATAACCTCAGCATTGTCAGAAAACTGGCCCACCGCGTGGCGGTAATGCAAAACGGTCGCtgtg     <  2:2801266/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
TTGATCGCGTTGGTATCCGCCAGGCGGCAAAACGGCTGACAGATTATCCGCATCAGCTCTCCGGCGGCGAACGGCAGCGGGTGATGATTGCGATGGCGCTGTTAACGCGACCGGAATTATTAATTGCCGATGAACCGACCACCGCACTGGACGTCTCTGTCCAGGCGCAGATTTTACAGCTGTTGCGCGAACTGCAAGGCGAGCTGAATATGGGCATGCTGTTTATTACTCATAACCTCAGCATTGTCAGAAAACTGGCCCACCGCGTGGCGGTAATGCAAAACGGTCGCTGTGTCG  >  CP009273/2270191‑2270487

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