Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,860,159 C→A 100% P33P (CCG→CCT ibpB ← heat shock chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,860,1590CA96.7% 1034.8 / 24.8 334P33P (CCG→CCTibpBheat shock chaperone
Reads supporting (aligned to +/- strand):  ref base C (10/1);  new base A (136/186);  total (147/187)
Fisher's exact test for biased strand distribution p-value = 1.49e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.96e-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 338 total aligned reads displayed.

TTTTTCTCTTCTTTTGGCTGCTCCGGCGTGCCTTTTACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGCGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCATAGTCATTTCTCCTTCTAAGAAGCGAGTAAGTACCTGCAAATCCGAAGA  >  CP009273/3860013‑3860305
                                                                                                                                                  |                                                                                                                                                  
tttttCTCTTCTTTTGGCTGCTCCGGCGTGCCTTTTACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTATCAATGTTGTACGGAggg                                                                                                                                                 >  2:866174/1‑150 (MQ=255)
 ttttCTCTTCTTTTGGCTGCTCCGGCGTGCCTTTTACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGa                                                                                                                                                <  2:1719607/150‑1 (MQ=255)
  tttCTCTTCTTTTGGCTGCTCCGGCGTGCCTTTTACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGaa                                                                                                                                               <  2:2430482/150‑1 (MQ=255)
    tctcTTCTTTTGGCTGCTCCGGCGTGCCTTTTACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAAAGTTGTACGGAGGGAAGc                                                                                                                                             >  2:931428/1‑150 (MQ=255)
     ctctTCTTTTGGCTGCTCCGGCGTGCCTTTTACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGct                                                                                                                                            <  2:1171977/150‑1 (MQ=255)
      tcttctTTTGGCTGCTCCGGCGTGCCTTTTACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGctc                                                                                                                                           <  2:2229719/150‑1 (MQ=255)
      tcttctTTTGGCTGCTCCGGCGTGCCTTTTACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTATCAATGTTGTTCGGAGGGCAGctc                                                                                                                                           >  1:1370782/1‑150 (MQ=255)
        ttcttTTGGCTGCTCCGGCGTGCCTTTTACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTg                                                                                                                                         <  2:4941696/150‑1 (MQ=255)
             ttGGCTGCTCCGGCGTGCCTTTTACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTAAATGTTGTACGGAGGGAAGCTCTGGCgtt                                                                                                                                    >  2:58235/1‑150 (MQ=255)
              tGGCTGCTCCGGCGTGCCTTTTACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTc                                                                                                                                   <  2:1036844/150‑1 (MQ=255)
              tGGCTGCTCCGGCGTGCCTTTTACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTc                                                                                                                                   <  1:382777/150‑1 (MQ=255)
               ggCTGCTCCGGCGTGCCTTTTACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCa                                                                                                                                  <  1:2810369/150‑1 (MQ=255)
                   gctCCGGCGTGCCTTTTACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCgg                                                                                                                              <  2:2994771/150‑1 (MQ=255)
                   gctCCGGCGTGCCTTTTACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTAACCgg                                                                                                                              >  1:853447/1‑150 (MQ=255)
                     tCCGGCGTGCCTTTTACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTTTTCGGAGGGTAGCTCTGGCTTTCACCGGCg                                                                                                                            >  1:3030615/1‑150 (MQ=255)
                       cGGCGTGCCTTTTACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTATCAATGTTGTACGGAGGGAAGCACAGGCTTTCACCGGCGtg                                                                                                                          >  1:3481498/1‑149 (MQ=255)
                        ggCGTGCCTTTTACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGtgt                                                                                                                         >  2:2502280/1‑148 (MQ=255)
                          cGTGCCTTTTACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTg                                                                                                                       >  2:4629602/1‑150 (MQ=255)
                            tGCCTTTTACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCa                                                                                                                     >  2:710334/1‑150 (MQ=255)
                             gCCTTTTACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAg                                                                                                                    >  1:2519799/1‑150 (MQ=255)
                             gCCTTTTACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAg                                                                                                                    >  2:2857393/1‑150 (MQ=255)
                              ccTTTTACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGt                                                                                                                   <  2:4137540/150‑1 (MQ=255)
                                ttttACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCACGGGTAATGCGGTAGTGGGTATCGTCGCTTTTCACAAAGTAGTACGGAGGGAAGCTCTGGCGTTAACCGGCGTTTTGCAGAGc                                                                                                                 >  2:4753079/1‑150 (MQ=255)
                                ttttACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGc                                                                                                                 >  2:3817494/1‑150 (MQ=255)
                                 tttACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCg                                                                                                                <  2:738140/150‑1 (MQ=255)
                                  ttACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGt                                                                                                               >  1:3356818/1‑150 (MQ=255)
                                   tACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGtt                                                                                                              <  2:5060007/150‑1 (MQ=255)
                                    aCGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTg                                                                                                             <  2:18671/150‑1 (MQ=255)
                                     cGCTCAGGCGCGTACCGTCCAGTTGAATCTCTAAATCTTCCTGACGGACACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCCATGGTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTgg                                                                                                            <  2:3220536/150‑1 (MQ=255)
                                      gCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAAGGTTGTACGGAGGGAAGCGCTGGCTTTCACCGGCGTTTTGCAGTGCGGTGGc                                                                                                           >  1:3366972/1‑150 (MQ=255)
                                         cAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAg                                                                                                        >  1:1084888/1‑150 (MQ=255)
                                            gcgcGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGttt                                                                                                     <  1:315701/150‑1 (MQ=255)
                                             cgcgTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTg                                                                                                    >  2:3793602/1‑150 (MQ=255)
                                             cgcgTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTg                                                                                                    <  1:1882237/150‑1 (MQ=255)
                                                   ccTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCaaaa                                                                                              >  2:5099676/1‑150 (MQ=255)
                                                   ccTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCaaaa                                                                                              >  2:1712188/1‑150 (MQ=255)
                                                     ttCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGCGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCCGTTTGTCAAAAcc                                                                                            >  2:1103653/1‑150 (MQ=255)
                                                     ttCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCCGCGTTTTGCAGTGCGGTGGCCCGTTTGTCaaaaac                                                                                            >  1:1229155/1‑148 (MQ=255)
                                                       ccAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGa                                                                                          <  1:1990978/150‑1 (MQ=255)
                                                        cAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATTCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTTCGTTGGCCAGTTTGTCAAAACCGAt                                                                                         <  1:4035106/150‑1 (MQ=255)
                                                         atttGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATc                                                                                        >  2:2170287/3‑150 (MQ=255)
                                                         aGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATc                                                                                        >  1:2338872/1‑150 (MQ=255)
                                                         aGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCGTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATc                                                                                        >  2:884796/1‑150 (MQ=255)
                                                         aGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTAAATGTTGTACGGAGGGAAGCTCTTGCTTTAACCGGCGTTTTTCAGTGCGGTGGCCCGTTTGTCAAAACCCATc                                                                                        >  1:1609821/1‑150 (MQ=255)
                                                           ttGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCa                                                                                      <  2:1324627/150‑1 (MQ=255)
                                                              aaTCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTg                                                                                   >  2:3704834/1‑150 (MQ=255)
                                                                tctcTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGAc                                                                                 <  1:589090/150‑1 (MQ=255)
                                                                   ctAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCa                                                                              <  2:1818521/150‑1 (MQ=255)
                                                                     aaaTCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATc                                                                            <  1:3937752/150‑1 (MQ=255)
                                                                        tCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGt                                                                         <  1:3111273/150‑1 (MQ=255)
                                                                          ttCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTgg                                                                       <  2:4175261/150‑1 (MQ=255)
                                                                          ttCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTgg                                                                       <  2:1410607/150‑1 (MQ=255)
                                                                            ccggACGGACACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCATTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTgggg                                                                     <  2:2913825/147‑1 (MQ=255)
                                                                              tGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGAt                                                                   <  1:4228427/150‑1 (MQ=255)
                                                                               gACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATa                                                                  <  2:3295700/150‑1 (MQ=255)
                                                                                aCGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTTTCAAAACCGATCCATTGACGCATCAGTGGGGATaa                                                                 >  1:2728153/1‑150 (MQ=255)
                                                                                  ggAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAAt                                                               <  1:1323656/150‑1 (MQ=255)
                                                                                      aCCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGaa                                                           <  1:2408210/150‑1 (MQ=255)
                                                                                         tGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGtt                                                        >  2:1685535/1‑150 (MQ=255)
                                                                                         tGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGtt                                                        <  2:1469131/150‑1 (MQ=255)
                                                                                         tGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGtt                                                        >  2:5056901/1‑150 (MQ=255)
                                                                                         tGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCAGCAGTGGGGATAAATCGAAGtt                                                        >  1:4696784/1‑150 (MQ=255)
                                                                                          gCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTa                                                       >  2:2603752/1‑150 (MQ=255)
                                                                                            cAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACg                                                     <  2:3768216/150‑1 (MQ=255)
                                                                                                gcAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCATa                                                 <  1:460242/150‑1 (MQ=255)
                                                                                                 cAAGGGTAATGCGGTTGTGGTTATCGTCGCTTTTCTCAATGTTGTCCGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCATAg                                                <  2:2647717/150‑1 (MQ=255)
                                                                                                 cAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCATAg                                                <  2:3555258/150‑1 (MQ=255)
                                                                                                 cAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAAAGTTGTACGGAGGGAAGCTCTGGCTTTAACCGGCGTTTTTCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCCTTGAAGCATCAGTGGGGCTAAATCGAAGTTACGCagag                                                >  1:4441918/1‑150 (MQ=255)
                                                                                                   aGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCATAGTc                                              <  2:3574971/150‑1 (MQ=255)
                                                                                                      gTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCATAGTCAtt                                           <  1:4814560/150‑1 (MQ=255)
                                                                                                      gTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCATAGTCAtt                                           <  1:502453/150‑1 (MQ=255)
                                                                                                        aaTGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCATAGTCATTtc                                         <  1:3660544/150‑1 (MQ=255)
                                                                                                           gCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGTGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGGGGGGATAAATCGAAGTTACGCATAGTCATTTCTcc                                      >  2:5047118/1‑150 (MQ=255)
                                                                                                             ggTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCATAGTCATTTCTCCtt                                    <  2:4400563/150‑1 (MQ=255)
                                                                                                              gTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCATAGTCATTTCTCCTTc                                   <  2:3441486/150‑1 (MQ=255)
                                                                                                               tAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCATAGTCATTTCTCCTTCt                                  <  2:4834338/150‑1 (MQ=255)
                                                                                                               tAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCATAGTCATTTCTCCTTCt                                  <  1:4282565/150‑1 (MQ=255)
                                                                                                                  tGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCATAGTCATTTCTCCTTCTaag                               <  1:2614443/150‑1 (MQ=255)
                                                                                                                   ggTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCATAGTCATTTCTCCTTCTaaga                              >  2:2953388/1‑150 (MQ=255)
                                                                                                                   ggTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCATAGTCATTTCTCCTTCTaaga                              <  2:193459/150‑1 (MQ=255)
                                                                                                                          gtcgCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCATAGTCATTTCTCCTTCTAAGAAGCGagt                       <  1:932215/150‑1 (MQ=255)
                                                                                                                           tcgCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCATAGTCATTTCTCCTTCTAAGAAGCGagta                      <  1:1833522/150‑1 (MQ=255)
                                                                                                                            cgCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCATAGTCATTTCTCCTTCTAAGAAGCGagtaa                     <  1:4163331/150‑1 (MQ=255)
                                                                                                                            cgCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCATAGTCATTTCTCCTTCTAAGAAGCGagtaa                     <  1:2193722/150‑1 (MQ=255)
                                                                                                                            cgCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTGTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCATAGTCATTTCTCCTTCTAAGAAGCGagtaa                     >  1:1662321/1‑150 (MQ=255)
                                                                                                                             gCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCATAGTCATTTCTCCTTCTAAGAAGCGagtaag                    <  1:1059314/150‑1 (MQ=255)
                                                                                                                              cTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCATAGTCATTTCTCCTTCTAAGAAGCGagtaagt                   <  1:4946068/150‑1 (MQ=255)
                                                                                                                                 ttCTCAATGTTGTAAGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAATTTACGCATAGTCATTTCTCCTTCTAAGAAGCGAGTAAGTAcc                <  1:3712671/150‑1 (MQ=255)
                                                                                                                                   ctcAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCATAGTCATTTCTCCTTCTAAGAAGCGAGTAAGTACCTg              <  2:2147060/150‑1 (MQ=255)
                                                                                                                                    tcAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCATAGTCATTTCTCCTTCTAAGAAGCGAGTAAGTACCTGc             <  2:3481556/150‑1 (MQ=255)
                                                                                                                                     cAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCATAGTCATTTCTCCTTCTAAGAAGCGAGTAAGTACCTGCa            <  1:360418/150‑1 (MQ=255)
                                                                                                                                       aTGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCATAGTCATTTCTCCTTCTAAGAAGCGAGTAAGTACCTGCaaa          <  1:2924319/150‑1 (MQ=255)
                                                                                                                                        tgttgtACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCATAGTCATTTCTCCTTCTAAGAAGCGAGTAAGTACCTGCAAAt         <  1:2301662/150‑1 (MQ=255)
                                                                                                                                        tgttgtACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCATAGTCATTTCTCCTTCTAAGAAGCGAGTAAGTACCTGCAAAt         <  2:1579801/150‑1 (MQ=255)
                                                                                                                                            gtACGGCGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCATAGTCATTTCTCCTTCTAAGAAGCGAGTAAGTACCTGCAAATCCGa     >  1:3025764/1‑150 (MQ=255)
                                                                                                                                            gtACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCATAGTCATTTCTCCTTCTAAGAAGCGAGTAAGTACCTGCAAATCCGa     >  1:3544841/1‑150 (MQ=255)
                                                                                                                                            gtACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCATAGTCATTTCTCCTTCTAAGAAGCGAGTAAGTACCTGCAAATCCGa     >  2:3458310/1‑150 (MQ=255)
                                                                                                                                            gtACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCATAGTCATTTCTCCTTCTAAGAAGCGAGTAAGTACCTGCAAATCCGa     >  1:1470617/1‑150 (MQ=255)
                                                                                                                                               cggaggGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCATAGTCATTTCTCCTTCTAAGAAGCGAGTAAGTACCTGCAAATCCGAAGa  >  1:3749053/1‑150 (MQ=255)
                                                                                                                                               cggaggGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCATAGTCATTTCTCCTTCTAAGAAGCGAGTAAGTACCTGCAAATCCGAAGa  <  1:4702038/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
TTTTTCTCTTCTTTTGGCTGCTCCGGCGTGCCTTTTACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGCGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCATAGTCATTTCTCCTTCTAAGAAGCGAGTAAGTACCTGCAAATCCGAAGA  >  CP009273/3860013‑3860305

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: