Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,941,334 T→A 7.5% intergenic (+60/+121) yhhJ → / ← PP_5722 translation‑like protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,941,3340TA 7.5% 195.7 / 5.2 67intergenic (+60/+121)yhhJ/PP_5722translation‑like protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (28/34);  new base A (2/3);  total (30/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.46e-01

GCCGGATATCCTGGCGCTGACGTTGATCGGCCTGGTGTTCTTCAGCGTGGCGCTGGCACGGTTTCGTCGGAGCCTGGCGTCTTGAAGGCATGCGCGAGCGTCCGCCGAAAATTGGAAAACTGGCACAAAGACAAAGGGCGCCTTTCGGCGCCCTACTCACAGACAGGTCGACTTGTCATCGCCATATACCTGCCCGGCTCTGTAATTGCTGGTCGAACCAGGATCCTCAGAGTCTTACAAGCCCCTTGCGGGAACAGGCTGATTATCACGCATCGAGCGTAACATTGACAA  >  NC_002947/5941190‑5941480
                                                                                                                                                |                                                                                                                                                  
gCCGGATATCCTGGCGCTGACGTTGATCGGCCTGGTGTTCTTCAGCGTGGCGCTGGCACGGTTTCGTCGGAGCCTGGCGTCTTGAAGGCATGCGCGAGCGTCCGCCGAAAATTGGAAAACTGGCACAAAGACAAAGGGCGCCTTTCGgc                                                                                                                                                >  4:140314/1‑149 (MQ=255)
     atatCCTGGCGCTGACGTTGATCGGCCTGGTGTTCTTCAGCGTGGCGCTGGCACGGTTTCGTCGGAGCCTGGCGTCTTGAAGGCATGCGCGAGCGTCCGCCGAAAATTGGAAAACTGGCACAAAGACAAAGGGCGCCTTTCGGCGCCCt                                                                                                                                           <  7:46131/149‑1 (MQ=255)
        tCCTGGCGCTGACGTTGATCGGCCTGGTGTTCTTCAGCGTGGCGCTGGCACGGTTTCGTCGGAGCCTGGCGTCTTGAAGGCATGCGCGAGCGTCCGCCGAAAATTGGAAAACTGGCACAAAGACAAAGGGCGCCGAAAGGCGCCCttt                                                                                                                                         <  1:56655/148‑3 (MQ=255)
          cTGGCGCTGACGTTGATCGGCCTGGTGTTCTTCAGCGTGGCGCTGGCACGGTTTCGTCGGAGCCTGGCGTCTTGAAGGCATGCGCGAGCGTCCGCCGAAAATTGGAAAACTGGCACAAAGACAAAGGGCGCCTTTCGGCGCCCTACTca                                                                                                                                      >  4:290417/1‑149 (MQ=255)
          cTGGCGCTGACGTTGATCGGCCTGGTGTTCTTCAGCGTGGCGCTGGCACGGTTTCGTCGGAGCCTGGCGTCTTGAAGGCATGCGCGAGCGTCCGCCGAAAATTGGAAAACTGGCACAAAGACAAAGGGCGCCTTTCGGCGCCCTACTca                                                                                                                                      <  3:191877/149‑1 (MQ=255)
          cTGGCGCTGACGTTGATCGGCCTGGTGTTCTTCAGCGTGGCGCTGGCACGGTTTCGTCGGAGCCTGGCGTCTTGAAGGCATGCGCAAGCGTCCGCCGAAAATTGGAAAACTGGCACAAAGACAAAGGGCGCCTTTCGGCGCCCTa                                                                                                                                          >  4:138093/1‑145 (MQ=255)
          cTGGCGCTGACGTTGATCGGCCTGGTGTTCTTCAGCGTGGCGCTGGCACGGTTTCGTCGGAGCCTGGCGTCTTGAAGGCATGCGCAAGCGTCCGCCGAAAATTGGAAAACTGGCACAAAGACAAAGGGCGCCTTTCGGCGCCCTa                                                                                                                                          <  3:138093/145‑1 (MQ=255)
             gcgcTGACGTTGATCGGCCTGGTGTTCTTCAGCGTGGCGCTGGCACGGTTTCGTCGGAGCCTGGCGTCTTGAAGGCATGCGCGAGCGTCCGCCGAAAATTGGAAAACTGGCACAAAGACAAAGGGCGCCTTTCGGCGCCCTACTcaca                                                                                                                                    <  5:229750/148‑1 (MQ=255)
               gcTGACGTTGATCGGCCTGGTGTACTTCAGCGTGGCGCTGGCACGGTTTCGTCGGAGCCTGGCGTCTTGAAGGCATGCGCGAGCGTCCGCCGAAAATTGGAAAACTGGCACAAAGACAAAGGGCGCCTTTCGGCGCCCTACTCacaga                                                                                                                                  >  7:254443/1‑148 (MQ=255)
                     gTTGATCGGCCTGGTGTTCTTCAGCGTGGCGCTGGCACGGTTTCGTCGGAGCCTGGCGTCTTGAAGGCATGCGCGAGCGTCCGCCGAAAATTGGAAAACTGGCACAAAGACAAAGGGCGCCGAAAGGCGCCCttt                                                                                                                                         <  7:241075/135‑3 (MQ=255)
                     gTTGATCGGCCTGGTGTTCTTCAGCGTGGCGCTGGCACGGTTTCGTCGGAGCCTGGCGTCTTGAAGGCATGCGCGAGCGTCCGCCGAAAATTGGAAAACTGGCACAAAGACAAAGGGCGCCGAAAGGCGCCCttt                                                                                                                                         >  8:241075/1‑133 (MQ=255)
                     gTTGATCGGCCTGGTGTTCTTCAGCGTGGCGCTGGCACGGTTTCGTCGGAGCCTGGCGTCTTGAAGGCATGCGCGAGCGTCCGCCGAAAATTGGAAAACTGGCACAAAGACAAAGGGCGCCGAAAGGCGCCCttt                                                                                                                                         >  4:51211/1‑133 (MQ=255)
                     gTTGATCGGCCTGGTGTTCTTCAGCGTGGCGCTGGCACGGTTTCGTCGGAGCCTGGCGTCTTGAAGGCATGCGCGAGCGTCCGCCGAAAATTGGAAAACTGGCACAAAGACAAAGGGCGCCGAAAGGCGCCCttt                                                                                                                                         <  3:51211/135‑3 (MQ=255)
                        gATCGGCCTGGTGTTCTTCAGCGTGGCGCTGGCACGGTTTCGTCCGAGCCTGGCGTCTAGAAGGCATGCGCGAGCGTCCGCCGAAAAATGGAAAACTGCCACAAAGACAAAGGGCGCCTTTCGCCGCCCAACTCACAGACAGGTCGAca                                                                                                                        >  8:182506/1‑148 (MQ=255)
                         aTCGGCCTGGTGTTCTTCAGCGTGGCGCTGGCACGGTTTCGTCGGAGCCTGGCGTCTTGAAGGCATGCGCGAGCGTCCGCCGAAAATTGGAAAACTGGCACAATGACAAAGGGCGGCTTTCGGGGACCTAGTCACAGACAGGGCGACtg                                                                                                                       >  4:266331/1‑148 (MQ=255)
                         aTCGGCCTGGTGTTCTTCAGCGTGGCGCTGGCACGGTTTCGTCGGAGCCTGGCGTCTTGAAGGCATGCGCGAGCGTCCGCCGAAAATTGGAAAACTGGCACAAAGACAAAGGGCGCCTTTCGGCGCCCTACTcaca                                                                                                                                    <  8:61926/136‑1 (MQ=255)
                         aTCGGCCTGGTGTTCTTCAGCGTGGCGCTGGCACGGTTTCGTCGGAGCCTGGCGTCTTGAAGGCATGCGCGAGCGTCCGCCGAAAATTGGAAAACTGGCACAAAGACAAAGGGCGCCTTTCGGCGCCCTACTCACAGACAGGTCGACtt                                                                                                                       <  3:290417/149‑1 (MQ=255)
                         aTCGGCCTGGTGTTCTTCAGCGGGGCGCTGGCACGGTTTCGTCGGAGCCTGGCGTCTTGAAGGCATGCGCGAGCGTCCGCCGAAAATTGGAAAACTGGCACAAAGACAAAGGGCGCCTTTCGGCGCCCTACTcaca                                                                                                                                    >  7:61926/1‑136 (MQ=255)
                         aTCGGCCTGTGTTCTTCAGCGTGGCGCTGGCACGGTTTCGTCGGAGCCTGGCGTCTTGAAGGCATGCGCGAGCGTCCGCCGAAAATTGGAAAACTGGCACAAAGACAAAGGGCGCCTTTCGGCGCCCTACTCACAGACAGGTCGACTTg                                                                                                                      >  8:138742/1‑149 (MQ=255)
                          tCGGCCTGGTGTTCTTCAGCGTGGCGCTGGCACGGTTTCGTCGGAGCCTGGCGTCTTGAAGGCATGCGCGAGCGTCCGCCGAAAATTGGAAAACTGGCACAAAGACAAAGGGCGCCTTTCGGCGCCCTACTCACAGACAGGTCGACTTg                                                                                                                      <  5:394411/149‑1 (MQ=255)
                          tCGGCCTGGTGTTCTTCAGCGTGGCGCTGGCACGGTTTCGTCGGAGCCTGGCGTCTTGAAGGCATGCGCGAGCGTCCGCCGAAAATTGGAAAACTGGCACAAAGACAAAGGGCGCCTTTCGGCGCCCTACTCACAGACAGGTCGACTTg                                                                                                                      >  3:359435/1‑149 (MQ=255)
                            ggCCTGGTGTTCTTCAGCGTGGCGCTGGCACGGTTTCGTCGGAGCCTGGCGTCTTGAAGGCATGCGCGAGCGTCCGCCGAAAATTGGAAAACTGGCACAAAGACAAAGGGCGCCTTTCGGCGCCCTACTCACAGACAGGTCGACTTGt                                                                                                                     <  2:100668/148‑1 (MQ=255)
                                tGGTGTACTTCAGCGTGCCGCTGGCACGGTTTCGTCGGAGCCTGGCGTCTAGAAGGAATGCGCGAGCTTCAGCCGAAAATTGGAAAACTGGCACAAAGACAAACGGCGCATTTCGGCGCCCTACTCACAGACAGGTCGACTTGTCATCg                                                                                                                >  4:137238/1‑149 (MQ=255)
                                  gtgtTCTTCAGCGTGGCGCTGGCACGGTTTCGTCGGAGCCTGGCGTCTTGAAGGCATGCGCGAGCGTCCGCCGAAAATTGGAAAACTGGCACAAAGACAAAGGGCGCCTTTCGGCGCCCTACTCACAGACAGGTCGACTTGTCATCGc                                                                                                               <  4:359435/148‑1 (MQ=255)
                                      tcttcAGCGTGGCGCTGGCACGGTTTCGTCGGAGCCTGGCGTCTTGAAGGCATGCGCGAGCGTCCGCCGAAAATTGGAAAACTGGCACAAAGACAAAGGGCGCCTTTCGGCGCCCTACTCACAGACAGGTCGACTTGTCATCGCCata                                                                                                           >  1:9060/1‑148 (MQ=255)
                                         tcAGCGTGGCGCTGGCACGGTTTCGTCGGAGCCTGGCGTCTTGAAGGCATGCGCGAGCGTCCGCCGAAAATTGGAAAACTGGCACAAAGACAAAGGGCGCCTTTCGGCGCCCTACTCACAGACAGGTCGACTTGTCATCGCCATATAcc                                                                                                       >  2:241113/1‑149 (MQ=255)
                                          cAGCGTGGCGCTGGCACGGTTTCGTCGGAGCCTGGCGTCTTGAAGGCATGCGCGAGCGTCCGCCGAAAATTGGAAAACTGGCACAAAGACAAAGGGCGCCTTTCGGCGCCCTACTCACAGACAGGTCGACTTGTCATCGCCATATACCt                                                                                                      >  7:29027/1‑149 (MQ=255)
                                           aGCGTGGCGCTGGCACGGTTTCGTCGGAGCCTGGCGTCTTGAAGGCATGCGCGAGCGTCCGCCGAAAATTGGAAAACTGGCACAAAGACAAAGGGCGCCTTTCGGCGCCCTACTCACAGACAGGTCGACTTGTCATCGCCATATACCTg                                                                                                     <  6:366041/149‑1 (MQ=255)
                                             cGTGGCGCTGGCACGGTTTCGTCGGAGCCTGGCGTCTTGAAGGCATGCGCGAGCGTCCGCCGAAAATTGGAAAACTGGCACAAAGACAAAGGGCGCCTTTCGGCGCCCTACTCACAGACAGGTCGACTTGTCATCGCCATATACCTGcc                                                                                                   >  5:186135/1‑149 (MQ=255)
                                              gTGGCGCTGGCACGGTTTCGTCGGAGCCTGGCGTCTTGAAGGCATGCGCGAGCGTCCGCCGAAAATTGGAAAACTGGCACAAAGACAAAGGGCGCCTTTCGGCGCCCTACTCACAGACAGGTCGACTTGTCATCGCCATATACCTGccc                                                                                                  <  8:29027/149‑1 (MQ=255)
                                                 gcgcTGGCACGGTTTCGTCGGAGCCTGGCGTCTTGAAGGCATGCGCGAGCGTCCGCCGAAAATTGGAAAACTGGCACAAAGACAAAGGGCGCCTTTCGGCGCCCTACTCACAGACAGGTCGACTTGTCATCGCCATATACCTGCCCGGc                                                                                               >  1:366169/1‑149 (MQ=255)
                                                       gCACGGTTTCGTCGGAGCCTGGCGTCTTGAAGGCATGCGCGAGCGTCCGCCGAAAATTGGAAAACTGGCACAAAGACAAAGGGCGCCTTTCGGCGCCCTACTCACAGACAGGTCGACTTGTCATCGCCATATACCTGCCCGGCTCTg                                                                                           <  3:140314/147‑1 (MQ=255)
                                                         aCGGTTTCGTCGGAGCCTGGCGTCTTGAAGGCATGCGCGAGCGTCCGCCGAAAATTGGAAAACTGGCACAAAGACAAAGGGCGCCTTTCGGCGCCCTACTCACAGACAGGTCGACTTGTCATCGCCATATACCTGCCCGGCTCTGTAAt                                                                                       >  4:67912/1‑149 (MQ=255)
                                                               tcgtcgGAGCCTGGCGTCTTGAAGGCATGCGCGAGCGTCAGCAGAAAATTGGAAAACTCGCACAAAGACAAAGGGCGCCTTTCGGCGCCCTACTCACAGACAGGTCGACTTGTCATAGCCATATACCTGCACGGCTCTGTAATTGCTgg                                                                                 >  3:362280/1‑149 (MQ=255)
                                                                   cgGAGCCTGGCGTCTTGAACGCATGCGCGAGCGTCCGCCGAAAATTGGAAAACTGGCACAAAGACAAAGGGCGCCTTTCGGCGCCCTACTCACAGACAGGTCGACTTGTCATCGCCATATACCTGCCCGGCTCTGTAATTGCTGGTCg                                                                              >  7:38386/1‑148 (MQ=255)
                                                                      aGCCTGGCGTCTTGAAGGCATGCGCGAGCGTCCGCCGAAAATTGGAAAACTGGCACAAAGACAAAGGGCGCCTTTCGGCGCCCTACTCACAGACAGGTCGACTTGTCATCACCATATACCTGCCCGGCTCTGTAATTGCTGGTCGAAcc                                                                          <  6:78244/149‑1 (MQ=255)
                                                                            gCGTCTTGAAGGCATGCGCGAGCGTCCGCCGAAAATTGGAAAACTGGCACAAAGACAAAGGGCGCCTTTCGGCGCCCTACTCACAGACAGGTCGACTTGTCATCGCCATATACCTGCCCGGCTCTGTAATTGCTGGTCGAACCAGGATc                                                                    >  5:367493/1‑149 (MQ=255)
                                                                            gCGTCTTGAAGGACTGCGCGAGCGTCCGCCGAAAATTGGAAAACTGGCACAAAGACATAGGGCGCCTTTCGGCGCCCTACTCACAGACAGGTCGACTTGTCATCGCCATATACCTGCCCGGCTCTGTAATTGCTGGTCGAACCAGGa                                                                      <  6:225274/147‑1 (MQ=255)
                                                                                 ttGAAGGCATGCGCGAGCGTCCGCCGAAAATTGGAAAACTGGCACAAAGACAAAGGGCGCCTTTCGGCGCCCTACTCACAGACAGGTCGACTTGTCATCGCCATATACCTGCCCGGCTCTGTAATTGCTGGTCGAACCAGGATCCt                                                                  >  2:120177/1‑146 (MQ=255)
                                                                                 ttGAAGGCATGCGCGAGCGTCCGCCGAAAATTGGAAAACTGGCACAAAGACAAAGGGCGCCTTTCGGCGCCCTACTCACAGACAGGTCGACTTGTCATCGCCATATACCTGCCCGGCTCTGTAATTGCTGGTCGAACCAGGATCCt                                                                  >  5:138804/1‑146 (MQ=255)
                                                                                   gAAGGCATGCGCGAGCGTCCGCCGAAAATTGGAAAACTGGCACAAAGACAAAGGGCGCCTTTCGGCGCCCTACTCACAGACAGGTCGACTTGTCATCGCCATATACCTGCCCGGCTCTGTAATTGCTGGTCGAACCAGGATCCTCa                                                                <  8:254443/146‑1 (MQ=255)
                                                                                           gcgcGAGCGTCCGCCGAAAATTGGAAAACTGGCACAAAGACAAAGGGCGCCTTTCGGCGCCCTACTCACAGACAGGTCGACTTGTCATCGCCATATACCTGCCCGGCTCTGTAATTGCTGGTCGAACCAGGATCCTCAGAGTCTTAc                                                       <  2:366169/147‑1 (MQ=255)
                                                                                             gcgAGCGTCCGCCGAAAATTGGAAAACTGGCACAAAGACAAAGGGCGCCTTTCGGCGCCCTACTCACAGACAGGTCGACTTGTCATCGCCATATACCTGCCCGGCTCTGTAATTGCTGGTCGAACCAGGATCCTCAGAGTCTTACAAGc                                                   >  6:100254/1‑149 (MQ=255)
                                                                                               gAGCGTCCGCCGAAAATTGGAAAACTGGCACAAAGACAAAGGGCGCCTTTCGGCGCCCTACTCACAGACAGGTCGACTTGTCATCGCCATATACCTGCCCGGCTCTGTAATTGCTGGTCGAACCAGGATCCTCAGAGTCTTACAAGccc                                                 <  2:54556/149‑1 (MQ=255)
                                                                                                 gCGTCCGCCGAAAATTGGAAAACTGGCACAAAGACAAAGGGCGCCTTTCGGCGCCCTACTCACAGACAGGTCGACTTGTCATCGCCATATACCTGCCCGGCTCTGTAATTGCTGGTCGAACCAGGATCCTCAGAGTCTTACAAGcccc                                                <  7:138742/148‑1 (MQ=255)
                                                                                                   gTCCGCCGAAAATTGGAAAACTGGCACAAAGACAAAGGGCGCCTTTCGGCGCCCTACTCACAGACAGGTCGACTTGTCATCGCCATATACCTGCCCGGCTCTGTAATTGCTGGTCGAACCAGGATCCTCAGAGTCTTACAAGCCCCTTg                                             >  2:156345/1‑149 (MQ=255)
                                                                                                    tCCGCCGAAAATTGGAAAACTGGCACAAAGACAAAGGGCGCCTTTCGGCGCCCTACTCACAGACAGGTCGACTTGTCATCGCCATATACCTGCCCGGCTCTGTAATTGCTGGTCGAACCAGGATCCTCAGAGTCTTACAAGCCCCTTGc                                            <  3:67912/149‑1 (MQ=255)
                                                                                                         cgAAAATTGGAAAACTGGCACAAAGACAAAGGGCGCCTTTCGGCGCCCTACTCACAGACAAGTCGAATTGTCATCGCCATATACCTGCCCGGCTCTGTAATTGCTGGTCGAACCAGGATCCTCAGAGTCTTACAAGCCCCTTGCGGGaa                                       >  7:12728/1‑149 (MQ=255)
                                                                                                          gAAAATTGGAAAACTGGCACAAAGACAAAGGGCGCCTTTCGGCGCCCTACTCACAGACAGGTCGACTTGTCATCGCCATATACCTGCCCGGCTCTGTAATTGCTGGTCGAACCAGGATCCTCAGAGTCTTACAAGCCCCTTGCGGGaa                                       >  3:310506/1‑148 (MQ=255)
                                                                                                          gAAAATTGGAAAACTGGCACAAAGACAAAGGGCGCCTTTCGGCGCCCTACTCACAGACAGGTCGACTTGTCATCGCCATATACCTGCCCGGCTCTGTAATTGCTGGTCGAACCAGGATCCTCAGAGTCTTACAAGCCCCTTGCGGGAAc                                      <  7:226691/149‑1 (MQ=255)
                                                                                                           aaaaTTGGAAAACTGGCACAAAGACAAAGGGCGCCTTTCGGCGCCCTACTCACAGACAGGTCGACTTGTCATCGCCATATACCTGCCCGGCTCTGTAATTGCTGGTCGAACCAGGATCCTCAGAGTCTTACAAGCCCCTTGCGGGAACa                                     <  1:241113/149‑1 (MQ=255)
                                                                                                             aaTTGTAAAAGTGGCTCAAAGACAAAGGGCGCGTTTCTGCGCCCTACTCACATACAAGTCGACTTGTCATCGCCATATACCTGCCCGGGTCTGTAATTGCTGTTCGAAGCTGGATCCTCAGAGTCTTACAAGCCCCTTGCGGGAAGAgg                                   <  8:12728/149‑1 (MQ=255)
                                                                                                              aTTGGAAAACTGGCACAAAGACAAAGGGCGCCTTTCGGCGCCCTACTCACAGACAGGTCGACTTGTCATCGCCATATACCTGCCCGGCTCTGTAATTGCTGGTCGAACCAGGATCCTCAGAGTCTTACAAGCCCCTTGCggg                                         >  2:26902/1‑142 (MQ=255)
                                                                                                              aTTGGAAAACTGGCACAAAGACAAAGGGCGCCTTTCGGCGCCCTACTCACAGACAGGTCGACTTGTCATCGCCATATACCTGCCCGGCTCTGTAATTGCTGGTCGAACCAGGATCCTCAGAGTCTTACAAGCCCCTTGCggg                                         <  1:26902/142‑1 (MQ=255)
                                                                                                              aTTGGAAAACTGGCACAAAGACAAAGGGCGCCTTTCGGCGCCCTACTCACAGACAGGTCGACTTGTCATCGCCATATACCTGCCCGGCTCTGTAATTGCTGGTCGAACCAGGATCCTCAGAGTCTTACAAGCCCCTTGCGGGAACAGGc                                  <  7:182506/149‑1 (MQ=255)
                                                                                                               ttGGAAAACTGGCACAAAGACAAAGGGCGCCTTTCGGCGCCCTACTCACAGACAGGTCGACTTGTCATCGCCATATACCTGCCCGGCTCTGTAATTGCTGGTCGAACCAGGATCCTCagag                                                             <  3:299723/121‑1 (MQ=255)
                                                                                                                  gAAAACTGGCACAAAGACAAAGGGCGCCTTTCGGCGCCCTACTCACAGACAGGTCGACTTGTCATCGCCATATACCTGCCCGGCTCTGTAATTGCTGGTCGAACCAGGATCCTCAGAGTCTTACAAGCCCCTTGCGGGAACAGGCTGAt                              >  2:11667/1‑149 (MQ=255)
                                                                                                                      aCTGGCACAAAGACAAAGGGCGCCTTTCGGCGCCCTACTCACAGACAGGTCGACTTGTCATCGCCATATACCTGCCCGGCTTTGTAATTGCTGGTCGAACCAGGATCCTCAGAGTCTTACAAGCCCCTTGCGGGAACAGGCTGATTATc                          <  1:120177/149‑1 (MQ=255)
                                                                                                                      aCTGGCACAAAGACAAAGGGCGCCTTTCGGCGCCCTACTCACAGACAGGTCGACTTGTCATCGCCATATACCTGCCCGGCTCTGTAATTGCTGGTCGAACCAGGATCCTCAGAGTCTTACAAGCCCCTTGCGGGAACAGGCTGATTATc                          <  8:38386/149‑1 (MQ=255)
                                                                                                                      aCTGGCACAAAGACAAAGGGCGCCTTTCGGCGCCCTACTCACAGACAGGTCGACTTGTCATCGCCATATACCTGCCCGGCTCTGTAATTGCTGGTCGAACCAGGATCCTCAGAGTCTTACAAGCCCCTTGCGGGAACAGGCTGATTATc                          <  6:138804/149‑1 (MQ=255)
                                                                                                                        tGGCACAAAGACAAAGGGCGCCTTTCGGCGCCCTACTCACAGACAGGTCGACTTGTCATCGCCATATACCTGCCCGGCTCTGTAATTGCTGGTCGAACCAGGATCCTCAGAGTCTTACAAGCCCCTTGCGGGAACAGGCTGATTATCa                         >  8:126233/1‑148 (MQ=255)
                                                                                                                          gCACAAAGACAAAGGGCGCCTTTCGGCGCCCTACTCACAGACAGGTCGACTTGTCATCGCCATATACCTGCCCGGCTCTGTAATTGCTGGTCGAACCAGGATCCTCAGAGTCTTACAAGCCCCTTGCGGGAACAGGCTGATTATCa                         <  7:126233/146‑1 (MQ=255)
                                                                                                                          gCACAAAGACAAAGGGCGCCTTTCGGCGCCCTACTCACAGACAGGTCGACTTGTCATCGCCATATACCTGCCCGGCTCTGTAATTGCTGGTCGAACCAGGATCCTCAGAGTCTTACAAGCCCCTTGCGGGAACAGGCTGATTATCAc                        <  2:115466/147‑1 (MQ=255)
                                                                                                                               tagagcaagGGCGCCTTTCGGCGCCCTACTCACAGAGAGGTCGACTTGTCCTCGCCATCTACCTGCCCGGCGCTGTAATTCCTGGTCGAACCAGCATCCTCAGAGGCTTACAAGCCCCTTGCGGGAACAGGCTGATTCTCCCGCATCGa                 <  4:362280/143‑1 (MQ=255)
                                                                                                                               aagacaaagGGCGCCTTTCGGCGCCCTACTCACAGACAGGTCGACTTGTCATCGCCATATACCTGCCCGGCTCTGTAATTGCTGGTCGAACCAGGATCCTCAGAGTCTTACAAGCCCCTTGCGGGAACAGGCTGATTATCACGCATCGa                 <  4:374513/149‑1 (MQ=255)
                                                                                                                                       gggCGCCTTTCGGCGCCCTACTCACAGACAGGTCGACTTGTCATCGCCATATACCTGCCCGGCTCTGTAATTGCTGGTCGAACCAGGATCCTCAGAGTCTTACAAGCCCCTTGCGGGAACAGGCTGATTATCACGCATCGAGCGTAACa         <  1:156345/149‑1 (MQ=255)
                                                                                                                                            ccTTTCGGCGCCCTACTCACAGACAGGTCGACTTGTCATCGCCATATACCTGCCCGGCTCTGTAATTGCTGGTCGAACCAGGATCCTCAGAGTCTTACAAGCCCCTTGCGGGAACAGGCTGATTATCACGCATCGAGCGTAACATTGac    <  6:186135/149‑1 (MQ=255)
                                                                                                                                              tttCGGCGCCCTACTCACAGACAGGTCGACTTGTCATCGCCATATACCTGCCCGGCTCTGTAATTGCTGGTCGAACCAGGATCCTCAGAGTCTTACAAGCCCCTTGCGGGAACAGGCTg                                >  2:260037/1‑119 (MQ=255)
                                                                                                                                              tttCGGCGCCCTACTCACAGACAGGTCGACTTGTCATCGCCATATACCTGCCCGGCTCTGTAATTGCTGGTCGAACCAGGATCCTCAGAGTCTTACAAGCCCCTTGCGGGAACAGGCTg                                <  1:260037/119‑1 (MQ=255)
                                                                                                                                              tttCGGCGCCCTACTCACAGACAGGTCGACTTGCCATCGCCATATACCTGCCCGGCTCTGTAATTGCTGGTCGAACCAGGATCCTCAGAGTCTTACAAGCCCCTTGCGGGAACAGGCTGATTATCACGCATCGAGCGTAACATTGacaa  >  6:331474/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                  
GCCGGATATCCTGGCGCTGACGTTGATCGGCCTGGTGTTCTTCAGCGTGGCGCTGGCACGGTTTCGTCGGAGCCTGGCGTCTTGAAGGCATGCGCGAGCGTCCGCCGAAAATTGGAAAACTGGCACAAAGACAAAGGGCGCCTTTCGGCGCCCTACTCACAGACAGGTCGACTTGTCATCGCCATATACCTGCCCGGCTCTGTAATTGCTGGTCGAACCAGGATCCTCAGAGTCTTACAAGCCCCTTGCGGGAACAGGCTGATTATCACGCATCGAGCGTAACATTGACAA  >  NC_002947/5941190‑5941480

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