Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 110,883 A→T 20.1% Y318F (TAC→TTC)  guaC → GMP reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273110,8830AT20.1% 702.8 / 171.7 289Y318F (TAC→TTC) guaCGMP reductase
Reads supporting (aligned to +/- strand):  ref base A (121/110);  new base T (20/38);  total (141/148)
Fisher's exact test for biased strand distribution p-value = 1.84e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.39e-01

Only 100 of 290 total aligned reads displayed.

CATGAGCTCCGAGTCTGCGATGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCAT  >  CP009273/110734‑111031
                                                                                                                                                     |                                                                                                                                                    
 aTGAGCTCCGAGTCTGCGATGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATAc                                                                                                                                                     >  2:5071766/1‑150 (MQ=255)
 aTGAGCTCCGAGTCTGCGATGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATAc                                                                                                                                                     >  2:5022693/1‑150 (MQ=255)
 aTGAGCTCCGAGTCTGCGATGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATAc                                                                                                                                                     <  2:1113492/150‑1 (MQ=255)
 aTGAGCTCCGAGTCTGCGATGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATAc                                                                                                                                                     >  1:2641091/1‑150 (MQ=255)
     gCTCCGAGTCTGCGATGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATTCGTTg                                                                                                                                                 <  2:203517/150‑1 (MQ=255)
       tCCGAGTCTGCGATGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTggg                                                                                                                                               >  1:4417388/1‑150 (MQ=255)
       tCCGAGTCTGCGATGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTggg                                                                                                                                               <  1:3422070/150‑1 (MQ=255)
            gTCTGCGATGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTc                                                                                                                                          <  2:1199122/150‑1 (MQ=255)
             tCTGCGATGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCa                                                                                                                                         >  2:263178/1‑150 (MQ=255)
             tCTGCGATGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCa                                                                                                                                         >  1:1811712/1‑150 (MQ=255)
              cTGCGATGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCAc                                                                                                                                        <  1:3867148/150‑1 (MQ=255)
               tGCGATGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACg                                                                                                                                       <  2:930127/150‑1 (MQ=255)
               tGCGATGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACg                                                                                                                                       <  2:3338785/150‑1 (MQ=255)
                 cGATGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGcc                                                                                                                                     >  2:2150319/1‑150 (MQ=255)
                  gATGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCt                                                                                                                                    <  1:4167896/150‑1 (MQ=255)
                  gATGAAACGTCACGGTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCt                                                                                                                                    <  2:3925594/150‑1 (MQ=255)
                   aTGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTg                                                                                                                                   <  1:2553062/150‑1 (MQ=255)
                    tGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATTCGTTGGGGCTTCACGCCTGa                                                                                                                                  <  1:2376857/150‑1 (MQ=255)
                    tGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGa                                                                                                                                  >  1:1007894/1‑150 (MQ=255)
                     gAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGaa                                                                                                                                 >  2:2859332/1‑150 (MQ=255)
                       aaCGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAag                                                                                                                               <  1:2009179/150‑1 (MQ=255)
                        aCGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTTCGTTCATCTTGTACATTCGTTGGGGCTTCACGCCTGAAaga                                                                                                                              <  1:4845734/150‑1 (MQ=255)
                        aCGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATTCGTTGGGGCTTCACGCCTGAAaga                                                                                                                              <  1:5092668/150‑1 (MQ=255)
                        aCGTCAAGTTTTCGGAGTTGCGTAATATCGCGCAGCAGATGGTAAAACCGTTTAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGAGATATATTTTTGGCGTACTGCTTTCAGCTTGTTCATTCGTTGGGGCTTCAAGCATGAAaga                                                                                                                              <  2:1614485/150‑1 (MQ=255)
                          gTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGc                                                                                                                            <  2:2388553/150‑1 (MQ=255)
                            cACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTg                                                                                                                          <  2:3759018/150‑1 (MQ=255)
                              cGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGAc                                                                                                                        >  1:3220083/1‑150 (MQ=255)
                               gTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGAcc                                                                                                                       >  2:4020466/1‑150 (MQ=255)
                                        gTTGCGGAATATCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATTCGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGcacc                                                                                                              >  1:4919109/1‑150 (MQ=255)
                                         ttGCGGAATATCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTACACGCCTGAAAGAGCTGACCAAGCGcacca                                                                                                             >  2:4143833/1‑150 (MQ=255)
                                          tGCGGAATATCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGcaccac                                                                                                            >  2:4690147/1‑150 (MQ=255)
                                          tGCGGAATATCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGAGATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAACGAGCTGACCAAGCGcaccac                                                                                                            >  2:3669242/1‑150 (MQ=255)
                                           gCGGAATATCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACg                                                                                                           <  1:5014468/150‑1 (MQ=255)
                                            cGGAATATCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGt                                                                                                          <  1:3103034/150‑1 (MQ=255)
                                             ggAATATCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGtt                                                                                                         <  2:4797047/150‑1 (MQ=255)
                                              gAATATCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGttt                                                                                                        >  2:4245323/1‑150 (MQ=255)
                                               aaTATCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTa                                                                                                       >  2:3468062/1‑150 (MQ=255)
                                                atatCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTAt                                                                                                      >  1:3446024/1‑150 (MQ=255)
                                                   tCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATTCGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCg                                                                                                   <  1:4720556/150‑1 (MQ=255)
                                                   tCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCg                                                                                                   <  2:4540642/150‑1 (MQ=255)
                                                    cgcgCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCgt                                                                                                  <  1:3998311/150‑1 (MQ=255)
                                                    cgcgCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCgt                                                                                                  <  2:802214/150‑1 (MQ=255)
                                                    cgcgCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCgt                                                                                                  <  2:1489766/150‑1 (MQ=255)
                                                      cgcAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCgtgt                                                                                                <  1:79859/150‑1 (MQ=255)
                                                        cagcagAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGc                                                                                              <  1:1929787/150‑1 (MQ=255)
                                                        cagcagAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGc                                                                                              >  1:1100819/1‑150 (MQ=255)
                                                         agcagAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCa                                                                                             >  2:4113686/1‑150 (MQ=255)
                                                           cagAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAgg                                                                                           >  1:2174766/1‑150 (MQ=255)
                                                           cagAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAgg                                                                                           <  1:2134946/150‑1 (MQ=255)
                                                              aaGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATTCGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAAc                                                                                        <  1:1370698/150‑1 (MQ=255)
                                                                ggTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATTCGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACaa                                                                                      <  1:4077937/150‑1 (MQ=255)
                                                                ggTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACaa                                                                                      >  2:1046037/1‑150 (MQ=255)
                                                                 gTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATTCGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAg                                                                                     <  1:1336478/150‑1 (MQ=255)
                                                                   aaaaCCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATTCGGTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGaa                                                                                   >  1:4765555/1‑150 (MQ=255)
                                                                      aCCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAAc                                                                                <  2:4281168/150‑1 (MQ=255)
                                                                       ccGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAAcc                                                                               >  2:197853/1‑150 (MQ=255)
                                                                        cGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCg                                                                              >  2:2654903/1‑150 (MQ=255)
                                                                         gTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATTCGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGc                                                                             >  1:1218629/1‑150 (MQ=255)
                                                                         gTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGc                                                                             >  1:600345/1‑150 (MQ=255)
                                                                             cgtcgCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCt                                                                         >  1:4199477/5‑150 (MQ=255)
                                                                             aGCTGCCGCTTCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCt                                                                         <  1:2835973/150‑1 (MQ=255)
                                                                             aGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCt                                                                         >  1:3577420/1‑150 (MQ=255)
                                                                                tGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTca                                                                      >  1:4386723/1‑150 (MQ=255)
                                                                                tGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTca                                                                      >  2:3881714/1‑150 (MQ=255)
                                                                                tGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTca                                                                      >  1:2913318/1‑150 (MQ=255)
                                                                                 gCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTcaa                                                                     >  1:4630838/1‑150 (MQ=255)
                                                                                  ccGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTcaac                                                                    >  1:351208/1‑150 (MQ=255)
                                                                                      tGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATTCGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAAcc                                                                >  1:4681301/1‑150 (MQ=255)
                                                                                          aGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTa                                                            <  2:1676365/150‑1 (MQ=255)
                                                                                             cccGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAAtc                                                         <  2:273899/150‑1 (MQ=255)
                                                                                               cGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAAtctc                                                       <  1:3787326/150‑1 (MQ=255)
                                                                                                 gTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTccc                                                     <  2:1869460/150‑1 (MQ=255)
                                                                                                  ttGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCa                                                    <  1:3299663/150‑1 (MQ=255)
                                                                                                     aaaaTACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACg                                                 >  2:2845404/1‑150 (MQ=255)
                                                                                                      aaaTACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGc                                                >  2:4375255/1‑150 (MQ=255)
                                                                                                       aaTACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCt                                               >  1:3106487/1‑150 (MQ=255)
                                                                                                           ccGCGCGTGATATTTTGGGCGGACTGAGTTCATCTTGTACATTCGTTGTGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCAAGTTTATTAGTGTGAAGGAACAAGAAAACAGCATCTTCAACATCATGTAATCTCCCAACGCTGGCg                                           <  1:3042781/150‑1 (MQ=255)
                                                                                                               gcgAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATTCGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGa                                       <  2:867785/150‑1 (MQ=255)
                                                                                                                cgAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGACAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAg                                      >  2:4534042/1‑150 (MQ=255)
                                                                                                                cgAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAg                                      <  1:4460988/150‑1 (MQ=255)
                                                                                                                 gagaTATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGc                                     >  2:4437853/1‑150 (MQ=255)
                                                                                                                   gaTATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCaa                                   <  1:3669147/150‑1 (MQ=255)
                                                                                                                     tatTTTGGGCGGCCTGCGTTCAGCTTGTACATTCGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAaca                                 <  1:27597/150‑1 (MQ=255)
                                                                                                                      atTTTGGGCGGCCTGCGTTCAGCTTGTACATTCGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAacac                                <  2:2204002/150‑1 (MQ=255)
                                                                                                                      atTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAacac                                >  2:1340409/1‑150 (MQ=255)
                                                                                                                         ttGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGcc                             <  2:485769/150‑1 (MQ=255)
                                                                                                                            ggcggcCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGc                          >  2:2995472/1‑150 (MQ=255)
                                                                                                                            ggcggcCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGc                          <  1:1814197/150‑1 (MQ=255)
                                                                                                                               ggcCTGCGTTCAGCTTGTACAGTCGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGtt                       <  2:731281/150‑1 (MQ=255)
                                                                                                                                      gTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTatcccat                <  2:767033/150‑1 (MQ=255)
                                                                                                                                         cAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTatcccatccc             >  2:2486303/1‑150 (MQ=255)
                                                                                                                                           gCTTGTACATTCGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCAc           <  2:1456968/150‑1 (MQ=255)
                                                                                                                                           gCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCAc           >  2:3859903/1‑150 (MQ=255)
                                                                                                                                              tGTACATTCGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTca        >  2:3152768/1‑150 (MQ=255)
                                                                                                                                                tACATACGTTGGGGCTTCACGCCTGAACGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTcatc      >  1:1721991/1‑150 (MQ=255)
                                                                                                                                                tACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACAGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTcatc      >  2:1042876/1‑150 (MQ=255)
                                                                                                                                                 aCATACGTTGGGGCTTCACGCCTGAACGAGCTGACCAAGCGCACCACGTTTATTCGTGAGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTcatcg     >  1:1551390/1‑150 (MQ=255)
                                                                                                                                                 aCATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTcatcg     <  2:2824352/150‑1 (MQ=255)
                                                                                                                                                  cATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTcatcgc    >  2:2302337/1‑150 (MQ=255)
                                                                                                                                                   aTACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTcatcgca   >  2:4083534/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CATGAGCTCCGAGTCTGCGATGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCAT  >  CP009273/110734‑111031

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: ‑