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,5000AT97.0% 730.5 / 9.0 235S22S (TCA→TCTyliEputative membrane‑anchored cyclic‑di‑GMP phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base A (2/5);  new base T (110/118);  total (112/123)
Fisher's exact test for biased strand distribution p-value = 4.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.78e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 240 total aligned reads displayed.

ATTTGGTTAAAAAGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCATTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGTTCGGCCCTGACATTAAAGGA  >  CP009273/868355‑868643
                                                                                                                                                 |                                                                                                                                               
 tttGGTTAAAAAGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGATGATAATTTTATTTTTATTGGCAGCACTGTCtttt                                                                                                                                              >  1:2216862/1‑148 (MQ=255)
  ttGGTTAAAAAGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTa                                                                                                                                             >  2:798438/1‑148 (MQ=255)
   tGGTTAAAAAGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTAt                                                                                                                                            <  1:3472221/148‑1 (MQ=255)
   tGGTTAAAAAGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTAt                                                                                                                                            <  1:1453855/148‑1 (MQ=255)
     gTTAAAAAGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTg                                                                                                                                          >  1:1021766/1‑148 (MQ=255)
       gaacaaGCCGCAGCAAGGGACAATTTTGGCCGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAAGATGCTGTGTTTATAAGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGt                                                                                                                                        <  2:1556079/144‑1 (MQ=255)
             gCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCttttt                                                                                                                                  >  1:730729/1‑148 (MQ=255)
             gCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCttttt                                                                                                                                  <  2:113844/148‑1 (MQ=255)
               cGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCttttttt                                                                                                                                >  2:3306899/1‑148 (MQ=255)
               cGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCttttttt                                                                                                                                >  2:120431/1‑148 (MQ=255)
                      aGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCAtt                                                                                                                         >  2:3553792/1‑148 (MQ=255)
                       gggACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTa                                                                                                                        >  2:937845/1‑148 (MQ=255)
                           cAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTa                                                                                                                    >  2:1312971/1‑148 (MQ=255)
                           cAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTa                                                                                                                    >  2:1314400/1‑148 (MQ=255)
                           cAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTa                                                                                                                    <  1:2456095/148‑1 (MQ=255)
                               ttttGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCaa                                                                                                                >  1:567399/1‑148 (MQ=255)
                               ttttGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCaa                                                                                                                <  2:221828/148‑1 (MQ=255)
                                  tGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTg                                                                                                             <  1:2954585/148‑1 (MQ=255)
                                    cAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCATTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGt                                                                                                           <  2:1830168/148‑1 (MQ=255)
                                     aGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTa                                                                                                          <  1:974921/148‑1 (MQ=255)
                                        ggCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCg                                                                                                       <  1:3153427/148‑1 (MQ=255)
                                          cacaGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGa                                                                                                     <  2:1205287/148‑1 (MQ=255)
                                              gCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAgcgc                                                                                                 >  1:2916354/1‑148 (MQ=255)
                                                 ttCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGcgcgcg                                                                                              <  1:777822/148‑1 (MQ=255)
                                                 ttCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGcgcgcg                                                                                              >  1:1590425/1‑148 (MQ=255)
                                                      aTAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTaaa                                                                                         <  2:2516960/148‑1 (MQ=255)
                                                      aTAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTaaa                                                                                         <  1:1188968/148‑1 (MQ=255)
                                                         gTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACg                                                                                      >  2:893209/1‑148 (MQ=255)
                                                          ttATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGt                                                                                     >  1:1367846/1‑148 (MQ=255)
                                                            aTTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGc                                                                                   >  1:3213943/1‑148 (MQ=255)
                                                             ttttGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGcc                                                                                  >  1:81568/1‑148 (MQ=255)
                                                               ttGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGa                                                                                <  1:1709708/148‑1 (MQ=255)
                                                                tGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGAt                                                                               <  2:2901903/148‑1 (MQ=255)
                                                                tGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAAATATCAACTGGTATCGGAGCGCGCGGTCAAACGTGCCGAt                                                                               >  1:2460009/1‑148 (MQ=255)
                                                                 gTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATa                                                                              >  2:1476323/1‑148 (MQ=255)
                                                                  ttAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAg                                                                             <  1:1892411/148‑1 (MQ=255)
                                                                  ttAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAg                                                                             >  1:2889390/1‑148 (MQ=255)
                                                                     aaTGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCg                                                                          >  1:533002/1‑148 (MQ=255)
                                                                       tGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATATTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCt                                                                        >  2:2166534/1‑148 (MQ=255)
                                                                         tATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCttt                                                                      >  1:1717298/1‑148 (MQ=255)
                                                                          aTTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTg                                                                     >  2:1282010/1‑148 (MQ=255)
                                                                          aTTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTg                                                                     >  1:106238/1‑148 (MQ=255)
                                                                           ttAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGa                                                                    >  1:3076834/1‑148 (MQ=255)
                                                                            tAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGaa                                                                   >  2:2027157/1‑148 (MQ=255)
                                                                             aaCATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAAc                                                                  <  1:471733/148‑1 (MQ=255)
                                                                               cATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACtt                                                                >  2:1533832/1‑148 (MQ=255)
                                                                                 tGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTAt                                                              <  2:656016/148‑1 (MQ=255)
                                                                                  gCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTAATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGAAAGCCGCTTTGAACTAAtt                                                             >  1:338070/1‑148 (MQ=255)
                                                                                     gAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAg                                                          >  2:344726/1‑148 (MQ=255)
                                                                                      aGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGa                                                         >  2:837463/1‑148 (MQ=255)
                                                                                       gTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGaa                                                        >  1:2162500/1‑148 (MQ=255)
                                                                                         ttATACTAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCATTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGaaaa                                                      <  2:1885925/148‑1 (MQ=255)
                                                                                            taCGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACg                                                   >  2:2198387/1‑148 (MQ=255)
                                                                                            taCGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACg                                                   >  2:404395/1‑148 (MQ=255)
                                                                                              cGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGtt                                                 <  1:1840093/148‑1 (MQ=255)
                                                                                               gAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTg                                                <  1:1860889/148‑1 (MQ=255)
                                                                                                 aaagataaagatAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGc                                              <  1:1122807/148‑1 (MQ=255)
                                                                                                  aagataaagataAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCt                                             >  1:3144704/1‑148 (MQ=255)
                                                                                                   agataaagataaGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTa                                            <  2:3005513/148‑1 (MQ=255)
                                                                                                    gataaagataaGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTAt                                           <  2:3223268/148‑1 (MQ=255)
                                                                                                     ataaagataaGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTAtt                                          >  1:2994326/1‑148 (MQ=255)
                                                                                                      taaagataaGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTAttt                                         <  1:637937/148‑1 (MQ=255)
                                                                                                      taaagataaGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTAttt                                         <  2:2083485/148‑1 (MQ=255)
                                                                                                       aaagataaGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTc                                        >  1:2030982/1‑148 (MQ=255)
                                                                                                          gataaGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCttt                                     <  1:62634/148‑1 (MQ=255)
                                                                                                           ataaGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTa                                    >  1:2287969/1‑148 (MQ=255)
                                                                                                            taaGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTaa                                   <  2:3170443/148‑1 (MQ=255)
                                                                                                             aaGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTaaa                                  >  1:2367348/1‑148 (MQ=255)
                                                                                                             aaGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTaaa                                  >  1:1647950/1‑148 (MQ=255)
                                                                                                              aGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAg                                 <  2:2268151/148‑1 (MQ=255)
                                                                                                               ggCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGa                                >  1:1765994/1‑148 (MQ=255)
                                                                                                                gCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGat                               >  1:3540175/1‑148 (MQ=255)
                                                                                                                 cTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGata                              <  2:1021766/148‑1 (MQ=255)
                                                                                                                 cTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGata                              <  1:3306899/148‑1 (MQ=255)
                                                                                                                 cTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGata                              <  1:120431/148‑1 (MQ=255)
                                                                                                                  ttaTAATTTTATTTTTATTGGCATCACTTTATTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGatat                             <  2:2324992/146‑1 (MQ=255)
                                                                                                                     tAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGa                          <  1:2669802/148‑1 (MQ=255)
                                                                                                                       attttattttTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAAc                        >  2:2379116/1‑148 (MQ=255)
                                                                                                                        ttttatttttatTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACg                       >  2:2426132/1‑148 (MQ=255)
                                                                                                                         tttatttttattGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCTATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGt                      <  1:3112332/148‑1 (MQ=255)
                                                                                                                         tttatttttattGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGt                      <  2:1910495/148‑1 (MQ=255)
                                                                                                                          ttatttttattGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGtt                     <  1:2857273/148‑1 (MQ=255)
                                                                                                                          ttatttttattGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGtt                     <  1:1394989/148‑1 (MQ=255)
                                                                                                                          ttatttttattGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGtt                     >  2:3091871/1‑148 (MQ=255)
                                                                                                                          ttatttttattGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGtt                     <  2:1712125/148‑1 (MQ=255)
                                                                                                                            atttttattGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGTTCg                   >  2:777222/1‑148 (MQ=255)
                                                                                                                               tttattGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGTTCGGcc                <  1:174926/148‑1 (MQ=255)
                                                                                                                                ttattGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGTTCGGccc               <  2:533002/148‑1 (MQ=255)
                                                                                                                                         gcaCTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGTTCGGCCCTGACATTaa      <  1:1994777/148‑1 (MQ=255)
                                                                                                                                          caCTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGTTCGGCCCTGACATTaaa     <  2:1690970/148‑1 (MQ=255)
                                                                                                                                          caCTGTCTTTTATTGGTCTTTATTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGTTCGGCCCTGACATTaaa     <  2:2232590/148‑1 (MQ=255)
                                                                                                                                           aCTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGTTCGGCCCTGACATTAAAg    <  2:1845730/148‑1 (MQ=255)
                                                                                                                                           aCTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGTTCGGCCCTGACATTAAAg    <  2:567399/148‑1 (MQ=255)
                                                                                                                                            cTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTTGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGTTCGGCCCTGACATTAAAgg   >  1:2234389/1‑148 (MQ=255)
                                                                                                                                             tGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGTCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGTTCGGCCCTGACATTAAAGGa  >  2:607959/1‑148 (MQ=255)
                                                                                                                                             tGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGTCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGTTCGGCCCTGACATTAAAGGa  >  2:606356/1‑148 (MQ=255)
                                                                                                                                             tGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGTTCGGCCCTGACATTAAAGGa  <  1:1554632/148‑1 (MQ=255)
                                                                                                                                             tGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGTTCGGCCCTGACATTAAAGGa  <  2:1964287/148‑1 (MQ=255)
                                                                                                                                             tGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGTTCGGCCCTGACATTAAAGGa  <  2:3365403/148‑1 (MQ=255)
                                                                                                                                             tGTCTTTTATTGGTCTTTATTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGTTCGGCCCTGACATTAAAGGa  <  2:2682745/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
ATTTGGTTAAAAAGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCATTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGTTCGGCCCTGACATTAAAGGA  >  CP009273/868355‑868643

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: