Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,684,192 A→G 5.7% intergenic (+136/+62) 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,1920AG 5.7% 289.2 / 5.0 87intergenic (+136/+62)yphH/glyAputative DNA‑binding transcriptional regulator/serine hydroxymethyltransferase
Reads supporting (aligned to +/- strand):  ref base A (37/45);  new base G (2/3);  total (39/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 = 8.39e-01

GACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTCGATAACGGCTTCATCATTGATGCTGTCC  >  NC_000913/2684046‑2684338
                                                                                                                                                  |                                                                                                                                                  
gACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAg                                                                                                                                                  <  2:783604/149‑1 (MQ=255)
 aCAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAg                                                                                                                                                  >  4:872825/1‑148 (MQ=255)
  cAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAgg                                                                                                                                                 >  4:940402/1‑148 (MQ=255)
    gATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATCTATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCg                                                                                                                                               <  3:888603/148‑1 (MQ=255)
      tttGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGtt                                                                                                                                             >  3:813586/1‑148 (MQ=255)
            ggCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCAcgccg                                                                                                                                      >  2:1058008/1‑149 (MQ=255)
                 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                                                                                                        >  1:89897/1‑132 (MQ=255)
                                                           aaaCGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACa                                                                                                        <  2:89897/132‑1 (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                             >  1:564245/1‑149 (MQ=255)
                                                                                                                     ttGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTCg                             >  2:96219/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:854856/1‑148 (MQ=255)
                                                                                                                                 aGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTCGATAACGGCTtc                  >  4:510178/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                                                                                        <  4:579064/71‑1 (MQ=255)
                                                                                                                                        tAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCg                                                                                        >  3:579064/1‑71 (MQ=255)
                                                                                                                                        tAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTCGATAACGGCTTCATCATTg           <  3:860307/148‑1 (MQ=255)
                                                                                                                                         acgCCGGATAAGGCGTTCACGCCGCATCCGGCATTAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTCGATAACGGCTTCATCATTGa          <  3:863303/146‑1 (MQ=255)
                                                                                                                                         aGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTCGATAACGGCTTCATCATTGa          <  3:872825/148‑1 (MQ=255)
                                                                                                                                          tgCCGGATGCGGCGTGAACGCCTTATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGAt                                    <  2:404374/120‑1 (MQ=38)
                                                                                                                                          tgCCGGATGCGGCGTGAACGCCTTATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGAt                                    >  1:404374/2‑121 (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)
                                                                                                                                                  |                                                                                                                                                  
GACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTCGATAACGGCTTCATCATTGATGCTGTCC  >  NC_000913/2684046‑2684338

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: