Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 885,619 G→T 6.8% A25A (GCG→GCT ybjM → inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273885,6190GT 6.8% 599.3 / 14.8 160A25A (GCG→GCTybjMinner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (74/75);  new base T (7/4);  total (81/79)
Fisher's exact test for biased strand distribution p-value = 5.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

Only 100 of 163 total aligned reads displayed.

AAGCGGCTTTATCTTACTGAACGCCACACTGGCGAAAAATGTGTTCGATAGACGCAGTGTCAGGAGGAACGAGTGAAACATAAACAACGTTGGGCGGGGGCAATCTGCTGTTTTGTCCTCTTCATTGTGGTGTGCCTTTTTCTGGCGACGCACATGAAAGGCGCTTTTCGGGCTGCCGGGCATCCTGAAATCGGCTTGCTATTTTTCATTCTTCCTGGAGCAGTCGCCAGCTTCTTTTCACAGCGTAGAGAAGTCCTGAAACCTCTGTTTGGCGCAATGCTGGCGGCACCCTGT  >  CP009273/885473‑885766
                                                                                                                                                  |                                                                                                                                                   
aaGCGGCTTTATCTTACTGAACGCCACACTGGCGAAAAATGTGTTCGATAGACGCAGTGTCAGGAGGAACGAGTGAAACATAAACAACGTTGGGCGGGGGCAATCTGCTGTTTTGTCCTCTTCATTGTGGTGTGCCTTTTTCTGGCGa                                                                                                                                                    <  1:374576/148‑1 (MQ=255)
 aGCGGCTTTATCTTACTGAACGCCACACTGGCGAAAAATGTGTTCGATAGACGCAGTGTCAGGAGGAACGAGTGAAACATAAACAACGTTGGGCGGGGGCAATCTGCTGTTTTGTCCTCTTCATTGTGGTGTGCCTTTTTCTGGCGAc                                                                                                                                                   >  1:1892107/1‑148 (MQ=255)
 aGCGGCTTTATCTTACTGAACGCCACACTGGCGAAAAATGTGTTCGATAGACGCAGTGTCAGGAGGAACGAGTGAAACATAAACAACGTTGGGCGGGGGCAATCTGCTGTTTTGTCCTCTTCATTGTGGTGTGCCTTTTTCTGGCGAc                                                                                                                                                   >  1:1891977/1‑148 (MQ=255)
  gCGGCTTTATCTTACTGAACGCCACACTGGCGAAAAATGTGTTCGATAGACGCAGTGTCAGGAGGAACGAGTGAAACATAAACAACGTTGGGCGGGGGCAATCTGCTGTTTTGTCCTCTTCATTGTGGTGTGCCTTTTTCTGGCGACg                                                                                                                                                  >  1:1893185/1‑148 (MQ=255)
   cGGCTTTATCTTACTGAACGCCACACTGGCGAAAAATGTGTTCGATAGACGCAGTGTCAGGAGGAACGAGTGAAACATAAACAACGTTGGGCGGGGGCAATCTGCTGTTTTGTCCTCTTCATTGTGGTGTGCCTTTTTCTGGCGACGc                                                                                                                                                 >  1:1361529/1‑148 (MQ=255)
     gCTTTATCTTACTGAACGCCACACTGGCGAAAAATGTGTTCGATAGACGCAGTGTCAGGAGGAACGAGTGAAACATAAACAACGTTGGGCGGGGGCAATCTGCTGTTTTGTCCTCTTCATTGTGGTGTGCCTTTTTCTGGCGACGcac                                                                                                                                               <  1:3242487/148‑1 (MQ=255)
      cTTTATCTTACTGAACGCCACACTGGCGAAAAATGTGTTCGATAGACGCAGTGTCAGGAGGAACGAGTGAAACATAAACAACGTTGGGCGGGGGCAATCTGCTGTTTTGTCCTCTTCATTGTGGTGTGCCTTTTTCTGGCTACGcaca                                                                                                                                              >  1:1141107/1‑148 (MQ=255)
           tCTTACTGAACGCCACACTGGCGAAAAATGTGTTCGATAGACGCAGTGTCAGGAGGAACGAGTGAAACATAAACAACGTTGGGCGGGGGCAATCTGCTGTTTTGTCCTCTTCATTGTGGTGTGCCTTTTTCTGGCGACGCACATGaaa                                                                                                                                         >  1:1195126/1‑148 (MQ=255)
              tACTGAACGCCACACTGGCGAAAAATGTGTTCGATAGACGCAGTGTCAGGAGGAACGAGTGAAACATAAACAACGTTGGGCGGGGGCAATCTGCTGTTTTGTCCTCTTCATTGTGGTGTGCCTTTTTCTGGCGACGCACATGAAAGgc                                                                                                                                      >  2:324302/1‑148 (MQ=255)
               aCTGAACGCCACACTGGCGAAAAATGTGTTCGATAGACGCAGTGTCAGGAGGAACGAGTGAAACATAAACAACGTTGGGCGGGGGCAATCTGCTGTTTTGTCCTCTTCATTGTGGTGTGCCTTTTTCTGGCGACGCACATGAAAGgcg                                                                                                                                     >  1:2084057/1‑148 (MQ=255)
                cTGAACGCCACACTGGCGAAAAATGTGTTCGATAGACGCAGTGTCAGGAGGAACGAGTGAAACATAAACAACGTTGGGCGGGGGCAATCTGCTGTTTTGTCCTCTTCATTGTGGTGTGCCTTTTTCTGGCGACGCACATGAAAGgcgc                                                                                                                                    <  1:2217929/148‑1 (MQ=255)
                 tGAACGCCACACTGGCGAAAAATGTGTTCGATAGACGCAGTGTCAGGAGGAACGAGTGAAACATAAACAACGTTGGGCGGGGGCAATCTGCTGTTTTGTCCTCTTCATTGTGGTGTGCCTTTTTCTGGCGACGCACATGAAAGGCGCt                                                                                                                                   >  1:1954764/1‑148 (MQ=255)
                   aaCGCCACACTGTCGAAATATGTGTTCTATAGACGCAGTGTCAGGAGGAACGAGTGAAACAATAACAACGTTGGGCGGGGGGAATCTGCTGTTTTGTCCTCTTCATTGTGGTGTGCCTTTTTCTGGCTACGCACATGAAATGCGCttt                                                                                                                                 <  1:2204910/148‑1 (MQ=255)
                     cGCCACACTGGCGAAAAATGTGTTCGATAGACGCAGTGTCAGGAGGAACGAGTGAAACATAAACAACGTTGGGCGGGGGCAATCTGCTGTTTTGTCCTCTTCATTGTGGTGTGCCTTTTTCTGGCGACGCACATGAAAGGCGCTTTTc                                                                                                                               <  1:2644856/148‑1 (MQ=255)
                      gCCACACTGGCGAAAAATGTGTTCGATAGACGCAGTGTCAGGAGGAACGAGTGAAACATAAACAACGTTGGGCGGGGGCAATCTGCTGTTTTGTCCTCTTCATTGTGGTGTGCCTTTTTCTGGCGACGCACATGAAAGGCGCTTTTCg                                                                                                                              <  1:2663767/148‑1 (MQ=255)
                          cacTGGCGAAAAATGTGTTCGATAGACGCAGTGTCAGGAGGAACGAGTGAAACATAAACAACGTTGGGCGGGGGCAATCTGCTGTTTTGTCCTCTTCATTGTGGTGTGCCTTTTTCTGGCGACGCACATGAAAGGCGCTTTTCGGGCt                                                                                                                          >  1:1353474/1‑148 (MQ=255)
                           acTGGCGAAAAATGTGTTCGATAGACGCAGTGTCAGGAGGAACGATTGATACATAAACAACGTTGGGCGGGGGCAATCTGCTGTTTTGTCCTCTTCATTGTGGTGTGCCTTTTTCTGGCGACGCACATGAAAGGCGCTTTTCGGGCTg                                                                                                                         <  2:1434651/148‑1 (MQ=255)
                            cTGGCGAAAAATGTGTTCGATAGACGCAGTGTCAGGAGGAACGAGTGAAACATAAACAACGTTGGGCGGGGGCAATCTGCTGTTTTGTCCTCTTCATTGTGGTGTGCCTTTTTCTGGCGACGCACATGAAAGGCGCTTTTCGGGCTGc                                                                                                                        <  1:2448056/148‑1 (MQ=255)
                                  aaaaaTGTGTTCGATAGACGCAGTGTCAGGAGGAACGAGTGAAACATTTTGAACGTTGGGCGGGGGCAATCTGCTGTTTTGTCCTCTTCATTGTGGTGTGCCTTTTTCTGGCGACGCACATGAAAGGCGCTTTTCGGGCTGCCGGGCa                                                                                                                  <  2:2663237/148‑1 (MQ=255)
                                  aaaaaTGTGTTCGATAGACGCAGTGTCAGGAGGAACGAGTGAAACATAAACAACGTTGGGCGGGGGCAATCTGCTGTTTTGTCCTCTTCATTGTGGTGTGCCTTTTTCTGGCGACGCACATGAAAGGCGCTTTTCGGGCTGCCGGGCa                                                                                                                  <  2:477890/148‑1 (MQ=255)
                                   aaaaTGTGTTCGATAGACGCAGTGTCAGGAGGAACGAGTGAAACATAAACAACGTTGGGCGGGGGCAATCTGCTGTTTTGTCCTCTTCATTGTGGTGTGCCTTTTTCTGGCGACGCACATGAAAGGCGCTTTTCGGGCTGCCGGGCAt                                                                                                                 >  2:643803/1‑148 (MQ=255)
                                     aaTGTGTTCGATAGACGCAGTGTCAGGAGGAACGAGTGAAACATAAACAACGTTGGGCGGGGGCAATCTGCTGTTTTGTCCTCTTCATTGTGGTGTGCCTTTTTCTGGCGACGCACATGAAAGGCGCTTTTCGGGCTGCCGGGCATcc                                                                                                               <  1:2360929/148‑1 (MQ=255)
                                     aaTGTGTTCGATAGACGCAGTGTCAGGAGGAACGAGTGAAACATAAACAACGTTGGGCGGGGGCAATCTGCTGTTTTGTCCTCTTCATTGTGGTGTGCCTTTTTCTGGCGACGCACATGAAAGGCGCTTTTCGGGCTGCCGGGCATcc                                                                                                               >  2:2354915/1‑148 (MQ=255)
                                     aaTGTGTTCGATAGACGCAGTGTCAGGAGGAACGAGTGAAACATAAACAACGTTGGGCGGGGGCAATCTGCTGTTTTGTCCTCTTCATTGTGGTGTGCCTTTTTCTGGCGACGCACATGAAAGGCGCTTTTCGGGCTGCCGGGCATcc                                                                                                               >  2:2354944/1‑148 (MQ=255)
                                           ttCGATAGACGCAGTGTCAGGAGGAACGAGTGAAACATAAACAACGTTGGGCGGGGGCAATCTGCTGTTTTGTCCTCTTCATTGTGGTGTGCCTTTTTCTGGCGACGCACATGAAAGGCGCTTTTCGGGCTGCCGGGCATCCTGAAAt                                                                                                         <  1:2591997/148‑1 (MQ=255)
                                           ttCGATAGACGCAGTGTCAGGAGGAACGAGTGAAACATAAACAACGTTGGGCGGGGGCAATCTGCTGTTTTGTCCTCTTCATTGTGGTGTGCCTTTTTCTGGCGACGCACATGAAAGGCGCTTTTCGGGCTGCCGGGCATCCTGAAAt                                                                                                         <  2:290926/148‑1 (MQ=255)
                                           ttCGATAGACGCAGTGTCAGGAGGAACGAGTGAAACATAAACAACGTTGGGCGGGGGCAATCTGCTGTTTTGTCCTCTTCATTGTGGTGTGCCTTTTTCTGGCGACGCACATGAAAGGCGCTTTTCGGGCTGCCGGGCATCCTGAAAt                                                                                                         <  2:3278510/148‑1 (MQ=255)
                                           ttCGATAGACGCAGTGTCAGGAGGAACGAGTGAAACATAAACAACGTTGGGCGGGGGCAATCTGCTGTTTTGTCCTCTTCATTGTGGTGTGCCTTTTTCTGGCGACGCACATGAAAGGCGCTTTTCGGGCTGCCGGGCATCCTGAAAt                                                                                                         <  2:290997/148‑1 (MQ=255)
                                               aTAGACGCAGTGTCAGGAGGAACGAGTGAAACATAAACAACGTTGGGCGGGGGCAATCTGCTGTTTTGTCCTCTTCATTGTGGTGTGCCTTTTTCTGGCGACGCACATGAAAGGCGCTTTTCGGGCTGCCGGGCATCCTGAAATCGGc                                                                                                     >  1:435368/1‑148 (MQ=255)
                                                   aCGCAGTGTCAGGAGGAACGAGTGAAACATAAACAACGTTGGGCGGGGGCAATCTGCTGTTTTGTCCTCTTCATTGTGGTGTGCCTTTTTCTGGCGACGCACATGAAAGGCGCTTTTCGGGCTGCCGGGCATCCTGAAATCGGCTTGc                                                                                                 >  2:3090281/1‑148 (MQ=255)
                                                         tgtCAGGAGGAACGAGTGAAACATAAACAACGTTGGGCGGGGGCAATCTGCTGTTTTGTCCTCTTCATTGTGGTGTGCCTTTTTCTGGCGACGCACATGAAAGGCGCTTTTCGGGCTGCCGGGCATCCTGAAATCGGCTTGCTAtttt                                                                                           >  1:2735189/1‑148 (MQ=255)
                                                          gtCAGGAGGAACGAGTGAAACATAAACAACGTTGGGCGGGGGCAATCTGCTGTTTTGTCCTCTTCATTGTGGTGTGCCTTTTTCTGGCGACGCACATGAAAGGCGCTTTTCGGGCTGCCGGGCATCCTGAAATCGGCTTGCTAttttt                                                                                          <  1:1838164/148‑1 (MQ=255)
                                                               gaggaACGAGTGAAACATAAACAACGTTGGGCGGGGGCAATCTGCTGTTTTGTCCTCTTCATTGTGGTGTGCCTTTTTCTGGCTACGCACATGAAAGGCGCTTTTCGGGCTGCCGGGCATCCTGAAATCGGCTTGCTATTTTTCAttc                                                                                     <  2:2023624/148‑1 (MQ=255)
                                                                 ggaACGAGTGAAACATAAACAACGTTGGGCGGGGGCAATCTGCTGTTTTGTCCTCTTCATTGTGGTGTGCCTTTTTCTGGCGACGCACATGAAAGGCGCTTTTCGGGCTGCCGGGCATCCTGAAATCGGCTTGCTATTTTTCAttctt                                                                                   <  1:2921213/148‑1 (MQ=255)
                                                                  gaACGAGTGAAACATAAACAACGTTGGGCGGGGGCAATCTGCTGTTTTGTCCTCTTCATTGTGGTGTGCCTTTTTCTGGCGACGCACATGAAAGGCGCTTTTCGGGCTGCCGGGCATCCTGAAATCGGCTTGCTATTTTTCAttcttc                                                                                  >  1:2079762/1‑148 (MQ=255)
                                                                  gaACGAGTGAAACATAAACAACGTTGGGCGGGGGCAATCTGCTGTTTTGTCCTCTTCATTGTGGTGTGCCTTTTTCTGGCGACGCACATGAAAGGCGCTTTTCGGGCTGCCGGGCATCCTGAAATCGGCTTGCTATTTTTCAttcttc                                                                                  <  1:65793/148‑1 (MQ=255)
                                                                    aCGAGTGAAACATAAACAACGTTGGGCGGGGGCAATCTGCTGTTTTGTCCTCTTCGTTGTGGTGTGCCTTTTTCTGGCGACGCACATGAAAGGCGCTTTTCGGGCTGCCGGGCATCCTGAAATCGGCTTGCTATTTTTCATTCTTCCt                                                                                >  2:2292395/1‑148 (MQ=255)
                                                                    aCGAGTGAAACATAAACAACGTTGGGCGGGGGCAATCTGCTGTTTTGTCCTCTTCATTGTGGTGTGCCTTTTTCTGGCGACGCACATGAAAGGCGCTTTTCGGGCTGCCGGGCATCCTGAAATCGGCTTGCTATTTTTCATTCTTCCt                                                                                <  2:2672291/148‑1 (MQ=255)
                                                                     cGAGTGAAACATAAACAACGTTGGGCGGGGGCAATCTGCTGTTTTGTCCTCTTCATTGTGGTGTGCCTTTTTCTGGCGACGCACATGAAAGGCGCTTTTCGGGCTGCCGGGCATCCTGAAATCGGCTTGCTATTTTTCATTCTTCCTg                                                                               >  1:1319796/1‑148 (MQ=255)
                                                                         tGAAACATAAACAACGTTGGGCGGGGGCAATCTGCTGTTTTGTCCTCTTCATTGTGGTGTGCCTTTTTCTGGCGACGCACATGAAAGGCGCTTTTCGGGCTGCCGGGCATCCTGAAATCGGCTTGCTATTTTTCATTCTTCCTGGAGc                                                                           <  1:35511/148‑1 (MQ=255)
                                                                         tGAAACATAAACAACGTTGGGCGGGGGCAATCTGCTGTTTTGTCCTCTTCATTGTGGTGTGCCTTTTTCTGGCGACGCACATGAAAGGCGCTTTTCGGGCTGCCGGGCATCCTGAAATCGGCTTGCTATTTTTCATTCTTCCTGGAGc                                                                           >  2:3264152/1‑148 (MQ=255)
                                                                          gAAACATAAACAACGTTGGGCGGGGGCAATCTGCTGTTTTGTCCTCTTCATTGTGGTGTGCCTTTTTCTGGCGACGCACATGAAAGGCGCTTTTCGGGCTGCCGGGCATCCTGAAATCGGCTTGCTATTTTTCATTCTTCCTGGAGCa                                                                          >  2:562721/1‑148 (MQ=255)
                                                                           aaaCATAAACAACGTTGGGCGGGGGCAATCTGCTGTTTTGTCCTCTTCATTGTGGTGTGCCTTTTTCTGGCGACGCACATGAAAGGCGCTTTTCGGGCTGCCGGGCATCCTGAAATCGGCTTGCTATTTTTCATTCTTCCTGGAGCAg                                                                         <  1:324302/148‑1 (MQ=255)
                                                                           aaaCATAAACAACGTTGGGCGGGGGCAATCTGCTGTTTTGTCCTCTTCATTGTGGTGTGCCTTTTTCTGGCGACGCACATGAAAGGCGCTTTTCGGGCTGCCGGGCATCCTGAAATCGGCTTGCTATTTTTCATTCTTCCTGGAGCAg                                                                         <  2:340566/148‑1 (MQ=255)
                                                                                   acaacGTTGGGCGGGGGCAATCTGCTGTTTTGTCCTCTTCATTGTGGTGTGCCTTTTTCTGGCGACGCACATGAAAGGCGCTTTTCGGGCTGCCGGGCATCCTGAAATCGGCTTGCTATTTTTCATTCTTCCTGGAGCAGTCGCCAGc                                                                 <  1:2923882/148‑1 (MQ=255)
                                                                                      acGTTGGGCGGGGGCAATCTGCTGTTTTGTCCTCTTCATTGTGGTGTGCCTTTTTCTGGCGACGCACATGAAAGGCGCTTTTCGGGCTGCCGGGCATCCTGAAATCGGCTTGCTATTTTTCATTCTTCCTGGAGCAGTCGCCAGcttc                                                              <  1:2783357/148‑1 (MQ=255)
                                                                                       cGTTGGGCGGGGGCAATCTGCTGTTTTGTCCTCTTCATTGTGGTGTGCCTTTTTCTGGCGACGCACATGAAAGGCGCTTTTCGGGCTGCCGGGCATCCTGAAATCGGCTTGCTATTTTTCATTCTTCCTGGAGCAGTCGCCAGcttct                                                             >  1:2105636/1‑148 (MQ=255)
                                                                                       cGTTGGGCGGGGGCAATCTGCTGTTTTGTCCTCTTCATTGTGGTGTGCCTTTTTCTGGCGACGCACATGAAAGGCGCTTTTCGGGCTGCCGGGCATCCTGAAATCGGCTTGCTATTTTTCATTCTTCCTGGAGCAGTCGCCAGcttct                                                             >  1:1988316/1‑148 (MQ=255)
                                                                                          tGGGCGGGGGCAATCTGCTGTTTTGTCCTCTTCATTGTGGTGTGCCTTTTTCTGGCGACGCACATGAAAGGCGCTTTTCGGGCTGCCGGGCATCCTGAAATCGGCTTGCTATTTTTCATTCTTCCTGGAGCAGTCGCCAGCTTCtttt                                                          <  2:2340246/148‑1 (MQ=255)
                                                                                          ggggCGGGGGCAATCTGCTGTTGTGTCCTCTTCATTGTGGTGTGCCTTTTTCTGGCGACGCACATGAAAGGCGCTTTTCGGGCTGCCGGGCATCCTGAAATCGGCTTGCTATTTTTCATTCTTCCTGGAGCAGTCGCCAGCTTCtttt                                                          <  1:1572129/147‑1 (MQ=255)
                                                                                           gggCGGGGGCAATCTGCTGTTTTGTCCTCTTCATTGTGGTGTGCCTTTTTCTGGCGACGCACATGAAAGGCGCTTTTCGGGCTGCCGGGCATCCTGAAATCGGCTTGCTATTTTTCATTCTTCCTGGAGCAGTCGCCAGCTTCTTTTc                                                         <  1:2024970/148‑1 (MQ=255)
                                                                                           gggCGGGGGCAATCTGCTGTTTTGTCCTCTTCATTGTGGTGTGCCTTTTTCTGGCGACGCACATGAAAGGCGCTTTTCGGGCTGCCGGGCATCCTGAAATCGGCTTGCTATTTTTCATTCTTCCTGGAGCAGTCGCCAGCTTCTTTTc                                                         <  1:2023483/148‑1 (MQ=255)
                                                                                           gggCGGGGGCAATCTGCTGTTTTGTCCTCTTCATTGTGGTGTGCCTTTTTCTGGCGACGCACATGAAAGGCGCTTTTCGGGCTGCCGGGCATCCTGAAATCGGCTTGCTATTTTTCATTCTTCCTGGAGCAGTCGCCAGCTTCTTTTc                                                         <  1:2023479/148‑1 (MQ=255)
                                                                                             gCGGGGGCAATCTGCTGTTTTGTCCTCTTCATTGTGGTGTGCCTTTTTCTGGCGACGCACATGAAAGGCGCTTTTCGGGCTGCCGGGCATCCTGAAATCGGCTTGCTATTTTTCATTCTTCCTGGAGCAGTCGCCAGCTTCTTTTcac                                                       <  2:2261952/148‑1 (MQ=255)
                                                                                                ggggCAATCTGCTGTTTTGTCCTCTTCATTGTGGTGTGCCTTTTTCTGGCGACGCACATGAAAGGCGCTTTTCGGGCTGCCGGGCATCCTGAAATCGGCTTGCTATTTTTCATTCTTCCTGGAGCAGTCGCCAGCTTCTTTTCACAGc                                                    >  2:506936/1‑148 (MQ=255)
                                                                                                 gggCAATCTGCTGTTTTGTCCTCTTCATTGTGGTGTGCCTTTTTCTGGCGACGCACATGAAAGGCGCTTTTCGGGCTGCCGGGCATCCTGAAATCGGCTTGCTATTTTTCATTCTTCCTGGAGCAGTCGCCAGCTTCTTTTCACAGCg                                                   <  2:1920595/148‑1 (MQ=255)
                                                                                                      aTCTGCTGTTTTGTCCTCTTCATTGTGGTGTGCCTTTTTCTGGCGACGCACATGAAAGGCGCTTTTCGGGCTGCCGGGCATCCTGAAATCGGCTTGCTATTTTTCATTCTTCCTGGAGCAGTCGCCAGCTTCTTTTCACAGCGTagag                                              >  2:1019241/1‑148 (MQ=255)
                                                                                                       tCTGCTGTTTTGTCCTCTTCATTGTGGTGTGCCTTTTTCTGGCGACGCACATGAAAGGCGCTTTTCGGGCTGCCGGGCATCCTGAAATCGGCTTGCTATTTTTCATTCTTCCTGGAGCAGTCGCCAGCTTCTTTTCACAGCGTAgaga                                             >  1:937453/1‑148 (MQ=255)
                                                                                                        ctgctgTTTTGTCCTCTTCATTGTGGTGTGCCTTTTTCTGGCGACGCACATGAAAGGCGCTTTTCGGGCTGCCGGGCATCCTGAAATCGGCTTGCTATTTTTCATTCTTCCTGGAGCAGTCGCCAGCTTCTTTTCACAGCGTAGAGaa                                            <  1:1719684/148‑1 (MQ=255)
                                                                                                         tgctgTTTTGTCCTCTTCATTGTGGTGTGCCTTTTTCTGGCTACGCACATGAAAGGCGCTTTTCGGGCTGCCGGGCATCCTGAAATCGGCTTGCTATTTTTCATTCTTCCTGGAGCAGTCGCCAGCTTCTTTTCACAGCGTAGAGAAg                                           <  2:1235031/148‑1 (MQ=255)
                                                                                                          gctgTTTTGTCCTCTTCATTGTGGTGTGCCTTTTTCTGGCGACGCACATGAAAGGCGCTTTTCGGGCTGCCGGGCATCCTGAAATCGGCTTGCTATTTTTCATTCTTCCTGGAGCAGTCGCCAGCTTCTTTTCACAGCGTAGAGAAGt                                          <  1:3268244/148‑1 (MQ=255)
                                                                                                          gctgTTTTGTCCTCTTCATTGTGGTGTGCCTTTTTCTGGCGACGCACATGAAAGGCGCTTTTCGGGCTGCCGGGCATCCTGAAATCGGCTTGCTATTTTTCATTCTTCCTGGAGCAGTCGCCAGCTTCTTTTCACAGCGTAGAGAAGt                                          >  1:1027641/1‑148 (MQ=255)
                                                                                                           ctgTTTTGTCCTCTTCATTGTGGTGTGCCTTTTTCTGGCGACGCACATGAAAGGCGCTTTTCGGGCTGCCGGGCATCCTGAAATCGGCTTGCTATTTTTCATTCTTCCTGGAGCAGTCGCCAGCTTCTTTTCACAGCGTAGAGAAGTc                                         <  1:460288/148‑1 (MQ=255)
                                                                                                            tgTTTTGTCCTCTTCATTGTGGTGTGCCTTTTTCTGGCGACGCACATGAAAGGCGCTTTTCGGGCTGCCGGGCATCCTGAAATCGGCTTGCTATTTTTCATTCTTCCTGGAGCAGTCGCCAGCTTCTTTTCACAGCGTCGAGAAGTcc                                        >  2:1998014/1‑148 (MQ=255)
                                                                                                                   tCCTCTTCATTGTGGTGTGCCTTTTTCTGGCGACGCACATGAAAGGCGCTTTTCGGGCTGCCGGGCATCCTGAAATCGGCTTGCTATTTTTCATTCTTCCTGGAGCAGTCGCCAGCTTCTTTTCACAGCGTAGAGAAGTCCTGAAAcc                                 >  2:2253869/1‑148 (MQ=255)
                                                                                                                    ccTCTTCATTGTGGTGTGCCTTTTTCTGGCGACGCACATGAAAGGCGCTTTTCGGGCTGCCGGGCATCCTGAAATCGGCTTGCTATTTTTCATTCTTCCTGGAGCAGTCGCCAGCTTCTTTTCACAGCGTAGAGAAGTCCTGAAACct                                <  1:1045933/148‑1 (MQ=255)
                                                                                                                     ctctTCATTGTGGTGTGCCTTTTTCTGGCTACGCACATGAAAGGCGCTTTTCGGGCTGCCGGGCATCCTGAAATCGGCTTGCTATTTTTCATTCTTCCTGGAGCAGTCGCCAGCTTCTTTTCACAGCGTAGAGAAGTCCTGAAACctc                               >  1:2300394/1‑148 (MQ=255)
                                                                                                                       ctTCATTGTGGTGTGCCTTTTTCTGGCGACGCACATGAAAGGCGCTTTTCGGGCTGCCGGGCATCCTGAAATCGGCTTGCTATTTTTCATTCTTCCTGGAGCAGTCGCCAGCTTCTTTTCACAGCGTAGAGAAGTCCTGAAACCTCTg                             >  2:758132/1‑148 (MQ=255)
                                                                                                                       ctTCATTGTGGTGTGCCTTTTTCTGGCGACGCACATGAAAGGCGCTTTTCGGGCTGCCGGGCATCCTGAAATCGGCTTGCTATTTTTCATTCTTCCTGGAGCAGTCGCCAGCTTCTTTTCACAGCGTAGAGAAGTCCTGAAACCTCTg                             >  1:738298/1‑148 (MQ=255)
                                                                                                                       ctTCATTGTGGTGTGCCTTTTTCTGGCGACGCACATGAAAGGCGCTTTTCGGGCTGCCGGGCATCCTGAAATCGGCTTGCTATTTTTCATTCTTCCTGGAGCAGTCGCCAGCTTCTTTTCACAGCGTAGAGAAGTCCTGAAACCTCTg                             >  1:738356/1‑148 (MQ=255)
                                                                                                                        ttCATTGTGGTGTGCCTTTTTCTGGCGACGCACATGAAAGGCGCTTTTCGGGCTGCCGGGCATCCTGAAATCGGCTTGCTATTTTTCATTCTTCCTGGAGCAGTCGCCAGCTTCTTTTCACAGCGTAGAGAAGTCCTGAAACCTCTGt                            <  2:2160563/148‑1 (MQ=255)
                                                                                                                          cATTGTGGTGTGCCTTTTTCTGGCTACGCACATGAAAGGCGCTTTTCGGGCTGCCGGGCATCCTGAAATCGGCTTGCTATTTTTCATTCTTCCTGGAGCAGTCGCCAGCTTCTTTTCACAGCGTAGAGAAGTCCTGAAACCTCTGttt                          >  2:1628344/1‑148 (MQ=255)
                                                                                                                          cATTGTGGTGTGCCTTTTTCTGGCTACGCACATGAAAGGCGCTTTTCGGGCTGCCGGGCATCCTGAAATCGGCTTGCTATTTTTCATTCTTCCTGGAGCAGTCGCCAGCTTCTTTTCACAGCGTAGAGAAGTCCTGAAACCTCTGttt                          >  1:2566772/1‑148 (MQ=255)
                                                                                                                          cATTGTGGTGTGCCTTTTTCTGGCTACGCACATGAAAGGCGCTTTTCGGGCTGCCGGGCATCCTGAAATCGGCTTGCTATTTTTCATTCTTCCTGGAGCAGTCGCCAGCTTCTTTTCACAGCGTAGAGAAGTCCTGAAACCTCTGttt                          >  2:3265969/1‑148 (MQ=255)
                                                                                                                          cATTGTGGTGTGCCTTTTTCTGGCGACGCACATGAAAGGCGCTTTTCGGGCTGCCGGGCATCCTGAAATCGGCTTGCTATTTTTCATTCTTCCTGGAGCAGTCGCCAGCTTCTTTTCACAGCGTAGAGAAGTCCTGAAACCTCTGttt                          <  2:1361529/148‑1 (MQ=255)
                                                                                                                            ttGTGGTGTGCCTTTTTCTGGCGACGCACATGAAAGGCGCTTTTCGGGCTGCCGGGCATCCTGAAATCGGCTTGCTATTTTTCATTCTTCCTGGAGCAGTCGCCAGCTTCTTTTCACAGCGTAGAGAAGTCCTGAAACCTCTGTTTgg                        >  2:1722314/1‑148 (MQ=255)
                                                                                                                            ttGTGGTGTGCCTTTTTCTGGCGACGCACATGAAAGGCGCTTTTCGGGCTGCCGGGCATCCTGAAATCGGCTTGCTATTTTTCATTCTTCCTGGAGCAGTCGCCAGCTTCTTTTCACAGCGTAGAGAAGTCCTGAAACCTCTGTTTgg                        >  2:3173600/1‑148 (MQ=255)
                                                                                                                             tgtgGTGTGCCTTTTTCTGGCGACGCACATGAAAGGCGCTTTTCGGGCTGCCGGGCATCCTGAAATCGGCTTGCTATTTTTCATTCTTCCTGGAGCAGTCGCCAGCTTCTTTTCACAGCGTAGAGAAGTCCTGAAACCTCTGTTTGgc                       >  1:3096759/1‑148 (MQ=255)
                                                                                                                             tgtgGTGTGCCTTTTTCTGGCGACGCACATGAAAGGCGCTTTTCGGGCTGCCGGGCATCCTGAAATCGGCTTGCTATTTTTCATTCTTCCTGGAGCAGTCGCCAGCTTCTTTTCACAGCGTAGAGAAGTCCTGAAACCTCTGTTTGgc                       >  2:274967/1‑148 (MQ=255)
                                                                                                                                ggtgTGCCTTTTTCTGGCGACGCACATGAAAGGCGCTTTTCGGGCTGCCGGGCATCCTGAAATCGGCTTGCTATTTTTCATTCTTCCTGGAGCAGTCGCCAGCTTCTTTTCACAGCGTAGAGAAGTCCTGAAACCTCTGTTTGGCGCa                    <  2:1805056/148‑1 (MQ=255)
                                                                                                                                ggtgTGCCTTTTTCTGGCGACGCACATGAAAGGCGCTTTTCGGGCTGCCGGGCATCCTGAAATCGGCTTGCTATTTTTCATTCTTCCTGGAGCAGTCGCCAGCTTCTTTTCACAGCGTAGAGAAGTCCTGAAACCTCTGTTTGGCGCa                    <  2:3202001/148‑1 (MQ=255)
                                                                                                                                ggtgTGCCTTTTTCTGGCGACGCACATGAAAGGCGCTTTTCGGGCTGCCGGGCATCCTGAAATCGGCTTGCTATTTTTCATTCTTCCTGGAGCAGTCGCCAGCTTCTTTTCACAGCGTAGAGAAGTCATGAAACCTCTGTTTGGCGCa                    <  1:684065/148‑1 (MQ=255)
                                                                                                                                 gtgtGCCTTTTTCTGGCGACGCACATGAAAGGCGCTTTTCGGGCTGCCGGGCATCCTGAAATCGGCTTGCTATTTTTCATTCTTCCTGGAGCAGTCGCCAGCTTCTTTTCACAGCGTAGAGAAGTCCTGAAACCTCTGTTTGGCGCaa                   <  1:1146652/148‑1 (MQ=255)
                                                                                                                                   gtgCCTTTTTCTGGCGACGCACATGAAAGGCGCTTTTCGGGCTGCCGGGCATCCTGAAATCGGCTTGCTATTTTTCATTCTTCCTGGAGCAGTCGCCAGCTTCTTTTCACAGCGTAGAGAAGTCCTGAAACCTCTGTTTGGCGCAATg                 <  1:1803183/148‑1 (MQ=255)
                                                                                                                                    tgCCTTTTTCTGGCGACGCACATGAAAGGCGCTTTTCGGGCTGCCGGGCATCCTGAAATCGGCTTGCTATTTTTCATTCTTCCTGGAGCAGTCGCCAGCTTCTTTTCACAGCGTAGAGAAGTCCTGAAACCTCTGTTTGGCGCAATga                >  1:1574253/1‑147 (MQ=255)
                                                                                                                                    tgCCTTTTTCTGGCGACGCACATGAAAGGCGCTTTTCGGGCTGCCGGGCATCCTGAAATCGGCTTGCTATTTTTCATTCTTCCTGGAGCAGTCGCCAGCTTCTTTTCACAGCGTAGAGAAGTCCTGAAACCTCTGTTTGGCGCAATGc                >  1:1574280/1‑148 (MQ=255)
                                                                                                                                     gCCTTTTTCTGGCGACGCACATGAAAGGCGCTTTTCGGGCTGCCGGGCATCCTGAAATCGGCTTGCTATTTTTCATTCTTCCTGGAGCAGTCGCCAGCTTCTTTTCACAGCGTAGAGAAGTCCTGAAACCTCTGTTTGGCGCAATGCt               <  2:435368/148‑1 (MQ=255)
                                                                                                                                     gCCTTTTTCTGGCGACGCACATGAAAGGCGCTTTTCGGGCTGCCGGGCATCCTGAAATCGGCTTGCTATTTTTCATTCTTCCTGGAGCAGTCGCCAGCTTCTTTTCACAGCGTAGAGAAGTCCTGAAACCTCTGTTTGGCGCAATGCt               >  2:7834/1‑148 (MQ=255)
                                                                                                                                     gCCTTTTTCTGGCGACGCACATGAAAGGCGCTTTTCGGGCTGCCGGGCATCCTGAAATCGGCTTGCTATTTTTCATTCTTCCTGGAGCAGTCGCCAGCTTCTTTTCACAGCGTAGAGAAGTCCTGAAACCTCTGTTTGGCGCAATGCt               >  2:2299175/1‑148 (MQ=255)
                                                                                                                                      ccTTTTTCTGGCGACGCACATGAAAGGCGCTTTTCGGGCTGCCGGGCATCCTGAAATCGGCTTGCTATTTTTCATTCTTCCTGGAGCAGTCGCCAGCTTCTTTTCACAGCGTAGAGAAGTCCTGAAACCTCTGTTTGGCGCAATGCTg              <  1:590659/148‑1 (MQ=255)
                                                                                                                                       cTTTTTCTGGCGACGCACATGAAAGGCGCTTTTCGGGCTGCCGGGCATCCTGAAATCGGCTTGCTATTTTTCATTCTTCCTGGAGCAGTCGCCAGCTTCTTTTCACAGCGTAGAGAAGTCCTGAAACCTCTGTTTGGCGCAATGCTgg             >  2:1609869/1‑148 (MQ=255)
                                                                                                                                        tttttCTGGCGACGCACATGAAAGGCGCTTTTCGGGCTGCCGGGCATCCTGAAATCGGCTTGCTATTTTTCATTCTTCCTGGAGCAGTCGCCAGCTTCTTTTCACAGCGTAGAGAAGTCCTGAAACCTCTGTTTGGCGCAATGCTggc            >  1:911082/1‑148 (MQ=255)
                                                                                                                                        tttttCTGGCGACGCACATGAAAGGCGCTTTTCGGGCTGCCGGGCATCCTGAAATCGGCTTGCTATTTTTCATTCTTCCTGGAGCAGTCGCCAGCTTCTTTTCACAGCGTAGAGAAGTCCTGAAACCTCTGTTTGGCGCAATGCTggc            >  1:2565062/1‑148 (MQ=255)
                                                                                                                                         ttttCTGGCGACGCACATGAAAGGCGCTTTTCGGGCTGCCGGGCATCCTGAAATCGGCTTGCTATTTTTCATTCTTCCTGGAGCAGTCGCCAGCTTCTTTTCACAGCGTAGAGACGTCCTGAAACCTCTGTTTGGCGCAATGCTggcg           >  1:2870533/1‑148 (MQ=255)
                                                                                                                                         ttttCTGGCGACGCACATGAAAGGCGCTTTTCGGGCTGCCGGGCATCCTGAAATCGGCTTGCTATTTTTCATTCTTCCTGGAGCAGTCGCCAGCTTCTTTTCACAGCGTAGAGAAGTCCTGAAACCTCTGTTTGGCGCAATGCTggcg           >  1:1654321/1‑148 (MQ=255)
                                                                                                                                           ttCTGGCGACGCACATGAAAGGCGCTTTTCGGGCTGCCGGGCATCCTGAAATCGGCTTGCTATTTTTCATTCTTCCTGGAGCAGTCGCCAGCTTCTTTTCACAGCGTAGAGAAGTCCTGAAACCTCTGTTTGGCGCAATGCTggcggc         >  2:1093374/1‑148 (MQ=255)
                                                                                                                                           ttCTGGCGACGCACATGAAAGGCGCTTTTCGGGCTGCCGGGCATCCTGAAATCGGCTTGCTATTTTTCATTCTTCCTGGAGCAGTCGCCAGCTTCTTTTCACAGCGTAGAGAAGTCCTGAAACCTCTGTTTGGCGCAATGCTggcggc         >  2:1093355/1‑148 (MQ=255)
                                                                                                                                                 cGACGCACATGAAAGGCGCTTTTCGGGCTGCCGGGCATCCTGAAATCGGCTTGCTATTTTTCATTCTTCCTGGAGCAGTCGCCAGCTTCTTTTCACAGCGTAGAGAAGTCCTGAAACCTCTGTTTGGCGCAATGCTGGCGGCACCCTg   >  2:1795935/1‑148 (MQ=255)
                                                                                                                                                  gACGCACATGAAAGGCGCTTTTCGGGCTGCCGGGCATCCTGAAATCGGCTTGCTATTTTTCATTCTTCCTGGAGCAGTCGCCAGCTTCTTTTCACAGCGTAGAGAAGTCCTGAAACCTCTGTTTGGCGCAATGCTGGCGGCACCCTGt  <  1:909246/148‑1 (MQ=255)
                                                                                                                                                  gACGCACATGAAAGGCGCTTTTCGGGCTGCCGGGCATCCTGAAATCGGCTTGCTATTTTTCATTCTTCCTGGAGCAGTCGCCAGCTTCTTTTCACAGCGTAGAGAAGTCCTGAAACCTCTGTTTGGCGCAATGCTGGCGGCACCCTGt  <  1:1531089/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
AAGCGGCTTTATCTTACTGAACGCCACACTGGCGAAAAATGTGTTCGATAGACGCAGTGTCAGGAGGAACGAGTGAAACATAAACAACGTTGGGCGGGGGCAATCTGCTGTTTTGTCCTCTTCATTGTGGTGTGCCTTTTTCTGGCGACGCACATGAAAGGCGCTTTTCGGGCTGCCGGGCATCCTGAAATCGGCTTGCTATTTTTCATTCTTCCTGGAGCAGTCGCCAGCTTCTTTTCACAGCGTAGAGAAGTCCTGAAACCTCTGTTTGGCGCAATGCTGGCGGCACCCTGT  >  CP009273/885473‑885766

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: