Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,908,486 G→C 8.7% A394G (GCC→GGC)  pyrG ← CTP synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,908,4860GC 8.7% 438.0 / 26.0 138A394G (GCC→GGC) pyrGCTP synthetase
Reads supporting (aligned to +/- strand):  ref base G (67/58);  new base C (8/4);  total (75/63)
Fisher's exact test for biased strand distribution p-value = 5.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.16e-01

Only 100 of 146 total aligned reads displayed.

GTGCGCCGAGACGCATGGTACCGCCGAGATCGCTCTTCTCGCTACGAACTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGATTGCGTCC  >  NC_000913/2908337‑2908635
                                                                                                                                                     |                                                                                                                                                     
gtgCGCCGAGACGCATGGTACCGCCGAGATCGCTCTTCTCGCTACGAACTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTgg                                                                                                                                                       >  2:1869544/1‑150 (MQ=255)
gtgCGCCGAGACGCATGGTACCGCCGAGATCGCTCTTCTCGCTACGAACTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTgg                                                                                                                                                       >  2:1799197/1‑150 (MQ=255)
gtgCGCCGAGACGCATGGTACCGCCGAGATCGCTCTTCTCGCTACGAACTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTgg                                                                                                                                                       <  2:2403142/150‑1 (MQ=255)
gtgCGCCGAGACGCATGGTACCGCCGAGATCGCTCTTCTCGCTACGAACTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTgg                                                                                                                                                       >  1:1680484/1‑150 (MQ=255)
gtgCGCCGAGACGCATGGTACCGCCGAGATCGCTCTTCTCGCTACGAACTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGAGTTCTCCATGTTgg                                                                                                                                                       <  1:1171763/150‑1 (MQ=255)
 tgCGCCGAGACGCATGGTACCGCCGAGATCGCTCTTCTCGCTACGAACTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGc                                                                                                                                                      <  1:2052127/150‑1 (MQ=255)
 tgCGCCGAGACGCATGGTACCGCCGAGATCGCTCTTCTCGCTACGAACTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGc                                                                                                                                                      >  1:1963913/1‑150 (MQ=255)
  gcgcCGAGACGCATGGTACCGCCGAGATCGCTCTTCTCGCTACGAACTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCa                                                                                                                                                     <  2:787215/150‑1 (MQ=255)
    gcCGAGACGCATGGTACCGCCGAGATCGCTCTTCTCGCTACGAACTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAAc                                                                                                                                                   <  2:463668/150‑1 (MQ=255)
        agaCGCATGGTACCGCCGAGATCGCTCTTCTCGCTACGAACTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATgg                                                                                                                                               <  2:1176549/150‑1 (MQ=255)
         gaCGCATGGTACCGCCGAGATCGCTCTTCTCGCTACGAACTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGGCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGc                                                                                                                                              <  2:1105729/150‑1 (MQ=255)
          aCGCATGGTACCGCCGAGATCGCTCTTCTCGCTACGAACTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCg                                                                                                                                             <  1:2137520/150‑1 (MQ=255)
           cGCATGGTACCGCCGAGATCGCTCTTCTCGCTACGAACTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGa                                                                                                                                            >  1:2427060/1‑150 (MQ=255)
            gCATGGTACCGCCGAGATCGCTCTTCTCGCTACGAACTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAg                                                                                                                                           <  2:1956704/150‑1 (MQ=255)
             cATGGTACCGCCGAGATCGCTCTTCTCGCTACGAACTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGc                                                                                                                                          >  1:1865315/1‑150 (MQ=255)
              aTGGTACCGCCGAGATCGCTCTTCTCGCTACGAACTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCg                                                                                                                                         <  1:1792493/150‑1 (MQ=255)
              aTGGTACCGCCGAGATCGCTCTTCTCGCTACGAACTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCAGGTTGCCAACATGGAGAACg                                                                                                                                         <  1:2289751/150‑1 (MQ=255)
               tGGTACCGCCGAGATCGCTCTTCTCGCTACGAACTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGCCAACATGGAGAACgc                                                                                                                                        <  1:1865221/150‑2 (MQ=255)
                  tACCGCCGAGATCGCTCTTCTCGCTACGAACTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAAt                                                                                                                                     >  1:678600/1‑150 (MQ=255)
                      gccgAGATCGCTCTTCTCGCTACGAACTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAAt                                                                                                                                 >  1:256331/1‑150 (MQ=255)
                      gccgAGATCGCTCTTCTCGCTACGAACTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAAt                                                                                                                                 <  1:2029525/150‑1 (MQ=255)
                       ccgAGATCGCTCTTCTCGCTACGAACTTCAACGTTGCCGTTTTCCTCGCGCCACTCCGTAATCAGCGCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCCTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTa                                                                                                                               >  1:43210/1‑150 (MQ=255)
                          agaTCGCTCTTCTCGCTACGAACTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAAc                                                                                                                             <  1:275707/150‑1 (MQ=255)
                           gaTCGCTCTTCTCGCTACGAACTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAAATAACg                                                                                                                            <  1:2132848/150‑1 (MQ=255)
                              cGCTCTTCTCGCTACGAACTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCa                                                                                                                         >  1:1421940/1‑150 (MQ=255)
                               gCTCTTCTCGCTACGAACTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCAc                                                                                                                        >  1:510843/1‑150 (MQ=255)
                                ctctTCTCGCTACGAACTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCAcc                                                                                                                       <  1:1577667/150‑1 (MQ=255)
                                    tctcGCTACGAACTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCa                                                                                                                   <  2:922864/150‑1 (MQ=255)
                                    tctcGCTACGAACTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCa                                                                                                                   <  2:1639646/150‑1 (MQ=255)
                                          tACGAACTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCa                                                                                                             >  2:2356606/1‑150 (MQ=255)
                                            cGAACTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAgg                                                                                                           >  1:2242132/1‑150 (MQ=255)
                                                cTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCaaa                                                                                                       >  2:1762414/1‑150 (MQ=255)
                                                  tCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATg                                                                                                     >  1:1922553/1‑150 (MQ=255)
                                                  tCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATg                                                                                                     <  2:847350/150‑1 (MQ=255)
                                                   cAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGTGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGc                                                                                                    <  1:1179099/150‑1 (MQ=255)
                                                   cAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGc                                                                                                    <  2:1344213/150‑1 (MQ=255)
                                                    aaCGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGGACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGGTCTCCATATTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGcc                                                                                                   <  2:2084075/150‑1 (MQ=255)
                                                      cGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCa                                                                                                 <  1:937500/150‑1 (MQ=255)
                                                      cGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCa                                                                                                 <  1:324221/150‑1 (MQ=255)
                                                         tGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGAt                                                                                              <  1:1320950/150‑1 (MQ=255)
                                                          gCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATa                                                                                             >  1:1247967/1‑150 (MQ=255)
                                                           ccGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCa                                                                                                 <  1:1433016/145‑1 (MQ=255)
                                                           ccGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCa                                                                                                 >  2:1433016/1‑145 (MQ=255)
                                                           ccGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATaa                                                                                            >  1:1593437/1‑150 (MQ=255)
                                                           ccGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATaa                                                                                            <  2:154597/150‑1 (MQ=255)
                                                               tttCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGaa                                                                                        <  1:237074/150‑1 (MQ=255)
                                                                 tCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAata                                                                                      >  2:2467463/1‑150 (MQ=255)
                                                                  cATCGCGCCACTGGGGAAGCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGTCGTTCTCCATGTTGGCAACCTGGCGAGCGTAATCAATTAAGGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAatat                                                                                     <  2:1754577/150‑1 (MQ=255)
                                                                  cATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAatat                                                                                     >  2:711470/1‑150 (MQ=255)
                                                                     cgcgCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAACTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGt                                                                                  >  2:1420977/1‑150 (MQ=255)
                                                                       cgcCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGGTAAGGAATATTGTtc                                                                                >  2:2137954/1‑150 (MQ=255)
                                                                         ccccTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTACCGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTtctc                                                                              >  1:1880169/4‑150 (MQ=255)
                                                                         ccACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTtctc                                                                              <  2:338507/150‑1 (MQ=255)
                                                                           aCTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCAc                                                                            >  1:746799/1‑150 (MQ=255)
                                                                            cTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGt                                                                                  >  1:929064/1‑143 (MQ=255)
                                                                            cTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGt                                                                                  <  2:929064/143‑1 (MQ=255)
                                                                            cTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCAcg                                                                           <  2:387880/150‑1 (MQ=255)
                                                                                 tAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCaa                                                                      >  1:1236991/1‑150 (MQ=255)
                                                                                   aTCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCaaaa                                                                    <  1:833145/150‑1 (MQ=255)
                                                                                      aGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACgcgcgttttgc                                                                 <  1:1799197/150‑8 (MQ=255)
                                                                                         gcCCCAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCACCATGGCGAGCGTAATCAGATACCGCCACCTGCATACACAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAAcgcgcg                                                              >  2:1326782/1‑150 (MQ=255)
                                                                                             caACCGGGTACTTACAGTCTGTCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCg                                                          <  1:1028633/150‑1 (MQ=255)
                                                                                                    gTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCAt                                                   <  2:1716874/150‑1 (MQ=255)
                                                                                                         tACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCtt                                              >  2:699655/1‑150 (MQ=255)
                                                                                                         tACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCtt                                              >  1:2281269/1‑150 (MQ=255)
                                                                                                             gTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTAc                                          >  1:154408/1‑150 (MQ=255)
                                                                                                             gTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTAc                                          >  2:891518/1‑150 (MQ=255)
                                                                                                             gTCTGGAACAAATTCCGTAGAGTTGGAGTTATCCATGTTGGCAACATCGCCAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTATCACGCGCAAAACGCGCGCTCGTAATCATGCCTTCTAc                                          >  2:58633/1‑150 (MQ=255)
                                                                                                               cTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGc                                        <  2:1615839/150‑1 (MQ=255)
                                                                                                                 ggCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCa                                      >  2:1162668/1‑150 (MQ=255)
                                                                                                                 ggCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCa                                      >  2:1516236/1‑150 (MQ=255)
                                                                                                                   cacaAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCAcc                                    >  1:1741190/1‑149 (MQ=255)
                                                                                                                   cacaAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACg                                    >  1:837960/1‑150 (MQ=255)
                                                                                                                        aTTCCGTAGAGTTGGCGTTCTCCATGTTGTCAACATGGCGAGCTTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACTCGCGGTCGTAATAATGACTTCTACGCCACGATAGc                               <  2:1751793/150‑1 (MQ=255)
                                                                                                                          tCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCg                             <  1:2356606/150‑1 (MQ=255)
                                                                                                                            cGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGaa                           >  2:1282526/1‑150 (MQ=255)
                                                                                                                             gTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGaaa                          >  1:966143/1‑150 (MQ=255)
                                                                                                                             gTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGaaa                          >  1:2190229/1‑150 (MQ=255)
                                                                                                                              tAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAAc                         >  2:1944575/1‑150 (MQ=255)
                                                                                                                                gagTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCg                       >  2:1509893/1‑150 (MQ=255)
                                                                                                                                 agTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGc                      >  1:1262135/1‑150 (MQ=255)
                                                                                                                                   ttGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCa                    <  2:65217/150‑1 (MQ=255)
                                                                                                                                    tGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAg                   <  2:215883/150‑1 (MQ=255)
                                                                                                                                     ggCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAgg                  >  1:1061731/1‑150 (MQ=255)
                                                                                                                                      gCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGt                 >  1:1880283/1‑150 (MQ=255)
                                                                                                                                       cGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTa                <  1:2467463/150‑1 (MQ=255)
                                                                                                                                       cGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTa                <  2:1906123/150‑1 (MQ=255)
                                                                                                                                        gctcgccATGTTGCCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCCCGATAGCCGAAACCGCCAGGTAc               >  2:2222803/6‑150 (MQ=255)
                                                                                                                                        gctcgccATGTTGCCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTAc               >  1:2085698/6‑150 (MQ=255)
                                                                                                                                         ctcgccATGTTGCCAACATGGCGCGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACg              >  2:2071708/5‑150 (MQ=255)
                                                                                                                                         ctcgccATGTTGCCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACg              >  2:1838516/5‑150 (MQ=255)
                                                                                                                                         ctcgccATGTTGCCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACg              >  1:1410409/5‑150 (MQ=255)
                                                                                                                                           ctcCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAg            <  1:1137241/150‑1 (MQ=255)
                                                                                                                                           cgcaaTGTTGCCAAAATGGAGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGACCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAg            >  2:79429/5‑150 (MQ=255)
                                                                                                                                            tcCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAgg           >  1:996470/1‑150 (MQ=255)
                                                                                                                                                tGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGATTg       <  2:1745522/150‑1 (MQ=255)
                                                                                                                                                 gTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGATTGc      >  2:1177737/1‑150 (MQ=255)
                                                                                                                                                   tGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGATTGCGt    >  2:239712/1‑150 (MQ=255)
                                                                                                                                                   tGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGATTGCGt    >  1:167652/1‑150 (MQ=255)
                                                                                                                                                     gCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGATTGCGTcc  >  2:2468210/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTGCGCCGAGACGCATGGTACCGCCGAGATCGCTCTTCTCGCTACGAACTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGATTGCGTCC  >  NC_000913/2908337‑2908635

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: