Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,487,141 C→T 15.9% intergenic (+77/‑70) purH → / → purD bifunctional phosphoribosylaminoimidazolecarboxamide formyltransferase/IMP cyclohydrolase/phosphoribosylamine‑‑glycine ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,487,1410CT15.9% 255.1 / 40.9 123intergenic (+77/‑70)purH/purDbifunctional phosphoribosylaminoimidazolecarboxamide formyltransferase/IMP cyclohydrolase/phosphoribosylamine‑‑glycine ligase
Reads supporting (aligned to +/- strand):  ref base C (42/57);  new base T (9/10);  total (52/71)
Fisher's exact test for biased strand distribution p-value = 8.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.52e-01

Only 100 of 132 total aligned reads displayed.

TGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGACCCACGCGTCGAGA  >  NC_002947/5486992‑5487289
                                                                                                                                                     |                                                                                                                                                    
tGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACcgccgc                                                                                                                                                      <  1:3192538/150‑1 (MQ=255)
tGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACcgccgc                                                                                                                                                      <  1:3192537/150‑1 (MQ=255)
    gTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTc                                                                                                                                                  <  2:182105/150‑1 (MQ=37)
      cATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTcgc                                                                                                                                                <  2:2027031/150‑1 (MQ=255)
      cATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATGACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTcgc                                                                                                                                                <  2:4816309/150‑1 (MQ=255)
       aTCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGGCGGTGAATCCACCATCGTATTcgcg                                                                                                                                               <  2:2752547/150‑1 (MQ=37)
       aTCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTcgcg                                                                                                                                               <  2:3576803/150‑1 (MQ=255)
       aTCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTcgcg                                                                                                                                               <  1:4073930/150‑1 (MQ=255)
        tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTCGCgg                                                                                                                                              <  1:2106152/150‑1 (MQ=37)
        tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTCGCgg                                                                                                                                              <  2:5094214/150‑1 (MQ=37)
        tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCgg                                                                                                                                              <  1:3705927/150‑1 (MQ=255)
        tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCgg                                                                                                                                              <  1:3751153/150‑1 (MQ=255)
        tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCgg                                                                                                                                              <  1:4686872/150‑1 (MQ=255)
        tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCgg                                                                                                                                              <  1:328869/150‑1 (MQ=255)
        tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACGAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCgg                                                                                                                                              <  1:5238062/150‑1 (MQ=255)
              ctgcCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGc                                                                                                                                        <  1:5334896/150‑1 (MQ=25)
              ctgcCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGc                                                                                                                                        <  2:3083621/150‑1 (MQ=25)
                   gacgaGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGCATACGCTGGTGGATTAACCGCCGCATTCGCGGGCAAGCCCGcg                                                                                                                                   >  2:4494959/1‑149 (MQ=255)
                   gacgaGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCAGTCGCGGGCAAGCCCGcg                                                                                                                                   >  2:4494547/1‑149 (MQ=255)
                    acgaGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCGTTCCCGCGCCTACGGTGGTGGCTTCACCGCCGCAGTCGCGGGCCAGCCCGCTc                                                                                                                                  >  2:5345698/1‑150 (MQ=255)
                            ggCATCGCGTGGGTCTTCACGGGCAGGCGCCACTGCCGCCACTAATTCGGCGGATCCCGCGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGGGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAgtg                                                                                                                          <  2:4909930/150‑1 (MQ=37)
                                         tCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCGTTCCCCCCAAAACCCTGGGGGATTCACCGCCCCAATCGCCGGCAAGCCCCCAACCACAGGGTTCGGCCCaacac                                                                                                             >  2:3110387/1‑147 (MQ=255)
                                            tCACTGGCTTGCGCAACTTCCGCCACTAATTACGGGGATCCCGCGGCACGGCTATCTCGTGGGGGCGGGGGCTGGCCGGCTAATCCTCTGGTGGATTCCCCGCCGCTTTCGGGGGCACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACg                                                                                                          <  2:2467792/150‑1 (MQ=255)
                                            tCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTCGCgg                                                                                                                                              >  1:1795562/1‑114 (MQ=18)
                                             cccTGGCAGGCGCCACTCCCGCCACTACTTACGCGGCTCCCGAGCCACGGCGCGAGCCTGTGGGTGGGGGGGTGCCGGCGACTCCGATGGGGGATCCCCCGCCGCATGCGGGGGCACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACgg                                                                                                         <  2:3732609/148‑1 (MQ=255)
                                              aCTGGCTTGCGCAACTTCCGCCACTACTGCCGCGGATCCCGAGGCAAGTCGTGATCCTGTGGGAGCGGGCTGGCCCGCGAATACGATGGTGGTTCCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACggg                                                                                                        <  2:3599185/150‑1 (MQ=37)
                                                   cttGCGCCCCACCCTCCACTAATTAGGCGGACCCGGCGGCAAGTCGATATCGTTTGGGAGCGGGGGTCCCCGCGAATACGCGGGTGGATACCCCGCCGATTTCGGGGGCACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCAc                                                                                                   <  2:1593359/148‑1 (MQ=255)
                                                             ttCCGCCACTAATTACGCGGATCCCGAGGCACGTCGAGTGCCGGGGGGCGCGGGCTTGCCCGCGAATACGTTGGTGGATTCACCGCCGCAGTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAg                                                                                         <  2:3757958/150‑1 (MQ=38)
                                                             tcccGCCACTAATTACGCGTATCCCGAGGCAAGTCGAGAGCCTGTGGGAGCGGGCTTGCCCGCGAATACGAGGGTGGTTTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAg                                                                                         <  2:3904382/148‑1 (MQ=255)
                                                                 gCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCGTGCCCGCGAAAACGATGGTGGGTTCACCGCCGCAATAGCGGGCAAGCCCCCTCCCACAGTGTTCGGCGCAAGCCACGGGATCTCATGAATCGAGGttt                                                                                     >  2:3131919/1‑150 (MQ=255)
                                                                     cTAATTACGCGGATCCCGCGGAATGTCGTGATCCTGGGGGCGCGGGGTTGCCCGCTAATAGGTTGGTGTATTCACCGCCGCATTCGCGGTCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGAc                                                                                 <  2:3314033/150‑1 (MQ=255)
                                                                           ccggggATCCCGCGGCACGTCGCGAGCCTGGGGGCGCGGGCTGGCCCGCGACTCCGAGGGTGGATTCACCGCCGCAGTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGaaa                                                                           <  2:4442429/146‑1 (MQ=255)
                                                                           aCGCGGATCCCGAGGCAAGTCGACATCCTGTGGGAGCGGGCTTGCCCGCTAATACGATGGTGCATACACCGCCCAATTCGCGGGAAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGaaa                                                                           <  2:3512712/150‑1 (MQ=38)
                                                                               ggATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCgg                                                                                                                                              >  2:127079/1‑79 (MQ=255)
                                                                               ggATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCgg                                                                                                                                              <  1:127079/79‑1 (MQ=255)
                                                                                gATCCCGCGGCAGGTCGAGAGCCTGTGGGCGCGGGCTTGCCCGCGAATACGTGGGTGTATACACCCCCGCTTACGGGGGCAACCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGtttt                                                                      <  1:3128261/150‑1 (MQ=255)
                                                                                 ttCCCGCGGCATGTCGATATCCTGGGGTTGCGGGTTTGCCGGCGAATACTTTGGTGTTTTCCCCGCGGCTTTGGGGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTg                                                                     <  2:2302376/149‑1 (MQ=255)
                                                                                  ccccGCGGCAATCCGAGATCGGGTGGGTGCGGGGTTGCCCGCCAACGCGCGGTTGGATTCACCCCCGCAGTGGCGGGCACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGa                                                                    <  2:40063/149‑1 (MQ=255)
                                                                                   cccGAGGCATGTCGATATCCTGTGGGTGCGGGCTTGCCCGCGAATACGATGCTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGat                                                                   <  2:339679/150‑1 (MQ=255)
                                                                                     cGAGGCAAGTCGAGATCCTGTGGGAGCGGGCGTGCCCGCGCAAAACATTGTGGAGTCACCCCCCCCATAGCGGGCAACCCCCCGCCCACAGTGTTCGGCGCAAGCCCCGGGCAAACATGAATACAGGGTTTGACAAGAAAGATTTTatca                                                                 >  2:3600664/1‑150 (MQ=255)
                                                                                                gagaTCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCgg                                                                                                                                              <  1:2250718/62‑1 (MQ=255)
                                                                                                gagaTCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCgg                                                                                                                                              >  2:2250606/1‑62 (MQ=255)
                                                                                                    tgcttggggTGCGGGCTTGCCCCTTAATCCTTTGTTTGATTCACCGCCGCTTTGGGGGGCACCCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGcc                                                  <  1:4876943/144‑1 (MQ=255)
                                                                                                     cgtgggggAGCGGCCGTGCCCGCGAATCCGAGGGTGGATCCACCGCCGCCTGCGCGGGCACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCg                                                 <  2:5206416/145‑1 (MQ=255)
                                                                                                      cTGTGGGCGGGGGCTTGCCCGCTATTACTATGGTGTATTCACCGCCGCATTGGCGGGAACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGt                                                <  2:233931/150‑1 (MQ=255)
                                                                                                        ttgGGAGGGGGCTTGCGCCGGAATACGATGTTGTTCTCACCGCCGCATTCGCGGTCAACCCCGCTCCCACAGTGTTCGGCGCAAGCCACTGGATCACTTGAATCGATGTTTTGACATGTAAGTTTTGATCATCGGCAGCGGTGGCCGTGa                                              <  1:2711809/149‑1 (MQ=255)
                                                                                                        gtgGGAGCGGGCTTGCCCGCGAATACGATGGTGTATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGa                                              <  2:3633192/150‑1 (MQ=255)
                                                                                                               cGGGCTTGCCCGCTAATACGAGGGTGGATTCACCGCCGAATTCGCGGGCACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGcc                                       <  1:3614634/150‑1 (MQ=255)
                                                                                                                  gCTGGCCCGCGAATCCTTTGGTGGTTACACCGCCGCATGCGCGGGCACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctg                                    <  2:4611896/150‑1 (MQ=255)
                                                                                                                   cTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctgg                                   >  1:1643022/1‑150 (MQ=255)
                                                                                                                    ttGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggc                                  >  1:2543740/1‑150 (MQ=25)
                                                                                                                     ggCCCGCGAATACGCTGGTGGATTCACCGCCGCATTCGCGGGCACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcc                                 <  2:3538219/149‑1 (MQ=255)
                                                                                                                      gCCCGCGATTACGTGGGTGGATCCCCCGCCGCATGCGCGGGCACCCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcct                                <  1:4881434/150‑1 (MQ=255)
                                                                                                                      gCCCGCGAATCCGATGGTGGATCCCCCGCCGCATTCGCGGGAAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcct                                <  2:4281716/150‑1 (MQ=255)
                                                                                                                       cccGCTAATCCGTGGGGGGTTCCACCGCCGCATGCGCGGGCACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctg                               <  1:363118/150‑1 (MQ=255)
                                                                                                                        ccGCGAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGTGCACGCcctggcctgg                              >  1:5261806/1‑150 (MQ=25)
                                                                                                                        ccGCGAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctgg                              >  1:1506545/1‑150 (MQ=37)
                                                                                                                        ccGCGAATGCGGCGGTGAATCCACCATCGAATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctgg                              >  1:1507231/1‑150 (MQ=37)
                                                                                                                        ccGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctgg                              >  2:2347118/1‑150 (MQ=255)
                                                                                                                        ccGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctgg                              >  2:1782602/1‑150 (MQ=255)
                                                                                                                        ccGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcatggcctgg                              >  2:3258108/1‑150 (MQ=255)
                                                                                                                         cgcgAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGa                                              >  1:1165183/1‑133 (MQ=25)
                                                                                                                         cgcgAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGa                                              <  2:1165183/133‑1 (MQ=25)
                                                                                                                         cgcgAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGCCCGCTCCCACAGGATCTCGACTTGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGaa                            >  2:2035673/1‑150 (MQ=17)
                                                                                                                         cgcgAATGCGGCGGTGAATCCACCATCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGCATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGa                                              >  2:1919944/1‑133 (MQ=25)
                                                                                                                         cgcgAATGCGGCGGTGAATCCACCATCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGa                                              <  1:1919944/133‑1 (MQ=25)
                                                                                                                         cgcgAATACTTTGGTGTTTACACCGCCCCCTTCGCGGGCACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGa                             <  1:446623/150‑1 (MQ=255)
                                                                                                                         cgcgAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGa                             <  2:4416167/150‑1 (MQ=255)
                                                                                                                          gcgAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGaa                            >  1:2998666/1‑150 (MQ=37)
                                                                                                                          gcgAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGaa                            >  2:1442009/1‑150 (MQ=37)
                                                                                                                           cgAATACGTTGGTGGATACACCGCCGCAGGCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGaaa                           <  1:3937737/150‑1 (MQ=255)
                                                                                                                           cgAATACGATGGTGTATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGaaa                           <  2:4530970/150‑1 (MQ=255)
                                                                                                                              ctccgtgggtggATTCCCCGCCGCATTCGCGGGCACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTc                        <  2:1125118/143‑1 (MQ=255)
                                                                                                                                aCGATGGTGTATTCCCCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGc                      <  2:2797194/150‑1 (MQ=255)
                                                                                                                                aCGATGGTGGATTCACCGCCGCTTTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGc                      <  2:5329312/150‑1 (MQ=255)
                                                                                                                                 cggcggtgaATCCACCATCGTATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGcc                     >  1:1792547/5‑150 (MQ=38)
                                                                                                                                  gATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATAATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGccc                    >  2:4858259/1‑150 (MQ=255)
                                                                                                                                    tggtggATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAg                  >  1:1261754/1‑150 (MQ=255)
                                                                                                                                    tggtggATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAg                  >  1:396691/1‑150 (MQ=255)
                                                                                                                                      gtggATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGa                >  2:1843971/1‑150 (MQ=255)
                                                                                                                                        ggATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGAcc              >  1:1682204/1‑150 (MQ=255)
                                                                                                                                           ttCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGACCCAc           <  2:2558187/150‑1 (MQ=255)
                                                                                                                                              aCCTCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCa                                                           >  2:1815885/1‑99 (MQ=255)
                                                                                                                                              aCCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGa                                                                                  <  1:4828824/76‑1 (MQ=255)
                                                                                                                                              aCCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTg                                                                     <  1:1817217/89‑1 (MQ=255)
                                                                                                                                              aCCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTg                                                                     >  2:1817217/1‑89 (MQ=255)
                                                                                                                                              aCCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGaaa                           >  2:2392990/1‑131 (MQ=255)
                                                                                                                                              aCCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGaaa                           <  1:2392990/131‑1 (MQ=255)
                                                                                                                                              aCCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGACCCAcgcg        >  2:3532489/1‑150 (MQ=255)
                                                                                                                                              aCCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGACCCAcgcg        >  1:3185474/1‑150 (MQ=255)
                                                                                                                                              aCCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGACCCAcgcg        >  2:2007913/1‑150 (MQ=255)
                                                                                                                                              aCCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGACCCAcgcg        >  2:2016885/1‑150 (MQ=255)
                                                                                                                                              aCCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTAGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGACCCAcgcg        >  1:3346739/1‑150 (MQ=255)
                                                                                                                                              aCCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTAGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGACCCAcgcg        >  1:3346183/1‑150 (MQ=255)
                                                                                                                                               ccgccgCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGaaa                                                                           <  2:2697520/82‑1 (MQ=255)
                                                                                                                                                cgccgcATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGACCCACGCGTc      >  2:3872041/1‑150 (MQ=255)
                                                                                                                                                cgccgcATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGACCCACGCGTc      >  1:1298021/1‑150 (MQ=255)
                                                                                                                                                 gccgcATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGACCCACGCGTCg     >  1:465888/1‑150 (MQ=255)
                                                                                                                                                  ccgcATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGACCCACGCGTCga    <  2:4040082/150‑1 (MQ=255)
                                                                                                                                                    gcATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGACAGTCGCCCAGGACCCACGCGTCgaga  >  2:4749354/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGACCCACGCGTCGAGA  >  NC_002947/5486992‑5487289

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