Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,490,135 T→G 11.2% intergenic (+44/‑103) 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,1350TG11.2% 332.0 / 4.4 124intergenic (+44/‑103)purH/purDbifunctional phosphoribosylaminoimidazolecarboxamide formyltransferase/IMP cyclohydrolase/phosphoribosylamine‑‑glycine ligase
Reads supporting (aligned to +/- strand):  ref base T (13/93);  new base G (16/2);  total (29/95)
Fisher's exact test for biased strand distribution p-value = 1.50e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 133 total aligned reads displayed.

GTGATCCAGCCGGGTGGTTCGATGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGC  >  NC_002947/5489987‑5490284
                                                                                                                                                      |                                                                                                                                                     
gTGATCCAGCCGGGTGGTTCGATGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCtt                                                                                                                                                      <  2:5043058/150‑1 (MQ=255)
gTGATCCAGCCGGGTGGTTCGATGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGACtt                                                                                                                                                      >  2:5472936/1‑150 (MQ=255)
gTGATCCAGCCGGGTGGTTCGATGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGACtt                                                                                                                                                      >  2:5472814/1‑150 (MQ=255)
  gATCCAGCCGGGTGGTTCGATGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGc                                                                                                                                                    >  2:1280972/1‑150 (MQ=255)
  gATCCAGCCGGGTGGTTCGATGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGc                                                                                                                                                    <  1:2909672/150‑1 (MQ=255)
  gATCCAGCCGGGTGGTTCGATGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGc                                                                                                                                                    <  2:4223717/150‑1 (MQ=255)
      cAGCCGGGTGGTTCGATGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGGGTGCCcgc                                                                                                                                                >  2:3853283/1‑150 (MQ=255)
       aGCCGGGTGGTTCGATGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCcgcg                                                                                                                                               <  2:1220395/150‑1 (MQ=255)
       aGCCGGGTGGTTCGATGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCcgcg                                                                                                                                               <  2:3120786/150‑1 (MQ=255)
       aGCCGGGTGGTTCGATGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCcgcg                                                                                                                                               <  2:3121153/150‑1 (MQ=255)
       aGCCGGGTGGTTCGATGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCGTGccccgc                                                                                                                                               >  2:5168358/1‑147 (MQ=255)
       aGCCGGGTGGTTCGATGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCGTGACcgcg                                                                                                                                               >  2:4688678/1‑150 (MQ=255)
         ccGGGTGGTTCGATGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGaa                                                                                                                                             <  2:1155538/150‑1 (MQ=255)
             gtggtTCGATGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATAACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCGTGCAcgcccaaaac                                                                                                                                         >  2:4693890/1‑143 (MQ=255)
               ggtTCGATGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGAt                                                                                                                                       <  2:2955053/150‑1 (MQ=255)
                   cGATGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCTCTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGAtggtg                                                                                                                                   <  1:1008870/150‑1 (MQ=255)
                   cGATGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGAtggtg                                                                                                                                   <  2:1167168/150‑1 (MQ=255)
                   cGATGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGAtggtg                                                                                                                                   <  1:626683/150‑1 (MQ=255)
                   cGATGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGAtggtg                                                                                                                                   <  1:4038857/150‑1 (MQ=255)
                   cGATGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGAtggtg                                                                                                                                   <  1:1008654/150‑1 (MQ=255)
                     aTGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGa                                                                                                                                 <  1:2566012/150‑1 (MQ=255)
                      tGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTACTTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGAt                                                                                                                                <  1:4218646/150‑1 (MQ=255)
                      tGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGAt                                                                                                                                <  1:4218854/150‑1 (MQ=255)
                      tGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGAt                                                                                                                                <  1:4218857/150‑1 (MQ=255)
                       gCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGTGGCTTGCGCCGAACACTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAAt                                                                                                                                <  1:4266119/150‑1 (MQ=21)
                       gCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGAtt                                                                                                                               <  2:1695915/150‑1 (MQ=255)
                        cGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTc                                                                                                                              <  1:1137803/150‑1 (MQ=255)
                          tGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCAc                                                                                                                            <  1:4986663/150‑1 (MQ=255)
                          tGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCAc                                                                                                                            <  1:4581793/150‑1 (MQ=255)
                          tGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCAc                                                                                                                            <  1:4582101/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:4152776/150‑4 (MQ=38)
                               cTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGTGGCTTGCGCCGAACACTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCAccatc                                                                                                                        <  1:4287692/150‑4 (MQ=17)
                                tGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGGCGGTGAATCCAccatcgt                                                                                                                      <  2:3094726/150‑6 (MQ=38)
                                tGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACcgccgc                                                                                                                      <  1:2642187/150‑1 (MQ=255)
                                tGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACcgccgc                                                                                                                      <  1:2642341/150‑1 (MQ=255)
                                tGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACcgccgc                                                                                                                      <  1:324358/150‑1 (MQ=255)
                                    gTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTc                                                                                                                  <  1:4813911/150‑1 (MQ=37)
                                    gTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTc                                                                                                                  <  1:1981917/150‑1 (MQ=37)
                                       aTCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTAACCGCCGCATTcgcg                                                                                                               >  2:4546723/1‑150 (MQ=255)
                                       aTCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCCCGAAAACGATGGGGGATTCACCGCCGCATTcgcg                                                                                                               >  2:2740296/1‑150 (MQ=255)
                                        tCGCTGCTGCCGACGAGGCCGTCATCGCTATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCgg                                                                                                              <  1:940249/150‑1 (MQ=255)
                                        tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTCGCgg                                                                                                              <  2:4064517/150‑1 (MQ=37)
                                        tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTCGCgg                                                                                                              <  2:5102386/150‑1 (MQ=37)
                                        tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTCGCgg                                                                                                              <  2:5546631/150‑1 (MQ=37)
                                        tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTCGCgg                                                                                                              <  2:1107338/150‑1 (MQ=37)
                                        tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTCGCgg                                                                                                              <  1:2268023/150‑1 (MQ=37)
                                        tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGTATTCGCgg                                                                                                              <  1:5555627/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                                                                                                              <  1:3929510/150‑1 (MQ=255)
                                        tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCgg                                                                                                              <  1:940267/150‑1 (MQ=255)
                                        tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCgg                                                                                                              <  1:5248897/150‑1 (MQ=255)
                                         cGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCggg                                                                                                             <  1:3323226/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)
                                                           cGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCGTTCCCGCGCAAACCATTGTGGATTCAACCCCCCCGTGGCGGGCAACCCCCCACCcacacg                                                                                           >  2:2558197/1‑148 (MQ=255)
                                                                tcgggttGGTCTTCCGGGGCAGGCGCCCCTCCCGCCCATAATTACCCGAATCCCGCCCATATTCGTGATCCGGGGGGAGCGGGCTGGCCCGCTATAACGTTGGTGGTTACACCGCCGCATTGGGGGCCAACCCCGCTCCCACAGTGTTCg                                                                                      <  2:137054/144‑1 (MQ=255)
                                                                      tGGTATACACTGGCTGGCGCCACTACCGCCACTAATTACGCGGACCCCGCGGCAGGTCGAGATCCTGTGGGCGCGGGCTGGCCCGCCAATCCTTTGGTGGATGCACCGCCGCATGCGCGGGCACGCCCGCTCCCACAGTGTTCGGCGCaa                                                                                <  2:1699375/150‑1 (MQ=255)
                                                                           ttCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGCTTCACCGCCGCATTCGCGGGCAAGCCCGCGCCCCCAGTGTTCGGCGCCAGCCAc                                                                           >  2:4622048/1‑150 (MQ=255)
                                                                           ttCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATACTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTCGCgg                                                                                                              <  2:3054645/115‑1 (MQ=18)
                                                                           ttCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATACTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTCGCgg                                                                                                              >  1:3054645/1‑115 (MQ=18)
                                                                              aCTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCGTTCCCCCCAAAACGCAGGTGGAGTACCCCCCGCCATCGCCGGCCAGCCCGCGCCCCCACGGTTGGGCGCAAGCCCCggg                                                                        >  2:4366481/1‑150 (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)
                                                                                   caggcgcCCCTCCCCCCACTATTTCCGCGCCCCCGGCGGCACGTCGTGAGCGGGGGGGGGGGGGCGTCCCGGCGAATCGGTTGGTGTATCCACCGCCGATTTCGCGGCCACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCAc                                                                   <  2:18990/147‑1 (MQ=255)
                                                                                         ccACTTCCCCCAATACTCCCGCGTATCCGGCGGCATGTGGTTTCCGTGTGGGAGCGTGCTTGCCCGCTAATAGGTTGGTGTATTCCCGGCCGAATGCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAAt                                                             <  2:2542714/150‑1 (MQ=255)
                                                                                                  ccACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGtttt                                                    <  2:4620336/150‑1 (MQ=255)
                                                                                                    aCTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGAGTGTGGTTTCACCCCCGCATTAGCGGGCAAGCCCGCTGCCACACGGTTCGGCGCAAGCCACGGGATCACTTGGATCGAGGgttata                                                  >  2:3681970/1‑147 (MQ=38)
                                                                                                           aCGGGGTTCCCGTGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCTAATACGATGGTGTATTCACCGCCGCATTCGGGGGAATGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGaaa                                           <  2:5418466/150‑1 (MQ=38)
                                                                                                           aCGCGTAGCCCGTGGAATGTCTAGATCGTGTGGGCGGGGGGTGGCCGGCGAATACGTTGGTTGATACACCGCCGCATGCGCGGGCACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGaaa                                           <  1:3624860/150‑1 (MQ=255)
                                                                                                               ggATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGTTTCACCGCCGCTTTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGttt                                       <  2:2332293/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)
                                                                                                                  tCCCGAGGCAAGTCGAGATCCTGTGGTAGCGGGCTTGCCCGCTAATACTATGGTGTATACACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGa                                    <  2:175004/150‑1 (MQ=39)
                                                                                                                                gaggTCCTGTGGGAGCGGGCTTGCCCGCGAATACGTTGGTGGATTCACCGCCGCATTCGCGGGCACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGt                      <  2:1278138/150‑1 (MQ=255)
                                                                                                                                gagaCCCGGTGGGAGCGGCCTTGCCCGCGAATCCGATGGTGTATACCCCGCCGCATCCGCGGGCACCCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGt                      <  1:3275924/150‑1 (MQ=255)
                                                                                                                                    tCCTGTGGGAGCGGGCTTGCCCGATAATACGTTGGTGGATTCACCGCCCCATGCGCGGGCAACCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGcc                  <  1:3635664/150‑1 (MQ=255)
                                                                                                                                    tCCTGTGGGAGCGGACTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGcc                  <  1:5472725/150‑1 (MQ=255)
                                                                                                                                    tCCTGTGGGAGCGGACTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGcc                  <  1:5472936/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)
                                                                                                                                                     cTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctgg   >  1:4161012/1‑150 (MQ=25)
                                                                                                                                                      ttGCCCGCGATTAGGTTGGTGTTTTCACCGCGGCTTGCGGGGGAAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggc  <  2:2520102/150‑1 (MQ=255)
                                                                                                                                                      ttGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggc  >  2:3197217/1‑150 (MQ=255)
                                                                                                                                                      ttGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggc  >  1:2495329/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
GTGATCCAGCCGGGTGGTTCGATGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGC  >  NC_002947/5489987‑5490284

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: