Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,975,167 A→T 9.2% intergenic (+88/+19) PP_3504 → / ← PP_3505 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,975,1670AT 9.2% 212.0 / 9.6 76intergenic (+88/+19)PP_3504/PP_3505hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (36/33);  new base T (2/5);  total (38/38)
Fisher's exact test for biased strand distribution p-value = 4.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.35e-01

AGTGCGCGCCGAATTGTGTGAAAAGGGCTACGACCTGCAACCGTGGTCGATGCCGAAGTAATATCGGCCAGCACAACAAGTGCGTGCATGGTCACGCACTTGTACAACTTCTCTGTTGCTGTTGTGGCAACTTTCCCAATTGCGACAAAGTGTCGCAACCTGCCAGTTCAATCCATCACCCTGAATTGTTGCAAGTGCTGGTAATCAGCCTGCAATTGTTCTGCCAGACGCTGCGCGCGTCCCAGGCGCACGGGCGCCAGTTCCATATCGACCAACAACGTGGCGCACGGCAT  >  NC_002947/3975019‑3975311
                                                                                                                                                    |                                                                                                                                                 
agTGCGCGCCGAATTGTGTGAAAAGGGCTACGACCTGCAACCGTGGTCGATGCCGAAGTAATATCGGCCAGCACAACAAGTGCGTGCATGGTCACGCACTTGTACAACTTCTCTGTTGCTGTTGTGGCAACTTTCCCAATTGCGACaaa                                                                                                                                                   <  6:301314/149‑1 (MQ=255)
   gcgcgcCGAATTGTGTGAAAAGGGCTACGACCTGCAACCGTGGTCGATGCCGAAGTAATATCGGCCAGCACAACAAGTGCGTGCATGGTCACGCACTTGTACAACTTCTCTGTTGCTGTTGTGGCAACTTTCCCAATTGCGACAAAgtg                                                                                                                                                >  1:290303/1‑149 (MQ=255)
     gcgcCGAATTGTGTGAAAAGGGCTACGACCTGCAACCGTGGTCGATGCCGAAGTAATATCGGCCAGCACAACAAGTGCGTGCATGGTCACGCACTTGTACAACTTCTCTGTTGCTGTTGTGGCAACTTTCCCAATTGCGACAAAGTGTc                                                                                                                                              <  8:239369/149‑1 (MQ=255)
         cGAATTGTGTGAAAAGGGCTACGACCTGCAACCGTGGTCGATGCCGAAGTAATATCGGCCAGCACAACAAGTGCGTGCATGGTCACGCACTTGTACAACTTCTCTGTTGCTGTTGTGGCAACTTTCCCAATTGCGACAAAGTGTCGCaa                                                                                                                                          >  3:122234/1‑149 (MQ=255)
         cGAATTGTGTGAAAAGGGCTACGACCTGCAACCCTGGTCGATGCCGAAGTAATATCGGCCAGCACAACAAGTGCGTGCATGGTCACGCACTTGTACAACTTCTCTGTTGCTGTTGTGGCAACTTTCCCAATTGCGACACTTTGTCGCaa                                                                                                                                          <  7:107202/149‑1 (MQ=255)
         cGAATTGTCTGAAAAGGGCTACGACCTGCAACCGTGGTCGATGCCGAAGTAATATCGGCCAGCACAACAAGTGCGTGCATGGTCACGCACTTGTACAACTTCTCTGTTGCTGTTGTGGCAACTTTCCCAATTGCGACACTTTGTCGCaa                                                                                                                                          <  7:112432/149‑1 (MQ=255)
           aaTTGTGTGAAAAGGGCTACGACCTGCAACCGTGGTCGATGCCGAAGTAATATCGGCCAGCACAACAAGTGCGTGCATGGTCACGCACTTGTACAACTTCTCTGATGCTGTTGTGGCAACTTTCCCAATTGCGACACTTTGTCGCaatt                                                                                                                                        <  8:69256/149‑3 (MQ=255)
            aTTGTGTGAAAAGGGCTACGACCTGCAACCGTGGTCGATGCCGAAGTAATATCGGCCAGCACAACAAGTGCGTGCATGGTCACGCACTTGTACAACTTCTCTGTTGCTGTTGTGGCAACTTTCCCAATTGCGACACTTTGTCGCaattg                                                                                                                                       <  2:275803/149‑4 (MQ=255)
               gtgtgtATAGGGCTAGGACCTGCAACCGTGGTCGATGCCGAAGTAATTTCGGCCAGCAGAACAAGTGCGTGCATGGTCACGCACTTGTACAACTTCTCTGTTGCTGTTGTGGCAACTTTCCCAATTGCGACAAAGTGTCGCAACCTGc                                                                                                                                     <  8:322008/148‑1 (MQ=255)
               gtgtgAAAAGGGCTACGACCTGCAACCGTGGTCGATGCCGAAGTAATATCGGCCAGCACAACAAGTGCGTGCATGGTCACGCACTTGTACAACTTCTCTGTTGCTGTTGTGGCAACTTTCCCAATTGCGACAAAGTGTCGCAACCTGc                                                                                                                                     <  4:87624/148‑1 (MQ=255)
               gtgtgAAAAGGGCTACGACCTGCAACCGTGGTCGATGCCGAAGTAATATCGGCCAGCACAACAAGTGCGTGCATGGTCACGCACTTGTACAACTTCTCTGTTGCTGTTGTGGCAACTTTCCCAATTGCGACAAAGTGTCGCAACCTGc                                                                                                                                     <  2:290303/148‑1 (MQ=255)
                tgtgAAAAGGGCTACGACCTGCAACCGTGGTCGATGCCGAAGTAATATCGGCCAGCACAACAAGTGCGTGCATGGTCACGCACTTGTACAACTTCTCTGTTGCTGTTGTGGCAACTTTCCCAATTGCGACAAAGTGTCGCAACCTGcc                                                                                                                                    <  6:176165/148‑1 (MQ=255)
                  tgAAAAGGGCTACGACCTGCAACCGTGGTCGATGCCGAAGTAATATCGGCCAGCACAACAAGTGCGTGCATGGTCACGCACTTGTACAACTTCTCTGTTGCTGTTGTGGCAACTTTCCCAATTGCGACAAAGTGTCGCAACCTGCCAGt                                                                                                                                 >  7:300637/1‑149 (MQ=255)
                     aaaGGGCTACGACCTGCAACCGTGGTCGATGCCGAAGTAATATCGGCCAGCACAACAAGTGCGTGCATGGTCACGCACTTGTACAACTTCTCTGTTGCTGTTGTGGAAAATTTCCCAATTGCGACAAAGTGTCGCAACCTGCCAGTTc                                                                                                                               >  5:222518/1‑148 (MQ=255)
                       aGGGCTACGACCTGCAACCGTGGTCGATGCCGAAGTAATATCGGCCAGCACAACAAGTGCGTGCATGGTCACGCACTTGTACAACTTCTCTGTTGCTGTTGTGGCAACTTTCCCAATTGCGACAAAGTGTCGCAACCTGCCAGTTCAAt                                                                                                                            <  1:184650/149‑1 (MQ=255)
                       aGGGCTACGACCTGCAACCGTGGTCGATGCCGAACTAATATCGGCCAGCACAACAAGTGCGTGCATGGTCACGCACTTGTACAACTTCTCTGTTGCTGTTGTGGCAACTTTCCCAATTGCGACAAAGTGTCGCAACCTGCCAGTTCAAt                                                                                                                            >  7:103805/1‑149 (MQ=255)
                         ggCTACGACCTGCAACCGTGGTCGATGCCGAAGTAATATCGGCCAGCACAACAAGTGCGTGCATGGTCACGCACTTGTACAACTTCTCTGTTGCTGTTGTGGCAACTTTCCCAATTGCGACAAAGTGTCGCAACCTGCCAGTTCAATc                                                                                                                           >  8:104333/1‑148 (MQ=255)
                           cTACGACCTGCAACCGTGGTCGATGCCGAAGTAATATCGGCCAGCACAACAAGTGCGTGCATGGTCACGCACTTGTACAACTTCTCTGTTGCTGTTGTGGCAACTTTCCCAATTGCGACACTTTGTCGCaattg                                                                                                                                       >  5:370884/1‑131 (MQ=255)
                               gACCTGCAACCGTGGTCGATGCCGAAGTAATATCGGCCAGCACAACAAGTGCGTGCATGGTCACGCACTTGTACAACTTCTCTGTTGCTGTTGTTGCAACTTTCCCAATTGCGACAAAGTGTCGCAACCTGCCAGTTCAATCCATCAcc                                                                                                                    >  4:23901/1‑149 (MQ=255)
                                 ccTGCAACCGTGGTCGATGCCGAAGTAATATCGGCCAGCACAACAAGTGCGTGCATGGTCACGCACTTGTACAACTTCTCTGTTGCTGTTGTGGCAACTTTCCCAATTGCGACAAAGTGTCGCAACCTGCCAGTTCAATCCATCACCCt                                                                                                                  <  8:73053/149‑1 (MQ=255)
                                 ccTGCAACCGTGGTCGATGCCGAAGTAATATCGGCCAGCACAACAAGTGCGTGCATGGTCACGCACTTGTACAACTTCTCTGTTGCTGTTGTGGCAACTTTCCCAATTGCGACAAAGTGTCGCAACCTGCCAGTTCAATCCATCACCCt                                                                                                                  <  2:342647/149‑1 (MQ=255)
                                    gCAACCGTGGTCGATGCCGAAGTAATATCGGCCAGCACAACAAGTGCGTGCATGGTCACGCACTTGTACAACTTCTCTGTTGCTGTTGTGGCAACTTTCCCAATTGCGACACTTTGTCGCaattgg                                                                                                                                      <  5:131772/126‑5 (MQ=255)
                                    gCAACCGTGGTCGATGCCGAAGTAATATCGGCCAGCACAACAAGTGCGTGCATGGTCACGCACTTGTACAACTTCTCTGTTGCTGTTGTGGCAACTTTCCCAATTGCGACACTTTGTCGCaattgg                                                                                                                                      >  6:131772/1‑122 (MQ=255)
                                     cAACCGTGGTCGATGCCGAAGTAATATCGGCCAGCACAACAAGTGCGTGCATGGTCACGCACTTGTACAACTTCTCTGTTGCTGTTGTGGCAACTTTCCCAATTGCGACAAAGTGTCGCAACCTGCCAGTTCAATCCATCACCCTGAAt                                                                                                             <  5:169246/149‑1 (MQ=255)
                                         cGTGGTCGATGCCGAAGTAATATCGGCCAGCACAACAAGTGCGTGCATGGTCACGCACTTGTACAACTTCTCTGTTGCTGTTGTGGCAACTTTCCCAATTGCGACAAAGTGTCGCAACCTGCCAGTTCAATCCATCACGCTGAAttgtt                                                                                                         <  6:362500/149‑1 (MQ=255)
                                             gTCGATGCCGAAGTAATATCGGCCAGCACAACAAGTGCGTGCATGGTCACGCACTTGTACAACTTCTCTGTTGCTGTTGTGGCAACTTTCCCAATTGCGACAAAGTGTCGCAACCTGCCAGTTCAATCCATCACCCTGAATTGTTGCaa                                                                                                     <  5:263446/149‑1 (MQ=255)
                                                gATGCCGAAGTAATATCGGCCAGCACAACAAGTGCGTGCATGGTCACGCACTTGTACAACTTCTCTGTTGCTGTTGTGGCAACTTTCCCAATTGCGACAAAGTGTCGCAACCTGCCAGTTCAATCCATCACCCTGAATTGTTGCAAGTg                                                                                                  >  7:200076/1‑149 (MQ=255)
                                                    ccGAAGTAATATCGGCCAGCACAACAAGTGCGTGCATGGTCACGCACTTGTACAACTTCTCTGTTGCTGTTGTGGCAACTTTCCCAATTGCGACAAAGTGTCGCAACCTGCCAGTTCAATCCATCACCCTGAATTGTTGCAAGTGCTgg                                                                                              <  7:104333/149‑1 (MQ=255)
                                                    ccGAAGTAATATCGGCCAGCACAACAAGTGCGTGCATGGTCACGCACTTGTACAACTTCTCTGTTGCTGTTGTGGCAACTTTCCCAATTGCGACAAAGTGTCGCAACCTGCCAGTTCAATCCATCACCCTGAATTGTTGCAAGTGCTg                                                                                               <  8:103805/148‑1 (MQ=255)
                                                      gAAGTAATATCGGCCAGCACAACAAGTGCGTGCATGGTCACGCACTTGTACAACTTCTCTGTTGCTGTTGTGGCAACTTTCCCAATTGCGACAAAGTGTCGCAACCTGCCAGTTc                                                                                                                               <  6:177306/115‑1 (MQ=255)
                                                      gAAGTAATATCGGCCAGCACAACAAGTGCGTGCATGGTCACGCACTTGTACAACTTCTCTGTTGCTGTTGTGGCAACTTTCCCAATTGCGACAAAGTGTCGCAACCTGCCAGTTc                                                                                                                               >  5:177306/1‑115 (MQ=255)
                                                          tAATATCGGCCAGCACAACAAGTGCGTGCATGGTCACGCACTTGTACAACTTCTCTGTTGCTGTTGTGGCAACTTTCCCAATTGCGACAAAGTGTCGCAACCTGCCAGTTCAATCCATCACCCTGAATTGTTGCAAGTGCTGGTAAt                                                                                          >  5:287327/1‑147 (MQ=255)
                                                            atatCGGCCAGCACAACAAGTGCGTGCATGGTCACGCACTTGTACAACTTCTCTGTTGCTGTTGTGGCAACTTTCCCAATTGCGACAAAGTGTCGCAACCTGCCAGTTCAATCCATCACCCTGAATTGTTGCAAGTGCTGGTAATCAGc                                                                                      <  2:156151/149‑1 (MQ=255)
                                                               tCGGCCAGCACAACAAGTGCGTGCATGGTCACGCACTTGTACAACTTCTCTGTTGCTGTTGTGGCAACTTTCCCAATTGCGACAAAGTGTCGCAACCTGCCAGTTCAATCCATCACCCTGAATTGTTGCAAGTGCTGGTAATCAGCCt                                                                                    <  6:287327/148‑1 (MQ=255)
                                                                  gCCAGCACAACAAGTGCGTGCATGGTCACGCACTTGTACAACTTCTCTGTTGCTGTTGTGGCAACTTTCCCAATTGCGACAAAGTGTCGCAACCTGCCAGTTCAATCCATGACCCTGAATTGTTGCAAGTGCTGGTAATCAGCCTGCaa                                                                                >  6:279134/1‑149 (MQ=255)
                                                                      gCACAACAAGTGCGTGCATGGTCACGCACTTGTACAACTTCTCTGTTGCTGTTGTGGCAACTTTCCCAATTGCGACAAAGTGTCGCAACCTGCCAGTTCAATCCATCACCCTGAATTGTTGCAAGTGCTGGTAATCAGCCTGCAATTGt                                                                            >  6:274055/1‑149 (MQ=255)
                                                                           acaaGTGCGTGCATGGTCACGCACTTGTACAACTTCTCTGTTGCTGTTGTTGCAACTTTCCCAATTGCGACAAAGTGTCGCAACCTGCCAGTTCAATCCATCACCCTGAATTGTTGCAAGTGCTGGTAATCAGCCTGCAATTGTTCTGc                                                                       <  3:23901/149‑1 (MQ=255)
                                                                             aaGTGCGTGCATGGTCACGCACTTGTACAACTTCTCTGTTGCTGTTGTGGCAACTTTCCCAATTGCGACAAAGTGTCGCAACCTGCCAGTTCAATCCATCACCCTGAATTGTTGCAAGTGCTGGTAATCAGCCTGCAATTGTTCTGc                                                                       >  4:240618/1‑147 (MQ=255)
                                                                               gtgcgtgcATGGTCACGCACTTGTACAACTTCTCTGTTGCTGTTCTGGCAACTTTCCCAATTGCGACAAAGTGTCGCAACCTGCCAGTTCAATCCATCACCCTGAATTGTTGCAAGTGCTGGTAATCAGCCTGCAATTGTTCTGCCAGa                                                                   >  2:45348/1‑149 (MQ=255)
                                                                                   gtgcATGGTCACGCACTTGTACAACTTCTCTGTTGCTGTTGTGGCAACTTTCCCAATTGCGACAAAGTGTCGCAACCTGCCAGTTCAATCCATCACCCTGAATTGTTGCAAGTGCTGGTAATCAGCCTGCAATTGTTCTGCCAGACGCt                                                               >  8:307321/1‑149 (MQ=255)
                                                                                   gtgcATGGTCACGCACTTGTACAACTTCTCTGTTGCTGTTGTGGCAACTTTCCCAATTGCGACAAAGTGTCGCAACCTGCCAGTTCAATCCATCACCCTGAATTGTTGAAAGTGCTGGTAATCAGCCTGCAATTGTTCTGCCAGACGCt                                                               >  7:349911/1‑149 (MQ=255)
                                                                                     gcATGGTCACGCACTTGTAGAACTTCTCTGTTGCTGTTGTGGCAACTTTCCCAATTGCGACAAAGTGTCGCAACCTGCCAGTTCAATCCATCACCCTGAATTGTTGCAAGTGCTGGTAATCAGCCTGCAATTGTTCTTCCAGACGCTgc                                                             >  4:307465/1‑149 (MQ=255)
                                                                                           tCACGCACTTGTACAACTTCTCTGTTGCTGTTGTGGCAACTTTCCCAATTGCGACAAAGTGTCGCAACCTGCCAGTTCAATCCATCACCCTGAATTGTTGCAAGTGCTGGTAATCAGCCTGCAATTGTTCTGCCAGACGCTGCGCGCGt                                                       <  1:11058/149‑1 (MQ=255)
                                                                                           tCACGCACTTGTACAACTTCTCTGTTGCTGTTGTGGCAACTTTCCCAATTGCGACAAAGTGTCGCAACCTGCCAGTTCAATCCATCACCCTGAATTGTTGCAAGTGCTGGTAATCAGCCTGCAATTGTTCTGCCAGACGCTGCGCGCGt                                                       <  1:193292/149‑1 (MQ=255)
                                                                                             aCGCACTTGTACAACTTCTCTGTTGCTGTTGTGGCAACTTTCCCAATTGCGACAAAGTGTCGCAACCTGCCAGTTCAATCCATCACCCTGAATTGTTGCAAGTGCTGGTAATCAGCCTGCAATTGTTCTGCCAGACGCTGCGCGCGTcc                                                     >  2:192596/1‑149 (MQ=255)
                                                                                              cGCACTTGTACAACTTCTCTGTTGCTGTTGTGGCAACTTTCCCAATTGCGACAAAGTGTCGCAACCTGCCAGTTCAATCCATCACCCTGAATTGTTGCAAGTGCTGGTAATCAGCCTGCAATTGTTCTGCCAGACGCTGCGCGCGTccc                                                    <  3:239220/149‑1 (MQ=255)
                                                                                              cGCACTTGTACAACTTCTCTGTTGCTGTTGTGGCAACTTTCCCAATTGCGACAAAGTGTCGCAACCTGCCAGTTCAATCCATCACCCTGAATTGTTGCAAGTGCTGGTAATCAGCCTGCAATTGTTCTGCCAGACGCTGCGCGCGTccc                                                    <  4:122234/149‑1 (MQ=255)
                                                                                               gCACTTGTACAACTTCTCTGTTGCTGTTGTGGCAACTTTCCCAATTGCGACAAAGTGTCGCAACCTGCCAGTTCAATCCATCACCCTGAATTGTTGCAAGTGCTGGTAATCAGCCTGCAATTGTTCTGCCAGACGCTGCGCGCGTcc                                                     >  5:134109/1‑147 (MQ=255)
                                                                                                 aCTTGTACAACTTCTCTGTTGCTGTTGTGGCAACTTTCCCAATTGCGACAAAGTGTCGCAACCTGCCAGTTCAATCCATCACCCTGAATTGTTGCAAGTGCTGGTAATCAGCCTGCAATTGTTCTGCCAGACGCTGCGCGCGTCCCAgg                                                 <  8:200076/149‑1 (MQ=255)
                                                                                                  cTTGTACAACTTCTCTGTTGCTGTTGTGGCAACTTTCCCAATTGCGACAAAGTGTCGCAACCTGCCAGTTCAATCCATCACCCTGAATTGTTGCAAGTGCTGGTAATCAGCCTGCAATTGTTCTGCCAGACGCTGCGCGCGTCCCAGgc                                                <  8:300637/149‑1 (MQ=255)
                                                                                                     gTACAACTTCTCTGTTGCTGTTGTGGCAACTTTCCCAATTGCGACAAAGTGTCGCAACCTGCCAGTTCAATCCATCACCCTGAATTGTTGCAAGTGCTGGTAATCAGCCTGCAATTGTTCTGCCAGACGCTGCGCGCGTCCCAGGCGCa                                             >  2:5072/1‑149 (MQ=255)
                                                                                                      tACAACTTCTCTGTTGCTGTTGTGGCAACTTTCCCAATTGCGACAAAGTGTCGCAACCTGCCAGTTCAATCCATCACCCTGAATTGTTGCAAGTGCTGGTAATCAGCCTGCAATTGTTCTGCCAGAc                                                                  >  8:284136/1‑127 (MQ=255)
                                                                                                      tACAACTTCTCTGTTGCTGTTGTGGCAACTTTCCCAATTGCGACAAAGTGTCGCAACCTGCCAGTTCAATCCATCACCCTGAATTGTTGCAAGTGCTGGTAATCAGCCTGCAATTGTTCTGCCAGAc                                                                  <  7:284136/127‑1 (MQ=255)
                                                                                                         aaCTTCTCTGTTGCTGTTGTGGCAACTTTCCCAATTGCGACAAAGTGTCGCAACCTGCCAGTTCAATCCATCACCCTGAATTGTTGCAAGTGCTGGTAATCAGCCTGCAATTGTTCt                                                                         <  5:15022/117‑1 (MQ=255)
                                                                                                         aaCTTCTCTGTTGCTGTTGTGGCAACTTTCCCAATTGCGACAAAGTGTCGCAACCTGCCAGTTCAATCCATCACCCTGAATTGTTGCAAGTGCTGGTAATCAGCCTGCAATTGTTCt                                                                         >  6:15022/1‑117 (MQ=255)
                                                                                                         aaCTTCTCTGTTGCTGTTGTGGCAACTTTCCCAATTGCGACAAAGTGTCGCAACCTGCCAGTTCAATCCATCACCCTGAATTGTTGCAAGTGCTGGTAATCAGCCTGCAATTGTTCTGCCAGACGCTGCGCGCGTCCCAGGCGCACggg                                         >  8:70327/1‑149 (MQ=255)
                                                                                                              ctctGTTGCTGTTGTGGCAACTTTCCCAATTGCGACAAAGTGTCGCAACCTGCCAGTTCAATCCATCACCCTGAATTGTTGCAAGTGCTGGTAATCAGCCTGCAATTGTTCTGCCAGACGCTGCGCGCGTCCCAGGCGCACGGgcgc                                      <  7:307321/147‑1 (MQ=255)
                                                                                                                 tgttgctgttgTGGCAACTTTCCCAATTGCGACAAAGTGTCGCAACCTGCCAGTTCAATCCATCACCCTGAATTGTTGCAAGTGCTGGTAATCAGCCTGCAATTGTTCTGCCAGACGCTGCGCGCGTCCCAGGCGCACGGGCGCCa                                    <  6:134109/146‑1 (MQ=255)
                                                                                                                 tgttgctgttgTGGCAACTTTCCCAATTGCGACAAAGTGTCGCAACCTGCCAGTTCAATCCATCACCCTGAATTGTTGCAAGTGCTGGTAATCAGCCTGCAATTGTTCTGCCAGACGCTGCGCGCGTCCCAGGCGCACGGGCGCCAGtt                                 >  4:124392/1‑149 (MQ=255)
                                                                                                                   ttgctgttgTGGCAACTTTCCCAATTGCGACAAAGTGTCGCAACCTGCCAGTTCAATCCATCACCCTGAATTGTTGCAAGTGCTGGTAATCAGCCTGCAATTGTTCTGCCAGACGCTGCGCGCGTCCCAGGCGCACGGGCGCCAGtt                                 <  3:240618/147‑1 (MQ=255)
                                                                                                                    tgctgttgTGGCAACTTTCCCAATTGCGACAAAGTGTCGCAACCTGCCAGTTCAATCCATCACCCTGAATTGTTGCAAGTGCTGGTAATCAGCCTGCAATTGTTCTGCCAGACGCTGCGCGCGTCCCAGGCGCACGGGCGCCAGTTcc                               <  1:45348/148‑1 (MQ=255)
                                                                                                                     gctgttgTGGCAACTTTCCCAATTGCGACAAAGTGTCGCAACCTGCCAGTTCAATCCATCACCCTGAATTGTTGCAAGTGCTGGTAATCAGCCTGCAATTGTTCTGCCAGACGCTGCGCGCGTCCCAGGCGCACGGGCGCCAGTTCCa                              >  2:240053/1‑148 (MQ=255)
                                                                                                                         ttgtGGCAACTTTCCCAATTGCGACAAAGTGTCGCAACCTGCCAGTTCAATCCATCACCCTGAATTGTTGCAAGTGCTGGTAATCAGCCTGCAATTGTTCTGCCAGACGCTGCGCGCGTCCCAGGCGCACGGGCGCCAGTTCCATATc                          <  7:186517/148‑1 (MQ=255)
                                                                                                                         ttgtGGCAACTTTCCCAATTGCGACAAAGTGTCGCAACCTGCCAGTTCAATCCATCACCCTGAATTGTTGCAAGTGCTGGTAATCAGCCTGCAATTGTTCTGCCAGACGCTGCGCGCGTCCCAGGCGCACGGGCGCCAGTTCCATATc                          <  1:240053/148‑1 (MQ=255)
                                                                                                                           gtgGCAACTTTCCCAATTGCGACAAAGTGTCGCAACCTGCCAGTTCAATCCATCACCCTGAATTGTTGCAAGTGCTGGTAATCAGCCTGCAATTGTTCTGCCAGACGCTGCGCGCGTCCCAGGCGCACGGGCGCCAGTTCCATATCGAc                       >  2:365745/1‑149 (MQ=255)
                                                                                                                            tgGCAACTTTCCCAATTGCGACAAAGTGTCGCAACCTGCCAGTTCAATCCATCACCCTGAATTGTTGCAAGTGCTGGTAATCAGCCTGCAATTGTTCTGCCAGACGCTGCGCGCGTCCCAGGCGCACGGGCGCCAGTTCCATATCGAc                       >  2:59956/1‑148 (MQ=255)
                                                                                                                            tgGCAACTTTCCCAATTGCGACAAAGTGTCGCAACCTGCCAGTTCAATCCATCACCCTGAATTGTTGCAAGTGCTGGTAATCAGCCTGCAATTGTTCTGCCAGACGCTGCGCGCGTCCCAGGCGCACGGGCGCCAGTTCCATATCGAc                       >  2:186200/1‑148 (MQ=255)
                                                                                                                               cAACTTTCCCAATTGCGACAAAGTGTCGCAACCTGCCAGTTCAATCCATCACCCTGAAATTGTTGCAAGTGCTGGTAATCAGCCTGCAATTGTTCTGCCAGACGCTGCGCGCGTCCCAGGCGCACGGGCGCCAGTTCCATATCGACcaa                    >  1:288896/1‑149 (MQ=255)
                                                                                                                               cAACTTTCCCAATTGCGACAAAGTGTCGCAACCTGCCAGTTCAATCCATCACCCTGAATTGTTGCAAGTGCTGGTAATCAGCCTGCAATTGTTCTGCCAGACGCTGCGCGCGTCCCAGGCGCACGGGCGCCAGTTCCATATCGACcaa                    >  4:53330/1‑148 (MQ=255)
                                                                                                                               cAACTTTCCCAATTGCGACAAAGTGTCGCAACCTGCCAGTTCAATCCATCACCCTGAATTGTTGCAAGTGCTGGTAATCAGCCTGCAATTGTTCTGCCAGACGCTGCGCGCGTCCCAGGCGCACGGGCGCCAGTTCCATATCGACcaa                    >  6:373278/1‑148 (MQ=255)
                                                                                                                                         aaTTGCGACAAAGTGTCGCAACCTGCCAGTTCAATCCATCACCCTGAATTGTTGCAAGTGCTGGTAATCAGCCTGCAATTGTTCTGCCAGACGCTGCGCGCGTCCCAGGCGCACGGGCGCCAGTTCCATATCGACCAACAACGTGgcg          >  6:5851/1‑148 (MQ=255)
                                                                                                                                          aTTGCGACAAAGTGTCGCAACCTGCCAGTTCAATCCATCACCCTGAATTGTTGCAAGTGCTGGTAATCAGCCt                                                                                    >  1:262895/1‑73 (MQ=255)
                                                                                                                                          aTTGCGACAAAGTGTCGCAACCTGCCAGTTCAATCCATCACCCTGAATTGTTGCAAGTGCTGGTAATCAGCCt                                                                                    <  2:262895/73‑1 (MQ=255)
                                                                                                                                               gACAAAGTGTCGCAACCTGCCAGTTCAATCCATCACCCTGAATTGTTGCAAGTGCTGGTAATCAGCCTGCAATTGTTCTGCCAGACGCTGCGCGCGTCCCAGGCGCACGGGCGCCAGTTCCATATCGACCAACAACGTg             <  7:48400/139‑1 (MQ=255)
                                                                                                                                               gACAAAGTGTCGCAACCTGCCAGTTCAATCCATCACCCTGAATTGTTGCAAGTGCTGGTAATCAGCCTGCAATTGTTCTGCCAGACGCTGCGCGCGTCCCAGGCGCACGGGCGCCAGTTCCATATCGACCAACAACGTg             >  8:48400/1‑139 (MQ=255)
                                                                                                                                                aCAAAGTGTCGCAACCTGCCAGTTCAATCCATCACCCTGAATTGTTGCAAGTGCTGGTAATCAGCCTGCAATTGTTCTGCCAGACGCTGCGCGCGTCCCAGGCGCACGGGCGCCAGTTCCATATCGACCAACAACGTGGCGCACGGCAt  >  6:134411/1‑149 (MQ=255)
                                                                                                                                                aCAAAGTGTCGCAACCTGCCAGTTCAATCCATCACCCTGAATTGTTGCAAGTGCTGGTAATCAGCCTGCAATTGTTCTGCCAGACGCTGCGCGCGTCCCAGGCGCACGGGCGCCAGTTCCATATCGACCAACAACGTGGCGCACGGCAt  >  7:336358/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
AGTGCGCGCCGAATTGTGTGAAAAGGGCTACGACCTGCAACCGTGGTCGATGCCGAAGTAATATCGGCCAGCACAACAAGTGCGTGCATGGTCACGCACTTGTACAACTTCTCTGTTGCTGTTGTGGCAACTTTCCCAATTGCGACAAAGTGTCGCAACCTGCCAGTTCAATCCATCACCCTGAATTGTTGCAAGTGCTGGTAATCAGCCTGCAATTGTTCTGCCAGACGCTGCGCGCGTCCCAGGCGCACGGGCGCCAGTTCCATATCGACCAACAACGTGGCGCACGGCAT  >  NC_002947/3975019‑3975311

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: