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,5000AT99.5% 681.0 / ‑4.6 209S22S (TCA→TCTyliEputative membrane‑anchored cyclic‑di‑GMP phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base A (2/0);  new base T (104/103);  total (106/103)
Fisher's exact test for biased strand distribution p-value = 4.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
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 211 total aligned reads displayed.

ATTTGGTTAAAAAGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCATTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGTTCGGCCCTGACATTAAAGGAC  >  CP009273/868355‑868644
                                                                                                                                                 |                                                                                                                                                
aTTTGGTTAAAAAGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCttt                                                                                                                                                <  1:776206/148‑1 (MQ=255)
  ttGGTTAAAAAGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTa                                                                                                                                              <  2:1250223/148‑1 (MQ=255)
     gTTAAAAAGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTg                                                                                                                                           <  2:1121718/148‑1 (MQ=255)
      ttAAAAAGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTgg                                                                                                                                          <  1:445096/148‑1 (MQ=255)
        aaaaaGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTc                                                                                                                                        <  1:461745/148‑1 (MQ=255)
          aaaGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCtt                                                                                                                                      <  2:637371/148‑1 (MQ=255)
          aaaGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCtt                                                                                                                                      <  2:1388356/148‑1 (MQ=255)
           aaGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCttt                                                                                                                                     <  2:2044408/148‑1 (MQ=255)
            aGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCtttt                                                                                                                                    <  1:1154727/148‑1 (MQ=255)
            aGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCtttt                                                                                                                                    >  1:3510648/1‑148 (MQ=255)
            aGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCtttt                                                                                                                                    >  1:2850920/1‑148 (MQ=255)
             gCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCttttt                                                                                                                                   >  1:965360/1‑148 (MQ=255)
                 cagcaAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTtca                                                                                                                               <  1:339474/148‑1 (MQ=255)
                    caAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTtcatca                                                                                                                            >  2:3443595/1‑148 (MQ=255)
                    caAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTAATTTTATTGGCAGCACTGTCATTTATTGGTCTTTTTTtcatca                                                                                                                            >  2:2542686/1‑148 (MQ=255)
                       gggACAATTTTTGCAGCGGCACAGCGTTCAGATATTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTa                                                                                                                         <  2:2782476/148‑1 (MQ=255)
                       gggACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTa                                                                                                                         <  2:2711794/148‑1 (MQ=255)
                          aCAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACt                                                                                                                      <  1:761217/148‑1 (MQ=255)
                          aCAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACt                                                                                                                      <  2:2458679/148‑1 (MQ=255)
                          aCAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACt                                                                                                                      <  2:1767254/148‑1 (MQ=255)
                           cAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTa                                                                                                                     >  1:1934278/1‑148 (MQ=255)
                           cAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTa                                                                                                                     <  2:701586/148‑1 (MQ=255)
                              tttttGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCa                                                                                                                  >  2:1845930/1‑148 (MQ=255)
                              tttttGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCa                                                                                                                  >  2:1228711/1‑148 (MQ=255)
                               ttttGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTTAGTTTATACTAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCaa                                                                                                                 <  2:2429530/148‑1 (MQ=255)
                               ttttGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCaa                                                                                                                 <  1:609662/148‑1 (MQ=255)
                               ttttGCAGCGGCACAGCGTTCAGATAGTTATATTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCaa                                                                                                                 <  2:2568114/148‑1 (MQ=255)
                                      gCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTAt                                                                                                          <  1:3356163/148‑1 (MQ=255)
                                       cGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATc                                                                                                         >  1:1923821/1‑148 (MQ=255)
                                          cacaGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGa                                                                                                      >  2:197982/1‑148 (MQ=255)
                                           acaGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAg                                                                                                     <  1:2214100/148‑1 (MQ=255)
                                           acaGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAg                                                                                                     <  1:2975654/148‑1 (MQ=255)
                                            caGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAgc                                                                                                    >  2:3374914/1‑148 (MQ=255)
                                                      aTAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTaaa                                                                                          <  1:3443642/148‑1 (MQ=255)
                                                      aTAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTaaa                                                                                          >  2:1741657/1‑148 (MQ=255)
                                                         gTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACg                                                                                       <  1:1461748/148‑1 (MQ=255)
                                                           tATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTg                                                                                     >  1:67999/1‑148 (MQ=255)
                                                           tATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTg                                                                                     >  1:1007147/1‑148 (MQ=255)
                                                           tATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTGTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTCAAACGTg                                                                                     >  1:1496113/1‑148 (MQ=255)
                                                                     aaTGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCg                                                                           <  2:1703324/148‑1 (MQ=255)
                                                                      aTGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGc                                                                          >  1:1165513/1‑148 (MQ=255)
                                                                       tGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCt                                                                         >  2:599158/1‑148 (MQ=255)
                                                                        gTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCtt                                                                        <  2:3187260/148‑1 (MQ=255)
                                                                          aTTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTg                                                                      >  1:2015946/1‑148 (MQ=255)
                                                                          aTTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTg                                                                      <  2:2307354/148‑1 (MQ=255)
                                                                          aTTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTg                                                                      >  1:2760796/1‑148 (MQ=255)
                                                                           ttAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGa                                                                     <  1:1891010/148‑1 (MQ=255)
                                                                            tAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGaa                                                                    <  1:3061231/148‑1 (MQ=255)
                                                                            tAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGaa                                                                    >  2:1165049/1‑148 (MQ=255)
                                                                             aaCATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAAc                                                                   <  2:2601626/148‑1 (MQ=255)
                                                                              aCATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCACTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACt                                                                  >  1:2310220/1‑148 (MQ=255)
                                                                               cATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACtt                                                                 >  2:856853/1‑148 (MQ=255)
                                                                                   cTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTc                                                             <  1:1110105/148‑1 (MQ=255)
                                                                                        tttATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGaaa                                                        >  1:1971698/1‑148 (MQ=255)
                                                                                         ttATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGaaaa                                                       >  2:2024844/1‑148 (MQ=255)
                                                                                         ttATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGaaaa                                                       <  2:279351/148‑1 (MQ=255)
                                                                                         ttATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGaaaa                                                       <  2:900798/148‑1 (MQ=255)
                                                                                           ataCGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAAc                                                     <  1:2395855/148‑1 (MQ=255)
                                                                                             aCGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGt                                                   >  1:2006505/1‑148 (MQ=255)
                                                                                               gAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTg                                                 >  2:3000363/1‑148 (MQ=255)
                                                                                                 aaagataaagatAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGc                                               >  2:1546730/1‑148 (MQ=255)
                                                                                                   agataaagataaGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTa                                             <  1:2993285/148‑1 (MQ=255)
                                                                                                     ataaagataaGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTAtt                                           <  1:3366462/148‑1 (MQ=255)
                                                                                                      taaagataaGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTAttt                                          <  2:2401154/148‑1 (MQ=255)
                                                                                                       aaagataaGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTc                                         <  1:3421265/148‑1 (MQ=255)
                                                                                                         agataaGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCtt                                       >  1:1289631/1‑148 (MQ=255)
                                                                                                          gataaGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCttt                                      >  2:1785408/1‑148 (MQ=255)
                                                                                                           ataaGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTa                                     >  1:2469489/1‑148 (MQ=255)
                                                                                                            taaGTCTGTTTATTTTATTTTTATTGGCTGCACTGTCTTTTATTGTTCTTTTTTTCATCATTAACTATCAACTTGTATCGGAGCTCGCGTTAAAACGTGCCGATAGCAGCTTTGAACTTATTCAGAAAAACTTTGGCTATTTCTTTaa                                    <  2:1065798/148‑1 (MQ=255)
                                                                                                            taaGGCTGATATTTTTATTTTTATTGGCATCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTTTCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTaa                                    <  2:3307200/148‑1 (MQ=255)
                                                                                                            taaGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTaa                                    >  2:1023925/1‑148 (MQ=255)
                                                                                                              aGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTAACAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAAAGTTGGCTATTTCTTTaaat                                  >  1:601264/1‑147 (MQ=255)
                                                                                                               ggCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGa                                 >  1:1719742/1‑148 (MQ=255)
                                                                                                               ggCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGa                                 >  1:2091618/1‑148 (MQ=255)
                                                                                                               ggCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGa                                 >  1:2759965/1‑148 (MQ=255)
                                                                                                               ggCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGa                                 >  2:792035/1‑148 (MQ=255)
                                                                                                               ggCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGa                                 <  2:1960652/148‑1 (MQ=255)
                                                                                                                gCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGat                                >  1:1828036/1‑148 (MQ=255)
                                                                                                                  tGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGatat                              >  1:3059972/1‑148 (MQ=255)
                                                                                                                   gATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATAtt                             >  1:2417307/1‑148 (MQ=255)
                                                                                                                    aTAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTg                            >  2:3093580/1‑148 (MQ=255)
                                                                                                                    aTAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTg                            >  1:2808549/1‑148 (MQ=255)
                                                                                                                    aTAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTg                            <  1:2449352/148‑1 (MQ=255)
                                                                                                                     tAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTGTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGa                           >  1:612530/1‑148 (MQ=255)
                                                                                                                      aaTTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGaa                          >  1:2924654/1‑148 (MQ=255)
                                                                                                                      aaTTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGaa                          <  2:3244762/148‑1 (MQ=255)
                                                                                                                        ttttatttttatTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACg                        <  2:12846/148‑1 (MQ=255)
                                                                                                                          ttatttttattGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGtt                      >  1:236841/1‑148 (MQ=255)
                                                                                                                          ttatttttattGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGtt                      <  2:3500319/148‑1 (MQ=255)
                                                                                                                          ttatttttattGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGtt                      <  1:759305/148‑1 (MQ=255)
                                                                                                                           tatttttattGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGTTc                     <  1:1741657/148‑1 (MQ=255)
                                                                                                                              ttttattGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGTTCGGc                  >  2:119625/1‑148 (MQ=255)
                                                                                                                                    tGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGTTCGGCCCTGAc            >  2:837932/1‑148 (MQ=255)
                                                                                                                                     ggCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGTTCGGCCCTGACa           <  1:2191375/148‑1 (MQ=255)
                                                                                                                                      gcagcaCTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGTTCGGCCCTGACAt          <  1:3207587/148‑1 (MQ=255)
                                                                                                                                       cagcaCTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGTTCGGCCCTGACAtt         >  1:1290177/1‑148 (MQ=255)
                                                                                                                                          caCTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGTTCGGCCCTGACATTaaa      >  2:2943250/1‑148 (MQ=255)
                                                                                                                                           aCTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATAAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAAAGTTCGGCCCTGACATTAAAg     >  1:131369/1‑148 (MQ=255)
                                                                                                                                            cTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGTTCGGCCCTGACATTAAAgg    >  2:22018/1‑148 (MQ=255)
                                                                                                                                             tGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGTTCGGCCCTGACATTAAAGGa   >  1:1143059/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
ATTTGGTTAAAAAGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCATTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGTTCGGCCCTGACATTAAAGGAC  >  CP009273/868355‑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: