Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,490,151 T→C 32.5% intergenic (+60/‑87) 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,490,1510TC32.5% 168.5 / 121.1 147intergenic (+60/‑87)purH/purDbifunctional phosphoribosylaminoimidazolecarboxamide formyltransferase/IMP cyclohydrolase/phosphoribosylamine‑‑glycine ligase
Reads supporting (aligned to +/- strand):  ref base T (29/63);  new base C (24/21);  total (58/89)
Fisher's exact test for biased strand distribution p-value = 1.61e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.09e-03

Only 100 of 158 total aligned reads displayed.

GGTTCGATGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAG  >  NC_002947/5490002‑5490300
                                                                                                                                                        |                                                                                                                                                     
  ttCGATGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCGTGCCCGCGCAAAAGAtgg                                                                                                                                                     >  1:3653394/1‑150 (MQ=255)
   tCGATGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGAtggt                                                                                                                                                    <  2:5144681/150‑1 (MQ=255)
    cGATGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCTCTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGAtggtg                                                                                                                                                   <  1:1008870/150‑1 (MQ=255)
    cGATGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGAtggtg                                                                                                                                                   <  1:4038857/150‑1 (MQ=255)
    cGATGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGAtggtg                                                                                                                                                   <  1:626683/150‑1 (MQ=255)
    cGATGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGAtggtg                                                                                                                                                   <  1:1008654/150‑1 (MQ=255)
      aTGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGa                                                                                                                                                 <  1:2566012/150‑1 (MQ=255)
      aTGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCGTGCCCGCGAAGAAGAtgggggg                                                                                                                                                 >  2:1319897/1‑149 (MQ=255)
       tGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTACTTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGAt                                                                                                                                                <  1:4218646/150‑1 (MQ=255)
       tGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGAt                                                                                                                                                <  1:4218857/150‑1 (MQ=255)
       tGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGAt                                                                                                                                                <  1:4218854/150‑1 (MQ=255)
        gCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGAtt                                                                                                                                               <  2:1695915/150‑1 (MQ=255)
           tGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCAc                                                                                                                                            <  1:4582101/150‑1 (MQ=255)
           tGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCAc                                                                                                                                            <  1:4581793/150‑1 (MQ=255)
           tGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCAc                                                                                                                                            <  1:4986663/150‑1 (MQ=255)
              tGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCAccgc                                                                                                                                         <  2:2285330/150‑1 (MQ=255)
               taccaaGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCAccatc                                                                                                                                        <  2:3574434/146‑4 (MQ=38)
               gCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCAccgcc                                                                                                                                        <  1:476091/150‑1 (MQ=255)
               gCTGAAGTCATCGCTGCTACCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCAccatc                                                                                                                                        <  2:4152749/150‑4 (MQ=38)
                cTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGTGGCTTGCGCCGAACACTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCAccatc                                                                                                                                        <  1:2650645/150‑4 (MQ=17)
                cTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGTGGCTTGCGCCGAACACTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCAccatc                                                                                                                                        <  1:4287692/150‑4 (MQ=17)
                cTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGTGGCTTGCGCCGAACACTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCAccatc                                                                                                                                        <  1:4287574/150‑4 (MQ=17)
                 tGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACcgccgc                                                                                                                                      <  2:1986429/150‑1 (MQ=255)
                 tGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACcgccgc                                                                                                                                      <  1:2642187/150‑1 (MQ=255)
                 tGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACcgccgc                                                                                                                                      <  1:324358/150‑1 (MQ=255)
                         tCGCTGCTGCCGACGAGGCCGTCATCGCTATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCgg                                                                                                                              <  1:940249/150‑1 (MQ=255)
                         tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTCGCgg                                                                                                                              <  1:2268023/150‑1 (MQ=37)
                         tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTCGCgg                                                                                                                              <  2:1107338/150‑1 (MQ=37)
                         tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTCGCgg                                                                                                                              <  1:4706771/150‑1 (MQ=37)
                         tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGTATTCGCgg                                                                                                                              <  1:5555627/150‑1 (MQ=255)
                         tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCgg                                                                                                                              <  1:5248897/150‑1 (MQ=255)
                         tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCgg                                                                                                                              <  2:1580376/150‑1 (MQ=255)
                         tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCgg                                                                                                                              <  1:940288/150‑1 (MQ=255)
                         tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCgg                                                                                                                              <  2:4173117/150‑1 (MQ=255)
                         tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCgg                                                                                                                              <  2:1102012/150‑1 (MQ=255)
                             tgctgcCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCGTGGGCCCCCAGACGCTGGTGGATTGAACGCCCCAGTagcgcggca                                                                                                                          >  1:3704838/1‑145 (MQ=255)
                                tgcCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCa                                                                                                                           <  2:5149829/146‑1 (MQ=25)
                                tgcCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCa                                                                                                                           >  1:5149829/1‑146 (MQ=25)
                                       gaGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCGTTCCCGCGAATAAGATTGTGGGTTCACCGCCGCCATCGCGGGCCAACCCGCGccc                                                                                                                >  2:2272638/1‑150 (MQ=255)
                                          gCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCAc                                                                                                                                            >  2:3286701/1‑119 (MQ=255)
                                          gCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCAc                                                                                                                                            <  1:3286701/119‑1 (MQ=255)
                                            cGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGGGTGCCCGCGAAAAAGAAGGTGGATTGACCCCCCCATTAGCGGGCCACCCCCCTCCcacact                                                                                                           >  2:2558193/1‑148 (MQ=255)
                                            cGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCGTTCCCGCGCAAACCATTGTGGATTCAACCCCCCCGTGGCGGGCAACCCCCCACCcacacg                                                                                                           >  2:2558197/1‑148 (MQ=255)
                                                       tGGTATACACTGGCTGGCGCCACTACCGCCACTAATTACGCGGACCCCGCGGCAGGTCGAGATCCTGTGGGCGCGGGCTGGCCCGCCAATCCTTTGGTGGATGCACCGCCGCATGCGCGGGCACGCCCGCTCCCACAGTGTTCGGCGCaa                                                                                                <  2:1699375/150‑1 (MQ=255)
                                                       gggTCTTCACGGGCCTGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCACGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAAGACGATGGTGGTTTCACCGCCGCATGCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCaa                                                                                                <  2:5551022/149‑1 (MQ=38)
                                                            ttCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGCTTCACCGCCGCATTCGCGGGCAAGCCCGCGCCCCCAGTGTTCGGCGCCAGCCAc                                                                                           >  2:4622048/1‑150 (MQ=255)
                                                            ttCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATACTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTCGCgg                                                                                                                              <  2:3054645/115‑1 (MQ=18)
                                                                 tGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGGGTGCCCCCGCAATACGGGGTTGGTTGACCCCCCCAATAGCGGGCCAAGCCGCGCCCCCACGGTTGGGCGCAAGcccccgggt                                                                                      >  2:1217694/1‑143 (MQ=255)
                                                                 tGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCGTGCCCCCGAAAACGAAGGTTGTTTCACCGCCCCAATAGCGGGCAAGCCCCCGCCCCCCGGGTTGGGCGCCAGCCCCggggg                                                                                      >  2:4430220/1‑148 (MQ=255)
                                                                  gtcATGCGCCCCTTCCGCCACTAATGACGCGGATCCGGCGGCAAGTGTAGATCGTGTGGGAGGGGGGTTGCCCGCGAATACGTTGGTGGATTCACCGCCGCATTCGGGGGAAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATc                                                                                     <  2:3617981/148‑1 (MQ=255)
                                                                  gtcATGCGCCACTTCCGCAAGTTTTTTCGGTGTTCCCTTGTCTTGTCGAGATACTGTGTGAGGGGGCTTGCCCTCTATTACGATTGGGCTTTCACGGGCGATTTCGCGGGCATGCCCTCTCCCACAGTGTTCGGCGCAAGCCACGGGATc                                                                                     <  2:1368153/148‑1 (MQ=255)
                                                                    cttGCGCCACTACCGCCACTACTTACGCGGAGCCCGCGGCAAGGGGAGACCGTGGGGGCGCGGGCTTGCCCGCCAATAGGTGGGTGTATTCCCCGCCGCATTGGCGGCAACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCAc                                                                                   <  2:18987/148‑1 (MQ=255)
                                                                    ccggcgcCACTGCCGCCACTATTTCGGGGTAGCCGGTGGCATGGGTTGTTCGTGGGGGCGGGGGCGTGCCCGCTATTCCTTGGGTGGATCCCCGGCCGCTGCCGGGGCCACCCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCAc                                                                                   <  2:18867/147‑1 (MQ=255)
                                                                                   ccACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGtttt                                                                    <  2:4620336/150‑1 (MQ=255)
                                                                                            aCGGGGTTCCCGTGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCTAATACGATGGTGTATTCACCGCCGCATTCGGGGGAATGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGaaa                                                           <  2:5418466/150‑1 (MQ=38)
                                                                                            aCGCGTAGCCCGTGGAATGTCTAGATCGTGTGGGCGGGGGGTGGCCGGCGAATACGTTGGTTGATACACCGCCGCATGCGCGGGCACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGaaa                                                           <  1:3624860/150‑1 (MQ=255)
                                                                                                 gATCCCGCGGCAAGTCGAGATCCTGGGGGAGCGGGCTTGCCCGCGAATACGATGGTGTATTCACCGCCGCTTTCGGGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGtttt                                                      <  2:453361/150‑1 (MQ=39)
                                                                                                  aTCCCGAGGCACGTCGAGATCCTGGGGGAGCGGGCTGGCCCGCGAATACGATGGTGTATTCACCGCCGCTTTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTg                                                     <  2:4979564/150‑1 (MQ=39)
                                                                                                  aTCCCGAGGCACGTCGAGATCGGTGGGTAGCGGGCTTGCCCGCGAATACGTTGGTGGTTTCACCGCCGCATTCGGGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTg                                                     <  2:4412739/150‑1 (MQ=38)
                                                                                                                     tCCTGTGGGAGCGGACTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGcc                                  <  1:5472725/150‑1 (MQ=255)
                                                                                                                     tCCTGTGGGAGCGGACTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGcc                                  <  1:5472814/150‑1 (MQ=255)
                                                                                                                          gtgGGCGGGGGCTGGCCCGCCAATACTTGGGTGTTTCCACCGCCGCTTTGGCGGGAACCCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGa                              <  1:1283598/150‑1 (MQ=255)
                                                                                                                                 gCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGc                        >  2:2487142/1‑150 (MQ=255)
                                                                                                                                     gCTTGCCGGCTAATACGATGGTGTATCCCCCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctg                    <  2:2855915/150‑1 (MQ=255)
                                                                                                                                     gCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctg                    >  1:3223810/1‑150 (MQ=255)
                                                                                                                                     gCTGGCCCGCCAATACGTTGGTGGATTCACCGCCGCATTCGGGGGAAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctg                    <  1:2778178/150‑1 (MQ=255)
                                                                                                                                      cTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctgg                   >  1:4160949/1‑150 (MQ=25)
                                                                                                                                       ttGCCCGCGATTAGGTTGGTGTTTTCACCGCGGCTTGCGGGGGAAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggc                  <  2:2520102/150‑1 (MQ=255)
                                                                                                                                       ttGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggc                  >  1:2495329/1‑150 (MQ=255)
                                                                                                                                       ttGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggc                  >  2:3197217/1‑150 (MQ=255)
                                                                                                                                         gCCTGCGAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcct                <  1:4855805/150‑1 (MQ=25)
                                                                                                                                          cccGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctg               >  2:1826825/1‑150 (MQ=255)
                                                                                                                                           ccGCGAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACGTGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctgg              >  2:4770857/1‑150 (MQ=25)
                                                                                                                                           ccGCGAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctgg              >  1:476039/1‑150 (MQ=37)
                                                                                                                                           ccGCGAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctgg              >  2:4771478/1‑150 (MQ=37)
                                                                                                                                           ccGCGAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctgg              >  2:4238273/1‑150 (MQ=37)
                                                                                                                                           ccGCGAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctgg              >  1:186700/1‑150 (MQ=37)
                                                                                                                                           ccGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGAc                                                                 <  2:4580628/99‑1 (MQ=255)
                                                                                                                                           ccGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAATGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctgg              >  2:117288/1‑150 (MQ=255)
                                                                                                                                           ccGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGa                                                    >  2:3412980/1‑112 (MQ=255)
                                                                                                                                           ccGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGa                                                    <  1:3412980/112‑1 (MQ=255)
                                                                                                                                           ccGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctgg              >  1:2431781/1‑150 (MQ=255)
                                                                                                                                           ccGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctgg              >  1:4010579/1‑150 (MQ=255)
                                                                                                                                           ccGCGAATACGATGGTGGATTCACCGCCGCAGTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctgg              >  2:3862261/1‑150 (MQ=255)
                                                                                                                                           ccGCGAAGGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCCCGGGACTCGAGGGGGTGACAGGAACGGTGTGATCATCGGCAGCGGTGGCCGGGAGCACGCccgggcctgg              >  2:2448158/1‑150 (MQ=12)
                                                                                                                                            cgcgAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGa             >  1:3048016/1‑150 (MQ=37)
                                                                                                                                            cgcgAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGa             >  1:810413/1‑150 (MQ=37)
                                                                                                                                            cgcgAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGa             >  1:1525425/1‑150 (MQ=37)
                                                                                                                                            cgcgAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGa             >  1:575708/1‑150 (MQ=37)
                                                                                                                                            cgcgAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGa             >  1:4476265/1‑150 (MQ=37)
                                                                                                                                             gcgAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGaa            >  2:4031170/1‑150 (MQ=37)
                                                                                                                                             gcgAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGaa            >  1:4183493/1‑150 (MQ=37)
                                                                                                                                             gcgAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGaa            >  1:2240138/1‑150 (MQ=255)
                                                                                                                                             gcgAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGaa            >  2:5387626/1‑150 (MQ=255)
                                                                                                                                             gcgAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGaa            >  2:1779980/1‑150 (MQ=255)
                                                                                                                                                  ctACGAGGGTGGATCCCCCGCCGCTTCCGCGGGCACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACAAGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTc        <  1:1332142/149‑1 (MQ=255)
                                                                                                                                                  aTACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTc        >  2:473748/1‑150 (MQ=255)
                                                                                                                                                  aTACGAGGGTGGTTTCCCCGCCCCATGCGGGGGCAACCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTc        <  2:2929998/150‑1 (MQ=255)
                                                                                                                                                     cgctggtggATCCACCGCCGCAGGCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGcc     <  2:3011665/147‑1 (MQ=255)
                                                                                                                                                        tggtggATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAg  >  1:1191996/1‑150 (MQ=255)
                                                                                                                                                        |                                                                                                                                                     
GGTTCGATGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAG  >  NC_002947/5490002‑5490300

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