Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,893,446 T→A 68.9% K63* (AAA→TAA)  ygcW ← putative dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,893,4460TA68.9% 379.9 / 244.3 254K63* (AAA→TAA) ygcWputative dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (31/48);  new base A (80/95);  total (111/143)
Fisher's exact test for biased strand distribution p-value = 3.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

Only 100 of 262 total aligned reads displayed.

GTCCAGCACCTTATTCAGCTTACAAATACCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAA  >  CP009273/2893297‑2893595
                                                                                                                                                     |                                                                                                                                                     
gTCCAGCACCTTATTCAGCTTACAAATACCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGttt                                                                                                                                                       >  2:3525100/1‑150 (MQ=255)
   cAGCACCTTATTCAGCTTACAAATACCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTc                                                                                                                                                    >  2:1871982/1‑150 (MQ=255)
   cAGCACCTTATTCAGCTTACAAATACCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTc                                                                                                                                                    >  2:4123924/1‑150 (MQ=255)
   cAGCACCTTATTCAGCTTACAAATACCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTc                                                                                                                                                    <  2:3828105/150‑1 (MQ=255)
     gCACCTTATTCAGCTTACAAATACCGGCATTGTTAACCAGAATATCAACTGTACCTAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCaa                                                                                                                                                  <  2:1555051/150‑1 (MQ=255)
      cACCTTATTCAGCTTACAAATACCGGCATTGTTAACCAGACTATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATCATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAAt                                                                                                                                                 >  2:1488323/1‑150 (MQ=255)
        cattATTCAGCTTACAAATTCCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGTCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCTAtga                                                                                                                                               <  2:3878298/148‑3 (MQ=255)
         cTTATTCAGCTTACAAATACCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCAt                                                                                                                                              >  2:2794924/1‑150 (MQ=255)
            aTTCAGCTTACAAATACCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTc                                                                                                                                           <  1:1334384/150‑1 (MQ=255)
            aTTCAGCTTACAAATACCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTc                                                                                                                                           <  2:3915060/150‑1 (MQ=255)
            aTTCAGCTTACAAATACCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTc                                                                                                                                           <  2:111212/150‑1 (MQ=255)
                 gCTTACAAATACCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTg                                                                                                                                      <  1:1117847/150‑1 (MQ=255)
                 gCTTACAAATACCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTg                                                                                                                                      >  1:3455476/1‑150 (MQ=255)
                  cTTACAAATACCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGt                                                                                                                                     <  1:1107274/150‑1 (MQ=255)
                    tACAAATACCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGttt                                                                                                                                   >  1:1428205/1‑150 (MQ=255)
                    tACAAATACCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGttt                                                                                                                                   <  1:762856/150‑1 (MQ=255)
                     aCAAATACCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTc                                                                                                                                  <  1:3245468/150‑1 (MQ=255)
                     aCAAATACCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTc                                                                                                                                  <  2:4290876/150‑1 (MQ=255)
                     aCAAATACCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTc                                                                                                                                  >  1:2229279/1‑150 (MQ=255)
                      cAAATACCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCg                                                                                                                                 <  2:542796/150‑1 (MQ=255)
                      cAAATACCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCg                                                                                                                                 <  2:846327/150‑1 (MQ=255)
                          tACCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGt                                                                                                                             <  1:4522357/150‑1 (MQ=255)
                           aCCGGCATTGTTAACCAGAATATCATCTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCGTCTGCGGTGATACCCACCTGAATGATGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGtt                                                                                                                            <  1:1037024/150‑1 (MQ=255)
                            ccGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTa                                                                                                                           >  1:4461591/1‑150 (MQ=255)
                               gCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCt                                                                                                                        >  1:3432252/1‑150 (MQ=255)
                               gCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCt                                                                                                                        >  1:1241821/1‑150 (MQ=255)
                                cATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCtt                                                                                                                       >  1:3917079/1‑150 (MQ=255)
                                   tGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGa                                                                                                                    >  2:2622163/1‑150 (MQ=255)
                                   tGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGa                                                                                                                    >  2:3176561/1‑150 (MQ=255)
                                     ttAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACg                                                                                                                  >  1:3043879/1‑150 (MQ=255)
                                     ttAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACg                                                                                                                  <  1:1908430/150‑1 (MQ=255)
                                     ttAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACg                                                                                                                  >  1:1539131/1‑150 (MQ=255)
                                       aaCCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGaa                                                                                                                <  1:2550973/150‑1 (MQ=255)
                                         ccAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAAc                                                                                                              <  1:3662334/150‑1 (MQ=255)
                                          cAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACt                                                                                                             <  2:1827365/150‑1 (MQ=255)
                                           aGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTa                                                                                                            >  1:3828803/1‑150 (MQ=255)
                                             aaTATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAgg                                                                                                          >  2:3601193/1‑150 (MQ=255)
                                             aaTATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAgg                                                                                                          >  2:1704188/1‑150 (MQ=255)
                                               tatCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGaa                                                                                                        <  1:1061420/150‑1 (MQ=255)
                                                 tCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATa                                                                                                      <  2:2945185/150‑1 (MQ=255)
                                                  cAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATat                                                                                                     <  1:2037596/150‑2 (MQ=255)
                                                   aaCTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATaaa                                                                                                    >  2:2188293/1‑150 (MQ=255)
                                                   aaCTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATaaa                                                                                                    >  2:2188085/1‑150 (MQ=255)
                                                   aaCTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATaaa                                                                                                    >  2:4221364/1‑150 (MQ=255)
                                                    aCTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACACTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAg                                                                                                   >  1:4439792/1‑150 (MQ=255)
                                                     cTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGa                                                                                                  >  1:4307193/1‑150 (MQ=255)
                                                      tGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCACCACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGat                                                                                                 <  2:408106/150‑1 (MQ=255)
                                                             aaaCGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTgcg                                                                                          >  2:4531298/1‑150 (MQ=255)
                                                               aCGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGcc                                                                                        >  1:1253852/1‑150 (MQ=255)
                                                                cGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCa                                                                                       <  1:54762/150‑1 (MQ=255)
                                                                 gCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAg                                                                                      >  1:2825038/1‑150 (MQ=255)
                                                                 gCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAg                                                                                      <  2:3455476/150‑1 (MQ=255)
                                                                   tCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCt                                                                                    <  2:4439792/150‑1 (MQ=255)
                                                                     acaGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCttt                                                                                  >  2:4477720/1‑150 (MQ=255)
                                                                     acaGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCttt                                                                                  <  1:1230558/150‑1 (MQ=255)
                                                                     acaGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCttt                                                                                  <  1:1232008/150‑1 (MQ=255)
                                                                      caGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTg                                                                                 >  2:1513090/1‑150 (MQ=255)
                                                                      caGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTg                                                                                 <  1:3135198/150‑1 (MQ=255)
                                                                        gcaagcaGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGc                                                                               <  1:1394472/150‑1 (MQ=255)
                                                                             cagcGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCCACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACg                                                                          <  2:2305534/150‑1 (MQ=255)
                                                                               gcGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGcc                                                                        <  1:2686387/150‑1 (MQ=255)
                                                                                 gATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCAt                                                                      >  1:4208457/1‑150 (MQ=255)
                                                                                  aTAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATg                                                                     <  2:2321736/150‑1 (MQ=255)
                                                                                   tAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATgg                                                                    >  2:2679005/1‑150 (MQ=255)
                                                                                      tcttctGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCaa                                                                 >  2:3094411/1‑150 (MQ=255)
                                                                                       cttctGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCaaa                                                                >  1:1688062/1‑150 (MQ=255)
                                                                                             cgncgCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCATg                                                          >  1:3504139/1‑150 (MQ=255)
                                                                                             cggcgCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTg                                                          <  2:1241790/150‑1 (MQ=255)
                                                                                                     cTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAATGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAAc                                                  >  2:2477676/1‑150 (MQ=255)
                                                                                                       tCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCg                                                <  1:2092582/150‑1 (MQ=255)
                                                                                                       tCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCg                                                <  1:2794924/150‑1 (MQ=255)
                                                                                                       tCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGAGCCAGCTGTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCg                                                >  2:2050960/1‑150 (MQ=255)
                                                                                                           cGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTAt                                            <  1:250263/150‑1 (MQ=255)
                                                                                                              tGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTcc                                         >  1:575336/1‑150 (MQ=255)
                                                                                                                  aCCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCAcc                                     <  1:3566875/150‑1 (MQ=255)
                                                                                                                  aCCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCAcc                                     <  2:3354025/150‑1 (MQ=255)
                                                                                                                      aCCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCACCGGTa                                 <  1:1188677/150‑1 (MQ=255)
                                                                                                                        cTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCACCGGTAAc                               <  2:1553091/150‑1 (MQ=255)
                                                                                                                          gCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACaa                             <  1:1507571/150‑1 (MQ=255)
                                                                                                                             tGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTg                          >  2:1980251/1‑150 (MQ=255)
                                                                                                                               aaGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCg                        <  1:2203876/150‑1 (MQ=255)
                                                                                                                                 gTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGt                      >  1:2059293/1‑150 (MQ=255)
                                                                                                                                   ccACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGttt                    >  2:1997704/1‑150 (MQ=255)
                                                                                                                                   ccACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGttt                    <  2:1428271/150‑1 (MQ=255)
                                                                                                                                    cACCTCAACACCCTGTTATTCAATCATTTCCCTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGtttt                   <  2:3917271/150‑1 (MQ=255)
                                                                                                                                     aCCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTa                  <  1:3354463/150‑1 (MQ=255)
                                                                                                                                     aCCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTa                  <  1:1853845/150‑1 (MQ=255)
                                                                                                                                      ccTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTAc                 >  2:1254309/1‑150 (MQ=255)
                                                                                                                                       cTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTagc                <  2:4391603/150‑3 (MQ=255)
                                                                                                                                       cTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTAcc                <  1:939054/150‑1 (MQ=255)
                                                                                                                                       cTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTAcc                <  1:4328680/150‑1 (MQ=255)
                                                                                                                                         cAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTACCtt              <  1:582754/150‑1 (MQ=255)
                                                                                                                                         cAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTACCtt              <  2:1724204/150‑1 (MQ=255)
                                                                                                                                         cAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTACCtt              <  1:3731211/150‑1 (MQ=255)
                                                                                                                                          aaCACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTACCttt             <  2:3278353/150‑1 (MQ=255)
                                                                                                                                           acacCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTc            >  2:89797/1‑150 (MQ=255)
                                                                                                                                           acacCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTc            >  1:2518140/1‑150 (MQ=255)
                                                                                                                                              cccTGTTTTTCAATCATTTCCATTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAgg         <  2:3457055/150‑1 (MQ=255)
                                                                                                                                                cTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGga       >  1:408450/1‑150 (MQ=255)
                                                                                                                                                cTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGga       >  2:1699641/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTCCAGCACCTTATTCAGCTTACAAATACCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAA  >  CP009273/2893297‑2893595

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: