Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 660,816 G→A 15.2% V318V (GTC→GTT mrdB ← cell wall shape‑determining protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273660,8160GA15.2% 966.8 / 133.0 341V318V (GTC→GTTmrdBcell wall shape‑determining protein
Reads supporting (aligned to +/- strand):  ref base G (116/173);  new base A (25/27);  total (141/200)
Fisher's exact test for biased strand distribution p-value = 2.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 341 total aligned reads displayed.

CGACAACATTTTCCTGTGGGTGTGGATTGACATTACAATCCCGAACCCAGCCATCAGCACAATTAGCGCCGATCCTCCATAACTGACCAGTGGGAGCGGAACCCCTACAACCGGCAGAATACCGCTTACCATACCAATATTTACGAAGACATAAACGAATAATATCAGCATTAAGCCGCCAGCCATGACGCGACCAAAGGTGGTTTGCGCTCTGGCGGCTATCCACAGCCCGCGCATGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCTCTTCCG  >  CP009273/660669‑660964
                                                                                                                                                   |                                                                                                                                                    
 gACAACATTTTCCTGTGGGTGTGGATTGACATTACAATCCCGAACCCAGCCATCAGCACAATTAGCGCCGATCCTCCATAACTGACCAGTGGGAGCGGAACCCCTACAACCGGCAGAATACCGCTTACCATACCAATATTTACGAAGACa                                                                                                                                                   <  2:4183400/150‑1 (MQ=255)
    aacaTTTTCCTGTGGGTGTGGATTGACATTACAATCCCGAACCCAGCCATCAGCACAATTAGCGCCGATCCTCCATAACTGACCAGTGGGAGCGGAACCCCTACAACCGGCAGAATACCGCTTACCATACCAATATTTACGAAGACCTaa                                                                                                                                                >  2:2688794/1‑150 (MQ=255)
     acaTTTTCCTGTGGGTGTGGATTGACATTACAATCCCGAACCCAGCCATCAGCACAATTAGCGCCGATCCTCCATAACTGACCAGTGGGAGCGGAACCCCTACAACCGGCAGAATACCGCTTACCATACCAATATTTACGAAGACATaaa                                                                                                                                               <  2:4738114/150‑1 (MQ=255)
     acaTTTTCCTGTGGGTGTGGATTGACATTACAATCCCGAACCCAGCCATCAGCACAATTAGCGCCGATCCTCCATAACTGACCAGTGGGAGCGGAACCCCTACAACCGGCAGAATACCGCTTACCATACCAATATTTACGAAGACATaaa                                                                                                                                               <  1:5105232/150‑1 (MQ=255)
            ccTGTGGGTGTGGATTGACATTACAATCCCGAACCCAGCCATCAGCACAATTAGCGCCGATCCTCCATAACTGACCAGTGGGAGCGGAACCCCTACAACCGGCAGAATACCGCTTACCATACCAATATTTACGAAGACATAAACGaataa                                                                                                                                        <  2:215446/150‑1 (MQ=255)
            ccTGTGGGTGTGGATTGACATTACAATCCCGAACCCAGCCATCAGCACAATTAGCGCCGATCCTCCATAACTGACCAGTGGGAGCGGAACCCCTACAACCGGCAGAATACCGCTTACCATACCAATATTTACGAAGACATAAACGaataa                                                                                                                                        <  1:1394590/150‑1 (MQ=255)
            ccTGTGGGTGTGGATTGACATTACAATCCCGAACCCAGCCATCAGCACAATTAGCGCCGATCCTCCATAACTGACCAGTGGGAGCGGAACCCCTACAACCGGCAGAATACCGCTTACCATACCAATATTTACGAAAACATAAACGaataa                                                                                                                                        <  1:3225852/150‑1 (MQ=255)
             cTGTGGGTGTGGATTGACATTACAATCCCGAACCCAGCCATCAGCACAATTAGCGCCGATCCTCCATAACTGACCAGTGGGAGCGGAACCCCTACAACCGGCAGAATACCGCTTACCATACCAATATTTACGAAGACATAAACGaataat                                                                                                                                       <  2:669408/150‑1 (MQ=255)
             cTGTGGGTGTGGATTGACATTACAATCCCGAACCCAGCCATCAGCACAATTAGCGCCGATCCTCCATAACTGACCAGTGGGAGCGGAACCCCTACAACCGGCAGAATACCGCTTACCATACCAATATTTACGAAAACATAAACGaataat                                                                                                                                       <  1:4631131/150‑1 (MQ=255)
             cTGTGGGTGTGGATTGACATTACAATCCCGAACCCAGCCATCAGCACAATTAGCGCCGATCCTCCATAACTGACCAGTGGGAGCGGAACCCCTACAACCGGCAGAATACCGCTTACCATACCAATATTTACGAAAACATAAACGaataat                                                                                                                                       <  1:5499792/150‑1 (MQ=255)
               gtgGGTGTGGATTGACATTACAATCCCGAACCCAGCCATCAGCACAATTAGCGCCGATCCTCCATAACTGACCAGTGGGAGCGGAACCCCTACAACCGGCAGAATACCGCTTACCATACCAATATTTACGAAGACATAAACGAATAatat                                                                                                                                     >  1:3296611/1‑150 (MQ=255)
               gtgGGTGTGGATTGACATTACAATCCCGAACCCAGCCATCAGCACAATTAGCGCCGATCCTCCATAACTGACCAGTGGGAGCGGAACCCCTACAACCGGCAGAATACCGCTTACCATACCAATATTTACGAAAACATAAACGAATAatat                                                                                                                                     >  1:1349695/1‑150 (MQ=255)
                tgGGTGTGGATTGACATTACAATCCCGAACCCAGCCATCAGCACAATTAGCGCCGATCCTCCATAACTGACCAGTGGGAGCGGAACCCCTACAACCGGCAGAATACCGCTTACCATACCAATATTTACGAAGACATAAACGAATAATATc                                                                                                                                    >  2:1826733/1‑150 (MQ=255)
                tgGGTGTGGATTGACATTACAATCCCGAACCCAGCCATCAGCACAATTAGCGCCGATCCTCCATAACTGACCAGTGGGAGCGGAACCCCTACAACCGGCAGAATACCGCTTACCATACCAATATTTACGAAGACATAAACGAATAATATc                                                                                                                                    >  2:754161/1‑150 (MQ=255)
                tgGGTGTGGATTGACATTACAATCCCGAACCCAGCCATCAGCACAATTAGCGCCGATCCTCCATAACTGACCAGTGGGAGCGGAACCCCTACAACCGGCAGAATACCGCTTACCATACCAATATTTACGAAGACATAAACGAATAATATc                                                                                                                                    <  1:2548431/150‑1 (MQ=255)
                  ggTGTGGATTGACATTACAATCCCGAACCCAGCCATCAGCACAATTAGCGCCGATCCTCCATAACTGACCAGTGGGAGCGGAACCCCTACAACCGGCAGAATACCGCTTACCATACCAATATTTACGAAGACATAAACGAATAATATCAg                                                                                                                                  <  2:3740042/150‑1 (MQ=255)
                      tgGATTGACATTACAATCCCGAACCCAGCCATCAGCACAATTAGCGCCGATCCTCCATAACTGACCAGTGGGAGCGGAACCCCTACAACCGGCAGAATACCGCTTACCATACCAATATTTACGAAGACATAAACGAATAATATCAGCAtt                                                                                                                              <  1:3834024/150‑1 (MQ=255)
                       ggATTGACATTACAATCCCGAACCCAGCCATCAGCACAATTAGCGCCGATCCTCCATAACTGACCAGTGGGAGCGGAACCCCTACAACCGGCAGAATACCGCTTACCATACCAATATTTACGAAAACATAAACGAATAATATCAGCATTa                                                                                                                             >  1:4106272/1‑150 (MQ=255)
                        gATTGACATTACAATCCCGAACCCAGCCATCAGCACAATTAGCGCCGATCCTCCATAACTGACCAGTGGGAGCGGAACCCCTACAACCGGCAGAATACCGCTTACCATACCAATATTTACGAAGACATAATCGAATAATATCAGCATTaa                                                                                                                            >  1:5613592/1‑150 (MQ=255)
                        gATTGACATTACAATCCCGAACCCAGCCATCAGCACAATTAGCGCCGATCCTCCATAACTGACCAGTGGGAGCGGAACCCCTACAACCGGCAGAATACCGCTTACCATACCAATATTTACGAAGACATAAACGAATAATATCAGCATTaa                                                                                                                            <  2:1494996/150‑1 (MQ=255)
                        gATTGACATTACAATCCCGAACCCAGCCATCAGCACAATTAGCGCCGATCCTCCATAACTGACCAGTGGGAGCGGAACCCCTACAACCGGCAGAATACCGCTTACCATACCAATATTTACGAAGACATAAACGAATAATATCAGCATTaa                                                                                                                            >  2:2394569/1‑150 (MQ=255)
                        gATTGACATTACAATCCCGAACCCAGCCATCAGCACAATTAGCGCCGATCCCCCATAACTGACCAGTGGGAGCGGAACCCCTACAACCGGCAGAATACCGCTTACCATACCAATATTTACGAAGACATAAACGAATAATATCAGCATTaa                                                                                                                            <  1:2031629/150‑1 (MQ=255)
                           tGACATTACAATCCCGAACCCAGCCATCAGCACAATTAGCGCCGATCCTCCATAACTGACCAGTGGGAGCGGAACCCCTACAACCGGCAGAATACCGCTTACCATACCAATATTTACGAAGACATAAACGAATAATATCAGCATTAAgcc                                                                                                                         <  1:1579813/150‑1 (MQ=255)
                           tGACATTACAATCCCGAACCCAGCCATCAGCACAATTAGCGCCGATCCTCCATAACTGACCAGTGGGAGCGGAACCCCTACAACCGGCAGAATACCGCTTACCATACCAATATTTACGAAGACATAAACGAATAATATCAGCATTAAgcc                                                                                                                         >  1:424823/1‑150 (MQ=255)
                            gACATTACAATCCCGAACCCAGCCATCAGCACAATTAGCGCCGATCCTCCATAACTGACCAGTGGGAGCGGAACCCCTACAACCGGCAGAATACCGCTTACCATACCAATATTTACGAAGACATAAACGAATAATATCAGCATTAAgccg                                                                                                                        >  2:2957356/1‑150 (MQ=255)
                             aCATTACAATCCCGAACCCAGCCATCAGCACAATTAGCGCCGATCCTCCATAACTGACCAGTGGGAGCGGAACCCCTACAACCGGCAGAATACCGCTTACCATACCAATATTTACGAAGACATAAACGAATAATATCAGCATTAAgccgc                                                                                                                       >  1:848325/1‑150 (MQ=255)
                                ttACAATCCCGAACCCAGCCATCAGCACAATTAGCGCCGATCCTCCATAACTGACCAGTGGGAGCGGAACCCCTACAACCGGCAGAATACCGCTTACCATACCAATATTTACGAAGACATAAACGAATAATATCAGCATTAAGCCgccag                                                                                                                    >  2:2990566/1‑150 (MQ=255)
                                    aaccccGAACCCACCCACCACCACAAAGACCGCCGAACCACCAACACCGACCAGAGGGACCGGAACCCCTCCCACCCGCACACTACCGCTTACCATACCAATATCTCCGACGACATACACGAATAATATCCGCATTAAGCCGCCAGCCAt                                                                                                                <  2:3837309/147‑1 (MQ=255)
                                      tCCCGAACCCAGCCATCAGCACAATTAGCGCCGATCCTCCATAACTGACCAGTGGGAGCGGAACCCCTACAACCGGCAGAATACCGCTTACCATACCAATATTTACGAAGACATAAACGAATAATATCAGCATTAAGCCGCCAGCCATGa                                                                                                              <  1:4059924/150‑1 (MQ=255)
                                        ccGAACCCAGCCATCAGCACAATTAGCGCCGATCCTCCATAACTGACCAGTGGGAGCGGAACCCCTACAACCGGCAGAATACCGCTTACCATACCAATATTTACGAAAACATAAACGAATAATATCAGCATTAAGCCGCCAGCCATGAcg                                                                                                            <  1:5332788/150‑1 (MQ=255)
                                           aaCCCAGCCATCAGCACAATTAGCGCCGATCCTCCATAACTGACCAGTGGGAGCGGAACCCCTCCAACCGGCAGAATACCGCTTACCATACCTAGATTTACGAAGACATAAACGTATAATATCAGCATTAAGCCGCCAGCCATGACGCGa                                                                                                         <  2:1296188/150‑1 (MQ=255)
                                            aCCCAGCCATCAGCACAATTAGCGCCGATCCTCCATAACTGACCAGTGGGAGCGGAACCCCTACAACCGGCAGAATACCGCTTACCATACCAATATTTACGAAGACATAAACGAATAATATCAGCATTAAGCCGCCAGCCATGACGCGAc                                                                                                        >  2:5098482/1‑150 (MQ=255)
                                            aCCCAGCCATCAGCACAATTAGCGCCGATCCTCCATAACTGACCAGTGGGAGCGGAACCCCTACAACCGGCAGAATACCGCTTACCATACCAATATTTACGAAGACATAAACGAATAATATCAGCATTAAGCCGCCAGCCATGACGCGAc                                                                                                        >  1:132110/1‑150 (MQ=255)
                                                aGCCATCAGCACAATTAGCGCCGATCCTCCATAACTGACCAGTGGGAGCGGAACCCCTACAACCGGCAGAATACCGCTTACCATACCAATATTTACGAAGACATAAACGAATAATATCAGCATTAAGCCGCCAGCCATGACGCGACCaaa                                                                                                    >  2:5624920/1‑150 (MQ=255)
                                                aGCCATCAGCACAATTAGCGCCGATCCTCCATAACTGACCAGTGGGAGCGGAACCCCTACAACCGGCAGAATACCGCTTACCATACCAATATTTACGAAGACATAAACGAATAATATCAGCATTAAGCCGCCAGCCATGACGCGACCaaa                                                                                                    >  1:307511/1‑150 (MQ=255)
                                                  ccATCAGCACAATTAGCGCCGATCCTCCATAACTGACCAGTGGGAGCGGAACCCCTACAACCGGCAGAATACCGCTTACCATACCAATATTTACGAAGACATAAACGAATAATATCAGCATTAAGCCGCCAGCCATGACGCGACCAAAgg                                                                                                  >  2:2726437/1‑150 (MQ=255)
                                                  ccATCAGCACAATTAGCGCCGATCCTCCATAACTGACCAGTGGGAGCGGAACCCCTACAACCGGCAGAATACCGCTTACCATACCAATATTTACGAAGACATAAACGAATAATATCAGCATTAAGCCGCCAGCCATGACGCGACCAAAgg                                                                                                  >  2:628166/1‑150 (MQ=255)
                                                    aTCAGCACTATTAGCGCCGATCCTCCATAACTGACCAGTGGGAGCGGAACCCCTACAACCGGCAGAATACCTCTTACCATACCAATATTTACGAAGACTTAAACGAATATTATCAGCATTAAGACGCCAGCCATGACGCGACCAAAggtg                                                                                                <  2:5451322/150‑1 (MQ=255)
                                                     tCAGCACAATTAGCGCCGATCCTCCATACCTGACCAGTGGGAGCGGAACCCCTCCAACCGGCAGAATACCGCTTACCATACCAATATTTACGAAGACATAAACGAATAATATCAGCATTAAGCCGCCAGCCATGACGCGACCAAAggtgg                                                                                               <  2:1625850/150‑1 (MQ=255)
                                                       aGCACAATTAGCGCCGATCCTCCATAACTGACCAGTGGGAGCGGAACCCCTACAACCGGCAGAATACCGCTTACCATACCAATATTTACGAAGACATAAACGAATAATATCAGCATTAAGCCGCCAGCCATGACGCGACCAAAGGTGGtt                                                                                             >  1:346318/1‑150 (MQ=255)
                                                        gCACAATTAGCGCCGATCCTCCATAACTGACCAGTGGGAGCGGAACCCCTACACCCGGCAGAATACCGCTTACCATACCAATATTTACGAAGACATAAACGAATAATATCAGCATTAAGCCGCCAGCCATGACGCGACCAAAGGTGGttt                                                                                            <  1:4213474/150‑1 (MQ=255)
                                                          acaATTAGCGCCGATCCTCCATAACTGACCAGTGGGAGCGGAACCCCTACAACCGGCAGAATACCGCTTACCATACCAATATTTACGAAGACATAAACGAATAATATCAGCATTAAGCCGCCAGCCATGACGCGACCAAAGGTGGTTTgc                                                                                          <  1:2457744/150‑1 (MQ=255)
                                                          acaATTAGCGCCGATCCTCCATAACTGACCAGTGGGAGCGGAACCCCTACAACCGGCAGAATACCGCTTACCATACCAATATTTACGAAGACATAAACGAATAATATCAGCATTAAGCCGCCAGCCATGACGCGACCAAAGGTGGTTTgc                                                                                          <  1:3154313/150‑1 (MQ=255)
                                                           caATTAGCGCCGATCCTCCATAACTGACCAGTGGGAGCGGAACCCCTACAACCGGCAGAATACCGCTTACCATACCAATATTTACGAAGACATAAACGAATAATATCAGCATTAAGCCGCCAGCCATGACGCGACCAAAGGTGGTTTgcg                                                                                         >  1:5224689/1‑150 (MQ=255)
                                                           caATTAGCGCCGATCCTCAATAACTGACCAGTGGGAGCGGAACCCCTACAACCGGCAGAATACCGCTTACCATACCAATATTTACGAAGACATAAACGAATAATATCAGCATTAAGCCGCCAGCCATGACGCGACCAAAGGTGGTTTgcg                                                                                         <  2:4756413/150‑1 (MQ=255)
                                                                aGCGCCGATCCTCCATAACTGACCAGTGGGAGCGGAACCCCTACAACCGGCAGAATACCGCTTACCATACCAATATTTACGAAGACATAAACGAATAATATCAGCATTAAGCCGCCAGCCATGACGCGACCAAAGGTGGTTTGCGCTCTg                                                                                    <  2:5658352/150‑1 (MQ=255)
                                                                aGCGCCGATCCTCCATAACTGACCAGTGGGAGCGGAACCCCTACAACCGGCAGAATACCGCTTACCATACCAATATTTACGAAGACATAAACGAATAATATCAGCATTAAGCCGCCAGCCATGACGCGACCAAAGGTGGTTTGCGCTCTg                                                                                    <  1:1009272/150‑1 (MQ=255)
                                                                 gcgcCGATCCTCCATAACTGACCAGTGGGAGCGGAACCCCTACAACCGGCAGAATACCGCTTACCATACCAATATTTACGAAGACATAAACGAATAATATCAGCATTAAGCCGCCAGCCATGACGCGACCAAAGGTGGTTTGCGCTCTgg                                                                                   <  1:1105571/150‑1 (MQ=255)
                                                                 gcgcCGATCCTCCATAACTGAACAGTGGGAGCGGAACCCCTTCAACCGGCAGAATACCGCTTACCTTACCAATATTTACGAAGACATAAACGAATAATATCAGCATTAAGCCGCCAGCCATGACGCGACCAAAGGTGGTTTGCGCTCTgg                                                                                   <  2:204731/150‑1 (MQ=255)
                                                                  cgcCGATCCTCCATAACTGACCAGTGGGAGCGGAACCCCTACAACCGGCAGAATACCGCTTACCATACCAATATTTACGAAAACATAAACGAATAATATCAGCATTAAGCCGCCAGCCATGACGCGACCAAAGGTGGTTTGCGCTCTggc                                                                                  >  2:590431/1‑150 (MQ=255)
                                                                   gcCGATCCTCCATAACTGACCAGTGGGAGCGGAACCCCTACAACCGGCAGAATACCGCTTACCATACCAATATTTACGAAGACATAAACGAATAATATCAGCATTAAGCCGCCAGCCATGACGCGACCAAAGGTGGTTTGCGCTCTggcg                                                                                 >  1:2377449/1‑150 (MQ=255)
                                                                      gATCCTCCATAACTGACCAGTGGGAGCGGAACCCCTACAACCGGCAGAATACCGCTTACCATACCAATATTTACGAAAACATAAACGAATAATATCAGCATTAAGCCGCCAGCCATGACGCGACCAAAGGTGGTTTGCGCTCTGGCGGCt                                                                              <  1:5051533/150‑1 (MQ=255)
                                                                       aTCCTCCATAACTGACCAGTGGGAGCGGAACCCCTACAACCGGCAGAATACCGCTTACCATACCAATATTTACGAAGACATAAACGAATAATATCAGCATTAAGCCGCCAGCCATGACGCGACCAAAGGTGGTTTGCGCTCTGGCGGCTa                                                                             <  2:4272740/150‑1 (MQ=255)
                                                                       aTCCTCCATAACTGACCAGTGGGAGCGGAACCCCTACAACCGGCAGAATACCGCTTACCATACCAATATTTACGAAGACATAAACGAATAATATCAGCATTAAGCCGCCAGCCATGACGCGACCAAAGGTGGTTTGCGCTCTGGCGGCTa                                                                             >  1:496952/1‑150 (MQ=255)
                                                                       aTCCTCCATAACTGACCAGTGGGAGCGGAACCCCTACAACCGGCAGAATACCGCTTACCATACCAATATTTACGAAGACATAAACGAATAATATCAGCATTAAGCCGCCAGCCATGACGCGACCAAAGGTGGTTTGCGCTCTGGCGGCTa                                                                             >  2:4673656/1‑150 (MQ=255)
                                                                        tcctccATAACTGACCAGTGGGAGCGGAACCCCTACAACCGGCAGAATACCGCTTACCATACCAATATTTACGAAGACATAAACGAATAATATCAGCATTAAGCCGCCAGCCATGACGCGACCAAAGGTGGTTTGCGCTCTGGCGGCTAt                                                                            >  1:4182935/1‑150 (MQ=255)
                                                                         cctccATAACTGACCAGTGGGAGCGGAACCCCTACAACCGGCAGAATACCGCTTACCATACCAATATTTACGAAGACATAAACGAATAATATCAGCATTAAGCCGCCAGCCATGACGCGACCAAAGGTGGTTTGCGCTCTGGCGGCTATc                                                                           <  1:3587314/150‑1 (MQ=255)
                                                                            ccATAACTGACCAGTGGGAGCGGAACCCCTACAACCGGCAGAATACCGCTTACCATACCAATATTTACGAAGACATAAACGAATAATATCAGCATTAAGCCGCCAGCCATGACGCGACCAAAGGTGGTTTGCGCTCTGGCGGCTATCcac                                                                        <  1:2238076/150‑1 (MQ=255)
                                                                             cATAACTGACCAGTGGGAGCGGAACCCCTACAACCGGCAGAATACCGCTTACCATACCAATATTTACGAAGACATAAACGAATAATATCAGCATTAAGCCGCCAGCCATGACGCGACCAAAGGTGGTTTGCGCTCTGGCGGCTATCcaca                                                                       <  2:2648244/150‑1 (MQ=255)
                                                                              aTAACTGACCAGTGGGAGCGGAACCCCTACAACCGGCAGAATACCGCTTACCATACCAATATTTACGAAGACATAAACGAATAATATCAGCATTAAGCCGCCAGCCATGACGCGACCAAAGGTGGTTTGCGCTCTGGCGGCTATCCACAg                                                                      <  1:3605648/150‑1 (MQ=255)
                                                                              aTAACTGACCAGTGGGAGCGGAACCCCTACAACCGGCAGAATACCGCTTACCATACCAATATTTACGAAGACATAAACGAATAATATCAGCATTAAGCCGCCAGCCATGACGCGACCAAAGGTGGTTTGCGCTCTGGCGGCTATCCACAg                                                                      >  2:1712305/1‑150 (MQ=255)
                                                                               tAACTGACCAGTGGGAGCGGAACCCCTACAACCGGCAGAATACCGCTTACCATACCAATATTTACGAAGACATAAACGAATAATATCAGCATTAAGCCGCCAGCCATGACGCGACCAAAGGTGGTTTGCGCTCTGGCGGCTATCCACAGc                                                                     <  2:433355/150‑1 (MQ=255)
                                                                                   tGACCAGTGGGAGCGGAACCCCTACAACCGGCAGAATACCGCTTACCATACCAATATTTACGAAGACATAAACGAATAATATCAGCATTAAGCCGCCAGCCATGACGCGACCAAAGGTGGTTTGCGCTCTGGCGGCTATCCACAGCCcgc                                                                 >  1:717013/1‑150 (MQ=255)
                                                                                       taGTGGGAGCGGAACCCCTACAACCGGCAGAATACCGCTTACCATACCAATATTTACGAAGACATAAACGAATAATATCAGCATTAAGCCGCCAGCCATGACGCGACCAAAGGTGGTTTGCGCTCTGGCGGCTATCCACAGCCCGCGCAt                                                             >  1:2421698/2‑150 (MQ=255)
                                                                                         gTGGGAGCGGAACCCCTACAACCGGCAGAATACCGCTTACCATACCAATATTTACGAAGACATAAACGAATAATATCAGCATTAAGCCGCCAGCCATGACGCGACCAAAGGTGGTTTGCGCTCTGGCGGCTATCCACAGCCCGCGCATGa                                                           <  1:2974545/150‑1 (MQ=255)
                                                                                          tGGGAGCGGAACCCCTACAACCGGCAGAATACCGCTTACCATACCAATATTTACGAAAACATAAACGAATAATATCAGCATTAAGCCGCCAGCCATGACGCGACCAAAGGTGGTTTGCGCTCTGGCGGCTATCCACAGCCCGCGCATGAt                                                          >  2:314347/1‑150 (MQ=255)
                                                                                            ggAGCGGAACCCCTACAACCGGCAGAATACCGCTTACCATACCAATATTTACGAAGACATAAACGAATAATATCAGCATTAAGCCGCCAGCCATGACGCGACCAAAGGTGGTTTGCGCTCTGGCGGCTATCCACAGCCCGCGCATGATca                                                        <  2:5076593/150‑1 (MQ=255)
                                                                                            ggAGCGGAACCCCTACAACCGGCAGAATACCGCTTACCATACCAATATTTACGAAAACATAAACGAATAATATCAGCATTAAGCCGCCAGCCATGACGCGACCAAAGGTGGTTTGCGCTCTGGCGGCTATCCACAGCCCGCGCATGATca                                                        <  2:4029671/150‑1 (MQ=255)
                                                                                                   aaCCCCTACAACCGGCAGAATACCGCTTACCATACCAATATTTACGAAGACATAAACGAATAATATCAGCATTAAGCCGCCAGCCATGACGCGACCAAAGGTGGTTTGCGCTCTGGCGGCTATCCACAGCCCGCGCATGATCAGCAGAAt                                                 <  1:3533755/150‑1 (MQ=255)
                                                                                                   aaCCCCTACAACCGGCAGAATACCGCTTACCAGACCACTATTTACGAAGACATAAACGAATAATATCAGCATTAAGCCGCCAGCCATGACGCGACCAAAGGTGGTTTGCGCTCTGGCGGCTATCCACAGCCCGCGCATGATCAGCAGAAt                                                 <  2:3550706/150‑1 (MQ=255)
                                                                                                        cTACAACCGGCAGAATACCGCTTACCATACCAATATTTACGAAAACATAAACGAATAATATCAGCATTAAGCCGCCAGCCATGACGCGACCAAAGGTGGTTTGCGCTCTGGCGGCTATCCACAGCCCGCGCATGATCAGCAGAATGTaga                                            >  2:4199777/1‑150 (MQ=255)
                                                                                                           cAACCGGCAGAATACCGCTTACCATACCAATATTTACGAAGACATAAACGAATAATATCAGCATTAAGCCGCCAGCCATGACGCGACCAAAGGTGGTTTGCGCTCTGGCGGCTATCCACAGCCCGCGCATGATCAGCAGAATGTagagag                                         <  2:5235859/150‑1 (MQ=255)
                                                                                                           cAACCGGCAGAATACCGCTTACCATACCAATATTTACGAAGACATAAACGAATAATATCAGCATTAAGCCGCCAGCCATGACGCGACCAAAGGTGGTTTGCGCTCTGGCGGCTATCCACAGCCCGCGCATGATCAGCAGAATGTagagag                                         <  1:1263807/150‑1 (MQ=255)
                                                                                                            aaCCGGCAGAATACCGCTTACCATACCAATATTTACGAAGACATAAACGAATAATATCAGCATTAAGCCGCCAGCCATGACGCGACCAAAGGTGGTTTGCGCTCTGGCGGCTATCCACAGCCCGCGCATGATCAGCAGAATGTAGAgagc                                        <  2:3266974/150‑1 (MQ=255)
                                                                                                              ccGGCAGAATACCGCTTACCATACCAATATTTACGAAGACATAAACGAATAATATCAGCATTAAGCCGCCAGCCATGACGCGACCAAAGGTGGTTTGCGCTCTGGCGGCTATCCACAGCCCGCGCATGATCAGCAGAATGTAGAgagcga                                      >  1:1472638/1‑150 (MQ=255)
                                                                                                              ccGGCAGAATACCGCTTACCATACCAATATTTACGAAAACATAAACGAATAATATCAGCATTAAGCCGCCAGCCATGACGCGACCAAAGGTGGTTTGCGCTCTGGCGGCTATCCACAGCCCGCGCATGATCAGCAGAATGAAGAgagcta                                      >  2:3391052/1‑148 (MQ=255)
                                                                                                               cGGCAGAATACCGCTTACCATACCAATATTTACGAAGACATAAACGAATAATATCAGCATTAAGCCGCCAGCCATGACGCGACCAAAGGTGGTTTGCGCTCTGGCGGCTATCCACAGCCCGCGCATGATCAGCAGAATGTAGAgagcgag                                     <  1:5625284/150‑1 (MQ=255)
                                                                                                                ggCAGAATACCGCTTACCATACCAATATTTACGAAAACATAAACGAATAATATCAGCATTAAGCCGCCAGCCATGACGCGACCAAAGGTGGTTTGCGCTCTGGCGGCTATCCACAGCCCGCGCATGATCAGCAGAATGTAGAgagcgagc                                    >  2:1326164/1‑150 (MQ=255)
                                                                                                                     aaTACCGCTTACCATACCAATATTTACGAAGACATAAACGAATAATATCAGCATTAAGCCGCCAGCCATGACGCGACCAAAGGTGGTTTGCGCTCTGGCGGCTATCCACAGCCCGCGCATGATCAGCAGAATGTAGAGAGCGAGcagaat                               <  2:58408/150‑1 (MQ=255)
                                                                                                                      aTACCGCTTACCATACCAATATTTACGAAGACATAAACGAATAATATCAGCATTAAGCCGCCAGCCATGACGCGACCAAAGGTGGTTTGCGCTCTGGCGGCTATCCACAGCCCGCGCATGATCAGCAGAATGTAGAGAGCGAGcagaatc                              >  2:116217/1‑150 (MQ=255)
                                                                                                                       tACCGCTTACCATACCAATATTTACGAAGACATAAACGAATAATATCAGCATTAAGCCGCCAGCCATGACGCGACCAAAGGTGGTTTGCGCTCTGGCGGCTATCCACAGCCCGCGCATGATCAGCAGAATGTAGAGAGCGAGcagaatca                             <  1:762460/150‑1 (MQ=255)
                                                                                                                        aCCGCTTACCATACCAATATTTACGAAGACATAAACGAATAATATCAGCATTAAGCCGCCAGCCATGACGCGACCAAAGGTGGTTTGCGCTCTGGCGGCTATCCACAGCCCGCGCATGATCAGCAGAATGTAGAGAGCGAGcagaatcag                            >  1:4792693/1‑150 (MQ=255)
                                                                                                                        aCCGCTTACCATACCAATATTTACGAAGACATAAACGAATAATATCAGCATTAAGCCGCCAGCCATGACGCGACCAAAGGTGGTTTGCGCTCTGGCGGCTATCCACAGCCCGCGCATGATCAGCAGAATGTAGAGAGCGAGcagaatcag                            >  1:5284640/1‑150 (MQ=255)
                                                                                                                         ccGCTTACCATACCAATATTTACGAAGACATAAACGAATAATATCAGCATTAAGCCGCCAGCCATGACGCGACCAAAGGTGGTTTGCGCTCTGGCGGCTATCCACAGCCCGCGCATGATCAGCAGAATGTAGAGAGCGAGcagaatcaga                           >  1:2749952/1‑150 (MQ=255)
                                                                                                                           gCTTACCATACCAATATTTACGAAGACATAAACGAATAATATCAGCATTAAGCCGCCAGCCATGACGCGACCAAAGGTGGTTTGCGCTCTGGCGGCTATCCACAGCCCGCGCATGATCAGCAGAATGTAGAGAGCGAGcagaatcagaat                         <  1:2176103/150‑1 (MQ=255)
                                                                                                                            cTTACCATACCAATATTTACGAAGACATAAACGAATAATATCAGCATTAAGCCGCCAGCCATGACGCGACCAAAGGTGGTTTGCGCTCTGGCGGCTATCCACAGCCCGCGCATGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATg                        <  1:635254/150‑1 (MQ=255)
                                                                                                                             ttACCATACCAATATTTACGAAAACATAAACGAATAATATCAGCATTAAGCCGCCAGCCATGACGCGACCAAAGGTGGTTTGCGCTCTGGCGGCTATCCACAGCCCGCGCATGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGc                       <  1:1797674/150‑1 (MQ=255)
                                                                                                                              taccataccaATATTTACGAAGACATAAACGAATAATATCAGCATTAAGCCGCCAGCCATGACGCGACCAAAGGTGGTTTGCGCTCTGGCGGCTATCCACAGCCCGCGCATGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGcc                      <  2:1550191/150‑1 (MQ=255)
                                                                                                                                     ccaATATTTACGAAGACATAAACGAATAATATCAGCATTAAGCCGCCAGCCATGACGCGACCAAAGGTGGTTTGCGCTCTGGCGGCTATCCACAGCCCGCGCATGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAAt               <  2:2245263/150‑1 (MQ=255)
                                                                                                                                     ccaATATTTACGAAGACATAAACGAATAATATCAGCATTAAGCCGCCAGCCATGACGCGACCAAAGGTGGTTTGCGCTCTGGCGGCTATCCACAGCCCGCGCATGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAAt               >  1:4377152/1‑150 (MQ=255)
                                                                                                                                         tatTTACGAAGACATAAACGAATAATATCAGCATTAAGCCGCCAGCCATGACGCGACCAAAGGTGGTTTGCGCTCTGGCGGCTATCCACAGCCCGCGCATGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCa           <  2:1117261/150‑1 (MQ=255)
                                                                                                                                         tatTTACGAAAACATAAACGAATAATATCAGCATTAAGCCGCCAGCCATGACGCGACCAAAGGTGGTTTGCGCTCTGGCGGCTATCCACAGCCCGCGCATGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCa           >  2:2072490/1‑150 (MQ=255)
                                                                                                                                          atTTACGAAGACATAAACGAATAATATCAGCATTAAGCCGCCAGCCATGACGCGACCAAAGGTGGTTTGCGCTCTGGCGGCTATCCACAGCCCGCGCATGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAg          >  1:2466369/1‑150 (MQ=255)
                                                                                                                                          atTTACGAAGACATAAACGAATAATATCAGCATTAAGCCGCCAGCCATGACGCGACCAAAGGTGGTTTGCGCTCTGGCGGCTATCCACAGCCCGCGCATGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAg          >  1:15939/1‑150 (MQ=255)
                                                                                                                                          atTTACGAAGACATAAACGAATAATATCAGCATTAAGCCACCAGCCATGACGCGACCAAAGGTGGTTTGCGCTCTGGCGGCTATCCACAGCCCGCGCATGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAg          <  2:2552439/150‑1 (MQ=255)
                                                                                                                                          atTTACGAAGACATAAACGAATAATATCAGCATTAAGCCACCAGCCATGACGCGACCAAAGGTGGTTTGCGCTCTGGCGGCTATCCACAGCCCGCGCATGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAg          <  2:5366301/150‑1 (MQ=255)
                                                                                                                                           tttACGAAGACATAAACGAATAATATCAGCATTAAGCCGCCAGCCATGACGCGACCAAAGGTGGTTTGCGCTCTGGCGGCTATCCACAGCCCGCGCATGATCAGCAGAATGTAGAgagcgagc                                    >  1:5543724/1‑123 (MQ=255)
                                                                                                                                           tttACGAAGACATAAACGAATAATATCAGCATTAAGCCGCCAGCCATGACGCGACCAAAGGTGGTTTGCGCTCTGGCGGCTATCCACAGCCCGCGCATGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGc         >  1:4485157/1‑150 (MQ=255)
                                                                                                                                            ttACGAAGACATAAACGAATAATATCAGCATTAAGCCGCCAGCCATGACGCGACCAAAGGTGGTTTGCGCTCTGGCGGCTATCCACAGCCCGCGCATGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGct        >  2:3501236/1‑150 (MQ=255)
                                                                                                                                               cGAAGACATAAACGAATAATATCAGCATTAAGCCGCCAGCCATGACGCGACCAAAGGTGGTTTGCGCTCTGGCGGCTATCCACAGCCCGCGCATGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCTCtt     <  1:2764965/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CGACAACATTTTCCTGTGGGTGTGGATTGACATTACAATCCCGAACCCAGCCATCAGCACAATTAGCGCCGATCCTCCATAACTGACCAGTGGGAGCGGAACCCCTACAACCGGCAGAATACCGCTTACCATACCAATATTTACGAAGACATAAACGAATAATATCAGCATTAAGCCGCCAGCCATGACGCGACCAAAGGTGGTTTGCGCTCTGGCGGCTATCCACAGCCCGCGCATGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCTCTTCCG  >  CP009273/660669‑660964

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: