Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,637,458 G→T 25.8% S14Y (TCT→TAT)  quuQ ← Qin prophage; putative antitermination protein Q

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,637,4580GT25.8% 620.5 / 163.7 283S14Y (TCT→TAT) quuQQin prophage; putative antitermination protein Q
Reads supporting (aligned to +/- strand):  ref base G (106/103);  new base T (35/39);  total (141/142)
Fisher's exact test for biased strand distribution p-value = 6.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 292 total aligned reads displayed.

CAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACGAAATATCAGCTCCAGTTGGGAGCCATACTTATCTTCAAATGCCACTGTATCCG  >  CP009273/1637309‑1637606
                                                                                                                                                     |                                                                                                                                                    
cAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGgag                                                                                                                                                      <  1:1893504/150‑1 (MQ=255)
 aaaCTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGgaga                                                                                                                                                     >  1:3990896/1‑150 (MQ=255)
  aaCTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGatat                                                                                                                                                    >  2:247549/1‑150 (MQ=255)
  aaCTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGatat                                                                                                                                                    >  1:2163412/1‑150 (MQ=255)
     tGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGATAtttt                                                                                                                                                 >  1:3454930/1‑150 (MQ=255)
     tGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGAtttt                                                                                                                                                 >  1:1398669/1‑150 (MQ=255)
     tGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGAtttt                                                                                                                                                 >  2:2539763/1‑150 (MQ=255)
       aGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTgg                                                                                                                                               <  1:392076/150‑1 (MQ=255)
       aGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTgg                                                                                                                                               <  2:3675167/150‑1 (MQ=255)
         ggTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGATATTTTgggg                                                                                                                                             <  1:1639973/150‑1 (MQ=255)
         ggTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTgggg                                                                                                                                             <  1:4277920/150‑1 (MQ=255)
          gTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGa                                                                                                                                            >  1:827422/1‑150 (MQ=255)
           tCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATAACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAg                                                                                                                                           >  1:4251373/1‑150 (MQ=255)
            ccTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGt                                                                                                                                          >  2:701818/1‑150 (MQ=255)
             cTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTa                                                                                                                                         >  2:2311188/1‑150 (MQ=255)
              tGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGATATTTTGGGGAGTaa                                                                                                                                        >  1:4540715/1‑150 (MQ=255)
              tGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTaa                                                                                                                                        <  2:858962/150‑1 (MQ=255)
              tGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTGTGGAGATTTTGGGGAGTaa                                                                                                                                        >  2:3116642/1‑150 (MQ=255)
               gAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAat                                                                                                                                       >  1:4080404/1‑150 (MQ=255)
                   acCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATAttt                                                                                                                                   >  1:1959332/1‑150 (MQ=255)
                      aaCTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTgg                                                                                                                                >  2:3938348/1‑150 (MQ=255)
                       aCTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGt                                                                                                                               >  2:3047302/1‑150 (MQ=255)
                          tttGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAg                                                                                                                            >  2:3485345/1‑150 (MQ=255)
                            tGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGATATTTTGGGGAGTAATATTTTGGTAAGGc                                                                                                                          >  1:93100/1‑150 (MQ=255)
                              ccAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGATATTTTGGGGAGTAATATTTTGGTAAGGCtt                                                                                                                        <  2:4020725/150‑1 (MQ=255)
                                aGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCa                                                                                                                      >  1:744623/1‑150 (MQ=255)
                                   aaTAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGATATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAgg                                                                                                                   >  2:662318/1‑150 (MQ=255)
                                   aaTAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGATATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAgg                                                                                                                   <  2:1783875/150‑1 (MQ=255)
                                      aaGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTc                                                                                                                >  1:1051229/1‑150 (MQ=255)
                                       aGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGATATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCa                                                                                                               <  2:1221637/150‑1 (MQ=255)
                                       aGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCa                                                                                                               <  1:3317173/150‑1 (MQ=255)
                                          ccAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGc                                                                                                            <  1:3669388/150‑1 (MQ=255)
                                           cAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCt                                                                                                           <  1:244449/150‑1 (MQ=255)
                                           cAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCt                                                                                                           <  1:3431645/150‑1 (MQ=255)
                                            aaCCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGATATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTc                                                                                                          <  2:1355316/150‑1 (MQ=255)
                                             aCCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCg                                                                                                         >  2:4885015/1‑150 (MQ=255)
                                                cAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTtc                                                                                                      <  1:1013274/150‑1 (MQ=255)
                                                 aaGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTtct                                                                                                     >  1:5152299/1‑150 (MQ=255)
                                                    tggtgCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTCAGGCTTCAAGGTTCATGCTCGTTCTCCa                                                                                                  >  2:2697931/1‑150 (MQ=255)
                                                         cTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGATATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTAc                                                                                             >  2:198577/1‑150 (MQ=255)
                                                          tttCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACg                                                                                            >  2:3808919/1‑150 (MQ=255)
                                                            tCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGATATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGcc                                                                                          >  2:1017775/1‑150 (MQ=255)
                                                            tCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGcc                                                                                          <  2:3407570/150‑1 (MQ=255)
                                                             cGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCa                                                                                         <  2:638965/150‑1 (MQ=255)
                                                              gACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAg                                                                                        >  1:3627188/1‑150 (MQ=255)
                                                               aCTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGt                                                                                       >  1:257673/1‑150 (MQ=255)
                                                                 tGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTAc                                                                                     <  1:328848/150‑1 (MQ=255)
                                                                     ccATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGATTCATGCTCGTTCTCCACTTACGCCAGTACGCCa                                                                                 <  1:4953676/150‑1 (MQ=255)
                                                                      cATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGATATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCaa                                                                                >  2:4003764/1‑150 (MQ=255)
                                                                       aTCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAAt                                                                               <  1:4886945/150‑1 (MQ=255)
                                                                        tCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAAtt                                                                              <  2:1673663/150‑1 (MQ=255)
                                                                        tCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAAtt                                                                              >  2:4795148/1‑150 (MQ=255)
                                                                          cccTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGc                                                                            >  2:497221/1‑150 (MQ=255)
                                                                            cTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCa                                                                          <  1:2654684/150‑1 (MQ=255)
                                                                             tGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAg                                                                         >  1:5108366/1‑150 (MQ=255)
                                                                              gcgcTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGATATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAgc                                                                        <  1:4988156/150‑1 (MQ=255)
                                                                               cgcTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAgcg                                                                       <  1:2206286/150‑1 (MQ=255)
                                                                                 cTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAgcgcg                                                                     <  2:3200882/150‑1 (MQ=255)
                                                                                    ccATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGcga                                                                  >  1:2197591/1‑150 (MQ=255)
                                                                                     cATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGATATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGcgat                                                                 <  1:4822825/150‑1 (MQ=255)
                                                                                      aTTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGATATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGcgatc                                                                <  2:93100/150‑1 (MQ=255)
                                                                                      aTTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGATATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGcgatc                                                                <  2:4738761/150‑1 (MQ=255)
                                                                                       ttACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGcgatcg                                                               <  1:890654/150‑1 (MQ=255)
                                                                                        tACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGcgatcga                                                              <  1:3756614/150‑1 (MQ=255)
                                                                                         aCATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGcgatcgat                                                             >  2:2042280/1‑150 (MQ=255)
                                                                                          taTCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATa                                                            >  2:1657532/2‑150 (MQ=255)
                                                                                          cATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATa                                                            <  2:4251373/150‑1 (MQ=255)
                                                                                           aTCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATaa                                                           >  1:2319742/1‑150 (MQ=255)
                                                                                             ccGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGATATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATaaaa                                                         <  2:4310225/150‑1 (MQ=255)
                                                                                                  aTTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGATATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACGaaa                                                    <  2:2163412/150‑1 (MQ=255)
                                                                                                    tGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACGAAata                                                  <  1:3109781/150‑1 (MQ=255)
                                                                                                     gTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACGAAatat                                                 <  1:4410750/150‑1 (MQ=255)
                                                                                                      ttAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACGAAATATc                                                <  1:1329121/150‑1 (MQ=255)
                                                                                                       tAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGATATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACGAAATATCa                                               <  2:1518461/150‑1 (MQ=255)
                                                                                                        aaaCAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACGAAATATCAg                                              >  2:3304459/1‑150 (MQ=255)
                                                                                                        aaaCAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACGAAATATCAg                                              <  2:2396105/150‑1 (MQ=255)
                                                                                                           caacCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACGAAATATCAGCTc                                           <  2:2760314/150‑1 (MQ=255)
                                                                                                           caacCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACGAAATATCAGCTc                                           >  2:4188127/1‑150 (MQ=255)
                                                                                                            aacCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGATATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACGAAATATCAGCTcc                                          >  2:3336257/1‑150 (MQ=255)
                                                                                                            aacCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCCGCGCGCGATCGATAAAACGAAATATCAGCTcc                                          >  1:2007852/1‑150 (MQ=255)
                                                                                                                ggTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGATATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACGAAATATCAGCTCCAGtt                                      <  2:3454725/150‑1 (MQ=255)
                                                                                                                ggTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACGAAATATCAGCTCCAGtt                                      >  1:4589782/1‑150 (MQ=255)
                                                                                                                 gTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACGAAATATCAGCTCCAGTTg                                     <  1:4149116/150‑1 (MQ=255)
                                                                                                                    ccTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACGAAATATCAGCTCCAGTTGGGa                                  <  2:5064549/150‑1 (MQ=255)
                                                                                                                     cTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACGAAATATCAGCTCCAGTTGGGAg                                 >  1:2616036/1‑150 (MQ=255)
                                                                                                                      tGTCGCTGGCGCGGCCTCGCTCAAGTTGGGCGATGTTGGGGCGGCATATTGTGGTAAGGCTTCCAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCACTTGCGAGCGCGCGAGCGATAAAACGTATTAGCCGCTCCTGGTGGGAGc                                >  2:2984816/1‑150 (MQ=255)
                                                                                                                        tCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACGAAATATCAGCTCCAGTTGGGAGCCa                              <  2:2850658/150‑1 (MQ=255)
                                                                                                                          gCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACGAAATATCAGCTCCAGTTGGGAGACATa                            >  2:3588431/1‑150 (MQ=255)
                                                                                                                            tGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACGAAATATCAGCTCCAGTTGGGAGCCATACt                          >  2:4143796/1‑150 (MQ=255)
                                                                                                                              gAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACGAAATATCAGCTCCAGTTGGGAGCCATACTTa                        <  1:4791588/150‑1 (MQ=255)
                                                                                                                                  gTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACGAAATATCAGCTCCAGTTGGGAGCCCTACTTATCtt                    >  2:1433560/1‑150 (MQ=255)
                                                                                                                                         cTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACGAAATATCAGCTCCAGTTGGGAGCCATACTTATCTTCAAATGc             <  1:1744027/150‑1 (MQ=255)
                                                                                                                                             aTTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACGAAATATCAGCTCCAGTTGGGAGCCATACTTATCTTCAAATGCCACt         <  1:1548805/150‑1 (MQ=255)
                                                                                                                                             aTTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACGAAATATCAGCTCCAGTTGGGAGCCATACTTATCTTCAAATGCCACt         <  2:2321654/150‑1 (MQ=255)
                                                                                                                                               tttGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACGAAATATCAGCTCCAGTTGGGAGCCATACTTATCTTCAAATGCCACTGt       >  2:111112/1‑150 (MQ=255)
                                                                                                                                               tttGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACGAAATATCAGCTCCAGTTGGGAGCCATACTTATCTTCAAATGCCACTGt       >  1:3374521/1‑150 (MQ=255)
                                                                                                                                                 tGGATATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACGAAATATCAGCTCCAGTTGGGAGCCATACTTATCTTCAAATGCCACTGTAt     <  2:1897414/150‑1 (MQ=255)
                                                                                                                                                 tGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACGAAATATCAGCTCCAGTTGGGAGCCATACTTATCTTCAAATGCCACTGTAt     >  2:3883190/1‑150 (MQ=255)
                                                                                                                                                  ggATATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACGAAATATCAGCTCCAGTTGGGAGCCATACTTATCTTCAAATGCCACTGTATc    >  1:4054228/1‑150 (MQ=255)
                                                                                                                                                   gagaTTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCCGCGCGCGATCGATAAAACGAAATATCAGCTCCAGTTGGGAGCCATACTTATCTTCAAATGCCACTGTATcc   >  1:2086479/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACGAAATATCAGCTCCAGTTGGGAGCCATACTTATCTTCAAATGCCACTGTATCCG  >  CP009273/1637309‑1637606

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: