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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273110,8830AT12.5% 619.6 / 65.9 207Y318F (TAC→TTC) guaCGMP reductase
Reads supporting (aligned to +/- strand):  ref base A (84/97);  new base T (7/19);  total (91/116)
Fisher's exact test for biased strand distribution p-value = 8.96e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.33e-01

Only 100 of 210 total aligned reads displayed.

CATGAGCTCCGAGTCTGCGATGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATC  >  CP009273/110734‑111032
                                                                                                                                                     |                                                                                                                                                     
 aTGAGCTCCGAGTCTGCGATGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATAc                                                                                                                                                      >  1:882262/1‑150 (MQ=255)
 aTGAGCTCCGAGTCTGCGATGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATAc                                                                                                                                                      >  2:65728/1‑150 (MQ=255)
   gAGCTCCGAGTCTGCGATGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGt                                                                                                                                                    <  2:3574669/150‑1 (MQ=255)
    aGCTCCGAGTCTGCGATGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGtt                                                                                                                                                   <  2:2383527/150‑1 (MQ=255)
     gCTCCGAGTCTGCGATGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTg                                                                                                                                                  <  2:1967052/150‑1 (MQ=255)
      cTCCGAGTCTGCGATGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTgg                                                                                                                                                 <  2:4011906/150‑1 (MQ=255)
      cTCCGAGTCTGCGATGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTgg                                                                                                                                                 <  2:1906494/150‑1 (MQ=255)
      cTCCGAGTCTGCGATGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTgg                                                                                                                                                 <  1:2502868/150‑1 (MQ=255)
        ccGAGTCTGCGATGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTgggg                                                                                                                                               <  2:4221483/150‑1 (MQ=255)
         cGAGTCTGCGATGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGc                                                                                                                                              <  1:4353877/150‑1 (MQ=255)
            gTCTGCGATGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTc                                                                                                                                           <  1:17017/150‑1 (MQ=255)
            gTCTGCGATGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTc                                                                                                                                           <  1:2539220/150‑1 (MQ=255)
               tGCGATGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACg                                                                                                                                        >  2:614522/1‑150 (MQ=255)
                 cGATGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATTCGTTGGGGCTTCACGcc                                                                                                                                      >  2:1892938/1‑150 (MQ=255)
                 cGATGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGcc                                                                                                                                      >  2:145838/1‑150 (MQ=255)
                 cGATGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGcc                                                                                                                                      >  2:3958061/1‑150 (MQ=255)
                      aaaCGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGaaa                                                                                                                                 >  2:696390/1‑150 (MQ=255)
                       aaCGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATTCGTTGGGGCTTCACGCCTGAAag                                                                                                                                <  2:4060020/150‑1 (MQ=255)
                       aaCGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAag                                                                                                                                >  2:4132252/1‑150 (MQ=255)
                           tCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCt                                                                                                                            <  2:894295/150‑1 (MQ=255)
                             aCGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGa                                                                                                                          <  2:3998722/150‑1 (MQ=255)
                                    cggcgTTGCGGAATATCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAgcg                                                                                                                   >  2:849300/1‑150 (MQ=255)
                                      gcgTTGCGGAATATCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGca                                                                                                                 <  1:1824883/150‑1 (MQ=255)
                                      gcgTTGCGGAATATCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGca                                                                                                                 <  1:3281628/150‑1 (MQ=255)
                                       cgTTGCGGAATATCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGcac                                                                                                                >  1:2083597/1‑150 (MQ=255)
                                       cgTAGCGGAATATCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGATGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGGTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAATAGCTGACCAAGCGcac                                                                                                                >  1:1627372/1‑150 (MQ=255)
                                             ggAATATCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGtt                                                                                                          >  1:416345/1‑150 (MQ=255)
                                                 tatCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTAtt                                                                                                      >  2:528065/1‑150 (MQ=255)
                                                 tatCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACAAAGCGCACCACGTTTAtt                                                                                                      >  2:3092864/1‑150 (MQ=255)
                                                  atCGCGCAGCAGAAGGTAACACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATTCGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTc                                                                                                     <  1:2734605/150‑1 (MQ=255)
                                                      cgcAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATTCGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCgtgt                                                                                                 <  2:204453/150‑1 (MQ=255)
                                                      cgcAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCgtgt                                                                                                 <  1:1152721/150‑1 (MQ=255)
                                                        cagcagAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGc                                                                                               >  2:138755/1‑150 (MQ=255)
                                                             gAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGaa                                                                                          <  2:764503/150‑1 (MQ=255)
                                                                  tAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGa                                                                                     >  2:1621918/1‑150 (MQ=255)
                                                                   aaaaCCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGaa                                                                                    >  2:3971055/1‑150 (MQ=255)
                                                                   aaaaCCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGaa                                                                                    <  2:1589457/150‑1 (MQ=255)
                                                                    aaaCCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGaaa                                                                                   <  1:1316675/150‑1 (MQ=255)
                                                                      aCCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATTCGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAAc                                                                                 <  1:4467862/150‑1 (MQ=255)
                                                                      aCCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATTCGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAAc                                                                                 <  1:2055066/150‑1 (MQ=255)
                                                                      aCCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAAc                                                                                 >  2:3850873/1‑150 (MQ=255)
                                                                      aCCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAAc                                                                                 <  2:12289/150‑1 (MQ=255)
                                                                       ccGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAANcc                                                                                <  1:1439638/150‑1 (MQ=255)
                                                                       ccGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAAcc                                                                                <  1:2310471/150‑1 (MQ=255)
                                                                        cGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATTCGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCg                                                                               <  1:690864/150‑1 (MQ=255)
                                                                        cGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCg                                                                               >  2:1881221/1‑150 (MQ=255)
                                                                          tttaGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCATCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGTAACAAGAAAACCGCa                                                                             <  1:1202480/147‑1 (MQ=255)
                                                                            aaGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATc                                                                           <  2:1749592/150‑1 (MQ=255)
                                                                            aaGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATc                                                                           <  2:281703/150‑1 (MQ=255)
                                                                             aGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCt                                                                          >  2:3157257/1‑150 (MQ=255)
                                                                               cTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTc                                                                        <  1:1692796/150‑1 (MQ=255)
                                                                                tGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATTCGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTca                                                                       <  1:4454607/150‑1 (MQ=255)
                                                                                 gCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATTCGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTcaa                                                                      <  1:3237608/150‑1 (MQ=255)
                                                                                   cGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTcaaca                                                                    <  2:2674352/150‑1 (MQ=255)
                                                                                   cGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTcaaca                                                                    <  2:2281041/150‑1 (MQ=255)
                                                                                    gCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCTCGTTTATTCGTGTGCAGGAACAAGAACACCGCATCTTcaacaa                                                                   >  1:347268/1‑150 (MQ=255)
                                                                                       gCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGACTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAATGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCt                                                                <  1:376281/150‑1 (MQ=255)
                                                                                        cGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTg                                                               >  1:372404/1‑150 (MQ=255)
                                                                                         gAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGt                                                              >  2:1318366/1‑150 (MQ=255)
                                                                                         gAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGt                                                              <  1:1748810/150‑1 (MQ=255)
                                                                                             cccGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATTCGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAAtc                                                          >  2:4357574/1‑150 (MQ=255)
                                                                                             cccGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAAtc                                                          >  2:1128864/1‑150 (MQ=255)
                                                                                             cccGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGAGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAAtc                                                          >  2:1128743/1‑150 (MQ=255)
                                                                                              ccGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAAtct                                                         <  2:198117/150‑1 (MQ=255)
                                                                                                     aaaaTACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATTCGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGATTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACg                                                  <  2:572744/150‑1 (MQ=255)
                                                                                                     aaaaTACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACg                                                  >  2:949571/1‑150 (MQ=255)
                                                                                                     aaaaTACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACg                                                  >  2:3272687/1‑150 (MQ=255)
                                                                                                       aaTACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATTCGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCt                                                <  1:4153908/150‑1 (MQ=255)
                                                                                                        aTACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTg                                               >  1:3948715/1‑150 (MQ=255)
                                                                                                        aTACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTg                                               <  2:1859900/150‑1 (MQ=255)
                                                                                                         tACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTgg                                              >  1:2149297/1‑150 (MQ=255)
                                                                                                          aCCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGc                                             <  2:1810956/150‑1 (MQ=255)
                                                                                                           ccGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCg                                            >  2:984078/1‑150 (MQ=255)
                                                                                                              cgcgAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTgg                                         <  2:328515/150‑1 (MQ=255)
                                                                                                               gcgAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGa                                        >  1:2960498/1‑150 (MQ=255)
                                                                                                                 gagaTATATTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGc                                      >  1:1163277/1‑150 (MQ=255)
                                                                                                                  agaTATTTTGGGCGGCCTGCGTTCAGCTTGTACATTCGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCa                                     >  1:3148592/1‑150 (MQ=255)
                                                                                                                     tatTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAaca                                  >  1:2186515/1‑150 (MQ=255)
                                                                                                                      atTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAacac                                 <  1:4415966/150‑1 (MQ=255)
                                                                                                                       ttttGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCCACGCTGGCGTGGAGCAACACg                                <  1:849300/150‑1 (MQ=255)
                                                                                                                          tGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCa                             <  2:3941798/150‑1 (MQ=255)
                                                                                                                           gggCGGCCTGCGTTCAGCTTGTACATTCGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAg                            >  2:2682938/1‑150 (MQ=255)
                                                                                                                           gggCGGCCTGCGTTCAGCTTGTACATTCGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAg                            >  1:2650860/1‑150 (MQ=255)
                                                                                                                           gggCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAg                            <  1:2016806/150‑1 (MQ=255)
                                                                                                                              cggcCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGt                         <  2:279319/150‑1 (MQ=255)
                                                                                                                                 ccTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTat                      <  2:1467336/150‑1 (MQ=255)
                                                                                                                                 ccTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTat                      <  1:3612423/150‑1 (MQ=255)
                                                                                                                                 ccTGCGGTCAGCTTGTACAGACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCCACGCTGGCGTGGAGCAACACGCCAGCGTTat                      <  2:3662008/150‑1 (MQ=255)
                                                                                                                                  cTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTACTCTCCCCACGCTGGCGTGGAGCAACACGCCAGCGTTatc                     <  2:1412371/150‑1 (MQ=255)
                                                                                                                                  cTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTatc                     >  2:3710997/1‑150 (MQ=255)
                                                                                                                                     cGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTatccca                  >  1:615149/1‑150 (MQ=255)
                                                                                                                                     cGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTatccca                  >  1:3269400/1‑150 (MQ=255)
                                                                                                                                          aGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTAtcccatccca             >  2:1089928/1‑150 (MQ=255)
                                                                                                                                          aGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTAtcccatccca             >  2:2808364/1‑150 (MQ=255)
                                                                                                                                            cTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACt           <  1:3598417/150‑1 (MQ=255)
                                                                                                                                                tACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTcatc       <  2:3693940/150‑1 (MQ=255)
                                                                                                                                                   aTACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTcatcgca    >  2:2946164/1‑150 (MQ=255)
                                                                                                                                                   aTACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTcatcgca    >  2:2340063/1‑150 (MQ=255)
                                                                                                                                                   aTACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACGcatcgca    >  2:1287973/1‑150 (MQ=255)
                                                                                                                                                     aCGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTcatcgcatc  <  1:3850873/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CATGAGCTCCGAGTCTGCGATGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATC  >  CP009273/110734‑111032

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