Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,490,151 T→C 18.9% 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,1510TC18.9% 381.4 / 82.9 197intergenic (+60/‑87)purH/purDbifunctional phosphoribosylaminoimidazolecarboxamide formyltransferase/IMP cyclohydrolase/phosphoribosylamine‑‑glycine ligase
Reads supporting (aligned to +/- strand):  ref base T (67/79);  new base C (14/21);  total (84/113)
Fisher's exact test for biased strand distribution p-value = 5.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.07e-03

Only 100 of 212 total aligned reads displayed.

GTTCGATGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAG  >  NC_002947/5490003‑5490300
                                                                                                                                                     |                                                                                                                                                      
gtTCGATGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGAtg                                                                                                                                                       <  1:669271/150‑1 (MQ=255)
 ttCGATGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCCCGAATAAGAtgg                                                                                                                                                      >  2:6952252/1‑150 (MQ=255)
  tCGATGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGAtggt                                                                                                                                                     <  1:5756452/150‑1 (MQ=255)
  tCGATGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGAtggt                                                                                                                                                     <  1:5756772/150‑1 (MQ=255)
   cGATGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGAtggtg                                                                                                                                                    <  1:1523534/150‑1 (MQ=255)
    gATGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGGTTGCCCGCGAATACGAtggtgg                                                                                                                                                   >  2:5393875/1‑150 (MQ=255)
          tGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCAc                                                                                                                                             >  2:6584005/1‑150 (MQ=255)
           gATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCAcc                                                                                                                                            <  2:6303503/150‑1 (MQ=255)
             tGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCTCTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGTGGCTTGCGCCGAACACTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCAcca                                                                                                                                           <  2:2062662/150‑2 (MQ=17)
             tGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCAccgc                                                                                                                                          <  1:7828503/150‑1 (MQ=255)
             tGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCAccgc                                                                                                                                          <  2:6243870/150‑1 (MQ=255)
              gCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCAccgcc                                                                                                                                         <  2:1617129/150‑1 (MQ=255)
              gCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCAccgcc                                                                                                                                         <  2:326679/150‑1 (MQ=255)
                     tCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTcg                                                                                                                                  <  1:4756537/150‑1 (MQ=255)
                      cATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTcgc                                                                                                                                 <  1:1413379/150‑1 (MQ=255)
                      cATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTcgc                                                                                                                                 <  1:5320083/150‑1 (MQ=255)
                       aTCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTcgcg                                                                                                                                <  1:1085852/150‑1 (MQ=255)
                       aTCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAAGTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTcgcg                                                                                                                                <  2:4059509/150‑1 (MQ=25)
                        tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGGTCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCgg                                                                                                                               <  2:2747012/150‑1 (MQ=255)
                        tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTCGCgg                                                                                                                               <  2:2329503/150‑1 (MQ=37)
                        tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTCGCgg                                                                                                                               <  1:3740925/150‑1 (MQ=37)
                        tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCgg                                                                                                                               <  1:8520336/150‑1 (MQ=255)
                        tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCgg                                                                                                                               <  1:7683758/150‑1 (MQ=255)
                            tgctgcCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCaa                                                                                                                          <  1:7312405/150‑1 (MQ=25)
                             gctgcCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCCTGCGCCACTTCCGCCACTAATTCCGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAg                                                                                                                         <  1:4568593/150‑1 (MQ=255)
                              ctgcCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGGGTGCCCGCGCAAACGCTGGTGGGTTCGCCGCCGCCTTCGCGGGCCAGc                                                                                                                        >  2:1268905/1‑150 (MQ=255)
                              ctgcCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGc                                                                                                                        <  1:7189579/150‑1 (MQ=255)
                              ctgcCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGc                                                                                                                        <  1:804776/150‑1 (MQ=255)
                              ctgcCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGc                                                                                                                        <  1:804520/150‑1 (MQ=255)
                              ctgcCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGc                                                                                                                        <  1:422128/150‑1 (MQ=255)
                               tgcCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGcc                                                                                                                       <  2:7934485/150‑1 (MQ=255)
                                         gCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTTCCCGCGAAAAAGATGGTGGCTTGAACGCCCCCATGGCGGGGAAGCCCCCGCCcaca                                                                                                             >  2:3417639/1‑150 (MQ=255)
                                         gCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAAAACGATGGTGGATTCACCGCCGCATTAGCGGGCAAGCCCGCACCcacc                                                                                                             >  2:5856302/1‑149 (MQ=255)
                                                    gATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGGTTGCCCGCGCAAACGATGGTGGGTTGACCGCCGCCTTCGCGGGCCAGCCCGCTCCCACCGGGTTCGgcgc                                                                                                  >  2:6948082/1‑150 (MQ=255)
                                                    gATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCGTGCCCGCGTCAAAGATGGTGGGTTGAACGGCGCAATCGCGGGCAAGCCCGCGCCCCCCGGGTTGGgcgc                                                                                                  >  2:3573741/1‑150 (MQ=255)
                                                     aTGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGGTTGGCCGCGGATAAGATGTTGGTTTAACCCCCCCAGTAGCCGGCCAGCCCCCGCCCACCGGGGTAGGCGCa                                                                                                 >  2:2576722/1‑150 (MQ=255)
                                                      tGGTCTTCCCGGGCAGGCGCCACTCCCCCCACTAATCCCGGGTAGCCCGCGGAACGTCTCTACCGTGTGGGAGCGGGCTGGCCCGCTACTACGTTGGTGGATACCCCGCCGCTTGCGGGGGAACGCCCGCTCCCACAGTGTTCGGCGCaa                                                                                                <  2:1838181/150‑1 (MQ=255)
                                                            tCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGTCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCGAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACg                                                                                          >  1:7078128/1‑150 (MQ=255)
                                                                tGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGAt                                                                                      <  2:1435255/150‑1 (MQ=255)
                                                                   cATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCGTTCCCGCGAAAAACATGGTGGATTCACCCCCGCATTAGCCGGCAAGCCCCCAACCACACTGTTACGCGCAAGCCACGGGATCAc                                                                                   >  2:3495736/1‑150 (MQ=255)
                                                                    aTGCGCCACTTCCGCCACTAATTACGCGGATCCCGCGGCAAGTCGAGATCCTGTGGGAGCGGGCTGGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACt                                                                                  <  2:7852589/150‑1 (MQ=255)
                                                                        gcCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTTCCCGCGAAAACGCTGGTGGATTCACCGCCGCAATAGCCGGCAAGCCCGCACCCACAGTGTTCGGCGCAAGCCACGGGATCACTTgta                                                                              >  2:8035834/1‑148 (MQ=255)
                                                                          cACTCCCGCCACTATTTACGGGGATCCGGAGGCAAGTCGAGTCCCTGTGGGCGGGGGCTTGCCCGCTAATAGGATGGTGTTTTCACCGCCGCTTTCGGGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATc                                                                            <  2:4642735/150‑1 (MQ=255)
                                                                            cTTCCGCCACTAATTCCGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCTTTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGa                                                                          <  2:5619020/150‑1 (MQ=255)
                                                                             ttCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCGTGCCCCCGAATACGGTGGTGGATTAACCGCCGCAATAGCGGGCAAGCCCGCGCCCACAGTGGTCGGCGCAAGCCAACGGATAACATGAAACgcg                                                                         >  2:3592407/1‑148 (MQ=255)
                                                                                cGCCACTAATGCCGCGTATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGtt                                                                      <  2:8551369/150‑1 (MQ=255)
                                                                                  ccACTATTGCCGCGGATCCCGCGGCAAGTCGTGAGCCTGTGGGCGCGGGCTTGCCCGCGAATACGTTGGTGTATACCCCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGtttt                                                                    <  1:8075429/150‑1 (MQ=37)
                                                                                  ccACTATTGACGCGGATCCGGCGGCACGTGGCGATCCGGGGGGCGGGGGCTTGCCCGCGAATACGTTGGGGGATTCACCGCCGCATTCGGGGGCACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGtttt                                                                    <  2:2704080/150‑1 (MQ=255)
                                                                                  ccACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGGTTTCCCGCCGAAACCATTGGGGGTTAAACGCCCCCATAGCGGGCAAGCCCGCTCCCACAGGGTTAGGCGCAACCCACCGGGTCAAGTGTAACGAGGtttt                                                                    >  2:2787024/1‑150 (MQ=255)
                                                                                       aaTTACGCGGATCCGGAGGCACGGCGAGATCCGGGGGGTGCGGGCTTGCCCGCGAATACGTTGGGGGATTCACCGCCGCATTCGGGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACAt                                                               <  2:8357454/150‑1 (MQ=255)
                                                                                         ttACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTTCCCCCGAAAAAGATGGTGGCTTCACCGCCGCATTGCCGGGCCAGCCCGCGCCCACCGGGTTCGGCCCAAGCCACCGGAACACATGAAACCAGGTTTTTGCATGa                                                             >  2:5765944/1‑150 (MQ=255)
                                                                                           ccggggAGCCCGCGGCATGTCGAGATCGTGGGGGCGGGGGCTGGCCCGCTAATACGTTGGTGGATACACCGCCGCTTTCGGGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGaaa                                                           <  1:5030052/146‑1 (MQ=255)
                                                                                             gcgGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAAAACGATGGTGGATTCACCCCCCCAATCGCGGGCCAGCCCCCTCCCACACTGTTCGGCGCAAGCCCCGGGCTAACTTTAAACGAGGGTTTTTCATGAAAGt                                                         >  2:7154583/1‑150 (MQ=255)
                                                                                               ggATCCCGAGGCAAGTCGAGATCCTGTGGGTGGGGGCTTGCCCGCTAATACGATGGTGTATCCACCGCCGCATTCGGGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGttt                                                       <  2:2870570/150‑1 (MQ=255)
                                                                                               ggATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGTTGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGttt                                                       <  2:5104382/150‑1 (MQ=255)
                                                                                                 ttCCCGCGGCTAGTCGATATCGTGTGGGAGCGGGGTGGCCCGCTAATACTTTGGTGTATTCACCGCCGCTTTCGGGGGCTAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTg                                                     <  2:5583774/149‑1 (MQ=255)
                                                                                                  tCCCGAGGCATTTCGAGATCCTGTGGGAGCGGGCTTGCCCGCTAATACGTGGGTGTATTCACCGCCGCATTGGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGa                                                    <  2:1666062/150‑1 (MQ=38)
                                                                                                                       gtgtgGGCGGGGGGTGGCCCGCCAATCCTTTGGTGGTTCCCCCGCCGAATTCGGGGGCACCCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGt                                <  1:6835746/149‑1 (MQ=255)
                                                                                                                         gtggtcgCGGGCTTGCCCGCGACTACGAGGGTGTATTCCCCGCCGATTACGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGa                              <  1:2782358/144‑1 (MQ=255)
                                                                                                                         gtgGGTGCGGGCTTGCCCGCTAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGa                              <  2:4865979/150‑1 (MQ=255)
                                                                                                                         gtgGGAGCGGGCTTGCCCGCTATTACGATGGTGGATTCCCCGCCGCTTTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGa                              <  2:7486615/150‑1 (MQ=255)
                                                                                                                         gtgGGAGCGGGCTTGCCCGCTAATACGTTGGTGTATACACCGCCGATTTCGCGGGCACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGa                              <  1:2782133/150‑1 (MQ=255)
                                                                                                                         gtgGGAGCGGGCTTGCCCGCGAATACGTGGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGa                              <  2:8223009/150‑1 (MQ=255)
                                                                                                                         gggggCGCGGGCTCGCCCGCTAATCGGTGGGTGTTTCCCCCGCCGATTGGGGGGCCACCCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGa                              <  1:2576722/148‑1 (MQ=255)
                                                                                                                             tcgCGGGCAAGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCAc                          >  2:2244089/3‑150 (MQ=255)
                                                                                                                             gggCGGGCTGGCCCGCGAATCCGTTGGTGGATTCACCGCCGCATTCGCGGGCACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCAc                          <  1:5784774/148‑1 (MQ=255)
                                                                                                                              cgCGGGCTTGCCCGCGAATACGAGGGTGGATCCACCGCCGAATTCGCGGGCACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACg                         <  2:3933045/149‑1 (MQ=255)
                                                                                                                                ggggCTTGCCCGCGATTACGTTGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGcc                       <  2:5033445/149‑1 (MQ=255)
                                                                                                                                cGGCCTTCCCCCCGACTACGTGGGGGGATCCACCGCCGCATCCGCGGCCACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGcc                       <  1:3441826/150‑1 (MQ=255)
                                                                                                                                     ttGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggc                  >  2:3118574/1‑150 (MQ=255)
                                                                                                                                     ttGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggc                  >  1:167628/1‑150 (MQ=255)
                                                                                                                                       gCCCGCTAATCCGCGGGTTCATCCCCCCCCGCTTTGGGGGGCACCCCCGCTCCCACAGTGTTCGGCACAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcct                <  1:8293649/150‑1 (MQ=255)
                                                                                                                                       gCCCGCTAATCCGCGGGTGTATTCACCGCCGCAGTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcct                <  2:315575/150‑1 (MQ=255)
                                                                                                                                       gCCCGCGAATCCGCTGGTGGATCCCCCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcct                <  2:1203231/150‑1 (MQ=255)
                                                                                                                                        cccGCGAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctg               >  1:8318618/1‑150 (MQ=37)
                                                                                                                                         ccGCGAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctgg              >  1:3041382/1‑150 (MQ=37)
                                                                                                                                         ccGCGAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctgg              >  2:3360556/1‑150 (MQ=37)
                                                                                                                                         ccGCGAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACAGGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctgg              >  1:279771/1‑150 (MQ=25)
                                                                                                                                         ccGCGAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGCCCGCTCCCACAGAGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctgg              >  2:3323686/1‑150 (MQ=25)
                                                                                                                                         ccGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctgg              >  2:8634906/1‑150 (MQ=255)
                                                                                                                                         ccGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctgg              >  1:1225699/1‑150 (MQ=255)
                                                                                                                                         ccGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctgg              >  1:562841/1‑150 (MQ=255)
                                                                                                                                         ccGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctgg              >  1:6160541/1‑150 (MQ=255)
                                                                                                                                         ccGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctgg              >  1:6408283/1‑150 (MQ=255)
                                                                                                                                          cgcgAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGa             >  2:7014717/1‑150 (MQ=37)
                                                                                                                                          cgcgAATGCGGCGGTGAATCCACCATCATATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGa             >  2:2885116/1‑150 (MQ=37)
                                                                                                                                           gcgAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGaa            >  1:4659428/1‑150 (MQ=37)
                                                                                                                                           gcgAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGaa            >  1:844193/1‑150 (MQ=37)
                                                                                                                                           gcgAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGaa            >  1:8409872/1‑150 (MQ=255)
                                                                                                                                            cgAATACTCTGGTGGATTCCCCGCCGCATGCGGGGCCACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGaaa           <  2:1688991/150‑1 (MQ=255)
                                                                                                                                            cgAATACGTTGGTGGATTCACCGCCGCATTCGCGGGCACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGaaa           <  2:3217386/150‑1 (MQ=255)
                                                                                                                                             gAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAg          >  1:5126534/1‑150 (MQ=255)
                                                                                                                                             gAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAg          >  2:8455356/1‑150 (MQ=255)
                                                                                                                                               aTACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTc        <  2:7115349/150‑1 (MQ=255)
                                                                                                                                                tACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCcaca                                                                                                             <  1:652912/48‑1 (MQ=255)
                                                                                                                                                tACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCcaca                                                                                                             >  2:652912/1‑48 (MQ=255)
                                                                                                                                                   gATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGccc    >  2:4418509/1‑150 (MQ=255)
                                                                                                                                                   gATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGccc    >  2:3005122/1‑150 (MQ=255)
                                                                                                                                                    ttggtggATTCACCGCCGCATTCGCGGGCACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCa   <  2:2528041/149‑1 (MQ=255)
                                                                                                                                                     tggtggATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAg  >  1:3527304/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
GTTCGATGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAG  >  NC_002947/5490003‑5490300

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: