Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,684,207 C→T 5.9% intergenic (+151/+47) yphH → / ← glyA putative DNA‑binding transcriptional regulator/serine hydroxymethyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,684,2070CT 5.9% 281.5 / 5.5 85intergenic (+151/+47)yphH/glyAputative DNA‑binding transcriptional regulator/serine hydroxymethyltransferase
Reads supporting (aligned to +/- strand):  ref base C (35/45);  new base T (2/3);  total (37/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

GATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTCGATAACGGCTTCATCATTGATGCTGTCCAGCACGT  >  NC_000913/2684063‑2684345
                                                                                                                                                |                                                                                                                                          
gATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATcc                                                                                                                                        >  2:758347/1‑149 (MQ=255)
gATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCATTCACGCCGCATcc                                                                                                                                        <  2:799117/149‑1 (MQ=255)
 aTAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCg                                                                                                                                       >  1:537032/1‑149 (MQ=255)
    aCGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATGCGGCGTGAACGCCTTATCCGGcc                                                                                                                                    <  1:615424/149‑2 (MQ=37)
     cGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATGCGGCGTGAACGCCTTATCCGGcct                                                                                                                                   <  1:95558/149‑3 (MQ=25)
         gAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGaac                                                                                                                               <  2:206064/149‑1 (MQ=255)
                         gaagCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACa                                                                                                               <  1:903936/149‑1 (MQ=255)
                         caagCCTTGATGCTCATAAACGAAAAATTGCCTGATGTGCTACGCTTATCAGGCCTACGAGGAGATTACAATATATTGAATTTGCAAAGTTTCGTAGGCCGGATAAAGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGAc                                                                                                                <  4:978549/147‑1 (MQ=21)
                              cTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCaa                                                                                                           >  3:335145/1‑148 (MQ=255)
                                tgatgaTCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAAt                                                                                                         >  3:913979/1‑148 (MQ=255)
                                tgatgaTCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAAt                                                                                                         <  4:310659/148‑1 (MQ=255)
                                       cATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGttt                                                                                                 <  1:200227/149‑1 (MQ=255)
                                        aTAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTc                                                                                                >  2:791917/1‑149 (MQ=255)
                                          aaaCGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACa                                                                                                               <  2:89897/132‑1 (MQ=255)
                                          aaaCGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACa                                                                                                               >  1:89897/1‑132 (MQ=255)
                                          aaaCGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGc                                                                                              >  2:262083/1‑149 (MQ=255)
                                           aaCGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCt                                                                                             >  1:787553/1‑149 (MQ=255)
                                              gAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGaac                                                                                                                               >  1:694737/1‑112 (MQ=255)
                                              gAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGaac                                                                                                                               <  2:694737/112‑1 (MQ=255)
                                              gAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTAt                                                                                          >  2:831443/1‑149 (MQ=255)
                                               aaaaaTTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTAt                                                                                          >  4:276907/1‑148 (MQ=255)
                                                aaaaTTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATg                                                                                         >  4:61410/1‑148 (MQ=255)
                                                aaaaTTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATg                                                                                         <  3:642928/148‑1 (MQ=255)
                                                aaaaTTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATg                                                                                         <  3:149221/148‑1 (MQ=255)
                                                aaaaTTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGc                                                                                        <  1:716526/149‑1 (MQ=255)
                                                          gATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCggg                                                                              <  2:387537/149‑1 (MQ=255)
                                                          gATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCgg                                                                               <  3:96112/148‑1 (MQ=255)
                                                           aTGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGt                                                                             <  1:685255/149‑1 (MQ=255)
                                                                 tACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGt                                                                        <  4:813586/148‑1 (MQ=255)
                                                                 tACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTg                                                                       >  1:782135/1‑149 (MQ=255)
                                                                  aCGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTg                                                                       >  4:580699/1‑148 (MQ=255)
                                                                     cTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTgcgc                                                                    >  4:235004/1‑148 (MQ=255)
                                                                     cTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCa                                                                   >  1:195506/1‑149 (MQ=255)
                                                                        aTCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGa                                                                 <  3:276907/148‑1 (MQ=255)
                                                                          cAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGt                                                              <  1:758347/149‑1 (MQ=255)
                                                                            ggCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTc                                                             <  4:335145/148‑1 (MQ=255)
                                                                             gCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCga                                                           <  2:707328/149‑1 (MQ=255)
                                                                             gCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCga                                                           <  2:559846/149‑1 (MQ=255)
                                                                                tACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGGGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCgaga                                                         >  4:287545/1‑148 (MQ=255)
                                                                                tACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTAt                                                                                          <  4:238674/115‑1 (MQ=255)
                                                                                tACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTAt                                                                                          >  3:238674/1‑115 (MQ=255)
                                                                                   acGGAGATTGCCATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACtt                                                     >  1:799239/1‑149 (MQ=255)
                                                                                    cGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACtt                                                     >  3:462950/1‑148 (MQ=255)
                                                                                      gagaTTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTAc                                                  <  1:791917/149‑1 (MQ=255)
                                                                                       agaTTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTAc                                                  <  4:938988/148‑1 (MQ=255)
                                                                                          ttGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCtt                                               <  3:507732/148‑1 (MQ=255)
                                                                                            gCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGGGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTg                                             <  3:287545/148‑1 (MQ=255)
                                                                                            gCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTg                                             <  3:647008/148‑1 (MQ=255)
                                                                                            gCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGa                                            <  2:787553/149‑1 (MQ=255)
                                                                                             cAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGa                                            <  4:630800/148‑1 (MQ=255)
                                                                                               atatatTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATgc                                         >  1:599104/1‑149 (MQ=255)
                                                                                                tatatTGGATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATgcg                                        >  2:580324/1‑149 (MQ=255)
                                                                                                  catcgAATTTGCAAAGTTTTGTAGGCCGGATAAGGAGTTCACGCCGCACCCGGCATGCTCAACGAGCAAAGTGACAGCAAATCACCGTTCCGCTTATGCGTAAACCGGGTAACGTGCGCAGATCTCGAGAACTTTACATGTGATgcga                                       <  4:764684/144‑2 (MQ=255)
                                                                                                   atTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCt                                      <  3:235004/148‑1 (MQ=255)
                                                                                                    ttGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTCg                                    >  2:96219/1‑149 (MQ=255)
                                                                                                    ttGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTCg                                    >  1:564245/1‑149 (MQ=255)
                                                                                                      gAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTCGAt                                  >  2:396340/1‑149 (MQ=255)
                                                                                                         tttGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTCGATaat                               <  1:831443/149‑2 (MQ=255)
                                                                                                             cAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTCGATAACGGc                            >  3:243997/1‑148 (MQ=255)
                                                                                                             cAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTCGATAACGGCt                           >  1:968668/1‑149 (MQ=255)
                                                                                                               aaGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTCGATAACGGCtt                          >  4:216852/1‑148 (MQ=255)
                                                                                                               aaGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTCGATAACGGCTtc                         <  2:479821/149‑1 (MQ=255)
                                                                                                                aGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTCGATAACGGCTtc                         >  4:510178/1‑148 (MQ=255)
                                                                                                                aGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTCGATAACGGCTtc                         >  4:854856/1‑148 (MQ=255)
                                                                                                                  ttttGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTCGATAACGGCTtcatc                      >  1:673442/1‑149 (MQ=255)
                                                                                                                   tttGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTCGATAACGGCTtcatca                     <  2:782135/149‑1 (MQ=255)
                                                                                                                    ttGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTCGATAACGGCTTcatcat                    <  2:536710/149‑1 (MQ=255)
                                                                                                                     tGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTCGATAACGGCTTcatcat                    <  3:510178/148‑1 (MQ=255)
                                                                                                                       tAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCg                                                                                               >  3:579064/1‑71 (MQ=255)
                                                                                                                       tAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCg                                                                                               <  4:579064/71‑1 (MQ=255)
                                                                                                                       tAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTCGATAACGGCTTCATCATTg                  <  3:860307/148‑1 (MQ=255)
                                                                                                                        acgCCGGATAAGGCGTTCACGCCGCATCCGGCATTAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTCGATAACGGCTTCATCATTGa                 <  3:863303/146‑1 (MQ=255)
                                                                                                                        aGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTCGATAACGGCTTCATCATTGa                 <  3:872825/148‑1 (MQ=255)
                                                                                                                         tgCCGGATGCGGCGTGAACGCCTTATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGAt                                           >  1:404374/2‑121 (MQ=38)
                                                                                                                         tgCCGGATGCGGCGTGAACGCCTTATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGAt                                           <  2:404374/120‑1 (MQ=38)
                                                                                                                         ggCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTCGATAACGGCTTCATCATTGATg               >  2:327261/1‑149 (MQ=255)
                                                                                                                          gCCGGATGCGGCGTGAACGCCTTATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGa                                            >  3:872878/1‑119 (MQ=38)
                                                                                                                           ccGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTCGATAACGGCTTCATCATTGATGc              <  3:216852/148‑1 (MQ=255)
                                                                                                                              gATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTCGATAACGGCTTCATCATTGATGCTGt           <  4:352210/148‑1 (MQ=255)
                                                                                                                              gATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTCGATAACGGCTTCATCATTGATGCTGTc          <  2:758066/149‑1 (MQ=255)
                                                                                                                               aTAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTCGATAACGGCTTCATCATTGATGCTGTcc         <  2:534482/149‑1 (MQ=255)
                                                                                                                                   ggCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTCGATAACGGCTTCATCATTGATGCTGTCCAGCa     <  2:537032/149‑1 (MQ=255)
                                                                                                                                      gTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTCGATAACGGCTTCATCATTGATGCTGTCCAGCACGt  >  2:282959/1‑149 (MQ=255)
                                                                                                                                      gTTCACGCCGCATCCGGCACGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTCGATAACGGCTTCATCAtt                   <  1:888465/132‑1 (MQ=255)
                                                                                                                                      gTTCACGCCGCATCCGGCACGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTCGATAACGGCTTCATCAtt                   >  2:888465/1‑132 (MQ=255)
                                                                                                                                                |                                                                                                                                          
GATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTCGATAACGGCTTCATCATTGATGCTGTCCAGCACGT  >  NC_000913/2684063‑2684345

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: