Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,061,559:1 +C 7.7% intergenic (‑136/+52) hupA ← / ← alkT DNA‑binding protein HU‑alpha/rubredoxin‑NAD(+) reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,061,5591.C 7.7% 309.9 / 19.6 78intergenic (‑136/+52)hupA/alkTDNA‑binding protein HU‑alpha/rubredoxin‑NAD(+) reductase
Reads supporting (aligned to +/- strand):  ref base . (39/33);  new base C (2/4);  total (41/37)
Fisher's exact test for biased strand distribution p-value = 4.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.26e-01

CTGGTTTGCGCATAGTGAAGCCTCTTTGACGGAATTCTTGTTGTTATGCCGTGCTGCCTCTGGGAGCAGCGCTCAAGGCACCGCAGGCTCTAATCTGCGGCAGACGGGAGTGAGAATGGCACGAGCTCAGGGGCCGCGCCAGTATCTGGGCGGCCATTGTGCAGGCAACAACAGGATAAATCCGACAGAACGCGCGCTATTTACGCCATAACCGGTGGCAACTGACGATTCAGCGCCAGCTTCTCCATCACCGCCGCTCCGGTCAGGGCGTAGCCCAGCAGCTTGCCGTCGGCAT  >  NC_002947/6061411‑6061705
                                                                                                                                                      |                                                                                                                                                  
cTGGTTTGCGCATAGTGAAGCCTCTTTGACGGAATTCTTGTTGTTATGCCGTGCTGCCTCTGGGAGCAGCGCTCAAGGCACCGCAGGCTCTAATCTGCGGCAGACGGGAGTGAGAATGGCACGAGCTCAGGGGCCGCGCCAGTATCTggg                                                                                                                                                   <  2:1253013/150‑1 (MQ=255)
 tGGTTTGCGCATAGTGAAGCCTCTTTGACGGAATTCTTGTTGTTATGCCGTGCTGCCTCTGGGAGCAGCGCTCAAGGCACCGCAGGCTCTAATCTGCGGCAGACGGGAGTGAGAATGGCACGAGCTCAGGGGCCGCGCCAGTATCTGggc                                                                                                                                                  <  1:1061623/150‑1 (MQ=255)
    tttGCGCATAGTGAAGCCTATTTGACGGAATTCTTGTTGTTATGCCGTGCTGCCTCTGGGAGCAGCGCTCAAGGCACCGCAGGCTCTAATCTGCGGCAGACGGGAGTGAGAATGGCACGAGCTCAGGGGCCGCGCCAGTATCTGggcggc                                                                                                                                               <  2:842994/150‑1 (MQ=255)
        ctgggctacgccctgaccggagcggcggtgatggagaagctggcgctgaatcgtcagttgccaccggttatggcgtaaatagcgcgcgttctgtcggatttatcctgttgttgcctgcacaatGGCCGCGCCAGTATCTGGGCGGCCAtt                                                                                                                                           >  2:809256‑M1/124‑150 (MQ=255)
         gcATAGTGAAGCCTCTTTGACGGAATTCTTGTTGTTATGCCGTGCTGCCTCTGGGAGCAGCGCTCAAGGCACCGCAGGCTCTAATCTGCGGCAGACGGGAGTGAGAATGGCACGAGCTCAGGGGCCGCGCCAGTATCTGGGCGGCCATtg                                                                                                                                          >  1:604063/1‑150 (MQ=255)
              gTGAAGCCTCTTTGACGGAATTCTTGTTGTTATGCCGTGCTGCCTCTGGGAGCAGCGCTCAAGGCACCGCAGGCTCTAATCTGCGGCAGACGGGAGTGAGAATGGCACGAGCTCAGGGGCCGCGCCAGTATCTGGGCGGCCATTGTGCAg                                                                                                                                     <  2:1194620/150‑1 (MQ=255)
                       ctTTGACGGAATTCTTGTTGTTATGCCGTGCTGCCTCTGGGGGCAGCGCTCAAGGCACCGCAGGCTCTAATCTGCGGCAGACGGGAGTGAGAATGGCACGAGCTCAGGGGCCGCGCCAGTATCTGGGCGGCCATTGTGCAGGCAacaaca                                                                                                                            <  2:607828/150‑1 (MQ=255)
                        tttGACGGAATTCTTGTTGTTATGCCGTGCTGCCTCTGGGAGCAGCGCTCAAGGCACCGCAGGCTCTAATCTGCGGCAGACGGGAGTGAGAATGGCACGAGCTCAGGGGCCGCGCCAGTATCTGGGCGGCCATTGTGCAGGCAACAACAg                                                                                                                           >  1:1000932/1‑150 (MQ=255)
                          tGACGGAATTCTTGTTGTTATGCCGTGCTGCCTCTGGGAGCAGCGCTCAAGGCACCGCAGGCTCTAATCTGCGGCAGACGCGAGTGAGAATGGCACGAGCTCAGGGGCCGCGCCAGTATCTGGGCGGCCATTGTGCAGGCAACAACAGGa                                                                                                                         <  2:535553/150‑1 (MQ=255)
                            aCGGAATTCTTGTTGTTATGCCGTGCTGCCTCTGGGAGCAGCGCTCAAGGCACCGCAGGCTCTAATCTGCGGCAGACGGGAGTGAGAATGGCACGAGCTCAGGGGCCGCGCCAGTATCTGGGCGGCCATTGTGCAGGCAACAACAGGATa                                                                                                                       >  2:429922/1‑150 (MQ=255)
                                     ttgttgTTATGCCGTGCTGCCTCTGGGAGCAGCGCTCAAGGCACCGCAGGCTCTAATCTGCGGCAGACGGGAGTGAGAATGGCACGAGCTCAGGGGCCGCGCCAGTATCTGGGCGGCCATTGTGCAGGCAACAACAGGATAAATCCGACa                                                                                                              <  2:604063/150‑1 (MQ=255)
                                     ttgttgTTATGCCGTGCTGCCTCTGGGAGCAGCGCTCAAGGCACCGCAGGCTCTAATCTGCGGCAGACGGGAGTGAGAATGGCACGAGCTCAGGGGCCGCGCCAGTATCTGGGCGGCCATTGTGCAGGCAACAACAGGATAAATCCGACa                                                                                                              >  2:1161934/1‑150 (MQ=255)
                                      tgttgtTATGCCGTGCTGCCTCTGGGAGCAGCGCTCAAGGCACCGCAGGCTCTAATCTGCGGCAGACGGGAGTGAGAATGGCACGAGCTCAGGGGCCGCGCCAGTATCTGGGCGGCCATTGTGCAGGCAACAACAGGATAAATCCGACAg                                                                                                             >  1:934698/1‑150 (MQ=255)
                                         tgttATGCCGTGCTGCCTCTGGGAGCAGCGCTCAAGGCACCGCAGGCTCTAATCTGCGGCAGACGGGAGTGAGAATGGCACGAGCTCAGGGGCCGCGCCAGTATCTGGGCGGCCATTGTGCAGGCAACAACAGGATAAATCCGACAGAAc                                                                                                          <  2:934698/150‑1 (MQ=255)
                                           ttATGCCGTGCTGCCTCTGGGAGCAGCGCTCAAGGCACCGCAGGCTCTAATCTGCGGCAGACGGGAGTGAGAATGGCACGAGCTCAGGGGCCGCGCCAGTATCTGGGCGGCCATTGTGCAGGCAACAACAGGATAAATCCGACAGAAcgc                                                                                                        >  1:99678/1‑150 (MQ=255)
                                             aTGCCGTGCTGCCTCTGGGAGCAGCGCTCAAGGCACCGCAGGCTCTAATCTGCGGCAGACGGGAGTGAGAATGGCACGAGCTCAGGGGCCGCGCCAGTATCTGGGCGGCCATTGTGCAGGCAACAACAGGATAAATCCGACAGAAcgcgc                                                                                                      <  2:1000932/150‑1 (MQ=255)
                                              tGCCGTGCTGCCTCTGGGAGCAGCGCTCAAGGCACCGCAGGCTCTAATCTGCGGCAGACGGGAGTGAGAATGGCACGAGCTCAGGGGCCGCGCCAGTATCTGGGCGGCGATTGTGCAGGCAACAGCGGGATAAATCCGCCAGAAcgcgca                                                                                                     >  2:271787/1‑149 (MQ=255)
                                              tGCCGTGCTGCCTCTGGGAGCAGCGCTCAAGCCACCGCAGGCTCTAATCTGCGGCAGACGGGAGTGAGAATGGCACGAGCTCAGGGGCCGCGCCAGTATCTGGGCGGCCATTGTGCAGGCAACAACAGGATAAATCCGACAGAAcgcgcg                                                                                                     >  2:790699/1‑150 (MQ=255)
                                               gCCGTGCTGCCTCTGGGAGCAGCGCTCAAGGCACCGCAGGCTCTAATCTGCGGCAGACGGGAGTGAGAATGGCACGAGCTCAGGGGCCGCGCCAGTATCTGGGCGGCCATTGTGCAGGCAACAACAGGATAAATCCGACAGAACgcgcgc                                                                                                    <  1:738249/150‑1 (MQ=255)
                                                     ctgcCTCTGGGAGCAGCGCTCAAGGCACCGCAGGCTCTAATCTGCGGCAGACGGGAGTGAGAATGGCACGAGCTCAGGGGCCGCGCCAGTATCTGGGCGGCCATTGTGCAGGCAACAACAGGATAAATCCGACAGAACGCGCGCTATTTa                                                                                              >  2:883793/1‑150 (MQ=255)
                                                       gcCTCTGGGAGCAGCGCTCAAGGCACCGCAGGCTCTAATCTGCGGCAGACGGGAGTGAGAATGGCACGAGCTCAGGGGCCGCGCCAGTATCTGGGCGGCCATTGTGCAGGCAACAACAGGATAAATCCGACAGAACGCGCGCTATTTACg                                                                                            >  1:427230/1‑150 (MQ=255)
                                                       gcCTCTGGGAGCAGCGCTCAAGGCACCGCAGGCTCTAATCTGCGGCAGACGGGAGTGAGAATGGCACGAGCTCAGGGGCCGCGCCAGTATCTGGGCGGCCATTGTGCAGGCAACAACAGGATAAATCCGACAGAACGCGCGCTATTTACg                                                                                            >  2:47062/1‑150 (MQ=255)
                                                          tctGGGAGCAGCGCTCAAGGCACCGCAGGCTCTAATCTGCGGCAGACGGGAGTGAGAATGGCACGAGCTCAGGGGCCGCGCCAGTATCTGGGCGGCCATTGTGCAGGCAACAACAGGATAAATCCGa                                                                                                                <  1:77763/127‑1 (MQ=255)
                                                          tctGGGAGCAGCGCTCAAGGCACCGCAGGCTCTAATCTGCGGCAGACGGGAGTGAGAATGGCACGAGCTCAGGGGCCGCGCCAGTATCTGGGCGGCCATTGTGCAGGCAACAACAGGATAAATCCGa                                                                                                                >  2:77763/1‑127 (MQ=255)
                                                             gggAGCAGCGCTCAAGGCACCGCAGGCTCTAATCTGCGGCAGACGGGAGTGAGAATGGCACGAGCTCAGGGGCCGCGCCAGTATCTGGGCGGCCATTGTGCAGGCAACAACAGGATAAATCCGACAGAACGCGCGCTATTTACGCCATaa                                                                                      >  1:127672/1‑150 (MQ=255)
                                                                 gcagcGCTCAAGGCACCGCAGGCTCTAATCTGCGGCAGACGGGAGTGAGAATGGCACGAGCTCAGGGGCCGCGCCAGTATCTGGGCGGCCATTGTGCAGGCAACAACAGGATAAATCCGACAGAACGCGCGCTATTTACGCCATAACCgg                                                                                  <  1:429922/150‑1 (MQ=255)
                                                                   agcGCTCAAGGCACCGCAGGCTCTAATCTGCGGCAGACGGGAGTGAGAATGGCACGAGCTCAGGGGCCGCGCCAGTATCTGGGCGGCCATTGTGCAGGCAACAACAGGATAAATCCGACAGAACGCGCGCTATTTACGCCATAACCGGTg                                                                                >  2:939251/1‑150 (MQ=255)
                                                                   agcGCTCAAGGCACCGCAGGCTCTAATCTGCGGCAGACGGGAGTGAGAATGGCACGAGCTCAGGGGCCGCGCCAGTATCTGGGCGGCCATTGTGCAGGCAACAACAGGATAAATCCGACAGAACGCGCGCTATTTACGCCATAACCGGTg                                                                                >  2:53669/1‑150 (MQ=255)
                                                                   agcGCTCAAGGCACCGCAGGCTCTAATCTGCGGCAGACGGGAGTGAGAATGGCACGAGCTCAGGGGCCGCGCCAGTATCTGGGCGGCCATTGTGCAGGCAACAACAGGATAAATCCGACAGAACGCGCGCTATTTACGCCATAACCGGTg                                                                                >  1:334231/1‑150 (MQ=255)
                                                                          aaGGCACCGCAGGCTCTAATCTGCGGCAGACGGGAGTGAGAATGGCACGAGCTCAGGGGCCGCGCCAGTATCTGGGCGGCCATTGTGCAGGCAACAACAGGATAAATCCGa                                                                                                                >  2:1090001/1‑111 (MQ=255)
                                                                          aaGGCACCGCAGGCTCTAATCTGCGGCAGACGGGAGTGAGAATGGCACGAGCTCAGGGGCCGCGCCAGTATCTGGGCGGCCATTGTGCAGGCAACAACAGGATAAATCCGa                                                                                                                <  1:1090001/111‑1 (MQ=255)
                                                                          aaGGCACCGCAGGCTCTAATCTGCGGCAGACGGGAGTGAGAATGGCACGAGCTCAGGGGCCGCGCCAGTATCTGGGCGGCCATTGTGCAGGCAACAACAGGATAAATCCGACAGAACGCGCGCTATTTACGCCATAACCGGTGGCAACTg                                                                         >  1:839699/1‑150 (MQ=255)
                                                                          aaGGCACCGCAGGCTCTAATCTGCGGCAGACGGGAGTGAGAATGGCACGAGCTCAGGGGCCGCGCCAGTATCTGGGCGGCCATTGTGCAGGCAACAACAGGATAAATCCGACAGAACGCGCGCTATTTACGCCATAACCGGTGGCAACTg                                                                         >  2:494129/1‑150 (MQ=255)
                                                                            ggCACCGCAGGCTCTAATCTGCGGCAGACGGGAGTGAGAATGGCACGAGCTCAGGGGCCGCGCCAGTATCTGGGCGGCCATTGTGCAGGCAACAACAGGATAAATCCGACAGAACGCGCGCTATTTACGCCATAACCGGTGGCAACTGAc                                                                       <  1:47062/150‑1 (MQ=255)
                                                                            ggCACCGCAGGCTCTAATCTGCGGCAGACGGGAGTGAGAATGGCACGAAATCATGTGCCGCGCCAGTATCTGGGCGGCCATTGTGCAGGCAACAACAGGATAAATCCGACAGAACGCGCGCTATTTACGCCATAACCGGTGGCAACTGAc                                                                       >  2:179857/1‑150 (MQ=255)
                                                                             gCACCGCAGGCTCTAATCTGCGGCAGACGGGAGTGAGAATGGCACGAGCTCAGGGGCCGCGCCAGTATCTGGGCGGCCATTGTGCAGGCAACAACAGGATAAATCCGACAGAACGCGCTCTATTTACGCCATAACCg                                                                                   >  2:1239912/1‑137 (MQ=255)
                                                                             gCACCGCAGGCTCTAATCTGCGGCAGACGGGAGTGAGAATGGCACGAGCTCAGGGGCCGCGCCAGTATCTGGGCGGCCATTGTGCAGGCAACAACAGGATAAATCCGACAGAACGCGCTCTATTTACGCCATAACCg                                                                                   <  1:1239912/137‑1 (MQ=255)
                                                                              cACCGCAGGCTCTAATCTGCGGCAGACGGGAGTGAGAATGGCACGAGCTCAGGGGCCGCGCCAGTATCTGGGCGGCCATTGTGCAGGCAACAACAGGATAAATCCGACAGAACGCGCGCTATTTACGCCATAACCGGTGGCAACTGACGa                                                                     <  2:1159446/150‑1 (MQ=255)
                                                                              cACCGCAGGCTCTAATCTGCGGCAGACGGGAGTGAGAATGGCACGAGCTCAGGGGCCGCGCCAGTATCTGGGCGGCCATTGTGCAGGCAACAACAGGATAAATCCGACAGAACGCGCGCTATTTACGCCATAACCGGTGGCAACTGACGa                                                                     <  1:53669/150‑1 (MQ=255)
                                                                                   cAGGCTCTAATCTGCGGCAGACGGGAGTGAGAATGGCACGAGCTCAGGGGCCGCGCCAGTATCTGGGCGGCCATTGTGCAGGCAACAACAGGAt                                                                                                                        >  2:533343/1‑94 (MQ=255)
                                                                                   cAGGCTCTAATCTGCGGCAGACGGGAGTGAGAATGGCACGAGCTCAGGGGCCGCGCCAGTATCTGGGCGGCCATTGTGCAGGCAACAACAGGAt                                                                                                                        <  1:533343/94‑1 (MQ=255)
                                                                                      gCTCTAATCTGCGGCAGACGGGAGTGAGAATGGCACGAGCTCAGGGGCCGCGCCAGTATCTGGGCGGCCATTGTGCAGGCAACAACAGGATAAATCCGACAGAACGCGCGCTATTTACGCCATAACCGGTGGCAACTGACGATTCAgcgc                                                             <  2:334231/150‑1 (MQ=255)
                                                                                        tctAATCTGCGGCAGACGGGAGTGAGAATGGCACGAGCTCAGGGGCCGCGCCAGTATCTGGGCGGCCATTGTGCAGGCAACAACAGGATAAATCCGACAGAACGCGCGCTATTTACGCCATAACCGGTGGCAACTGACGATTCAGCGCCa                                                           <  2:839699/150‑1 (MQ=255)
                                                                                        tctAATCTGCGGCAGACGGGAGTGAGAATGGCACGAGCTCAGGGGCCGCGCCAGTATCTGGGCGGCCATTGTGCAGGCAACAACAGGATAAATCCGACAGAACGCGCGCTATTTACGCCATAACCGGTGGCAACTGACGATTCAGCGCCa                                                           >  2:138025/1‑150 (MQ=255)
                                                                                          tAATCTGCGGCAGACGGGAGTGAGAATGGCACGAGCTCAGGGGCCGCGCCAGTATCTGGGCGGCCATTGTGCAGGCAACAACAGGATAAATCCGACAGAACGCGCGCTATTTACGCCATAACCGGTGGCAACTGACGATTCAGCGCCAGc                                                         <  1:494129/150‑1 (MQ=255)
                                                                                            aTCTGCGGCAGACGGGAGTGAGAATGGCACGAGCTCAGGGGCCGCGCCAGTATCTGGGCGGCCATTGTGCAGGCAACAACAGGATAAATCCGACAGAACGCGCGCTATTTACGCCATAACCGGTGGCAACTGa                                                                        >  2:1044717/1‑133 (MQ=255)
                                                                                            aTCTGCGGCAGACGGGAGTGAGAATGGCACGAGCTCAGGGGCCGCGCCAGTATCTGGGCGGCCATTGTGCAGGCAACAACAGGATAAATCCGACAGAACGCGCGCTATTTACGCCATAACCGGTGGCAACTGa                                                                        <  1:1044717/133‑1 (MQ=255)
                                                                                               tGCGGCAGACGGGAGTGAGAATGGCACGAGCTCAGGGGCCGCGCCAGTATCTGGGCGGCCATTGTGCAGGCAACAACAGGATAAATCCGACAGAACGCGCGCTATTTACGCCATAAc                                                                                     <  1:886628/117‑1 (MQ=255)
                                                                                               tGCGGCAGACGGGAGTGAGAATGGCACGAGCTCAGGGGCCGCGCCAGTATCTGGGCGGCCATTGTGCAGGCAACAACAGGATAAATCCGACAGAACGCGCGCTATTTACGCCATAAc                                                                                     >  2:886628/1‑117 (MQ=255)
                                                                                                  ggCAGACGGGAGTGAGAATGGCACGAGCTCAGGGGCCGCGCCAGTATCTGGGCGGCCATTGTGCAGGCAACAACAGGATAAATCCGACAGAACGCGCGCTATTTACGCCATAACCGGTGGCAACTGACGATTCAGCGCCAGCTTCTCCAt                                                 <  1:138025/150‑1 (MQ=255)
                                                                                                   gCAGACGGGAGTGAGAATGGCACGAGCTCAGGGGCCGCGCCAGTATCTGGGCGGCCATTGTGCAGGCAACAACAGGATAAATCCGACAGAACGCGCGCTATTTACGCCATAACCGGTGGCAACTGACGATTCAGCGCCAGCTTCTCCATc                                                >  1:1169888/1‑150 (MQ=255)
                                                                                                     aGACGGGAGTGAGAATGGCACGAGCTCAGGGGCCGCGCCAGTATCTGGGCGGCCATTGTGCAGGCAACAACAGGATAAATCCGACAGAACGCGCGCTATTTACGCCATAACCGGTGGCAACTGACGATTCAGCGCCAGCTTCTCCATCAc                                              <  2:508970/150‑1 (MQ=255)
                                                                                                      gACGGGAGTGAGAATGGCACGAGCTCAGGGGCCGCGCCAGTATCTGGGCGGCCATTGTGCAGGCAACAACAGGATAAATCCGACAGAACGCGCGCTATTTACGCCATAACCGGTGGCAAc                                                                           >  2:865471/1‑120 (MQ=255)
                                                                                                      gACGGGAGTGAGAATGGCACGAGCTCAGGGGCCGCGCCAGTATCTGGGCGGCCATTGTGCAGGCAACAACAGGATAAATCCGACAGAACGCGCGCTATTTACGCCATAACCGGTGGCAAc                                                                           <  1:865471/120‑1 (MQ=255)
                                                                                                       aCGGGAGTGAGAATGGCACGAGCTCAGGGGCTGCGCCAGGATCTGGGCGGCCATTGTGCAGGCAACAACAGGATAAATCCGACAGAACGCGCGCTATTTACGCCATAACCGGTGGCAACTGACGATTCAGCGCCAGCTTCTCCATCAccg                                            >  1:102536/1‑150 (MQ=255)
                                                                                                       aCGGGAGTGAGAATGGCACGAGCTCAGGGGCCGCGCCAGTATCTGGGCGGCCATTGTGCAGGCAACAACAGGATAAATCCGACAGAACGCGCGCTATTTACGCCATAACCGGTGGCAACTGACGATTCAGCGCCAGCTTCTCCATCAccg                                            >  2:485410/1‑150 (MQ=255)
                                                                                                           ttatcctgttgttgcctgcacaatGGCCGCCCAGATACTGGCGCGGCCATTGTGCAGGCAACAACAGGATAAATCCGACAGAACGCGCGCTATTTACGCCATAACCGGTGGCAACTGACGATTCAGCGCCAGCTTCTCCATCACcgccgc                                        <  1:794480‑M1/126‑1 (MQ=255)
                                                                                                            aGTGAGAATGGCACGAGCTCAGGGGCCGCGCCAGTATCTGGGCGGCCATTGTGCAGGCAACAACAGGATAAATCCGACAGAACGCGCGCTATTTACGCCATAACCGGTGGCAACTGACGATTCAGCGCCAGCtt                                                       >  2:370329/1‑134 (MQ=255)
                                                                                                            aGTGAGAATGGCACGAGCTCAGGGGCCGCGCCAGTATCTGGGCGGCCATTGTGCAGGCAACAACAGGATAAATCCGACAGAACGCGCGCTATTTACGCCATAACCGGTGGCAACTGACGATTCAGCGCCAGCtt                                                       <  1:370329/134‑1 (MQ=255)
                                                                                                              tGAGAATGGCACGAGCTCAGGGGCCGCGCCAGTATCTGGGCGGCCATTGTGCAGGCAACAACAGGATAAAt                                                                                                                    >  2:466080/1‑71 (MQ=255)
                                                                                                              tGAGAATGGCACGAGCTCAGGGGCCGCGCCAGTATCTGGGCGGCCATTGTGCAGGCAACAACAGGATAAAt                                                                                                                    <  1:466080/71‑1 (MQ=255)
                                                                                                                  aaTGGCACGAGCTCAGGGGCCGCGCCAGTATCTGGGCGGCCATTGTGCAGGCAACAACAGGATAAATCCGACAGAACGCGCGCTATTTACGCCATAACCGGTGGCAACTGACGATTCAGCGCCAGCTTCTCCATCACCGCCGCTCCGGTc                                 <  1:485410/150‑1 (MQ=255)
                                                                                                                      gCACGAGCTCAGGGGCCGCGCCAGTATCTGGGCGGCCATTGTGCAGGCAACAACAGGATAAATCCGACAGAACGCGCGCTATTTACGCCATAACCGGTGGCAACTGACGATTCAGCGCCAGCTTCTCCATCACCGCCGCTCCGGTCAggg                             >  1:1273889/1‑150 (MQ=255)
                                                                                                                       cACGAGCTCAGGGGCCGCGCCAGTATCTGGGCGGCCATTGTGCAGGCAACAACAGGATAAATCCGACAGAACGCGCGCTATTTACGCCATAACCGGTGGCAACTGACGATTCAGCGCCAGCTTCTCCATCACCGCCGCTCCGGTCAGGGc                            >  1:615750/1‑150 (MQ=255)
                                                                                                                          gAGCTCAGGGGCCGCGCCAGTATCTGGGCGGCCATTGTGCAGGCAACAACAGGATAAATCCGACAGAACGCGCGCTATTTACGCCATAACCGGTGGCAACTGACGATTCAGCGCCAGCTTCTCCATCACCGCCGCTCCGGTCAGGGCGTa                         <  2:1273889/150‑1 (MQ=255)
                                                                                                                          ctgcacaatGGCCGCCCAGATACTGGCGCGGCCATTGTGCAGGCAACAACAGGATAAATCCGACAGAACGCGCGCTATTTACGCCATAACCGGTGGCAACTGACGATTCAGCGCCAGCTTCTCCATCACCGCCGCTCCGGTCAGGGCGTa                         <  1:462855‑M1/141‑1 (MQ=255)
                                                                                                                           tgcacaatGGCCGCCCAGATACTGGCGCGGCCATTGTGCAGGCAACAACAGGATAAATCCGACAGAACGCGCGCTATTTACGCCATAACCGGTGGCAACTGACGATTCAGCGCCAGCTTCTCCATCACCGCCGCTCCGGTCa                                <  1:770349‑M1/134‑1 (MQ=255)
                                                                                                                           tgcacaatGGCCGCCCAGATACTGGCGCGGCCATTGTGCAGGCAACAACAGGATAAATCCGACAGAACGCGCGCTATTTACGCCATAACCGGTGGCAACTGACGATTCAGCGCCAGCTTCTCCATCACCGCCGCTCCGGTCa                                >  2:770349‑M1/9‑142 (MQ=255)
                                                                                                                           aGCTCAGGGGCCGCGCCAGTATCTGGGCGGCCATTGTGCAGGCAACAACAGGATAAATCCGACAGAACGCGCGCTATTTACGCCATAACCGGTGGCAACTGACGATTCAGCGCCAGCTTCTCCATCACCGCCGCTCCGGTCAGGGCGTAg                        <  1:1161934/150‑1 (MQ=255)
                                                                                                                            gcacaatGGCCGCCCAGATACTGGCGCGGCCATTGTGCAGGCAACAACAGGATAAATCCGACAGAACGCGCGCTATTTACGCCATAACCGGTgg                                                                               <  2:610616‑M1/87‑1 (MQ=255)
                                                                                                                            gcacaatGGCCGCCCAGATACTGGCGCGGCCATTGTGCAGGCAACAACAGGATAAATCCGACAGAACGCGCGCTATTTACGCCATAACCGGTgg                                                                               >  1:610616‑M1/8‑94 (MQ=255)
                                                                                                                             cTCAGGGGCCGCGCCAGTATCTGGGCGGCCATTGTGCAGGCAACAACAGGATAAATCCGACAGAACGCGCGCTATTTACGCCATAACCGGTGGCAACTGACGATTCAGCGCCAGCTTCTCCATCACCGCCGCTCCGGTCAGGGCGTAGcc                      >  1:663354/1‑150 (MQ=255)
                                                                                                                               cAGGGGCCGCGCCAGTATCTGGGCGGCCATTGTGCAGGCAACAACAGGATAAATCCGACAGAACGCGCGCTATTTACGCCATAACCGGTGGCAACTGACGATTCAGCGCCAGCTTCTCCATCACCGCCGCTCCGGTCAGGGCGTAGCCca                    <  1:939251/150‑1 (MQ=255)
                                                                                                                                aGGGGCCGCGCCAGTATCTGGGCGGCCATTGTGCAGGCAACAACAGGATAAATCCGACAGAACGCGCGCTATTTACGCCATAACCGGTGGCAACTGa                                                                        >  2:17337/1‑97 (MQ=255)
                                                                                                                                aGGGGCCGCGCCAGTATCTGGGCGGCCATTGTGCAGGCAACAACAGGATAAATCCGACAGAACGCGCGCTATTTACGCCATAACCGGTGGCAACTGa                                                                        <  1:17337/97‑1 (MQ=255)
                                                                                                                                         gcCAGTATCTGGGCGGCCATTGTGCAGGCAACAACAGGATAAATCCGACAGAACGCGCGCTAtt                                                                                                <  1:344014/64‑1 (MQ=255)
                                                                                                                                         gcCAGTATCTGGGCGGCCATTGTGCAGGCAACAACAGGATAAATCCGACAGAACGCGCGCTAtt                                                                                                <  1:343925/64‑1 (MQ=255)
                                                                                                                                         gcCAGTATCTGGGCGGCCATTGTGCAGGCAACAACAGGATAAATCCGACAGAACGCGCGCTAtt                                                                                                >  2:344014/1‑64 (MQ=255)
                                                                                                                                         gcCAGTATCTGGGCGGCCATTGTGCAGGCAACAACAGGATAAATCCGACAGAACGCGCGCTAtt                                                                                                >  2:343925/1‑64 (MQ=255)
                                                                                                                                             gTATCTGGGCGGCCATTGTGCAGGCAACACCAGGATGAATCCGACAGAACGCGCGCTATTTACGCCATAACCGGTGGCAACTGACGATTCAGCGCCAGCTTCTCCATCACCGCCGCTCCGGTCAGGGCGTAGCCCAGCAGCTTGCCGTCg      <  2:102536/150‑1 (MQ=255)
                                                                                                                                                  cTGGGCGGCCATTGTGCAGGCAACAACAGGATAAATCCGACAGAACGCGCGCTATTTACGCCATAACCGGTGGCAACTGACGATTCAGCGCCAGCTTCTCCATCACCGCCGCTCCGGTCAGGGCGTAGCCCAGCAGCTTGCCGTCGGCAt  >  2:1286363/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                  
CTGGTTTGCGCATAGTGAAGCCTCTTTGACGGAATTCTTGTTGTTATGCCGTGCTGCCTCTGGGAGCAGCGCTCAAGGCACCGCAGGCTCTAATCTGCGGCAGACGGGAGTGAGAATGGCACGAGCTCAGGGGCCGCGCCAGTATCTGGGCGGCCATTGTGCAGGCAACAACAGGATAAATCCGACAGAACGCGCGCTATTTACGCCATAACCGGTGGCAACTGACGATTCAGCGCCAGCTTCTCCATCACCGCCGCTCCGGTCAGGGCGTAGCCCAGCAGCTTGCCGTCGGCAT  >  NC_002947/6061411‑6061705

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