Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,209,162 C→A 100% I9I (ATC→ATA yohO → putative membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,209,1620CA98.2% 1106.4 / 6.5 334I9I (ATC→ATAyohOputative membrane protein
Reads supporting (aligned to +/- strand):  ref base C (3/3);  new base A (208/120);  total (211/123)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 334 total aligned reads displayed.

ACAGCAGCAACTCACCAGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATCGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACG  >  CP009273/2209016‑2209306
                                                                                                                                                  |                                                                                                                                                
aCAGCAGCAACTCACCAGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGc                                                                                                                                                <  1:3420579/149‑1 (MQ=255)
  agcagcAACTCACCAGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGccc                                                                                                                                              <  2:736158/149‑1 (MQ=255)
  agcagcAACTCACCAGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGccc                                                                                                                                              <  1:326425/149‑1 (MQ=255)
     agcaACTCACCAGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCcctgt                                                                                                                                           >  2:3299861/1‑149 (MQ=255)
     agcaACTCACCAGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACAGACGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCcctgt                                                                                                                                           >  2:1330053/1‑149 (MQ=255)
      gcaACTCACCAGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCcctgtt                                                                                                                                          <  2:1563792/149‑1 (MQ=255)
       caACTCACCAGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCcctgttc                                                                                                                                         >  2:1104635/1‑149 (MQ=255)
           tCACCAGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCcctgttcctgt                                                                                                                                     >  2:2777244/1‑149 (MQ=255)
             aCCAGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCCTTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGttt                                                                                                                                   >  2:1023155/1‑149 (MQ=255)
              ccAGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTa                                                                                                                                  <  1:3344165/149‑1 (MQ=255)
               cAGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTAt                                                                                                                                 <  1:4666608/149‑1 (MQ=255)
               cAGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTAt                                                                                                                                 <  2:4493446/149‑1 (MQ=255)
               cAGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTAt                                                                                                                                 >  2:2831116/1‑149 (MQ=255)
                aGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATg                                                                                                                                >  2:3142996/1‑149 (MQ=255)
                aGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATg                                                                                                                                >  2:1785463/1‑149 (MQ=255)
                aGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATg                                                                                                                                >  2:3375453/1‑149 (MQ=255)
                aGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTGAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATg                                                                                                                                >  1:3500738/1‑149 (MQ=255)
                 gTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATgg                                                                                                                               >  1:1693026/1‑149 (MQ=255)
                  ttGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGc                                                                                                                              <  1:3140044/149‑1 (MQ=255)
                        aGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCCTTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAgg                                                                                                                        >  2:4230334/1‑149 (MQ=255)
                              aGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGact                                                                                                                  >  1:984045/1‑147 (MQ=255)
                               ggCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAAtt                                                                                                                 >  1:423712/1‑149 (MQ=255)
                                gCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATtg                                                                                                                >  1:2972291/1‑149 (MQ=255)
                                gCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATtg                                                                                                                >  2:3075357/1‑149 (MQ=255)
                                 ccATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGACTtgg                                                                                                               >  2:1694393/1‑149 (MQ=255)
                                  cATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATtggt                                                                                                              >  2:3866695/1‑149 (MQ=255)
                                   aTGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATtggtg                                                                                                             >  2:1699185/1‑149 (MQ=255)
                                   aTGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATtggtg                                                                                                             >  1:159975/1‑149 (MQ=255)
                                    tGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATtggtgg                                                                                                            <  1:2866882/149‑1 (MQ=255)
                                    tGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATtggtgg                                                                                                            >  2:2016374/1‑149 (MQ=255)
                                    tGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATtggtgg                                                                                                            >  2:898363/1‑149 (MQ=255)
                                     gATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGc                                                                                                           <  2:3470268/149‑1 (MQ=255)
                                        atTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTc                                                                                                        >  1:3149020/1‑149 (MQ=255)
                                         ttttCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCa                                                                                                       >  2:1421115/1‑149 (MQ=255)
                                          tttCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCAt                                                                                                      <  2:2667111/149‑1 (MQ=255)
                                             cccACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCt                                                                                                   >  2:1932639/1‑149 (MQ=255)
                                               cACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCg                                                                                                 >  1:1221411/1‑149 (MQ=255)
                                                    ggCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGt                                                                                            <  1:249257/149‑1 (MQ=255)
                                                     gCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGtt                                                                                           >  2:3187819/1‑149 (MQ=255)
                                                       aTTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTtat                                                                                         <  1:2733389/149‑1 (MQ=255)
                                                        tttAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGACTTGGTGGCGTCATGCTCGCAGGTtata                                                                                        >  2:1297053/1‑149 (MQ=255)
                                                          tAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTATTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTAATAGCCCTGTTCCTGTTAATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATAcc                                                                                      >  1:1415381/1‑149 (MQ=255)
                                                           aaTGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCt                                                                                     >  2:1875778/1‑149 (MQ=255)
                                                              gATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTa                                                                                  >  1:2525271/1‑149 (MQ=255)
                                                               ttgccatAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAACTTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTAt                                                                                 <  2:4566221/143‑1 (MQ=255)
                                                                  cGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTAtttt                                                                              >  2:1966887/1‑149 (MQ=255)
                                                                  cGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTAtttt                                                                              <  2:436904/149‑1 (MQ=255)
                                                                  cGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTAtttt                                                                              <  2:1896804/149‑1 (MQ=255)
                                                                     tAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCg                                                                           >  2:3797153/1‑149 (MQ=255)
                                                                       aCAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTg                                                                         >  1:4597143/1‑149 (MQ=255)
                                                                        cAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGc                                                                        <  2:2315871/149‑1 (MQ=255)
                                                                           aTGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTgg                                                                     >  1:2399984/1‑149 (MQ=255)
                                                                              cgcTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTa                                                                  >  1:3989078/1‑149 (MQ=255)
                                                                               gcTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTaa                                                                 >  2:2929744/1‑149 (MQ=255)
                                                                               gcTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTaa                                                                 <  1:788312/149‑1 (MQ=255)
                                                                                 tttACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAgc                                                               >  2:4733718/1‑149 (MQ=255)
                                                                                   tACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAgcgc                                                             <  1:3499476/149‑1 (MQ=255)
                                                                                     cTGTACAATCCTTTCGTTAACAGAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCt                                                           <  2:2187716/149‑1 (MQ=255)
                                                                                       gTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGAGTGCTGGCTAAGCGACTGc                                                         >  1:590798/1‑149 (MQ=255)
                                                                                        tACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCa                                                        >  2:4759443/1‑149 (MQ=255)
                                                                                        tACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCa                                                        <  1:3800587/149‑1 (MQ=255)
                                                                                             tCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGc                                                   >  1:3734224/1‑149 (MQ=255)
                                                                                             tCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGc                                                   >  1:2292428/1‑149 (MQ=255)
                                                                                              ccTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGcc                                                  <  2:3625530/149‑1 (MQ=255)
                                                                                                tttCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCtt                                                >  1:3099511/1‑149 (MQ=255)
                                                                                                 ttCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCttt                                               >  2:689920/1‑149 (MQ=255)
                                                                                                       aaCATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCTCCAGCCTTTCAAACa                                         <  1:721490/149‑1 (MQ=255)
                                                                                                       aaCATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACa                                         >  2:782786/1‑149 (MQ=255)
                                                                                                       aaCATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACa                                         <  1:1933979/149‑1 (MQ=255)
                                                                                                       aaCATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACa                                         >  2:2989792/1‑149 (MQ=255)
                                                                                                           tAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCg                                     >  1:1543153/1‑149 (MQ=255)
                                                                                                           tAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCg                                     >  2:4159346/1‑149 (MQ=255)
                                                                                                           tAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCg                                     >  1:2586277/1‑149 (MQ=255)
                                                                                                            aaGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCgg                                    <  2:2889271/149‑1 (MQ=255)
                                                                                                             aGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGt                                   >  1:1860403/1‑149 (MQ=255)
                                                                                                              ggAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTc                                  <  1:2129045/149‑1 (MQ=255)
                                                                                                              ggAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTc                                  <  2:1613150/149‑1 (MQ=255)
                                                                                                               gAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCt                                 >  1:1952106/1‑149 (MQ=255)
                                                                                                               gAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCt                                 >  1:4088694/1‑149 (MQ=255)
                                                                                                               gAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCt                                 >  1:1005903/1‑149 (MQ=255)
                                                                                                                aGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTg                                <  1:4197080/149‑1 (MQ=255)
                                                                                                                aGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTg                                >  2:465749/1‑149 (MQ=255)
                                                                                                                 gTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGc                               <  1:4477628/149‑1 (MQ=255)
                                                                                                                   gCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCga                             >  1:1521889/1‑149 (MQ=255)
                                                                                                                   gCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCga                             >  1:3233394/1‑149 (MQ=255)
                                                                                                                    cATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCgat                            <  2:2455834/149‑1 (MQ=255)
                                                                                                                       tATGCGCATAGCTAAAATTGGGGTCATCGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCgatgat                         <  1:91296/149‑1 (MQ=255)
                                                                                                                            gcATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATcgccg                    >  2:1484272/1‑149 (MQ=255)
                                                                                                                              aTAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCgccgcc                  >  1:511521/1‑149 (MQ=255)
                                                                                                                               tAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCa                 >  2:2231173/1‑149 (MQ=255)
                                                                                                                               tAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCa                 <  1:413055/149‑1 (MQ=255)
                                                                                                                                 gCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGt               >  2:585515/1‑149 (MQ=255)
                                                                                                                                  cTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTg              >  1:360327/1‑149 (MQ=255)
                                                                                                                                   tAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGc             <  1:429503/149‑1 (MQ=255)
                                                                                                                                    aaaaTTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGcc            >  1:1742555/1‑149 (MQ=255)
                                                                                                                                    aaaaTTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCAGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGcc            <  1:2512740/149‑1 (MQ=255)
                                                                                                                                      aaTTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGccac          >  2:4091331/1‑149 (MQ=255)
                                                                                                                                         tGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAg       <  2:1988666/149‑1 (MQ=255)
                                                                                                                                         tGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAg       <  2:756707/149‑1 (MQ=255)
                                                                                                                                            ggTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTaa    <  2:4018686/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
ACAGCAGCAACTCACCAGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATCGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACG  >  CP009273/2209016‑2209306

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: