Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,860,653 C→A 100% E44* (GAA→TAA)  ibpA ← heat shock chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,860,6530CA97.1% 630.8 / 8.9 204E44* (GAA→TAA) ibpAheat shock chaperone
Reads supporting (aligned to +/- strand):  ref base C (1/5);  new base A (102/96);  total (103/101)
Fisher's exact test for biased strand distribution p-value = 1.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 204 total aligned reads displayed.

TCAGCGATGCCCTGGTACAGATAGGTGCGCTCTTTTTGTTCGTCGGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTCGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCC  >  CP009273/3860508‑3860795
                                                                                                                                                 |                                                                                                                                              
tCAGCGATGCCCTGGTACAGATAGGTGCGCTCTTTTTGTTCGTCGGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGt                                                                                                                                              >  2:2209123/1‑148 (MQ=255)
  aGCGATGCCCTGGTACAGATAGGTGCGCTCTTTTTGTTCGTCGGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCt                                                                                                                                            >  2:2001045/1‑148 (MQ=255)
   gCGATGCCCTGGTACAGATAGGTGCGCTCTTTTTGGTCGTCGGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTa                                                                                                                                           <  1:1015425/148‑1 (MQ=255)
      aTGCCCTGGTACAGATAGGTGCGCTCTTTTTGTTCGTCGGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCa                                                                                                                                        >  2:529173/1‑148 (MQ=255)
        gCCCTGGTACAGATAGGTGCGCTCTTTTTGTTCGTCGGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGt                                                                                                                                      >  1:3281927/1‑148 (MQ=255)
         cccTGGTACAGATAGGTGCGCTCTTTTTGTTCGTCGGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGtt                                                                                                                                     >  1:2122281/1‑148 (MQ=255)
           cTGGTACAGATAGGTGCGCTCTTTTTGTTCGTCGGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCa                                                                                                                                   >  2:1279007/1‑148 (MQ=255)
            tGGTACAGATAGGTGCGCTCTTTTTGTTCGTCGGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCaa                                                                                                                                  >  1:765369/1‑148 (MQ=255)
             ggTACAGATAGGTGCGCTCTTTTTGTTCGTCGGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAAc                                                                                                                                 >  1:3527929/1‑148 (MQ=255)
             ggTACAGATAGGTGCGCTCTTTTTGTTCGTCGGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAAc                                                                                                                                 >  1:3527946/1‑148 (MQ=255)
               tACAGATAGGTGCGCTCTTTTTGTTCGTCGGCGTGAGCACCTTTCACCACCAGCAGATTCTCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTCGTCTACCAGTTCAACGt                                                                                                                               <  1:2217566/148‑1 (MQ=255)
                 cAGATAGGTGCGCTCTTTTTGTTCGTCGGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTta                                                                                                                             <  2:3541676/148‑1 (MQ=255)
                  aGATAGGTGCGCTCTTTTTGTTCGTCGGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTtat                                                                                                                            <  1:1023600/148‑1 (MQ=255)
                  aGATAGGTGCGCTCTTTTTGTTCGTCGGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTtat                                                                                                                            <  1:713810/148‑1 (MQ=255)
                  aGATAGGTGCGCTCTTTTTGTTCGTCGGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTtat                                                                                                                            >  2:299522/1‑148 (MQ=255)
                          gcgcTCTTTTTGTTCGTCGGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAggg                                                                                                                    >  2:2241304/1‑148 (MQ=255)
                          gcgcTCTTTTTGTTCGTCGGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAggg                                                                                                                    <  1:290132/148‑1 (MQ=255)
                             ctctTTTTGTTCGTCGGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTCGTCTACCAGTTCAACGTTATACGGAGGGTAg                                                                                                                 <  2:174978/148‑1 (MQ=255)
                               ctTTTTGTTCGTCGGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAgcc                                                                                                               >  2:2043815/1‑148 (MQ=255)
                               ctTTTTGTTCGTCGGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAgcc                                                                                                               >  2:2043851/1‑148 (MQ=255)
                               ctTTTTGTTCGTCGGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAgcc                                                                                                               >  2:3467691/1‑148 (MQ=255)
                               ctTTTTGTTCGTCGGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAgcc                                                                                                               >  1:469212/1‑148 (MQ=255)
                                tttttGTTCGTCGGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAgccg                                                                                                              <  1:3475265/148‑1 (MQ=255)
                                 ttttGTTCGTCGGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAgccgc                                                                                                             <  2:1192571/148‑1 (MQ=255)
                                      ttggtcgGCTTGAGCACCTTTCACCACAAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTa                                                                                                        <  1:946436/145‑1 (MQ=255)
                                            ggCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTActctgg                                                                                                  >  1:2614339/1‑148 (MQ=255)
                                            ggCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTActctgg                                                                                                  >  1:828398/1‑148 (MQ=255)
                                              cGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTActctggct                                                                                                >  2:3163152/1‑148 (MQ=255)
                                                tGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTActctggctct                                                                                              >  2:1181280/1‑148 (MQ=255)
                                                  aGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTActctggctctgg                                                                                            <  2:229008/148‑1 (MQ=255)
                                                   gCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGgt                                                                                           <  1:2709883/148‑1 (MQ=255)
                                                    cACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGgtt                                                                                          <  2:1766585/148‑1 (MQ=255)
                                                     aCCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTCGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGgttg                                                                                         <  2:890560/148‑1 (MQ=255)
                                                     aCCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGgttg                                                                                         <  2:1461975/148‑1 (MQ=255)
                                                         ttCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGtttt                                                                                     >  2:2554003/1‑148 (MQ=255)
                                                          tCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTc                                                                                    <  1:2043832/148‑1 (MQ=255)
                                                            accaccAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTa                                                                                  <  1:321422/148‑1 (MQ=255)
                                                             ccaccaGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTaa                                                                                 >  1:360841/1‑148 (MQ=255)
                                                              caccaGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAg                                                                                >  1:530534/1‑148 (MQ=255)
                                                                ccaGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTg                                                                              >  2:2666952/1‑148 (MQ=255)
                                                                 cagcagATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTgg                                                                             <  2:3376996/148‑1 (MQ=255)
                                                                 cagcagATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTgg                                                                             <  1:1329955/148‑1 (MQ=255)
                                                                 cagcagATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTgg                                                                             <  2:2180526/148‑1 (MQ=255)
                                                                  agcagATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGt                                                                            <  2:1498748/148‑1 (MQ=255)
                                                                   gcagATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGtt                                                                           >  1:3403556/1‑148 (MQ=255)
                                                                       aTTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTaaac                                                                       >  1:1555223/1‑148 (MQ=255)
                                                                        ttATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTaaaca                                                                      >  2:171430/1‑148 (MQ=255)
                                                                        ttATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTaaaca                                                                      >  2:1393083/1‑148 (MQ=255)
                                                                         tATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTaaacaa                                                                     >  2:36674/1‑148 (MQ=255)
                                                                          aTCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTaaacaaa                                                                    <  1:1143390/148‑1 (MQ=255)
                                                                           tCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTaaacaaac                                                                   >  1:1345835/1‑148 (MQ=255)
                                                                               gggCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTc                                                               >  1:2033878/1‑148 (MQ=255)
                                                                                  cGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTAAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCaaa                                                            >  2:3085749/1‑148 (MQ=255)
                                                                                    gTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATc                                                          >  2:2201156/1‑148 (MQ=255)
                                                                                       aTTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCaa                                                       <  2:1710979/148‑1 (MQ=255)
                                                                                       aTTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTGTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCaa                                                       >  2:2614444/1‑148 (MQ=255)
                                                                                          tCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAg                                                    <  2:363094/148‑1 (MQ=255)
                                                                                           ccAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGc                                                   >  1:2106084/1‑148 (MQ=255)
                                                                                            cAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCa                                                  <  2:3514455/148‑1 (MQ=255)
                                                                                              gTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGa                                                <  2:2716677/148‑1 (MQ=255)
                                                                                              gTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGa                                                <  1:3539905/148‑1 (MQ=255)
                                                                                                tCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAAc                                              >  2:2326160/1‑148 (MQ=255)
                                                                                                 cGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGCACg                                             >  2:2820199/1‑148 (MQ=255)
                                                                                                  gCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACgg                                            <  2:2955773/148‑1 (MQ=255)
                                                                                                   ctctCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGt                                           <  2:1537282/148‑1 (MQ=255)
                                                                                                        aGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGc                                      <  1:2554003/148‑1 (MQ=255)
                                                                                                            aaaCCAGCCACAGCGATAGCAATGCGGTAATGGTTTTCGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCgggg                                  >  1:1767352/1‑148 (MQ=255)
                                                                                                            aaaCCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCgggg                                  <  2:3184343/148‑1 (MQ=255)
                                                                                                             aaCCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGa                                 <  2:1858740/148‑1 (MQ=255)
                                                                                                              aCCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGTCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGAt                                <  1:523242/148‑1 (MQ=255)
                                                                                                              aCCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGAt                                <  1:3227819/148‑1 (MQ=255)
                                                                                                              aCCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGAt                                <  1:1455144/148‑1 (MQ=255)
                                                                                                               ccAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATa                               >  1:1189616/1‑148 (MQ=255)
                                                                                                               ccAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATa                               >  1:3131516/1‑148 (MQ=255)
                                                                                                                 aGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATaaa                             <  2:108765/148‑1 (MQ=255)
                                                                                                                  gCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAAt                            >  1:2365945/1‑148 (MQ=255)
                                                                                                                   ccACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATc                           <  2:758944/148‑1 (MQ=255)
                                                                                                                   ccACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATc                           >  1:1798318/1‑148 (MQ=255)
                                                                                                                    cacaGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCa                          <  2:1023540/148‑1 (MQ=255)
                                                                                                                        gCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGt                      >  1:1157535/1‑148 (MQ=255)
                                                                                                                        gCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGt                      <  2:2927391/148‑1 (MQ=255)
                                                                                                                         cGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGtt                     >  1:576222/1‑148 (MQ=255)
                                                                                                                         cGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGtt                     >  1:2009715/1‑148 (MQ=255)
                                                                                                                            tAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACg                  >  1:706578/1‑148 (MQ=255)
                                                                                                                            tAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACg                  >  2:352457/1‑148 (MQ=255)
                                                                                                                            tAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACg                  >  1:2292053/1‑148 (MQ=255)
                                                                                                                             aGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGc                 <  1:784777/148‑1 (MQ=255)
                                                                                                                              gCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCa                <  1:3356105/148‑1 (MQ=255)
                                                                                                                                aaTGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATa              <  2:3403556/148‑1 (MQ=255)
                                                                                                                                aaTGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATa              <  1:1981429/148‑1 (MQ=255)
                                                                                                                                aaTGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATa              <  2:866368/148‑1 (MQ=255)
                                                                                                                                  tGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAAt            <  1:2161350/148‑1 (MQ=255)
                                                                                                                                  tGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAAt            >  1:435671/1‑148 (MQ=255)
                                                                                                                                    cGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCa          <  2:84064/148‑1 (MQ=255)
                                                                                                                                    cGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCa          <  2:469212/148‑1 (MQ=255)
                                                                                                                                    cGGTAATGGTTTTAGTCTACCAGTTAAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCa          <  1:3085749/148‑1 (MQ=255)
                                                                                                                                     ggTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCaa         <  1:3269239/148‑1 (MQ=255)
                                                                                                                                      gTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAAt        <  1:1519260/148‑1 (MQ=255)
                                                                                                                                           ggTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTc   <  2:2409252/148‑1 (MQ=255)
                                                                                                                                            gTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTcc  <  2:2738805/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                              
TCAGCGATGCCCTGGTACAGATAGGTGCGCTCTTTTTGTTCGTCGGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTCGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCC  >  CP009273/3860508‑3860795

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: