Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,984,987 C→A 100% T38T (ACG→ACT yecA ← UPF0149 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,984,9870CA97.2% 664.4 / 8.5 214T38T (ACG→ACTyecAUPF0149 family protein
Reads supporting (aligned to +/- strand):  ref base C (0/6);  new base A (111/97);  total (111/103)
Fisher's exact test for biased strand distribution p-value = 1.15e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 214 total aligned reads displayed.

CGGTATCGGCCATATGTTGAAAAGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCCGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGATCTCAACTCGC  >  CP009273/1984842‑1985128
                                                                                                                                                 |                                                                                                                                             
  gTATCGGCCATATGTTGAAAAGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCa                                                                                                                                           <  2:424984/148‑1 (MQ=255)
  gTATCGGCCATATGTTGAAAAGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCa                                                                                                                                           <  2:1519267/148‑1 (MQ=255)
       ggCCATATGTTGAAAAGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATaaa                                                                                                                                      <  2:1637585/148‑1 (MQ=255)
            tatGTTGAAAAGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTc                                                                                                                                 <  2:1679148/148‑1 (MQ=255)
            tatGTTGAAAAGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTc                                                                                                                                 >  1:767393/1‑148 (MQ=255)
             atGTTGAAAAGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTcc                                                                                                                                >  2:1893147/1‑148 (MQ=255)
                  gAAAAGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTc                                                                                                                           >  1:1407827/1‑148 (MQ=255)
                   aaaaGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTcc                                                                                                                          >  1:3408532/1‑148 (MQ=255)
                    aaaGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCg                                                                                                                         >  1:1960535/1‑148 (MQ=255)
                         cAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCcaca                                                                                                                    <  2:48156/148‑1 (MQ=255)
                         cAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCcaca                                                                                                                    >  1:590092/1‑148 (MQ=255)
                            aTTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGACCTCAACACCGCAGTCAATAAACCGTCCCGCTCCGCCACATCa                                                                                                                 >  2:2457590/1‑148 (MQ=255)
                            aTTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCa                                                                                                                 >  2:2853211/1‑148 (MQ=255)
                             ttCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCaa                                                                                                                <  1:1538784/148‑1 (MQ=255)
                               cATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAgg                                                                                                              >  1:1498104/1‑148 (MQ=255)
                                aTAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGa                                                                                                             <  2:1328659/148‑1 (MQ=255)
                                  aaaGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTAAATCTCTTGCGGAGAACTCAACACCGAAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATg                                                                                                           <  1:1987405/148‑1 (MQ=255)
                                   aaGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATgg                                                                                                          <  2:1018187/148‑1 (MQ=255)
                                    aGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGAt                                                                                                            <  2:2853919/145‑1 (MQ=255)
                                    aGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGAt                                                                                                            >  1:2853919/1‑145 (MQ=255)
                                    aGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGc                                                                                                         >  2:907656/1‑148 (MQ=255)
                                    aGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGc                                                                                                         >  1:3320439/1‑148 (MQ=255)
                                    aGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGc                                                                                                         >  2:3066467/1‑148 (MQ=255)
                                    aGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGc                                                                                                         >  1:1609991/1‑148 (MQ=255)
                                          tCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCCGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTc                                                                                                   <  1:1723004/148‑1 (MQ=255)
                                           cATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCa                                                                                                  >  1:3344641/1‑148 (MQ=255)
                                            aTCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCGTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAg                                                                                                 <  1:1829507/148‑1 (MQ=255)
                                                ctTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGtgtt                                                                                             >  1:2153837/1‑148 (MQ=255)
                                                  ttCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGtgttgt                                                                                           >  1:2671708/1‑148 (MQ=255)
                                                  ttCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGtgttgt                                                                                           >  2:64713/1‑148 (MQ=255)
                                                  ttCTCTGACGCACAGCGCTGCACATATTCAGCCCCACCCCACACTGCAAAGAGCCACTGTTCCTTTTCAATCTCTTGCTGTGAACTCTTCACGGACTTCTATAAACCGTCCAGCTCCGCCACATCAAGGATTGCTTGTTCAGtgttgt                                                                                           <  2:3091695/148‑1 (MQ=255)
                                                    ctctGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCTGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTAt                                                                                         <  2:2659936/148‑1 (MQ=255)
                                                     tctGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTAtt                                                                                        >  2:1194126/1‑148 (MQ=255)
                                                         aCGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGt                                                                                    >  2:2104080/1‑148 (MQ=255)
                                                          cGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTc                                                                                   <  1:3433541/148‑1 (MQ=255)
                                                            cccAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAg                                                                                 >  1:3374324/1‑148 (MQ=255)
                                                               aGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAat                                                                              <  2:795626/148‑1 (MQ=255)
                                                               aGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAat                                                                              >  2:428453/1‑148 (MQ=255)
                                                                gcgcGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAAAACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAata                                                                             >  2:2230501/1‑148 (MQ=255)
                                                                 cgcgGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCCGTGTTGTATTTGGTCAGAatat                                                                            >  2:2594054/1‑148 (MQ=255)
                                                                 cgcgGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAatat                                                                            >  2:1063244/1‑148 (MQ=255)
                                                                   cgGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCg                                                                          >  1:1701384/1‑148 (MQ=255)
                                                                    ggCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGt                                                                         <  1:3540157/148‑1 (MQ=255)
                                                                     gCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTc                                                                        <  2:1498104/148‑1 (MQ=255)
                                                                      cacaTAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTcc                                                                       >  2:1824587/1‑148 (MQ=255)
                                                                        caTAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAg                                                                     >  2:1353568/1‑148 (MQ=255)
                                                                        caTAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAg                                                                     >  1:1426342/1‑148 (MQ=255)
                                                                         aTAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGc                                                                    <  2:2241788/148‑1 (MQ=255)
                                                                         aTAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGc                                                                    <  2:1319394/148‑1 (MQ=255)
                                                                          tAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGcc                                                                   >  2:2741160/1‑148 (MQ=255)
                                                                           aGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCa                                                                  >  2:2945730/1‑148 (MQ=255)
                                                                           aGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCa                                                                  >  2:1166875/1‑148 (MQ=255)
                                                                           aGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCa                                                                  <  1:2970093/148‑1 (MQ=255)
                                                                            gTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCAt                                                                 >  1:855068/1‑148 (MQ=255)
                                                                            gTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCAt                                                                 >  1:2641016/1‑148 (MQ=255)
                                                                              cAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTc                                                               <  2:681629/148‑1 (MQ=255)
                                                                                   ccaccccacACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACt                                                          >  1:1871920/1‑148 (MQ=255)
                                                                                     accccacACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCa                                                        >  1:1626425/1‑148 (MQ=255)
                                                                                       cccacACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACt                                                      >  2:2004956/1‑148 (MQ=255)
                                                                                         cacaCGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCCGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACttt                                                    <  1:2766318/148‑1 (MQ=255)
                                                                                         cacaCGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACttt                                                    <  1:1036691/148‑1 (MQ=255)
                                                                                         cacaCGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGAGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACttt                                                    >  2:1732039/1‑148 (MQ=255)
                                                                                           cacGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCGAGCCATTCCAACTCACTTTCg                                                  <  1:2404355/148‑1 (MQ=255)
                                                                                           cacGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCg                                                  >  1:3240807/1‑148 (MQ=255)
                                                                                           cacGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCg                                                  >  2:242063/1‑148 (MQ=255)
                                                                                            acGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCCGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGt                                                 <  1:641019/148‑1 (MQ=255)
                                                                                            acGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGt                                                 >  1:3495704/1‑148 (MQ=255)
                                                                                            acGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGt                                                 >  1:3492599/1‑148 (MQ=255)
                                                                                               gCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTa                                              <  1:38573/148‑1 (MQ=255)
                                                                                                cAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTaa                                             <  1:1270771/148‑1 (MQ=255)
                                                                                                 aaCCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAAc                                            >  1:2219897/1‑148 (MQ=255)
                                                                                                  aCCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACg                                           >  1:10840/1‑148 (MQ=255)
                                                                                                  aCCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTGTAACg                                           >  1:11244/1‑148 (MQ=255)
                                                                                                     aGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTc                                        <  2:630134/148‑1 (MQ=255)
                                                                                                     aGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTc                                        <  1:64713/148‑1 (MQ=255)
                                                                                                     aGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTc                                        <  1:558033/148‑1 (MQ=255)
                                                                                                         aCTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGt                                    >  2:815331/1‑148 (MQ=255)
                                                                                                          cTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGtt                                   >  2:1204226/1‑148 (MQ=255)
                                                                                                           tGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGttt                                  <  2:350009/148‑1 (MQ=255)
                                                                                                             ttCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTc                                <  2:2340959/148‑1 (MQ=255)
                                                                                                               ccGGTTCAATCTCGTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCAt                              <  1:1824587/148‑1 (MQ=255)
                                                                                                                    tCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCt                         <  2:3361005/148‑1 (MQ=255)
                                                                                                                        tctcTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTc                     <  2:837913/148‑1 (MQ=255)
                                                                                                                        tctcTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTc                     <  2:3489862/148‑1 (MQ=255)
                                                                                                                         ctctTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTcc                    <  2:70860/148‑1 (MQ=255)
                                                                                                                         ctctTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTcc                    <  2:3212545/148‑1 (MQ=255)
                                                                                                                           ctTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCtt                  >  2:2856166/1‑148 (MQ=255)
                                                                                                                              gCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTtgtg               <  1:907656/148‑1 (MQ=255)
                                                                                                                              gCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTtgtg               <  1:3066467/148‑1 (MQ=255)
                                                                                                                               cGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTgg              <  2:3415239/148‑1 (MQ=255)
                                                                                                                                ggAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGa             <  1:160235/148‑1 (MQ=255)
                                                                                                                                 gagaACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGAt            <  1:1199362/148‑1 (MQ=255)
                                                                                                                                  agaTCTCAACACCGCAGTCAATACACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTACCGTTTTCATACGCTTTTCCTTGTGGAtc           <  2:781767/148‑1 (MQ=255)
                                                                                                                                  agaACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGAtc           <  1:2945730/148‑1 (MQ=255)
                                                                                                                                  agaACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGAtc           <  1:1166875/148‑1 (MQ=255)
                                                                                                                                   gaACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGAtct          <  1:3026919/148‑1 (MQ=255)
                                                                                                                                   gaACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGAtct          <  2:1606581/148‑1 (MQ=255)
                                                                                                                                    aaCTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGAtctc         >  2:19418/1‑148 (MQ=255)
                                                                                                                                     aCTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGATCTCa        <  1:1252473/148‑1 (MQ=255)
                                                                                                                                       tCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGATCTCAAc      >  1:1242608/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                             
CGGTATCGGCCATATGTTGAAAAGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCCGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGATCTCAACTCGC  >  CP009273/1984842‑1985128

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: