Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 868,500 A→T 100% S22S (TCA→TCT yliE → putative membrane‑anchored cyclic‑di‑GMP phosphodiesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273868,5000AT98.5% 876.8 / 1.7 271S22S (TCA→TCTyliEputative membrane‑anchored cyclic‑di‑GMP phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base A (1/3);  new base T (158/109);  total (159/112)
Fisher's exact test for biased strand distribution p-value = 3.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.11e-01
Rejected as polymorphism: E-value score below prediction cutoff.
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 274 total aligned reads displayed.

AATTTGGTTAAAAAGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCATTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGTTCGGCCCTGACATTAAAGGAC  >  CP009273/868354‑868644
                                                                                                                                                  |                                                                                                                                                
 aTTTGGTTAAAAAGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCttt                                                                                                                                                >  1:1045256/1‑148 (MQ=255)
  tttGGTTAAAAAGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCtttt                                                                                                                                               <  1:2088635/148‑1 (MQ=255)
   ttGGTTAAAAAGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTa                                                                                                                                              <  1:1653136/148‑1 (MQ=255)
    tGGTTAAAAAGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTAt                                                                                                                                             <  2:2406974/148‑1 (MQ=255)
      gTTAAAAAGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTg                                                                                                                                           <  1:915224/148‑1 (MQ=255)
       ttAAAAAGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTgg                                                                                                                                          >  2:3635888/1‑148 (MQ=255)
       ttAAAAAGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTgg                                                                                                                                          <  2:1335985/148‑1 (MQ=255)
        tAAAAAGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGt                                                                                                                                         <  1:1366033/148‑1 (MQ=255)
          aaaaGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCt                                                                                                                                       <  2:3279078/148‑1 (MQ=255)
             aGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCtttt                                                                                                                                    <  1:2650970/148‑1 (MQ=255)
              gCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGATTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCttttt                                                                                                                                   <  2:938161/148‑1 (MQ=255)
              gCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCttttt                                                                                                                                   >  1:2741439/1‑148 (MQ=255)
                     caAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTtcatca                                                                                                                            <  2:762548/148‑1 (MQ=255)
                     caAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTtcatca                                                                                                                            >  1:1428082/1‑148 (MQ=255)
                       aGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCAtt                                                                                                                          >  2:1743988/1‑148 (MQ=255)
                       aGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCAtt                                                                                                                          >  2:1746707/1‑148 (MQ=255)
                        gggACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTa                                                                                                                         <  1:2576477/148‑1 (MQ=255)
                          gACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAAc                                                                                                                       >  1:1615059/1‑148 (MQ=255)
                          gACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAAc                                                                                                                       <  1:1813989/148‑1 (MQ=255)
                           aCAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGATTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACt                                                                                                                      <  2:2169694/148‑1 (MQ=255)
                           aCAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACt                                                                                                                      >  1:2498656/1‑148 (MQ=255)
                              aTTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATc                                                                                                                   >  1:660360/1‑148 (MQ=255)
                              aTTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATc                                                                                                                   >  1:2762812/1‑148 (MQ=255)
                                ttttGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCaa                                                                                                                 >  2:106057/1‑148 (MQ=255)
                                ttttGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCaa                                                                                                                 >  2:1167839/1‑148 (MQ=255)
                                ttttGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCaa                                                                                                                 >  2:2036392/1‑148 (MQ=255)
                                   tGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTg                                                                                                              <  2:2115936/148‑1 (MQ=255)
                                    gCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTgg                                                                                                             <  2:1505651/148‑1 (MQ=255)
                                     cAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGt                                                                                                            >  2:3758651/1‑148 (MQ=255)
                                      aGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTa                                                                                                           <  2:893634/148‑1 (MQ=255)
                                          gCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCgg                                                                                                       >  1:19543/1‑148 (MQ=255)
                                             caGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAgc                                                                                                    >  2:3453443/1‑148 (MQ=255)
                                              aGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAgcg                                                                                                   <  2:3170717/148‑1 (MQ=255)
                                              aGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAgcg                                                                                                   >  1:1486411/1‑148 (MQ=255)
                                              aGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAgcg                                                                                                   <  1:2891099/148‑1 (MQ=255)
                                                cGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAgcgcg                                                                                                 >  2:2263225/1‑148 (MQ=255)
                                                cGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAgcgcg                                                                                                 >  1:3204702/1‑148 (MQ=255)
                                                 gTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAgcgcgc                                                                                                <  1:106057/148‑1 (MQ=255)
                                                   tCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCgg                                                                                              >  1:1791503/1‑148 (MQ=255)
                                                        tAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTaaaa                                                                                         <  1:2667227/148‑1 (MQ=255)
                                                            tATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTg                                                                                     <  1:1348509/148‑1 (MQ=255)
                                                              ttttGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGcc                                                                                   >  2:1907360/1‑148 (MQ=255)
                                                               tttGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCg                                                                                  >  1:3163846/1‑148 (MQ=255)
                                                                ttGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCATCACTGTCTTTTATTGGTCTTTTTTTCATCATAAACTATCAAATGGTATCGGAGAGAGCGGTAAAAAGTGCCGa                                                                                 >  1:2465979/1‑148 (MQ=255)
                                                                   ttAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAg                                                                              >  2:2702565/1‑148 (MQ=255)
                                                                    tAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGc                                                                             >  2:3478922/1‑148 (MQ=255)
                                                                      aaTGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCg                                                                           >  2:2452067/1‑148 (MQ=255)
                                                                       aTGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGc                                                                          <  1:2177244/148‑1 (MQ=255)
                                                                       aTGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTGTATTGGTCTTTTTTTCATCATTAACTATAAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGc                                                                          >  1:3012029/1‑148 (MQ=255)
                                                                        tGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATAATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCt                                                                         >  1:100489/1‑148 (MQ=255)
                                                                         gTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCtt                                                                        <  2:726201/148‑1 (MQ=255)
                                                                            ttAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGa                                                                     >  1:2935735/1‑148 (MQ=255)
                                                                            ttAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGa                                                                     <  1:2563691/148‑1 (MQ=255)
                                                                            ttAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGa                                                                     >  2:2582419/1‑148 (MQ=255)
                                                                               aCATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACt                                                                  >  1:3057694/1‑148 (MQ=255)
                                                                                cATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACtt                                                                 <  2:2974952/148‑1 (MQ=255)
                                                                                cATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACtt                                                                 <  1:3453016/148‑1 (MQ=255)
                                                                                cATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACtt                                                                 >  2:3393710/1‑148 (MQ=255)
                                                                                   gCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTAtt                                                              >  2:1682247/1‑148 (MQ=255)
                                                                                    cTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTc                                                             <  1:115361/148‑1 (MQ=255)
                                                                                       aGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGa                                                          >  1:3767503/1‑148 (MQ=255)
                                                                                         tttATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGaaa                                                        >  2:561051/1‑148 (MQ=255)
                                                                                               cGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGtt                                                  <  2:277110/148‑1 (MQ=255)
                                                                                                gAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTg                                                 >  2:3515585/1‑148 (MQ=255)
                                                                                                 aaaaGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCATTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTgg                                                <  1:2031400/148‑1 (MQ=255)
                                                                                                   aagataaagataAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCt                                              >  2:2926569/1‑148 (MQ=255)
                                                                                                   aagataaagataAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCt                                              >  2:2072165/1‑148 (MQ=255)
                                                                                                   aagataaagataAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCATTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCt                                              <  1:310197/148‑1 (MQ=255)
                                                                                                         aagataaGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCAGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCt                                        <  1:1478305/148‑1 (MQ=255)
                                                                                                          agataaGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCtt                                       <  1:1746691/148‑1 (MQ=255)
                                                                                                          agataaGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCtt                                       >  1:966829/1‑148 (MQ=255)
                                                                                                           gataaGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCttt                                      >  1:582110/1‑148 (MQ=255)
                                                                                                            ataaGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGATCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTa                                     >  1:2014779/1‑148 (MQ=255)
                                                                                                            ataaGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTa                                     >  1:1890984/1‑148 (MQ=255)
                                                                                                              aaGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTaaa                                   >  2:242445/1‑148 (MQ=255)
                                                                                                              aaGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTaaa                                   >  2:1575608/1‑148 (MQ=255)
                                                                                                               aGGATGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAg                                  <  2:3778620/148‑1 (MQ=255)
                                                                                                                ggCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGa                                 >  1:196423/1‑148 (MQ=255)
                                                                                                                  cTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGata                               <  2:183905/148‑1 (MQ=255)
                                                                                                                    gATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATAtt                             >  2:407651/1‑148 (MQ=255)
                                                                                                                       aaTTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGaa                          >  2:1030764/1‑148 (MQ=255)
                                                                                                                         ttttatttttatTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACg                        >  1:2807427/1‑148 (MQ=255)
                                                                                                                         ttttatttttatTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACg                        >  2:2526160/1‑148 (MQ=255)
                                                                                                                         ttttatttttatTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACg                        >  1:2807485/1‑148 (MQ=255)
                                                                                                                         ttttatttttatTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACg                        >  2:3146467/1‑148 (MQ=255)
                                                                                                                          tttatttttattGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGt                       >  1:770525/1‑148 (MQ=255)
                                                                                                                          tttatttttattGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGt                       >  1:2230984/1‑148 (MQ=255)
                                                                                                                             atttttattGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCTCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGTTCg                    >  1:510650/1‑148 (MQ=255)
                                                                                                                                tttattGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGTTCGGcc                 >  2:1454400/1‑148 (MQ=255)
                                                                                                                                 ttattGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGTTCGGccc                >  1:481960/1‑148 (MQ=255)
                                                                                                                                    ttGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGTTCGGCCCTGa             <  1:3509114/148‑1 (MQ=255)
                                                                                                                                         agcaCTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGTTCGGCCCTGACATTa        <  1:3635888/148‑1 (MQ=255)
                                                                                                                                          gcaCTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGTTCGGCCCTGACATTaa       <  1:1234611/148‑1 (MQ=255)
                                                                                                                                           caCTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAATATATTTAAATTTCGGCCCTGACATTaaa      >  1:3276452/1‑148 (MQ=255)
                                                                                                                                           caCTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGTTCGGCCCTGACATTaaa      >  1:2394487/1‑148 (MQ=255)
                                                                                                                                            aCTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGTTCGGCCCTGACATTAAAg     <  2:3085175/148‑1 (MQ=255)
                                                                                                                                             cTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGTTCGGCCCTGACATTAAAgg    >  1:1151646/1‑148 (MQ=255)
                                                                                                                                              tGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGTTCGGCCCTGACATTAAAGGa   >  2:1962161/1‑148 (MQ=255)
                                                                                                                                              tGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGTTCGGCCCTGACATTAAAGGa   <  1:1697867/148‑1 (MQ=255)
                                                                                                                                               gTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGTTCGGCCCTGACATTAAAGGAc  <  2:2835565/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
AATTTGGTTAAAAAGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCATTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGTTCGGCCCTGACATTAAAGGAC  >  CP009273/868354‑868644

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: