Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,637,458 G→T 11.7% 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,4580GT11.7% 1051.1 / 68.3 318S14Y (TCT→TAT) quuQQin prophage; putative antitermination protein Q
Reads supporting (aligned to +/- strand):  ref base G (132/148);  new base T (18/20);  total (150/168)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 331 total aligned reads displayed.

CAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACGAAATATCAGCTCCAGTTGGGAGCCATACTTATCTTCAAATGCCACTGTATCCGT  >  CP009273/1637309‑1637607
                                                                                                                                                     |                                                                                                                                                     
     tGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGAtttt                                                                                                                                                  >  2:5031057/1‑150 (MQ=255)
       aGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTgg                                                                                                                                                <  1:1837419/150‑1 (MQ=255)
          gTCCTGAACACCAACTTTTGCCAGAAATAATGCCAACCCAAGTGGTGCTTTCGACTGCACCAACCCCTGCGATGCCATTACATCCGTAATTGTTAAACAACAGGAGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGa                                                                                                                                             >  2:5230351/1‑150 (MQ=255)
          gTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGa                                                                                                                                             >  2:1426271/1‑150 (MQ=255)
              tGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTaa                                                                                                                                         >  1:2223172/1‑150 (MQ=255)
               gAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAat                                                                                                                                        <  2:46112/150‑1 (MQ=255)
                  cacCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATAtt                                                                                                                                     <  1:3188063/150‑1 (MQ=255)
                     cAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTg                                                                                                                                  <  1:673006/150‑1 (MQ=255)
                      aaCTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTgg                                                                                                                                 <  2:4844315/150‑1 (MQ=255)
                        cTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTa                                                                                                                               >  2:783725/1‑150 (MQ=255)
                        cTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTa                                                                                                                               >  1:5021060/1‑150 (MQ=255)
                          tttGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGATATTTTGGGGAGTAATATTTTGGTAAg                                                                                                                             <  1:3500214/150‑1 (MQ=255)
                          tttGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGATATTTTGGGGAGTAATATTTTGGTAAg                                                                                                                             <  1:4187149/150‑1 (MQ=255)
                          tttGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAg                                                                                                                             <  1:3115078/150‑1 (MQ=255)
                             gCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGATATTTTGGGGAGTAATATTTTGGTAAGGCt                                                                                                                          <  2:4291885/150‑1 (MQ=255)
                             gCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGATATTTTGGGGAGTAATATTTTGGTAAGGCt                                                                                                                          <  2:1307978/150‑1 (MQ=255)
                             gCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCt                                                                                                                          >  2:2571676/1‑150 (MQ=255)
                              ccAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCtt                                                                                                                         <  1:3228615/150‑1 (MQ=255)
                               cAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGATATTTTGGGGAGTAATATTTTGGTAAGGCTTc                                                                                                                        <  2:1395735/150‑1 (MQ=255)
                                aGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCa                                                                                                                       <  1:3470726/150‑1 (MQ=255)
                                  aaaTAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAg                                                                                                                     >  1:1505376/1‑150 (MQ=255)
                                  aaaTAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAg                                                                                                                     >  1:2882719/1‑150 (MQ=255)
                                   aaTAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAgg                                                                                                                    <  1:2280446/150‑1 (MQ=255)
                                        ggCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCAt                                                                                                               <  2:1143377/150‑1 (MQ=255)
                                         gCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCAtt                                                                                                              >  2:1639724/1‑149 (MQ=255)
                                           cAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCt                                                                                                            >  2:848/1‑150 (MQ=255)
                                            aaCCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTc                                                                                                           >  2:4698664/1‑150 (MQ=255)
                                            aaCCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTc                                                                                                           <  1:4083520/150‑1 (MQ=255)
                                             aCCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCg                                                                                                          <  2:2384647/150‑1 (MQ=255)
                                                  aGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTtctc                                                                                                     >  2:3526009/1‑150 (MQ=255)
                                                    tggtgCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCa                                                                                                   <  1:2899626/150‑1 (MQ=255)
                                                     ggtgCTTTCGANTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAAACGGTGCCTGTCGCTGGAGCGTCAACGCTCAATTTTGGAGATTTTTGGGATTAATATTTTGGTAAGGCTTCAAGGTTCAAGCTAGTTCTCCAc                                                                                                  >  1:2762879/1‑150 (MQ=255)
                                                     ggtgCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCAc                                                                                                  <  2:1951588/150‑1 (MQ=255)
                                                     ggtgCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCAc                                                                                                  >  1:1232015/1‑150 (MQ=255)
                                                     ggtgCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCAc                                                                                                  >  1:5141777/1‑150 (MQ=255)
                                                     ggtgCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCAc                                                                                                  >  2:359165/1‑150 (MQ=255)
                                                      gtgCTTTCGTCTGCACCATCCCCTGCGCTGCCATTACTTCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCATGGTTCATGCTCGTTCTCCACt                                                                                                 <  2:1452396/150‑1 (MQ=255)
                                                         cTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTAc                                                                                              >  1:4894420/1‑150 (MQ=255)
                                                             cGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCa                                                                                          >  1:2107186/1‑150 (MQ=255)
                                                              gACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAg                                                                                         <  1:3661173/150‑1 (MQ=255)
                                                                cTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTa                                                                                       <  1:4744230/150‑1 (MQ=255)
                                                                    aCCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGcc                                                                                   >  1:727671/1‑150 (MQ=255)
                                                                    aCCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGcc                                                                                   <  2:4463733/150‑1 (MQ=255)
                                                                     ccATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCa                                                                                  <  1:5259927/150‑1 (MQ=255)
                                                                       aTCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAAt                                                                                <  2:2805488/150‑1 (MQ=255)
                                                                       aTCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAAt                                                                                <  2:3869033/150‑1 (MQ=255)
                                                                                 cTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGATATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAgcgcg                                                                      >  1:2006951/1‑150 (MQ=255)
                                                                                      aTTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTAAGCCAATTGCCAGCGCGcgatc                                                                 >  1:1108305/1‑150 (MQ=255)
                                                                                           aTCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGATATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATaa                                                            >  2:3426005/1‑150 (MQ=255)
                                                                                           aTCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGATATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATaa                                                            <  1:3004127/150‑1 (MQ=255)
                                                                                               gTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACg                                                        >  2:649838/1‑150 (MQ=255)
                                                                                                 aaTTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACGaa                                                      >  2:1085335/1‑150 (MQ=255)
                                                                                                   ttGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACGAAat                                                    <  1:1172960/150‑1 (MQ=255)
                                                                                                   ttGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACGAAat                                                    <  1:1890164/150‑1 (MQ=255)
                                                                                                     gTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACGAAatat                                                  >  2:4947861/1‑150 (MQ=255)
                                                                                                      ttAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACGAAATATc                                                 <  2:4168110/150‑1 (MQ=255)
                                                                                                      ttAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACGAAATATc                                                 <  1:2571676/150‑1 (MQ=255)
                                                                                                       tAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACGAAATATCa                                                <  2:2037610/150‑1 (MQ=255)
                                                                                                        aaaCAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTAGGGGAGTAATATTTTTGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACATACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACTAAATAACAg                                               >  2:4513629/1‑150 (MQ=255)
                                                                                                         aacaacCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACGAAATATCAGc                                              >  2:1773844/1‑150 (MQ=255)
                                                                                                          acaacCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACGAAATATCAGCt                                             <  2:5288603/150‑1 (MQ=255)
                                                                                                           caacCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACGAAATATCAGCTc                                            <  2:4537501/150‑1 (MQ=255)
                                                                                                            aacCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTTGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCATTACGCCAATTGCCAGCGCGCGATCGATAAAACGAAATATCAGCTcc                                           >  1:4559934/1‑150 (MQ=255)
                                                                                                            aacCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACGAAATATCAGCTcc                                           <  2:231748/150‑1 (MQ=255)
                                                                                                            aacCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACGAAATATCAGCTcc                                           >  2:3289968/1‑150 (MQ=255)
                                                                                                             acCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACGAAATATCAGCTCCa                                          <  1:1600223/150‑1 (MQ=255)
                                                                                                             acCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACGAAATATCAGCTCCa                                          <  1:4056507/150‑1 (MQ=255)
                                                                                                             acCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACGAAATATCAGCTCCa                                          <  1:359152/150‑1 (MQ=255)
                                                                                                               cGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACGAAATATCAGCTCCAGt                                        >  2:1561345/1‑150 (MQ=255)
                                                                                                               cGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACGAAATATCAGCTCCAGt                                        <  2:1898598/150‑1 (MQ=255)
                                                                                                                 gTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACGAAATATCAGCTCCAGTTg                                      >  2:3926410/1‑150 (MQ=255)
                                                                                                                  tGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACGAAATATCAGCTCCAGTTgg                                     <  1:4306629/150‑1 (MQ=255)
                                                                                                                    ccTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGGAAGGCTTCAATGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAAAGAAATATCAGCTCCAGTTGGGa                                   >  2:4909194/1‑150 (MQ=255)
                                                                                                                     cTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACGAAATATCAGCTCCAGTTGGGAg                                  <  2:4767195/150‑1 (MQ=255)
                                                                                                                     cTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACGAAATATCAGCTCCAGTTGGGAg                                  >  1:237285/1‑150 (MQ=255)
                                                                                                                        tCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACGAAATATCAGCTCCAGTTGGGAGCCa                               <  2:3721854/150‑1 (MQ=255)
                                                                                                                         cGCTGGAGCGTCATCGCTCAATTTTGGATATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACGAAATATCAGCTCCAGTTGGGAGCCAt                              >  1:1552471/1‑150 (MQ=255)
                                                                                                                          gCTGGAGCGTCATCGCTCAATTTTGGATATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACGAAATATCAGCTCCAGTTGGGAGCCATa                             <  2:434278/150‑1 (MQ=255)
                                                                                                                            tGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACGACATATCCGCTCCAGTTGGGAGCCATACt                           >  2:3043655/1‑150 (MQ=255)
                                                                                                                              gAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACGAAATATCAGCTCCAGTTGGGAGCCATACTTa                         <  1:1426252/150‑1 (MQ=255)
                                                                                                                                 cGTCATCGCTCAATTTTGGATATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACGAAATATCAGCTCCAGTTGGGAGCCATACTTATCt                      <  2:1020240/150‑1 (MQ=255)
                                                                                                                                 cGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACGAAATATCAGCTCCAGTTGGGAGCCATACTTATCt                      <  1:1139642/150‑1 (MQ=255)
                                                                                                                                  gTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACGAAATATCAGCTCCAGTTGGGAGCCATACTTATCtt                     >  1:5162163/1‑150 (MQ=255)
                                                                                                                                  gTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACGAAATATCAGCTCCAGTTGGGAGCCAGACTTATCtt                     >  2:379224/1‑150 (MQ=255)
                                                                                                                                    cATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACGAAATATCAGCTCCAGTTGGGAGCCATACTTATCTTCa                   <  1:3784529/150‑1 (MQ=255)
                                                                                                                                    cATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACGAAATATCAGCTCCAGTTGGGAGCCATACTTATCTTCa                   >  2:2697065/1‑150 (MQ=255)
                                                                                                                                     aTCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACGAAATATCAGCTCCAGTTGGGAGCCATACTTATCTTCaa                  >  1:1021993/1‑150 (MQ=255)
                                                                                                                                       cGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACGAAATATCAGCTCCAGTTGGGAGCCATACTTATCTTCAAAt                <  1:1801123/150‑1 (MQ=255)
                                                                                                                                          tCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACGAAATATCAGCTCCAGTTGGGAGCCATACTTATCTTCAAATGcc             >  2:3212338/1‑150 (MQ=255)
                                                                                                                                            aaTTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACGAAATATCAGCTCCAGTTGGGAGCCATACTTATCTTCAAATGCCAc           >  1:4484592/1‑150 (MQ=255)
                                                                                                                                            aaTTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACGAAATATCAGCTCCAGTTGGGAGCCATACTTATCTTCAAATGCCAc           <  2:2105345/150‑1 (MQ=255)
                                                                                                                                             aTTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACGAAATATCAGCTCCAGTTGGGAGCCATACTTATCTTCAAATGCCACt          <  1:3903407/150‑1 (MQ=255)
                                                                                                                                               tttGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACGAAATATCAGCTCCAGTTGGGAGCCATACTTATCTTCAAATGCCACTGt        <  2:3976396/150‑1 (MQ=255)
                                                                                                                                                 tGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACGAAATATCAGCTCCAGTTGGGAGCCATACTTATCTTCAAATGCCACTGTAt      <  1:2632274/150‑1 (MQ=255)
                                                                                                                                                  ggAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACGAAATATCAGCTCCAGTTGGGAGCCATACTTATCTTCAAATGCCACTGTATc     >  1:5101622/1‑150 (MQ=255)
                                                                                                                                                  ggAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACGAAATATCAGCTCCAGTTGGGAGCCATACTTATCTTCAAATGCCACTGTATc     >  1:2743296/1‑150 (MQ=255)
                                                                                                                                                    agaTTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACGAAATATCAGCTCCAGTTGGGAGCCATACTTATCTTCAAATGCCACTGTATCCg   >  1:3321974/1‑150 (MQ=255)
                                                                                                                                                     gaTTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATTCTCGTTCTCCACTTACTCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACGAAATATCAGCTCCAGTTGGGAGCCATACTTATCTTCAAATGCCACTGTATCCGt  >  2:3503973/1‑150 (MQ=255)
                                                                                                                                                     gaTTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACGAAATATCAGCTCCAGTTGGGAGCCATACTTATCTTCAAATGCCACTGTATCCGt  <  2:2223172/150‑1 (MQ=255)
                                                                                                                                                     gaTTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACGAAATATCAGCTCCAGTTGGGAGCCATACTTATCTTCAAATGCCACTGTATCCGt  <  2:962765/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACGAAATATCAGCTCCAGTTGGGAGCCATACTTATCTTCAAATGCCACTGTATCCGT  >  CP009273/1637309‑1637607

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: