Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,918,445 T→A 7.2% intergenic (+79/+16) ygiF → / ← PP_5188 thiamine triphosphatase/AsnC/Lrp family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,918,4450TA 7.2% 193.0 / 6.0 71intergenic (+79/+16)ygiF/PP_5188thiamine triphosphatase/AsnC/Lrp family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (34/31);  new base A (3/2);  total (37/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.29e-01

GGAAGCGCGTATCGAACAGGCTCGGGCGGTGGACCAGAGCCGCGGCTGGAAGCACCTGCTCAAGGCGTGATACCGCGTCGCCTGCTTCGCGGGTAAACCCGCTCCTACAAAGGTTCGCGCCAGGCGTGAGGCCGGTGTAGGCGCGGGTTTACCCGCGAATGGCCTCACCGCGGTAGCGGCAGGCTAGTGGTGGACTTGATCTCGGAAAGCGCCACAATCGAGTTCACCTCCTGAATCCCCGGCACGTTCGACAGTTTTTCGAAGAAGAACCGCTCGTAGGCCTCGATGTCCGACGTGA  >  NC_002947/5918297‑5918594
                                                                                                                                                    |                                                                                                                                                     
ggAAGCGCGTATCGAACAGGCTCGGGCGGTGGACCAGAGCCGCGGCTGGAAGCACCTGCTCAAGGCGTGATACCGCGTCGCCTGCTTCGCGGGTAAACCCGCTCCTACAAAGGTTCGCGGGGGGCGTGAGGCCGGTGTAGGCGCGGGgtt                                                                                                                                                      >  1:480243/1‑150 (MQ=255)
     cgcgTATCGAACAGGCTCGGGCGGTGGACCAGAGCCGCGGCTGGAAGCACCTGCTCAAGGCGTGATACCGCGTCGCCTGCTTCGCGGGTAAACCCGCTCCTACAAAGGTTCGCGCCAGGCGTGAGGCCGGTGTAGGCGCGGGTTTACCcg                                                                                                                                                 <  2:1246558/150‑1 (MQ=255)
      gcgTATCGAACAGGCTCGGGCGGTGGACCAGAGCCGCGGCTGGAAGCACCTGCTCAAGGCGTGATACCGCGTCGCCTGCTTCGCGGGTAAACCCGCTCCTACAAAGGTTCGCGCCAGGCGTGAGGCCGGTGTAGGCGCGGGTTTACCcgc                                                                                                                                                >  2:954603/1‑150 (MQ=255)
             gAACAGGCTCGGGCGGTGGACCAGAGCCGCGGCTGGAAGCACCTGCTCAAGGCGTGATACCGCGTCGCCTGCTTCGCGGGTAAACCCGCTCCTACAAAGGTTCGCGCCAGGCGTGAGGCCGGTGTAGGCGCGGGTTTACCCGCGAATGGc                                                                                                                                         <  2:1837635/150‑1 (MQ=255)
              aaCAGGCTCGGGCGGTGGACCAGAGCCGCGGCTGGAAGCACCTGCTCAAGGCGTGATACCGCGTCGCCTGCTTCGCGGGTAAACCCGCTCCTACAAAGGTTCGCGCCAGGCGTGAGGCCGGTGTAGGCGCGGGTTTACCCGCGAATGGcc                                                                                                                                        <  1:658429/150‑1 (MQ=255)
                                aCCAGAGCCGCGGCTGGAAGCACCTGCTCAAGGCGTGATACCGCGTCGCCTGCTTCGCGGGTAAACCCGCTCCTACAAAGGTTCGCGCCAGGCGTGAGGCCGGTGTAGGCGCGGGTTTACCCGCGAATGGCCTCACCGCGGTAGCGGCAg                                                                                                                      <  2:1038304/150‑1 (MQ=255)
                                      gCCGCGGCTGGAAGCACCTGCTCAAGGCGTGATACCGCGTCGCCTGCTTCGCGGGTAAACCCGCTCCTACAAAGGTTCGCGCCAGGCGTGAGGCCGGTGTAGGCGCGGGTTTACCCGCGAATGGCCTCACCGCGGTAGCGGCAGGCTAgt                                                                                                                <  1:1014276/150‑1 (MQ=255)
                                       ccGCGGCTGGAAGCACCTGCTCAAGGCGTGATACCGCGTCGCCTGCTTCGCGGGTAAACCCGCTCCTACAAAGGTTCGCGCCAGGCGTGAGGCCGGTGTAGGCGCGGGTTTACCCGCGAATGGCCt                                                                                                                                       <  1:685925/126‑1 (MQ=255)
                                       ccGCGGCTGGAAGCACCTGCTCAAGGCGTGATACCGCGTCGCCTGCTTCGCGGGTAAACCCGCTCCTACAAAGGTTCGCGCCAGGCGTGAGGCCGGTGTAGGCGCGGGTTTACCCGCGAATGGCCt                                                                                                                                       >  2:685925/1‑126 (MQ=255)
                                       ccGCGGCTGGAAGCACCTGCTCAAGGCGTGATACCGCGTCGCCTGCTTCGCGGGTAAACCCGCTCCTACAAAGGTTCGCGCCAGGCGTGAGGCCGGTGTAGGCGCGGGTTTACCCGCGAATGGCCTCACCGCGGTAGCGGCAGGCTAgtg                                                                                                               <  1:1115161/150‑1 (MQ=255)
                                       ccGCGGCTGGAAGCACCTGCTCAAGGCGTGATACCGCGTCGCCTGCTTCGCGGGTAAACCCGCTCCTACAAAGGTTCGCGCCAGGCGTGAGGCCGGTGTAGGCGCGGGTTTACCCGCGAATGGCCTCACCGCGGTAGCGGCAGGCTAgtg                                                                                                               >  1:459111/1‑150 (MQ=255)
                                         gcgGCTGGAAGCACCTGCTCAAGGCGTGATACCGCGTCGCCTGCTTCGCGGGTAAACCCGCTCCTACAAAGGTTCGCGCCAGGCGTGAGGCCGGTGTAGGCGCGGGTTTACCCGCGAATGGCCTCACCGCGGTAGCGGCAGGCTAgtggt                                                                                                             <  1:508256/150‑1 (MQ=255)
                                          cgGCTGGAAGCACCTGCTCAAGGCGTGATACCGCGTCGCCTGCTTCGCGAGTAAACCCGCTCCTACAAAGGTTCGCGCCAGGCGTGAGGCCGGTGTAGGCGCGGGTTTACCCGCGAATGGCCTCACCGCGGTAGCGGCAGGCTAgtggtg                                                                                                            <  1:213560/150‑1 (MQ=255)
                                            gCTGGAAGCACCTGCTCAAGGCGTGATACCGCGTCGCCTGCTTCGCGGGTAAACCCGCTCCTACAAAGGTTCGCGCCAGGCGTGAGGCCGGTGTAGGCGCGGGTTTACCCGCGAATGGCCTCACCGCGGTAGCGGCAGTCTAGTGGTGGa                                                                                                          <  1:807853/150‑1 (MQ=255)
                                            gCTGGAAGCACCTGCTCAAGGCGTGATACCGCGTCGCCTGCTTCGCGGGTAAACCCGCTCCTACAAAGGTTCGCGCCAGGCGTGAGGCCGGTGTAGGCGCGGGTTTACCCGCGAATGGCCTCACCGCGGTAGCGGCAGGCTAGTGGTGGa                                                                                                          <  2:235644/150‑1 (MQ=255)
                                            gCTGGAAGCACCTGCTCAAGGCGTGATACCGCGTCGCCTGCTTCGCGGGTAAACCCGCTCCTACAAAGGTTCGCGCCAGGCGTGAGGCCGGTGTAGGCGCGGGTTTACCCGCGAATGGCCTCACCGCGGTAGCGGCAGGCTAGTGGTGGa                                                                                                          <  1:573240/150‑1 (MQ=255)
                                             atgggggCACCGGCGGAAGGCGTGGGACCGCGTCGCCTGCTTCGCGGGTAAACCCGCTCCTACAAAGGTTCGCGCCAGGCGTGAGGCCGGTGTAGGCGCGGGTTTACCCGCGAATGGCCTCACCGCGGTAGCGGCAGGCTAGTGGTGGAc                                                                                                         <  2:1849336/144‑1 (MQ=255)
                                                gAAGCACCTGCTCAAGGCGTGATACCGCGTCGCCTGCTTCGCGGGTAAACCCGCTCCTACAAAGGTTCGCGCCAGGCGTGAGGCCGGTGTAGGCGCGGGTTTACCCGCGAATGGCCTCACCGCGGTAGCGGCAGGCCAGTGGGGGGCttt                                                                                                      >  2:1361046/1‑149 (MQ=255)
                                                 aaGCACCTGCTCAAGGCGTGATACCGCGTCGCCTGCTTCGCGGGTAAACCCGCTCCTACAAAGGTTCGCGCCAGGCGTGAGGCCGGTGTAGGCGCGGGTTTACCCGCGAATGGCCTCACCGCGGTAGCGGCAGGCTAGTGGTGGACTTGa                                                                                                     <  1:954603/150‑1 (MQ=255)
                                                   gCACCTGCTCAAGGCGTGATACCGCGTCGCCTGCTTCGCGGGTAAACCCGCTCCTACAAAGGTTCGCGCCAGGCGTGAGGCCGGTGTAGGCGCGGGTTTACCCGCGAATGGCCTCACCGCGGTAGCGGCAGGCTAGTGGTGGACTTGAtc                                                                                                   >  1:1408295/1‑150 (MQ=255)
                                                    cACCTGCTCAAGGCGTGACACCGCGTCGCCTGCTTCGCGGGTAAACCCGCTCCTACAAAGGTTCGCGCCAGGCGTGAGGCCGGTGTAGGCGCGGGCTCACCCGCGAATGGCCTCACCGCGGTAGCGGCAGGCTAGTGGTGGACTTGAtct                                                                                                  <  1:137560/150‑1 (MQ=255)
                                                     aCCTGCTCAAGGCGTGATACCGCGTCGCCTGCTTCGCGGGTAAACCCGCTCCTACAAAGGTTCGCGCCAGGCGTGAGGCCGGTGTAGGCGCGGGTTTACCCGCGAATGGCCTCACCGCGGTAGCGGCAGGCTAGTGGTGGACTTGAtctc                                                                                                 <  2:1141010/150‑1 (MQ=255)
                                                      ccTGCTCAAGGCGTGATACCGCGTCGCCTGCTTCGCGGGTAAACCCGCTCCTACAAAGGTTCGCGCCAGGCGTGAGGCCGGTGTAGGCGCGGGTTTACCCGCGAATGGCCTCACCGCGGTAGCGGCAGGCTAGTGGTGGACTTGATCTCg                                                                                                >  1:1906719/1‑150 (MQ=255)
                                                      ccTGCTCAAGGCGTGATACCGCGTCGCCTGCTTCGCGGGTAAACCCGCTCCTACAAAGGTTCGCGCCAGGCGTGAGGCCGGTGTAGGCGCGGGTTTACCCGCGAATGGCCTCACCGCGGTAGCGGCAGGCTAGTGGTGGACTTGATCTCg                                                                                                <  1:668895/150‑1 (MQ=255)
                                                      ccTGCTCAAGGCGTGATACCGCGTCGCCTGCTTCGCGGGTAAACCCGCTCCTACAAAGGTTCGCGCCAGGCGTGAGGCCGGTGTAGGCGCGGGTTTACCCGCGAATGGCCTCACCGCGGTAGCGGCAGGCTAGTGGTGGACTTGATCTCg                                                                                                <  1:1361046/150‑1 (MQ=255)
                                                      ccTGCTCAAGGCGTGATACCGCGTCGCCTGCTTCGCGGGTAAACCCGCTCCTACAAAGGTTCGCGCCAGGCGTGAGGCCGGTAGCGCAACGCTGGTTCCCGCGAATGGCCTGACCGCGGTAGCGGCAGGCTAGTGGTGGACGTGATCTCg                                                                                                <  2:544047/150‑1 (MQ=37)
                                                       cTGCTCAAGGCGTGATACCGCGTCGCCTGCTTCGCGGGTAAACCCGCTCCTACAAAGGTTCGCGCCAGGCGTGAGGCCGGTGTAGGCGCGGGTTTACCCGCGAATGGCCTCACCGCGGTAGCGGCAGGCTAGTGGTGGACTTGATCTCgg                                                                                               >  1:1709412/1‑150 (MQ=255)
                                                        tGCTCAAGGCGTGATACCGCGTCGCCTGCTTCGCGGGTAAACCCGCTCCTACAAAGGTTCGCGCCAGGCGTGAGGCCGGTGTAGGCGCGGGTTTACCCGCGAATGGCCTCACCGCGGTAGCGGCAGGCTAGTGGTGGGCTTGATCTCGGa                                                                                              >  1:590405/1‑150 (MQ=255)
                                                        tGCTCAAGGCGTGATACCGCGTCGCCTGCTTCGCGGGTAAACCCGCTCCTACAAAGGTTCGCGCCAGGCGTGAGGCCGGGGTAGGCGCGGGTTTACCCGCGAATGGCCTCACCGCGGTAGCGGCAGGCTAGTGGTGGACTTGATCTCGGa                                                                                              >  2:640495/1‑150 (MQ=255)
                                                             aaGGCGTGATACCGCGTCGCCTGCTTCGCGGGTAAACCCGCTCCTACAAAGGTTCGCGCCAGGCGTGAGGCCGGTGTAGGCGCGGGTTTACCCGCGAATGGCCTCACCGCGGTAGCGGCAGGCTAgtggtg                                                                                                            >  1:613991/1‑131 (MQ=255)
                                                             aaGGCGTGATACCGCGTCGCCTGCTTCGCGGGTAAACCCGCTCCTACAAAGGTTCGCGCCAGGCGTGAGGCCGGTGTAGGCGCGGGTTTACCCGCGAATGGCCTCACCGCGGTAGCGGCAGGCTAgtggtg                                                                                                            <  2:613991/131‑1 (MQ=255)
                                                             aaGGCGTGATACCGCGTCGCCTGCTTCGCGGGTAAACCCGCTCCTACAAAGGTTCGCGCCAGGCGTGAGGCCGGTGTAGGCGCGGGTTTACCCGCGAATGGCCTCACCGCGGTAGCGGCAGGCTAGTGGTGGACTTGATCTCGGAAAgcg                                                                                         >  1:872931/1‑150 (MQ=255)
                                                             aaGGCGTGATACCGCGTCGCCTGCTTCGCGGGTAAACCCGCTCCTACAAAGGTTCGCGCCAGGCGTGAGGCCGGTGTAGGCGCGGGTTTACCCGCGAATGGCCTCACCGCGGTAGCGGCAGGCTAGTGGTGGACTTGATCTCGGAAAgcg                                                                                         >  1:881174/1‑150 (MQ=255)
                                                                     aTACCGCGTCGCCTGCTTCGCGGGTAAACCCGCTCCTACAAAGGTTCGCGCCAGGCGTGAGGCCGGTGTAGGCGCGGGTTTACCCGCGAATGGCCTCACCGCGGTAGCGGCAGGCTAGTGGTGGACTTGATCTCGGAAAGCGCCACAATc                                                                                 <  1:781501/150‑1 (MQ=255)
                                                                       aCCGCGTCGCCTGCTTCGCGGGTAAACCCGCTCCTACAAAGGTTCGCGCCAGGCGTGAGGCCGGTGTAGGCGCGGGTTTACCCGCGAATGGCCTCACCGCGGTAGCGGCAGGCTAGTGGTGGACTTGATCTCGGAAAGCGCCACAATCGa                                                                               <  1:249592/150‑1 (MQ=255)
                                                                       aCCGCGTCGCCTGCTTCGCGGGTAAACCCGCTCCTACAAAGGTTCGCGCCAGGCGTGAGGCCGGTGTAGGCGCGGGTTTACCCGCGAATGGCCTCACCGCGGTAGCGGCAGGCTAGTGGTGGACTTGATCTCGGAAAGCGCCACAATCGa                                                                               >  2:249592/1‑150 (MQ=255)
                                                                       aCCGCGTCGCCTGCTTCGCGGGTAAACCCGCTCCTACAAAGGTTCGCGCCAGGCGTGAGGCCGGTGTAGGCGCGGGTTTACCCGCGAATGGCCTCACCGCGGTAGCGGCAGGCTAGTGGTGGACTTGATCTCGGAAAGCGCCACAATCGa                                                                               >  2:1626249/1‑150 (MQ=255)
                                                                           cgTCGCCTGCTTCGCGGGTAAACCCGCTCCTCCAAAGGTTCGCGCCAGGCGTGAGGCCGGTGTAGGCGCGGGTTTACCCGCGAATGGCCTCACCGCGGTAGCGGCAGGCTAGTGGTGGACTTGATCTCGGAAAGCGCCACAATCGAGTTc                                                                           <  1:123288/150‑1 (MQ=255)
                                                                                 cTGCTTCGCGGGTAAACCCGCTCCTACAAAGGTTCGCGCCAGGCGTGAGGCCGGTGTAGGCGCGGGTTTACCCGCGAATGGCCTCACCGCGGTAGCGGCAGGCTAGTGGTGGACTTGATCTCGGAAAGCGCCACAATCGAGTTCAcctcc                                                                     >  1:715470/1‑150 (MQ=255)
                                                                                     ttCGCGGGTAAACCCGCTCCTACAAAGGTTCGCGCCAGGCGTGAGGCCGGTGTAGGCGCGGGTTTACCCGCGAATGGCCTCACCGCGGTAGCGGCAGGCTAGTGGTGGACTTGATCTCGGAAAgcg                                                                                         >  2:1617194/1‑126 (MQ=255)
                                                                                     ttCGCGGGTAAACCCGCTCCTACAAAGGTTCGCGCCAGGCGTGAGGCCGGTGTAGGCGCGGGTTTACCCGCGAATGGCCTCACCGCGGTAGCGGCAGGCTAGTGGTGGACTTGATCTCGGAAAgcg                                                                                         <  1:1617194/126‑1 (MQ=255)
                                                                                        gcgGGTAAACCCGCTCCTACAAAGGTTCGCGCCAGGCGTGAGGCCGGTGTAGGCGCGGGTTTACCCGCGAATGGCCTCACCGCGGTAGCGGCAGGCTAGTGGTGGACTTGATCTCGGAAAGCGCCACAATCGAGTTCACCTCCTGAATcc                                                              >  2:1177945/1‑150 (MQ=255)
                                                                                          gggTAAACCCGCTCCTACAAAGGTTCGCGCCAGGCGTGAGGCCGGTGTAGGCGCGGGTTTACCCGCGAATGGCCTCACCGCGGTAGCGGCAGGCTAGTGGTGGACTTGATCTCGGAAAGCGCCACAATCGAGTTCACCTCCTGAATcccc                                                            <  1:381783/150‑1 (MQ=255)
                                                                                           ggTAAACCCGCTCCTACAAAGGTTCGCGCCAGGCGTGAGGCCGGTGTAGGCGCGGGTTTACCCGCGAATGGCCTCACCGCGGTAGCGGCAGGCTAGTGGTGGACTTGATCTCGGAAAGCGCCACAATCGAGTTCACCTCCTGAATCCCCg                                                           >  1:1402394/1‑150 (MQ=255)
                                                                                            gtcaaCCCGCCCCTGCAAAGGTTCGCGCCAGACGTGAGGCCGGTGTAGGCGCGGGTTTACACGCGAAGGGCCGCACCGCGGTAGCGACAGGCTAGTGGTGGACTTGATCTCGGAAGGCTCCACAACCGAGTGCACATCCTGAATCCCCgg                                                          >  2:607218/4‑150 (MQ=255)
                                                                                               aaCCCGCTCCTACAAAGGTTCGCGCCAGGCGTGAGGCCGGTGTAGGCGCGGGTTTACCCGCGAATGGCCTCACCGCGGTAGCGGCAGGCTAGTGGTGGACTTGATCTCGGAAAGCGCCACAATCGAGTTCACCTCCTGAATCCCCGGCAc                                                       >  2:1266552/1‑150 (MQ=255)
                                                                                                 cccGCTCCTACAAAGGTTCGCGCCAGGCGTGAGGCCGGTGTAGGCGCGGGTTTACCCGCGAATGGCCTCACCGCGGTAGCGGCAGGCTAGTGGTGGCATGGCACTCGGAAAGCGCCACAATCGAGTTCACCTCCTGAATCCCCGGCACGt                                                     >  1:593541/1‑150 (MQ=255)
                                                                                                 cccGCTCCTACAAAGGTTCGCGCCAGGCGTGAGGCCGGTGTAGGCGCGGGTTTACCCGCGAATGGCCTCACCGCGGTAGCGGCAGGCTAGTGGTGGACTTGATCTCGGAAAGCGCCACAATCGCGTTCACCTCCTGAATCCCCCGCCCGt                                                     >  2:255203/1‑150 (MQ=255)
                                                                                                 cccGCTCCTACAAAGGTTCGCGCCAGGCGTGAGGCCGGGGTAGGCGCGGGTTTACCCGCGAATGGCCTCACTGCGGTAGCGGCAGGCTAGTGGTGGACTTGATCTCGGAAAGCGCAACAATCGAGTTAACATCCTGAATCCCCGGCACGt                                                     >  2:461146/1‑150 (MQ=255)
                                                                                                   cGCTCCTACAAAGGTTCGCGCCAGGCGTGAGGCCGGTGTAGGCGCGGGTTTACCCGCGAATGGCCTCACCGCGGTAGCGGCAGGCTAGTGGTGGACTTGATCTCGGAAAGCGCCACAATCGAGTTCACCTCCTGAATCCCCGGCACGTTc                                                   <  1:1111823/150‑1 (MQ=255)
                                                                                                         tACAAAGGTTCGCGCCAGGCGTGAGGCCGGTGTAGGCGCGGGTTTACCCGCGAATGGCCTCACcg                                                                                                                                  >  2:1793801/1‑65 (MQ=255)
                                                                                                         tACAAAGGTTCGCGCCAGGCGTGAGGCCGGTGTAGGCGCGGGTTTACCCGCGAATGGCCTCACcg                                                                                                                                  <  1:1793801/65‑1 (MQ=255)
                                                                                                         tACAAAGGTTCGCGCCAGGCGTGAGGCCGGTGTAGGCGCGGGTTTACCCGCGAATGGCCTCACCGCGGTAGCGGCAGGCTAGTGGTGGACTTGATCTCGGAAAGCGCCACAATCGAGTTCACCTCCTGAATCCCCGGCACGTTCGACAGt                                             >  1:918201/1‑150 (MQ=255)
                                                                                                            aaaGGTTCGCGCCAGGCGTGAGGCCGGTGTAGGCGCGGGTTTACCCGCGAATGGCCTCACCGCGGTAGCGGCAGGCTAGTGGTGGACTTGATCTCGGAAAGCGCCACAATCGAGTTCACCTCCTGAATCCCCGGCACGTTCGACAGtttt                                          >  2:158785/1‑150 (MQ=255)
                                                                                                                 ttCGCGCCAGGCGGGAGGCCGGTGTAGGCGCGGGTTTACCCGCGAATGGCCTCACCGCGGTAGCGGCAGGCTAGTGGTGGACTTGATCTCGGAAAGCGCCACAATCGAGTTCACTTCCTGAATCCCCGGCACGTTCGACAGCTTCTCgac                                     <  2:1014613/150‑2 (MQ=255)
                                                                                                                        cAGGCGTGAGGCCGGTGTAGGCGCGGGTTTACCCGCGAATGGCCTCACCGCGGTAGCGGCAGGCTAGTGGTGGACTTGATCTCGGAAAGCGCCACAATCGAGTTCACCTCCTGAATCCCCGGCACGTTCGACAGTTTTTCGAAGAAGAAc                              <  2:715470/150‑1 (MQ=255)
                                                                                                                              tGAGGCCGGTGTAGGCGCGGGTTTACCCGCGAATGGCCTCACCGCGGTAGCGGCAGGCTAGTGGTGGACTTGATCTCGGAAAGCGCCACAATCGAGTTCAcctcct                                                                    <  1:627511/106‑1 (MQ=255)
                                                                                                                              tGAGGCCGGTGTAGGCGCGGGTTTACCCGCGAATGGCCTCACCGCGGTAGCGGCAGGCTAGTGGTGGACTTGATCTCGGAAAGCGCCACAATCGAGTTCAcctcct                                                                    >  2:627511/1‑106 (MQ=255)
                                                                                                                                aGGCCGGTGTAGGCGCGGGTTTACCCGCGAATGGCCTCACCGCGGTAGCGGCAGGCTAGTGGTGGACTTGATCTCGGAAAGCGCCACAATCGAGTTCACCTCCTGAATCCCCGGCACGTTTGACAGTTTTTCGAAGAAGAACCGCTCGTa                      <  1:923467/150‑1 (MQ=255)
                                                                                                                                  gCCGGTGTAGGCGCGGGTTTACCCGCGAATGGCCTCACCGCGGTAGCGGCAGGCTAGTGGTGGACTTGATCTCGGAAAGCGCCACAATCGAGTTCACCTCCTGAATCCCCGGCACGTTCGACAGTTTTTCGAAGAAGAACCGCTCGTAgg                    >  2:1904306/1‑150 (MQ=255)
                                                                                                                                  gCCGGTGTAGGCGCGGGTTTACCCGCGAATGGCCTCACCGCGGTAGCGGCAGGCTAGTGGTGGACTTGATCTCGGAAAGCGCCACAATCGAGTTCACCTCCTGAATCCCCGGCACGTTCGACAGTTTTTCGAAGAAGAACCGCTCGTAgg                    >  2:225941/1‑150 (MQ=255)
                                                                                                                                   ccGGTGTAGGCGCGGGTTTACCCGCGAATGGCCTCACCGCGGTAGCGGCAGGCTAGTGGTGGACTTGATCTCGGAAAGCGCCACAATCGAGTTCACCTCCTGAATCCCCGGCACGTTCGACAGTTTTTCGAAGAAGAACCGCTCGTAGGc                   >  1:841960/1‑150 (MQ=255)
                                                                                                                                         tAGGCGCGGGTTTACCCGCGAATGGCCTCACCGCGGTAGCGGCAGGCTAGTGGTGGACTTGATCTCGGAAAGCGCCACAATCGAGTTCACCTCCTGAATCCCCGGCACGTTCGACAGTTTTTCGAAGAAGAACCGCTCGTAGGCCTCGAt             <  1:640495/150‑1 (MQ=255)
                                                                                                                                         cattcgcgGGTAAACCCGCGAATGGCCTCACCGCGGTAGCGGCAGGCTAGTGGTGGACTTGATCTCGGAAAGCGCCACAATCGAGTTCACCTCCTGAATCCCCGGCACGTTCGACAGTTTTTCGAAGAAGAACCGCTCGTAGGCCTCGAt             >  1:1364347/5‑150 (MQ=255)
                                                                                                                                          attcgcgGGTAAACCCGCGAATGGCCTCACCGCGGTAGCGGCAGGCTAGTGGTGGACTTGATCTCGGAAAGCGCCACAATCGAGTTCACCTCCTGAATCCCCGGCACGTTc                                                   <  2:507435/108‑1 (MQ=255)
                                                                                                                                          attcgcgGGTAAACCCGCGAATGGCCTCACCGCGGTAGCGGCAGGCTAGTGGTGGACTTGATCTCGGAAAGCGCCACAATCGAGTTCACCTCCTGAATCCCCGGCACGTTc                                                   >  1:507435/4‑111 (MQ=255)
                                                                                                                                          attcgcgGGTAAACCCGCGAATGGCCTCACCGCGGTAGCGGCAGGCTAGTGGTGGACTTGATCTCGGAAAGCGCCACAATCGAGTTCACCTCCTGAATCCCCGGCACGTTCGACAGTTTTTCga                                      >  2:965352/4‑124 (MQ=255)
                                                                                                                                          attcgcgGGTAAACCCGCGAATGGCCTCACCGCGGTAGCGGCAGGCTAGTGGTGGACTTGATCTCGGAAAGCGCCACAATCGAGTTCACCTCCTGAATCCCCGCCACGTTCGACAGTTTTTCga                                      <  1:965352/121‑1 (MQ=255)
                                                                                                                                               cgGGTTTACCCGCGAATGGCCTCACCGCGGTAGCGGCAGGCTAGTGGTGGATTTTAGCTCGGCCAGCGCCACCATTGCACTGCCTTCCAGAATCTCCGTCACGTTCGACAGTTTTTCGAAGAAGAACCGCTCGTAGGCCGCGAGGTCCGa       >  1:1382899/1‑150 (MQ=255)
                                                                                                                                                gggTTTACCCGCGAATGGCCTCACCGCGGTAGCGGCAGGCTAGTGGTGGACTTGATCTCGGAAAGCGCCACAATCGAGTTCACCTCCTGAATCCCCGGCACGTTCGACAGTTTTTCGAAGAAGAACCGCTCGTAGGCCTCGATGTCCGAc      >  1:911148/1‑150 (MQ=255)
                                                                                                                                                 ggTTTACCCGCGAATGGCCTCACCGCGGTAGCGGCAGGCTAGTGGTGGACTTGATCTCGGAAAGCGCCACAATCGAGTTCACCTCCTGAATCCCCGGCACGTTCGACAGTTTTTCGAAGAAGAACCGCTCGTAGGCCTCGATGTCCGACg     >  2:245727/1‑150 (MQ=255)
                                                                                                                                                  gTTTACCCGCGAATGGCCTCACCGCGGTAGCGGCAGGCTAGTGGTGGACTTGATCTCGGAAAGCGCCACAATCGAGTTCACCTCCTGAATCCCCGGCACGTTCGACAGTTTTTCGAAGAAGAACCGCTCGTAGGCCTCGATGTCCGACGt    >  1:867611/1‑150 (MQ=255)
                                                                                                                                                   tttACCCGCGAATGGCCTCACCGCGGTAGCGGCAGGCTAGTGGTGGACTTGATCTCGGAAAGCGCCACAATCGAGTTCACCTCCTGAATCCCCGGCACGTTc                                                   >  1:767619/1‑102 (MQ=255)
                                                                                                                                                   tttACCCGCGAATGGCCTCACCGCGGTAGCGGCAGGCTAGTGGTGGACTTGATCTCGGAAAGCGCCACAATCGAGTTCACCTCCTGAATCCCCGGCACGTTc                                                   <  2:767619/102‑1 (MQ=255)
                                                                                                                                                    ttACCCGCGAATGGCCTCACCGCGGTAGCGGCAGGCTAGTGGTGGACTTGATCTCGGAAAGCGCCACAATCGAGTTCACCTCCTGAATCCCCGGCACGTTCGACAGTTTTTCGAAGAAGAACCGCTCGTAGGCCTCGATGTCCGACGTGa  <  1:245727/150‑1 (MQ=255)
                                                                                                                                                    ttACCCGCGAATGGCCTCACCGCGGTAGCGGCAGGCTAGTGGTGGACTTGATCTCGGAAAGCGCCACAATCGAGTTCACCTCCTGAATCCCCGGCACGTTCGACAGTTTTTCGAAGAAGAACCGCTCGTAGGCCTCGATGTCCGACGTGa  <  1:1626249/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GGAAGCGCGTATCGAACAGGCTCGGGCGGTGGACCAGAGCCGCGGCTGGAAGCACCTGCTCAAGGCGTGATACCGCGTCGCCTGCTTCGCGGGTAAACCCGCTCCTACAAAGGTTCGCGCCAGGCGTGAGGCCGGTGTAGGCGCGGGTTTACCCGCGAATGGCCTCACCGCGGTAGCGGCAGGCTAGTGGTGGACTTGATCTCGGAAAGCGCCACAATCGAGTTCACCTCCTGAATCCCCGGCACGTTCGACAGTTTTTCGAAGAAGAACCGCTCGTAGGCCTCGATGTCCGACGTGA  >  NC_002947/5918297‑5918594

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