Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 161,307 G→T 6.1% D31Y (GAC→TAC)  mrcB → fused glycosyl transferase and transpeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273161,3070GT 6.1% 808.1 / 17.4 213D31Y (GAC→TAC) mrcBfused glycosyl transferase and transpeptidase
Reads supporting (aligned to +/- strand):  ref base G (124/76);  new base T (9/4);  total (133/80)
Fisher's exact test for biased strand distribution p-value = 7.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 214 total aligned reads displayed.

TGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTACGAAGATGACGATGATTACGACGATTATGATGACTATGAGGATGAAGAACCGATGCCGCGCAAAGGTAAGGGCAAAGGCAAAGGGCGTAAGCCTCGTGGCAAACGCGGCTGGCTATGGCTACTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGA  >  CP009273/161161‑161454
                                                                                                                                                  |                                                                                                                                                   
tGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTACGCCGTCGTCGTTACGAAGATGa                                                                                                                                                    >  2:425003/1‑148 (MQ=255)
       gaagaaCAGAAACTCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCTTGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTACGAAGATGACgatgat                                                                                                                                             <  1:2930509/148‑1 (MQ=255)
        aagaaCAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTACGAAGATGACGATGAtt                                                                                                                                            >  1:1467191/1‑148 (MQ=255)
          gaaCAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTACGAAGATGACGATGATTac                                                                                                                                          >  2:1470097/1‑148 (MQ=255)
          gaaCAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTACGAAGATGACGATGATTac                                                                                                                                          >  2:2168419/1‑148 (MQ=255)
           aaCAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTACGAAGATTACGATGATTacg                                                                                                                                         >  2:1507595/1‑148 (MQ=255)
            aCAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTACGAAGATGACGATGATTacga                                                                                                                                        >  2:1603549/1‑148 (MQ=255)
              cgAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTACGAAGATGACGATGATTacgacg                                                                                                                                      <  1:3552344/147‑1 (MQ=255)
               gAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTACGAAGATGACGATGATTAcgacga                                                                                                                                     <  1:1022782/148‑1 (MQ=255)
               gAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTACGAAGATGACGATGATTAcgacga                                                                                                                                     <  1:210291/148‑1 (MQ=255)
                aaaaTCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTACGAAGATGACGATGATTACGACGAt                                                                                                                                    >  1:43615/1‑148 (MQ=255)
                aaaaTCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTACGAAGATGACGATGATTACGACGAt                                                                                                                                    >  1:1975034/1‑148 (MQ=255)
                aaaaTCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTACGAAGATGACGATGATTACGACGAt                                                                                                                                    >  1:2031576/1‑148 (MQ=255)
                aaaaTCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAAGGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTACGAAGATGACGATGATTACGACGAt                                                                                                                                    <  1:269309/148‑1 (MQ=255)
                    tCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTACGAAGATGACGATGATTACGACGATTatg                                                                                                                                >  2:2829997/1‑148 (MQ=255)
                        gCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTACGAAGATGACGATGATTACGACGATTAtgatga                                                                                                                            <  2:535090/148‑1 (MQ=255)
                             tGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTACGAAGATGACGATGATTACGACGATTATGATGACTATg                                                                                                                       <  2:584540/148‑1 (MQ=255)
                                   tGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTACGAAGATGACGATGATTACGACGATTATGATGACTATGAGGATg                                                                                                                 >  2:573232/1‑148 (MQ=255)
                                    gAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTACGAAGATGACGATGATTACGACGATTATGATGACTATGAGGATga                                                                                                                <  1:3464535/148‑1 (MQ=255)
                                       tatTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTACGAAGATGACGATGATTACGACGATTATGATGACTATGAGGATgaaga                                                                                                             <  2:3276384/148‑1 (MQ=255)
                                       tatTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTACGAAGATGACGATGATTACGACGATTATGATGACTATGAGGATgaaga                                                                                                             >  1:272290/1‑148 (MQ=255)
                                       tatTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTACGAAGATGACGATGATTACGACGATTATGATGACTATGAGGATgaaga                                                                                                             <  1:2341255/148‑1 (MQ=255)
                                        atTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCACAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTACGAAGATGACGATGATTACGACGATTATGATGACTATGAGGATgaagaa                                                                                                            >  2:2112794/1‑148 (MQ=255)
                                        atTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTACGAAGATGACGATGATTACGACGATTATGATGACTATGAGGATgaagaa                                                                                                            <  2:2281486/148‑1 (MQ=255)
                                        atTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTACGAAGATGACGATGATTACGACGATTATGATGACTATGAGGATgaagaa                                                                                                            >  2:1780772/1‑148 (MQ=255)
                                          tGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTACGAAGATGACGATGATTACGACGATTATGATGACTATGAGGATGAAGAAcc                                                                                                          >  2:3497507/1‑148 (MQ=255)
                                           gCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTACGAAGATGACGATGATTACGACGATTATGATGACTATGAGGATGAAGAACCg                                                                                                         <  1:1396164/148‑1 (MQ=255)
                                           gCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTACGAAGATGACGATGATTACGACGATTATGATGACTATGAGGATGAAGAACCg                                                                                                         <  1:1666819/148‑1 (MQ=255)
                                              gagaAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTACGAAGATGACGATGATTACGACGATTATGATGACTATGAGGATGAAGAACCGATg                                                                                                      >  1:2459533/1‑148 (MQ=255)
                                               agaAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTACGAAGATGACGATGATTACGACGATTATGATGACTATGAGGATGAAGAACCGATGc                                                                                                     >  2:2308490/1‑148 (MQ=255)
                                                gaAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTACGAAGATGACGATGATTACGACGATTATGATGACTATGAGGATGAAGAACCGATGcc                                                                                                    >  1:1649168/1‑148 (MQ=255)
                                                   aaaGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAGAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTACGAAGATGACGATGATTACGACGATTATGATGACTATGAGGATGAAGAACCGATGCcgcg                                                                                                 >  1:3649036/1‑148 (MQ=255)
                                                   aaaGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTACGAAGATGACGATGATTACGACGATTATGATGACTATGAGGATGAAGAACCGATGCcgcg                                                                                                 >  1:3649325/1‑148 (MQ=255)
                                                    aaGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTACGAAGATGACGATGATTACGACGATTATGATGACTATGAGGATGAAGAACCGATGCCgcgc                                                                                                >  1:2468246/1‑148 (MQ=255)
                                                     aGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTACGAAGATGACGATGATTACGACGATTATGATGACTATGAGGATGAAGAACCGATGCCGCGCa                                                                                               >  1:3475513/1‑148 (MQ=255)
                                                      gCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTACGAAGATGACGATGATTACGACGATTATGATGACTATGAGGATGAAGAACCGATGCCGCGCaa                                                                                              >  1:275757/1‑148 (MQ=255)
                                                      gCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTACGAAGATGACGATGATTACGACGATTATGATGACTATGAGGATGAAGAACCGATGCCGCGCaa                                                                                              >  1:275890/1‑148 (MQ=255)
                                                         tGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTACGAAGATGACGATGATTACGACGATTATGATGACTATGAGGATGAAGAACCGATGCCGCGCAAAgg                                                                                           <  1:3816774/148‑1 (MQ=255)
                                                          ggCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTACGAAGATGACGATGATTACGACGATTATGATGACTATGAGGATGAAGAACCGATGCCGCGCAAAGGt                                                                                          >  2:447080/1‑148 (MQ=255)
                                                          ggCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTACGAAGATGACGATGATTACGACGATTATGATGACTATGAGGATGAAGAACCGATGCCGCGCAAAGGt                                                                                          >  2:3268821/1‑148 (MQ=255)
                                                           gCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTACGAAGATGACGATGATTACGACGATTATGATGACTATGAGGATGAAGAACCGATGCCGCGCAAAGGTa                                                                                         >  1:3776011/1‑148 (MQ=255)
                                                            cctggAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTACGAAGATGACGATGATTACGACGATTATGATGACTATGAGGATGAAGAACCGATGCCGCGCAAAGGTaa                                                                                        <  2:798323/145‑1 (MQ=255)
                                                               ggAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTACGAAGATGACGATGATTACGACGATTATGATGACTATGAGGATGAAGAACCGATGCCGCGCAAAGGTAAggg                                                                                     >  1:3736798/1‑148 (MQ=255)
                                                                gAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTACGAAGATGACGATGATTACGACGATTATGATGACTATGAGGATGAAGAACCGATGCCGCGCAAAGGTAAGggc                                                                                    >  2:2533005/1‑148 (MQ=255)
                                                                   tGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTACGAAGATGACGATGATTACGACGATTATGATGACTATGAGGATGAAGAACCGATGCCGCGCAAAGGTAAGggcaaa                                                                                 >  1:3290504/1‑148 (MQ=255)
                                                                            gCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTACGAAGATGACGATGATTACGACGATTATGATGACTATGAGGATGAAGAACCGATGCCGCGCAAAGGTAAGGGCAAAGGCAAAggg                                                                        <  2:2614832/148‑1 (MQ=255)
                                                                            gCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTACGAAGATGACGATGATTACGACGATTATGATGACTATGAGGATGAAGAACCGATGCCGCGCAAAGGTAAGGGCAAAGGCAAAggg                                                                        >  2:1473568/1‑148 (MQ=255)
                                                                                aTTGGACGCAAAGGGAAACCGCCGCGTCCGGTCAAACCAAAGGGAAGCCGTCGTCGTTACGAAGATGACGATGATTACGACGATTATGATGACTATGAGGATGAAGAACCGATGCCGCGCAAAGGTAAGGGCAAAGGCAAAGGGCGTa                                                                    <  2:520021/148‑1 (MQ=255)
                                                                                aTTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTACGAAGATGACGATGATTACGACGATTATGATGACTATGAGGATGAAGAACCGATGCCGCGCAAAGGTAAGGGCAAAGGCAAAGGGCGTa                                                                    >  2:2077345/1‑148 (MQ=255)
                                                                                aTTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTACGAAGATGACGATGATTACGACGATTATGATGACTATGAGGATGAAGAACCGATGCCGCGCAAAGGTAAGGGCAAAGGCAAAGGGCGTa                                                                    >  1:1317902/1‑148 (MQ=255)
                                                                                aTTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTACGAAGATGACGATGATTACGACGATTATGATGACTATGAGGATGAAGAACCGATGCCGCGCAAAGGTAAGGGCAAAGGCAAAGGGCGTa                                                                    >  1:1930176/1‑148 (MQ=255)
                                                                                aTTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTACGAAGATGACGATGATTACGACGATTATGATGACTATGAGGATGAAGAACCGATGCCGCGCAAAGGTAAGGGCAAAGGCAAAGGGCGTa                                                                    >  2:1105247/1‑148 (MQ=255)
                                                                                 ttGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTACGAAGATTACGATGATTACGACGATTATGATGACTATGAGGATGAAGAACCGATGCCGCGCAAAGGTAAGGGCAAAGGCAAAGGGCGTaa                                                                   >  1:864194/1‑148 (MQ=255)
                                                                                 ttGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTACGAAGATTACGATGATTACGACGATTATGATGACTATGAGGATGAAGAACCGATGCCGCGCAAAGGTAAGGGCAAAGGCAAAGGGCGTaa                                                                   >  1:2778509/1‑148 (MQ=255)
                                                                                 ttGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTACGAAGATGACGATGATTACGACGATTATGATGACTCTGAGGATGAAGAACCGATGCCGCGCAAAGGTAAGGGCAAAGGCAAAGGGCGTaa                                                                   >  2:384525/1‑148 (MQ=255)
                                                                                 ttGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTACGAAGATGACGATGATTACGACGATTATGATGACTATGAGGATGAAGAACCGATGCCGCGCAAAGGTAAGGGCAAAGGCAAAGGGCGTaa                                                                   >  1:295078/1‑148 (MQ=255)
                                                                                 ttGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTACGAAGATGACGATGATTACGACGATTATGATGACTATGAGGATGAAGAACCGATGCCGCGCAAAGGTAAGGGCAAAGGCAAAGGGCGTaa                                                                   >  1:594340/1‑148 (MQ=255)
                                                                                   ggACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTACGAAGATGACGATGATTACGACGATTATGATGACTATGAGGATGAAGAACCGATGCCGCGCAAAGGTAAGGGCAAAGGCAAAGGGCGTAAGc                                                                 <  1:2546432/148‑1 (MQ=255)
                                                                                   ggACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTACGAAGATGACGATGATTACGACGATTATGATGACTATGAGGATGAAGAACCGATGCCGCGCAAAGGTAAGGGCAAAGGCAAAGGGCGTAAGc                                                                 <  1:3175513/148‑1 (MQ=255)
                                                                                          aaGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTACGAAGATGACGATGATTACGACGATTATGATGACTATGAGGATGAAGAACCGATGCCGCGCAAAGGTAAGGGCAAAGGCAAAGGGCGTAAGCCTCGTgg                                                          >  2:2284993/1‑148 (MQ=255)
                                                                                           aGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTACGAAGATGACGATGATTACGACGATTATGATGACTATGAGGATGAAGAACCGATGCCGCGCAAAGGTAAGGGCAAAGGCAAAGGGCGTAAGCCTCGTGGc                                                         <  2:584731/148‑1 (MQ=255)
                                                                                           aGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTACGAAGATGACGATGATTACGACGATTATGATGACTATGAGGATGAAGAACCGATGCCGCGCAAAGGTAAGGGCAAAGGCAAAGGGCGTAAGCCTCGTGGc                                                         >  2:373166/1‑148 (MQ=255)
                                                                                             ggAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTACGAAGATGACGATGATTACGACGATTATGATGACTATGAGGATGAAGAACCGATGCCGCGCAAAGGTAAGGGCAAAGGCAAAGGGCGTAAGCCTCGTGGCaa                                                       >  2:1643457/1‑148 (MQ=255)
                                                                                                aaCCGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTACGAAGATGACGATGATTACGACGATTATGATGACTATGAGGATGAAGAACCGATGCCGCGCAAAGGTAAGGGCAAAGGCAAAGGGCGTAAGCCTCGTGGCAAAcg                                                    <  1:2168419/148‑1 (MQ=255)
                                                                                                 aCCGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTACGAAGATGACGATGATTACGACGATTATGATGACTATGAGGATGAAGAACCGATGCCGCGCAAAGGTAAGGGCAAAGGCAAAGGGCGTAAGCCTCGTGGCAAAcgc                                                   <  2:2458128/148‑1 (MQ=255)
                                                                                                  ccGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTACGAAGATGACGATGATTACGACGATTATGATGACTATGAGGATGAAGAACCGATGCCGCGCAAAGGTAAGGGCAAAGGCAAAGGGCGTAAGCCTCGTGGCAAAcgcg                                                  >  2:2219733/1‑148 (MQ=255)
                                                                                                   cGACGCGTCCGGTCAAACAAAACGTAAGCCGTCGTCGTTACGAAGATGACGATGATTACGACGATTATGATGACTATGAGGATGAAGAACCGATGCCGCGCAAAAGTAAGGGCAAAGGCAAAGGGCGTAAGCCTCGTGGCAAACGCgg                                                 <  2:856209/148‑1 (MQ=255)
                                                                                                     aCGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTACGAAGATGACGATGATTACGACGATTATGATGACTATGAGGATGAAGAACCGATGCCGCGCAAAGGTAAGGGCAAAGGCAAAGGGCGTAAGCCTCGTGGCAAACGCggct                                               >  1:1119371/1‑148 (MQ=255)
                                                                                                      cgcgTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTACGAAGATGACGATGATTACGACGATTATGATGACTATGAGGATGAAGAACCGATGCCGCGCAAAGGTAAGGGCAAAGGCAAAGGGCGTAAGCCTCGTGGCAAACGCggctg                                              <  2:1709279/148‑1 (MQ=255)
                                                                                                      cgcgTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTACGAAGATGACGATGATTACGACGATTATGATGACTATGAGGATGAAGAACCGATGCCGCGCAAAGGTAAGGGCAAAGGCAAAGGGCGTAAGCCTCGTGGCAAACGCggctg                                              >  2:968778/1‑148 (MQ=255)
                                                                                                           ccGGTCAAACAAAAGGTAAGACGTCGTCGTTACGAAGATGACGATGATTACGACGATTATGATGACTATGAGGATGAAGAACCGATGCCGCGCAAAGGTAAGGGCAAAGGCAAAGGGCGTAAGCCTCGTGGCAAACGCGGCtggctat                                         >  1:1796353/1‑148 (MQ=255)
                                                                                                                     aaaaGGTAAGCCGTCGTCGTTACGAAGATGACGATGATTACGACGATTATGATGACTATGAGGATGAAGAACCGATGCCGCGCAAAGGTAAGGGCAAAGGCAAAGGGCGTAAGCCTCGTGGCAAACGCGGCTGGCTATGGCTACTGCt                               >  1:391343/1‑148 (MQ=255)
                                                                                                                     aaaaGGTAAGCCGTCGTCGTTACGAAGATGACGATGATTACGACGATTATGATGACTATGAGGATGAAGAACCGATGCCGCGCAAAGGTAAGGGCAAAGGCAAAGGGCGTAAGCCTCGTGGCAAACGCGGCTGGCTATGGCTACTGCt                               <  2:2915557/148‑1 (MQ=255)
                                                                                                                      aaaGGTAAGCCGTCGTCGTTACGAAGATGACGATGATTACGACGATTATGATGACTATGAGGATGAAGAACCGATGCCGCGCAAAGGTAAGGGCAAAGGCAAAGGGCGTAAGCCTCGTGGCAAACGCGGCTGGCTATGGCTACTGCTa                              >  2:743384/1‑148 (MQ=255)
                                                                                                                         ggTAAGCCGTCGTCGTTACGAAGATGACGATGATTACGACGATTATGATGACTATGAGGCTGAAGAACCGATGCCGCGCAAAGGTAAGGGCAAAGGCAAAGGGCGTAAGCCTCGTGGCAAACGCGGCTGGCTATGGCTACTGCTaaaa                           <  1:1525177/148‑1 (MQ=255)
                                                                                                                         ggTAAGCCGTCGTCGTTACGAAGATGACGATGATTACGACGATTATGATGACTATGAGGATGAAGAACCGATGCCGCGCAAAGGTAAGGGCAAAGGCAAAGGGCGTAAGCCTCGTGGCAAACGCGGCTGGCTATGGCTACTGCTaaaa                           <  2:3290504/148‑1 (MQ=255)
                                                                                                                         ggTAAGCCGTCGTCGTTACGAAGATGACGATGATTACGACGATTATGATGACTATGAGGATGAAGAACCGATGCCGCGCAAAGGTAAGGGCAAAGGCAAAGGGCGTAAGCCTCGTGGCAAACGCGGCTGGCTATGGCTACTGCTaaaa                           >  1:1591678/1‑148 (MQ=255)
                                                                                                                         ggTAAGCCGTCGTCGTTACGAAGATGACGATGATTACGACGATTATGATGACTATGAGGATGAAGAACCGATGCCGCGCAAAGGTAAGGGCAAAGGCAAAGGGCGTAAGCCTCGTGGCAAACGCGGCTGGCTATGGCTACTGCTaaaa                           <  1:1442280/148‑1 (MQ=255)
                                                                                                                            aaGCCGTCGTCGTTACGAAGATTACGATGATTACGACGATTATGATGACTATGAGGATGAAGAACCGATGCCGCGCAAAGGTAAGGGCAAAGGCAAAGGGCGTAAGCCTCGTGGCAAACGCGGCTGGCTATGGCTACTGCTAAAACTg                        <  1:480638/148‑1 (MQ=255)
                                                                                                                             aGCCGTCGTCGTTACGAAGATGACGATGATTACGACGATTATGATGACTATGAGGATGAAGAACCGATGCCGCGCAAAGGTAAGGGCAAAGGCAAAGGGCGTAAGCCTCGTGGCAAACGCGGCTGGCTATGGCTACTGCTAAAACTgg                       <  1:917085/148‑1 (MQ=255)
                                                                                                                              gCCGTCGTCGTTACGAAGATGACGATGATTACGACGATTATGATGACTATGAGGATGAAGAACCGATGCCGCGCAAAGGTAAGGGCAAAGGCAAAGGGCGTAAGCCTCGTGGCAAACGCGGCTGGCTATGGCTACTGCTAAAACTGGc                      <  2:3541690/148‑1 (MQ=255)
                                                                                                                              gCCGTCGTCGTTACGAAGATGACGATGATTACGACGATTATGATGACTATGAGGATGAAGAACCGATGCCGCGCAAAGGTAAGGGCAAAGGCAAAGGGCGTAAGCCTCGTGGCAAACGCGGCTGGCTATGGCTACTGCTAAAACTGGc                      >  1:913732/1‑148 (MQ=255)
                                                                                                                              gCCGTCGTCGTTACGAAGATGACGATGATTACGACGATTATGATGACTATGAGGATGAAGAACCGATGCCGCGCAAAGGTAAGGGCAAAGGCAAAGGGCGTAAGCCTCGTGGCAAACGCGGCTGGCTATGGCTACTGCTAAAACTGGc                      <  2:1140280/148‑1 (MQ=255)
                                                                                                                                cgtcgtcgtTACGAAGATGACGATGATTACGACGATTATGATGACTATGAGGATGAAGAACCGATGCCGCGCAAAGGTAAGGGCAAAGGCAAAGGGCGTAAGCCTCGTGGCAAACGCGGCTGGCTATGGCTACTGCTAAAACTGGCTa                    >  1:1638847/1‑148 (MQ=255)
                                                                                                                                 gtcgtcgtTACGAAGATGACGATGATTACGACGATTATGATGACTATGAGGATGAAGAACCGATGCCGCGCAAAGGTAAGGGCAAAGGCAAAGGGCGTAAGCCTCGTGGCAAACGCGGCTGGCTATGGCTACTGCTAAAACTGGCTAt                   <  1:1109724/148‑1 (MQ=255)
                                                                                                                                    gtcgtTACGAAGATGACGATGATTACGACGATTATGATGACTATGAGGATGAAGAACCGATGCCGCGCAAAGGTAAGGGCAAAGGCAAAGGGCGTAAGCCTCGTGGCAAACGCGGCTGGCTATGGCTACTGCTAAAACTGGCTATCGt                <  1:2219733/148‑1 (MQ=255)
                                                                                                                                     tcgtTACGAAGATGACGATGATTACGACGATTATGATGACTATGAGGATGAAGAACCGATGCCGCGCAAAGGTAAGGGCAAAGGCAAAGGGCGTAAGCCTCGTGGCAAACGCGGCTGGCTATGGCTACTGCTAAAACTGGCTATCGtt               <  2:3791239/148‑1 (MQ=255)
                                                                                                                                     tcgtTACGAAGATGACGATGATTACGACGATTATGATGACTATGAGGATGAAGAACCGATGCCGCGCAAAGGTAAGGGCAAAGGCAAAGGGCGTAAGCCTCGTGGCAAACGCGGCTGGCTATGGCTACTGCTAAAACTGGCTATCGtt               >  1:37946/1‑148 (MQ=255)
                                                                                                                                     tcgtTACGAAGATGACGATGATTACGACGATTATGATGACTATGAGGATGAAGAACCGATGCCGCGCAAAGGTAAGGGCAAAGGCAAAGGGCGTAAGCCTCGTGGCAAACGCGGCTGGCTATGGCTACTGCTAAAACTGGCTATCGtt               >  1:3263530/1‑148 (MQ=255)
                                                                                                                                      cgtTACGAAGATTACGATGATTACGACGATTATGATGACTATGAGGATGAAGAACCGATGCCGCGCAAAGGTAAGGGCAAAGGCAAAGGGCGTAAGCCTCGTGGCAAACGCGGCTGGCTATGGCTACTGCTAAAACTGGCTATCGttt              <  1:1507537/148‑1 (MQ=255)
                                                                                                                                      cgtTACGAAGATGACGATGATTACGACGATTATGATGACTATGAGGATGAAGAACCGATGCCGCGCAAAGGTAAGGGCAAAGGCAAAGGGCGTAAGCCTCGTGGCAAACGCGGCTGGCTATGGCTACTGCTAAAACTGGCTATCGttt              >  2:324772/1‑148 (MQ=255)
                                                                                                                                       gtTACGAAGATGACGATGATTACGACGATTATGATGACTATGAGGATGAAGAACCGATGCCGCGCAAAGGTAAGGGCAAAGGCAAAGGGCGTAAGCCTCGTGGCAAACGCGGCTGGCTATGGCTACTGCTAAAACTGGCTATCGtttt             <  1:1780772/148‑1 (MQ=255)
                                                                                                                                         tACGAAGATGACGATGATTACGACGATTATGATGACTATGAGGATGAAGAACCGATGCCGCGCAAAGGTAAGGGCAAAGGCAAAGGGCGTAAGCCTCGTGGCAAACGCGGCTGGCTATGGCTACTGCTAAAACTGGCTATCGTTTTTg           >  2:685274/1‑148 (MQ=255)
                                                                                                                                         tACGAAGATGACGATGATTACGACGAGTATGATGACGATGAGGCAGAAGCACCGATGCCGCGCAAAGGTAAGGGCAAAGGCAAAGGGCGTAAGCCTCGTGCCAAACGCGGCTGGCTATGGCTACTCCTAAAACTGGCTATCGTTTTTg           <  2:1836498/148‑1 (MQ=255)
                                                                                                                                              aGATGACGATGATTACGACGATTATGATGACTATGAGGATGAAGAACCGATGCCGCGCAAAGGTAAGGGCAAAGGCAAAGGGCGTAAGCCTCGTGGCAAACGCGGCTGGCTATGGCTACTGCTAAAACTGGCTATCGTTTTTGCCGTg      >  1:1368/1‑148 (MQ=255)
                                                                                                                                               gATGACGATGATTACGACGATTATGATGACTATGATGATGAATAACCGATGCCGCGCAAAGGTAAGGGCAAAGGCAAAGGTCGTAAGCCTCGTGGCAAACTCGGCTGGCTATGGCTACTGCTAAAACTGGCTATCGTTTTTGCCGTGc     >  1:795577/1‑148 (MQ=255)
                                                                                                                                               gATGACGATGATTACGACGATTATGATGACTATGAGGATGAAGAACCGATGCCGCGCAAAGGTAAGGGCAAAGGCAAAGGGCGTAAGCCTCGTGGCAAACGCGGCTGGCTATGGCTACTGCTAAAACTGGCTATCGTTTTTGCCGTGc     >  2:302854/1‑148 (MQ=255)
                                                                                                                                               gATGACGATGATTACGACGATTATGATGACTATGAGGATGAAGAACCGATGCCGCGCAAAGGTAAGGGCAAAGGCAAAGGGCGTAAGCCTCGTGGCAAACGCGGCTGGCTATGGCTACTGCTAAAACTGGCTATCGTTTTTGCCGTGc     <  2:2712951/148‑1 (MQ=255)
                                                                                                                                               gATGACGATGATTACGACGATTATGATGACTATGAGGATGAAGAACAGATGCCGCGCAAAGGTAAGGGCAAAGGCAAATTGCGTAAGCCTCGTGGCAAACGCGGCTGGCTATGGCTACTGCTAAATCTGGCTATCGTTTTTGCAGTGc     >  2:2867258/1‑148 (MQ=255)
                                                                                                                                                aTGACGATGATTACGACGATTATGATGACTATGAGGATGAAGAACCGATGCCGCGCAAAGGTAAGGGCAAAGGCAAAGGGCGTAAGCCTCGTGGCAAACGCGGCTGGCTATGGCTACTGCTAAAACTGGCTATCGTTTTTGCCGTGCt    <  1:3727621/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
TGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTACGAAGATGACGATGATTACGACGATTATGATGACTATGAGGATGAAGAACCGATGCCGCGCAAAGGTAAGGGCAAAGGCAAAGGGCGTAAGCCTCGTGGCAAACGCGGCTGGCTATGGCTACTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGA  >  CP009273/161161‑161454

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: