Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,037,577 A→T 6.6% intergenic (+171/+28) PP_1810 → / ← PP_5478 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,037,5770AT 6.6% 173.1 / 2.8 61intergenic (+171/+28)PP_1810/PP_5478hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (30/27);  new base T (2/2);  total (32/29)
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.94e-01

CATGGGATGGATGGCACCGGCGGTGCCGGTGTTCGCGGGCGCGCCCGCTCCCACAGTGTGGGGTGCTTGGCCAATGCAGTGGTTTGCCCTCGCCACGTCTGGCCTGATGAAGCCCATTTCTTCTGACTGCAAGCGAGGCATAGCCTCGCCTGGTCTCGCGTCATTCCCTTCAGTGGTGCCTGTTTTTCGCAATCTCGATGTCGTTCATCAGGGCTTTTGCGGTGGCCAGCAGGTAATCCGCTGCGCCCAGCATCTTGTGGTCGTCTTCCATGTCACCTTCGCGTACCA  >  NC_002947/2037436‑2037723
                                                                                                                                             |                                                                                                                                                  
cATGGGATGGATGGCACCGGCGGTGCCGGTGTTCGCGGGCGCGCCCGCTCCCACAGTGTGGGGTGCTTGGCCAATGCAGTGGTTTGCCCTCGCCACGTCTGGCCTGATGAAGCCCATTTCTTCTGACTGCAAGCGAGGCATAGCCTCGc                                                                                                                                             <  4:52607/149‑1 (MQ=255)
 aTGGGATGGATGGCACCGGCGGTGCCGGTGTTCGCGGGCGCGCCCGCTCCCACAGTGTGGGGTGCTTGGCCAATGCAGTGGTTTGCCCTCGCCACGTCTGGCCTGATGAAGCCCATTTCTTCTGACTGCAAGCGAGGCATAGCCTCGcc                                                                                                                                            >  7:82808/1‑149 (MQ=255)
     gatggatgGCACCGGCGGTGCCGGTGTTCGCGGGCGCGCCCGCTCCCACAGTGTGGGGTGCTTGGCCAATGCAGTGGTTTGCCCTCGCCACGTCTGGCCTGATGAAGCCCATTTCTTCTGACTGCAAGCGAGGCTATGCCTCGCTTgca                                                                                                                                        <  1:211560/149‑3 (MQ=255)
      atggatggCACCGGCGGTGCCGGTGTTCGCGGGCGCGCCCGCTCCCACAGTGTGGGGTGCTTTGCCAATGCAGTGGTTTGCCCTCGCCACGTCTGGCCTGATGAAGCCCATTTCTTCTGACTGCAAGCGAGGCATAGCCTCGCCTGGta                                                                                                                                       >  3:206032/1‑148 (MQ=255)
        ggatggCACCGGCGGTGCCGGTGTTCGCGGGCGCGCCCGCTCCCACAGTGTGGGGTGCTTGGCCAATGCAGTGGTTTGCCCTCGCCACGTCTGGCCTGATGAAGCCCATTTCTTCTGACTGCAAGCGAGGCATAGCCTCGCCTGGtctc                                                                                                                                     >  2:294614/1‑149 (MQ=255)
        ggatggCACCGGCGGTGCCGGTGTTAGCGGGCGCGCCCGCTCCCAAATTGAGGGGTGATTGGCCAATGCAGTGGTTTGCCCTCGCCACGTCTGGCCTGATGAAGCCCATTTCTTCTGACTGCAAGCGAGGCATAGCCTCGCCTGGtctc                                                                                                                                     >  1:120133/1‑149 (MQ=255)
          atggCACCGGCGGTGCCGGTGTTCGCGGGCGCGCCCGCTCCCACAGTGTGGGGTGCTTGGCCAATGCAGTGGTTTGCCCTCGCCACGTCTGGCCTGATGAAGCCCATTTCTTCTGACTGCAAGCGAGGCATAGCCTCGCCTGGTCTcgc                                                                                                                                   >  8:14543/1‑149 (MQ=255)
              cACCGGCGGTGCCGGTGTTCGCGGGCGCGCCCGCTCCCACAGTGTGGGGTGCTTGGCCAATGCAGTGGTTTGCCCTCGCCACGTCTGGCCTGATGAAGCCCATTTCTTCTGACTGCAAGCGAGGCATAGCCTCGCCTGGTCTCGCGt                                                                                                                                 >  7:123327/1‑147 (MQ=255)
              cACCGGCGGTGCCGGTGTTCGCGGGCGCGCCCGCTCCCACAGTGTGGGGTGCTTGGCCAATGCAGTGGTTTGCCCTCGCCACGTCTGGCCTGATGAAGCCCATTTCTTCTGACTGCAAGCGAGGCATAGCCTCGCCTGGTCTCGCGt                                                                                                                                 <  8:123327/147‑1 (MQ=255)
                cggggggTGCCGGTGTTCGCGGGCGCGCCCTCTCGCACAGTGGGGGGTGCTTGGCCAATGCAGTTGTTTGCCCTCGCCACGTCTGGCCTGATGAAGCCCATTTCTTCTGACTGCAAGCGAGGCATAGCCTCGCCTGGTCTCGCGTCAtt                                                                                                                             <  7:222363/144‑1 (MQ=255)
                 cggcggTGCCGGTGTTCGCGGGCGCGCCCGCTCCCACAGTGTGGGGTGCTTGGCCAATGCAGTGGTTTGCCCTCGCCACGTCTGGCCTGATGAAGCCCATTTCTTCTGACTGCAAGCGAGGCATAGCCTCGCCTGGTCTCGCGTCATTc                                                                                                                            >  5:188011/1‑149 (MQ=255)
                     ggTGCCGGTGTTCGCGGGCGCGCCCGCTCCCACAGTGTGGGGTGCTTGGCCAATGCAGTGGTTTGCCCTCGCCACGTCTGGCCTGATGAAGCCCATTTCTTCTGACTGCAAGCGAGGCATAGCCTCGCCTGGTCTCGCGTCATTCCCtt                                                                                                                        <  6:229657/149‑1 (MQ=255)
                          cGGTGTTCGCGGGCGCGCCCGCTCCCACAGTGTGGGGTGCTTGGCCAATGCAGTGGTTTGCCCTCGCCACGTCTGGCCTGATGAAGCCCATTTCTTCTGACTGCAAGCGAGGCATAGCCTCGCCTGGTCTCGCGTCATTCCCTTCAgtg                                                                                                                   >  8:143146/1‑149 (MQ=255)
                           ggTGTTCGCGGGCGCGCCCGCTCCCACAGTGTGGGGTGCTTGGCCAATGCAGTGGTTTGCCCTCGCCACGTCTGGCCTGATGAAGCCCATTTCTTCTGACTGCAAGCGAGGCATAGCCTCGCCTGGTCTCGCGTCATTCCCTTCAgtgg                                                                                                                  >  2:147229/1‑149 (MQ=255)
                            gtgtTCGCGGGCGCGCCCGCTCCCACAGTGTGGGGTGCTTGGCCAATGCAGTGGTTTGCCCTCGCCACGTCTGGCCTGATGAAGCCCATTTCTTCTGACTGCAAGCGAGGCATAGCCTCGCCTGGTCTCGCGTCATTCCCTTCAgtggt                                                                                                                 >  1:279618/1‑149 (MQ=255)
                                 cgcgGGCGCGCCCGCTCCCACAGTGTGGGGTGCTTGGCCAATGCAGTGGTTTGCCCTCGCCACGTCTGGCCTGATGAAGCCCATTTCTTCTGACTGCAAGCGAGGCATAGCCTCGCCTGGTCTCGCGTCATTCCCTTCAGTGGTGCCt                                                                                                             <  2:279618/148‑1 (MQ=255)
                                 cgcgGGCGCGCCCGCTCCCACAGTGTGGGGTGCTTGGCCAATGCAGTGGTTTGCCCTCGCCACGTCTGGCCTGATGAAGCCCATTTCTTCTGACTGCAAGCGAGGCATAGCCTCGCCTGGTCTCGCGTCATTCCCTTCAGTGGTGCCt                                                                                                             <  1:193859/148‑1 (MQ=255)
                                     ggCGCGCCCGCTCCCACAGTGTGGGGTGCTTGGCCAATGCAGTGGTTTGCCCTCGCCACGTCTGGCCTGATGAAGCCCATTTCTTCTGACTGCAAGCGAGGCTATGCCTCGCTTgca                                                                                                                                        >  6:250666/1‑115 (MQ=255)
                                     ggCGCGCCCGCTCCCACAGTGTGGGGTGCTTGGCCAATGCAGTGGTTTGCCCTCGCCACGTCTGGCCTGATGAAGCCCATTTCTTCTGACTGCAAGCGAGGCTATGCCTCGCTTgca                                                                                                                                        <  5:250666/117‑3 (MQ=255)
                                       cgcgcCCGCTCCCACAGTGTGGGGTGCTTGGCCAATGCAGTGGTTTGCCCTCGCCACGTCTGGCCTGATGAAGCCCATTTCTTCTGACTGCAAGCGAGGCATAGCCTCGCCTGGTCTCGCGTCATTCCCTTCAGTGGTGCCTGTTTTTc                                                                                                      <  7:14543/149‑1 (MQ=255)
                                              gccccccCAGTGCGGGGTGCTTGGCCCATCCAGTGCCTTGCCCTCGCCACGTCTGGCCTCCCCACCCCCACTTCTTCAGACTGCAAGCGACCCATAGCCCCGCCTGGTCTCGCGTCATTCCCTTCAGTGGTGCCTGCTTTTCGCAAtct                                                                                               <  8:140593/146‑1 (MQ=255)
                                              gCTCCCACAGTGTGGGGTGCTTGGCCAATGCAGTGGTTTGCCCTCGCCACGTCTGGCCTGATGAAGCCCATTTCTTCTGACTGCAAGCGAGGCATAGCCTCGCCTGGTCTCGCGTCATTCCCTTCAGTGGTGCCTGTTTTTCGCAAtct                                                                                               <  1:147229/149‑1 (MQ=255)
                                                  ccACAGTGTGGGGTGCTTGGCCAATGCAGTGGTTTGCCCTCGCCACGTCTGGCCTGATGAAGCCCATTTCTTCTGACTGCAAGCGAGGCATAGCCTCGCCTGGTCTCGCGTCATTCCCTTCAGTGGTGCCTGTTTTTCGCAATCTCGAt                                                                                           <  3:16248/149‑1 (MQ=255)
                                                  ccACAGTGTGGGGTGCTTGGCCAATGCAGTGGTTTGCCCTCGCCACGTCTGGCCTGATGAAGCCCATTTCTTCTGACTGCAAGCGAGGCATAGCCTCGCCTGGTCTCGCGTCATTCCCTTCAGTGGTGCCTGTTTTTCGCAATCTCGAt                                                                                           >  5:324668/1‑149 (MQ=255)
                                                   cacaGTGTGGGGTGCTTGGCCAATGCAGTGGTTTGCCCTCGCCACGTCTGGCCTGATGAAGCCCATTTCTTCTGACTGCAAGCGAGGCATAGCCTCGCCTGGTCTCGCGTCATTCCCTTCAGTGGTGCCTGTTTTTCGCAATCTCGATg                                                                                          >  7:92544/1‑149 (MQ=255)
                                                          tgGGGTGCTTGGCCAATGCAGTGGTTTGCCCTCGCCACGTCTGGCCTGATGAAGCCCATTTCTTCTGACTGCAAGCGAGGCATAGCCTCGCCTGGTCTCGCGTCATTCCCTTCAGTGGTGCCTGTTTTTCg                                                                                                     <  2:6851/131‑1 (MQ=255)
                                                          tgGGGTGCTTGGCCAATGCAGTGGTTTGCCCTCGCCACGTCTGGCCTGATGAAGCCCATTTCTTCTGACTGCAAGCGAGGCATAGCCTCGCCTGGTCTCGCGTCATTCCCTTCAGTGGTGCCTGTTTTTCg                                                                                                     >  1:6851/1‑131 (MQ=255)
                                                          tgGGGTGCTTGGCCAATGCAGTGGTTTGCCCTCGCCACGTCTGGCCTGATGAAGCCCATTTCTTCTGACTGCAAGCGAGGCATAGCCTCGCCTGGTCTCGCGTCATTCCCTTCAGTGGTGCCTGTTTTTCGCAATCTCGATGTCGTtc                                                                                    <  6:188011/148‑1 (MQ=255)
                                                              gTGCTTGGCCAATGCAGTGGTTTGCCCTCGCCACGTCTGGCCTGATGAAACCCATTTCTTCTGACTGCAAGCGAGGCATAGCCTCGCCTGGTCTCGCGTCATTCCCTTCAGTGGTGCCTGTTTTTCGCAATCTCGATGTCGTTCATCAg                                                                               >  1:82318/1‑149 (MQ=255)
                                                              gTGCTTGGCCAATGCAGTGGGTTGCCCTCGCCACGTCTGGCCTGATGAAGCCCATTTCTTCTGACTGCAAGCGAGGCATAGCCTCGCCTGGTCTCGCGTCATTCCCTTCAGTGGTGCCTGTTTTTCGCAATCTCGATGTCGTTCATCAg                                                                               <  8:82808/149‑1 (MQ=255)
                                                                   tGGCCAATGCAGTGGTTTGCCCTCGCCACGTCTGGCCTGATGAAGCCCATTTCTTCTGACTGCAAGCGAGGCATAGCCTCGCCTGGTCTCGCGTCATTCCCTTCAGTGGTGCCTGTTTTTCGCAATCTCGATGTCGTTCATCAGGGCtt                                                                          >  4:143007/1‑149 (MQ=255)
                                                                     gCCAATGCAGTGGTTTGCCCTCGCCACGTCTGGCCTGATGAAGCCCATTTCTTCTGACTGCAAGCGAGGCATAGCCTCGCCTGGTCTCGCGTCATTCCCTTCAGTGGTGCCTGTTTTTCGCAATCTCGATGTCGTTCATCAGGGCttt                                                                         <  6:324668/148‑1 (MQ=255)
                                                                       cAATGCAGTGGTTTGCCCTCGCCACGTCTGGCCTGATGAAGCCCATTTCTTCTGACTGCAAGCGAGGCATAGCCTCGCCTGGTCTCGCGTCATTCCCTTCAGTGGTGCCTGTTTTTCGCAATCTCGATGTCGTTCATCAGGGCTTTTGc                                                                      >  7:305802/1‑149 (MQ=255)
                                                                       cAATGCAGTGGTTTGCCCTCGCCACGTCTGGCCTGATGAAGCCCATTTCTTCTGACTGCAAGCGAGGCATAGCCTCGCCTGGTCTCGCGTCATTCCCTTCAGTGGTGCCTGTTTTTCGCAATCTCGATGTCGTTCATCAGGGCTTTTGc                                                                      >  1:240516/1‑149 (MQ=255)
                                                                            cAGTGGTTTGCCCTCGCGACGTCTGGCCTGATGAAGCCCATTTCTTCTGACTGCAAGCGAGGCATAGCCTCGCCTGGTCTCGCGTCATTCCCTTCAGTGGTGCCTGTTTTTCGCAATCTCGATGTCGTTCATCAGGGCTTTTGCGGTgg                                                                 >  3:184315/1‑149 (MQ=255)
                                                                                    tGCCCTCGCCACGTCTGGCCTGATGAAGCCCATTTCTTCTGACTGCAAGCGAGGCATAGCCTCGCCTGGTCTCGCGTCATTCCCTTCAGTGGTGCCTGTTTTTCGCAATCTCGATGTCGTTCATCAGGGCTTTTGCGGTGGCcagcag                                                          <  1:294614/148‑1 (MQ=255)
                                                                                        cTCGCCACGTCTGGCCTGATGAAGCCCATTTCTTCTGACTGCAAGCGAGGCATAGCCTCGCCTGGTCTCGCGTCATTCCCTTCAGTGGTGCCTGTTTTTCGCAATCTCGATGTCGTTCATCAGGGCTTTTGCGGTGGCCAGCAGGTAAt                                                     <  7:201237/149‑1 (MQ=255)
                                                                                         tCGCCACGTCTGGCCTGATGAAGCCCATTTCTTCTGACTGCAAGCGAGGCATAGCCTCGCCTGGTCTCGCGTCATTCCCTTCACTGGTGCCTGTTTTTCGCAATCTCGATGTCGTTCATCAGGGCTTTTGCGGTGGCCAGCAGGTAAt                                                     >  5:209248/1‑148 (MQ=255)
                                                                                           gCCTCGTCTGTCCTGATGATGCCCATTTCTTCTGACTGCAAGCGACGCATAGCCTAGCCTGGTCTCGCGTCATTCCCTTCAGTGGTGCCTGTTTTTCGCAATCTCGTTGTCGTTCATCAGGGCTTTTGCGGTGGCCAGCAGGTAATc                                                    <  6:209248/147‑1 (MQ=255)
                                                                                              aCGTCTGGCCTGATGAAGCCCATTTCTTCTGACTGCAAGCGAGGCATAGCCTCGCCTGGTCTCGCGTCATTCCCTTCAGTGGTGCCTGTTTTTCGCAATCTCGATGTCGTTCATCAGGGCTTTTGCGGTGGCCAGCAGGTAATCCGCTg                                               >  7:210914/1‑149 (MQ=255)
                                                                                                   tGGCCTGATGAAGCCCATTTCTTCTGACTGCAAGCGAGGCATAGCCTCGCCTGGTCTCGCGTCATTCCCTTCAGTGGTGCCTGTTTTTCGCAATCTCGATGTCGTTCATCAGGGCTTTTGCGGTGGCCAGCAGGTAATCCGCTGCGcc                                           >  6:257352/1‑148 (MQ=255)
                                                                                                       cTGATGAAGCCCATTTCTTCTGACTGCAAGCGAGGCATAGCCTCGCCTGGTCTCGCGTCATTCCCTTCAGTGGTGCCTGTTTTTCGCAATCTCGATGTCGTTCATCAGGGCTTTTg                                                                       >  5:134473/1‑116 (MQ=255)
                                                                                                       cTGATGAAGCCCATTTCTTCTGACTGCAAGCGAGGCATAGCCTCGCCTGGTCTCGCGTCATTCCCTTCAGTGGTGCCTGTTTTTCGCAATCTCGATGTCGTTCATCAGGGCTTTTg                                                                       <  6:134473/116‑1 (MQ=255)
                                                                                                        tgatgaAGCCCATTTCTTCTGACTGCAAGCGAGGCATAGCCTCGCCTGGTCTCGCGTCATTCCCTTCAGTGGTGCCTGTTTTTCGCAATCTCGATCTCGTTCATCAGGGCTTTTGCGGTGGCCAGCAGGTAATCCGCTGCGCCCAGCa                                      >  8:266562/1‑148 (MQ=255)
                                                                                                         gatgaAGCCCCTTTCTTCTGACTGCAAACAAGACATAGCCTCGCCTGGTCTCGAGTCATTACATTCAGTGGTGCATGTTTTTCGCAATCTCGATGTCGTTCATCAGGTCTTTTGCGGTGGCCAGCAGGTAATACGCTGCGCACAGCATc                                    >  8:48665/1‑149 (MQ=255)
                                                                                                          atgaAGCCCATTTCTTCTGACTGCAAGCGAGGCATAGCCTCGCCTGGTCTCGCGTCATTCCCTTCAGTGGTGCCTGTTTTTCGCAATCTCGATGTCGTTCATCAGGGCTTTTGCGGTGGCCAGCAGGTAATCCGCTGCGCCCAGCATCt                                   >  2:154183/1‑149 (MQ=255)
                                                                                                           tgaAGCCCATTTCTTCTGACTGCAAGCGAGGCATAGCCTCGCCTGGTCTCGCGTCATTCCCTTCAGTGGTGCCTGTTTTTCGCAATCTCGATGTCGTTCATCAGGGCTTTTGCGGTGGCCAGCAGGTAATCCGCTGCGCCCAGCATCtt                                  <  8:210914/149‑1 (MQ=255)
                                                                                                           ggaAGCCCATTTCTTCTGACTGCAAGCGAGGCATAGCCTCGCCTGGTCTCGCGTCATTCCCTTCAGTGGTGCCTGTTTTTCGCAATCTCGATGTCGTTCATCAGGGCTTTTGCGGTGGCCAGCAGGTAATCCGCTGCGCCCAGCATCtt                                  <  5:99710/148‑1 (MQ=255)
                                                                                                              aGCCCATTTCTTCTGACTGCAAGCGAGGCATAGCCTCGCCTGGTCTCGCGTCATTCCCTTCAGTGGTGCCTGTTTTTCGCAATCTCGATGTCGTTCATCAGGGCTTTTGCGGTGGCCAGCAGGTAATCCGCTGCGCCCAGCATCTtgtg                               <  4:184315/149‑1 (MQ=255)
                                                                                                                 ccATTTCTTCTGACTGCAAGCGAGGCATAGCCTCGCCTGGTCTCGCGTCATTCCCTTCAGTGGTGCCTGTTTTTCGCAATCTCGATGTCGTTCATCAGGGCTTTTGCGGTGGCCAGCAGGTAATCCGCTGCGCCCAGCATCTTGTGgtc                            <  2:240516/149‑1 (MQ=255)
                                                                                                                   aTTTCTTCTGACTGCAAGCGAGGCATAGCCTCGCCTGGTCTCGCGTCATTCCCTTCAGTGGTGCCTGTTTTTCGCAATCTCGATGTCGTTCATCAGGGCTTTTGCGGTGGCCAGCAGGTAATCCGCTGCGCCCAGCATCTTGTGgtcgt                          >  6:208216/1‑149 (MQ=255)
                                                                                                                   aTTTCTTCTGACTGCAAGCGAGGCATAGCCTCGCCTGGTCTCGCGTCATTCCCTTCAGTGGTGCCTGTTTTTCGCAATCTCGATGTCGTTCATCAGGGCTTTTGCGGTGGCCAGCAGGTAATCCGCTGCGCCCAGCATCTTGTGgtcgt                          >  3:235131/1‑149 (MQ=255)
                                                                                                                        ttctGACTGCAAGCGAGGCATAGCCTCGCCTGGTCTCGCGTCATTCCCTTCAGTGGTGCCTGTTTTTCGCAATCTCGATGTCGTTCATAAGGGCTTTTGCGGTGGCCAGCAGGTACTCCGCTGCGCCCAGCATCTTGTGGTCGTCTTc                      >  5:85714/1‑148 (MQ=255)
                                                                                                                          ctGACTGCAAGCGAGGCATAGCCTCGCCTGGTCTCGCGTCATTCCCTTCAGTGGTGCCTGTTTTTCGCAATCTCGATGTCGTTCATCAGGGCTTTTGCGGTGGCCAGCAGGTAATCCGCTGCGCCCAGCATCTTGTGGTCGTCTTc                      <  1:154183/146‑1 (MQ=255)
                                                                                                                          ctGACTGCAAGCGAGGCATAGCCTCGCCTGGTCTCGCGTCATTCCCTTCAGTGGTGCCTGTTTTTCGCAATCTCGATGTCGTTCATCAGGGCTTTTGCGGTGGCCAGCAGGTAATCCGCTGCGCCCAGCATCTTGTGGTCGTCTTCCa                    <  5:257352/148‑1 (MQ=255)
                                                                                                                           tGACTGCAAGCGAGGCATAGCCTCGCCTGGTCTCGCGTCATTCCCTTCAGTGGTGCCTGTTTTTCGCAATCTCGATGTCGTTCATCAGGGCTTTTGCGGTGGCCAGCAGGTAATCCGCTGCGCCCAGCATCTTGTGGTCGTCTTCCATg                  >  5:143590/1‑149 (MQ=255)
                                                                                                                                 ctaGCGAGGCATAGCCTCGCCTGGTCTCGCGTCATTCCCTTCAGTGGTGCCTGTTTTTCGCAATCTCGATGTCGTTCATCAGGGCTTTTGCGGTGGCCAGCAGGTAATCCGCTGCGCCCAGCATCTTGTGGTCGTCTTCCATGTCACCt            <  2:120133/147‑1 (MQ=255)
                                                                                                                                    gCGAGGCTATGCCTCGCCTGGTCTCGCGTCATTCCCTTCAGTGGTGCCTGTTTTTCGCAATCTCGATGTCGTTCATCAGGGCTTTTGCGGTGGCCAGCAGGTAATCCGCTGCGCCCAGCATCTTGTGGTCGTCTTCCATGTCACCTTcg         >  8:10417/1‑149 (MQ=255)
                                                                                                                                    gCGAGGCATAGCCTCGCCTGGTCTCGCGTCATTCCCTTCAGTGGTGCCTGTTTTTCGCAATCTCGATGTCGTTCATCAGGGCTTTTGCGGTGGCCAGCAGGTAATCCGCTGCGCCCAGCATCTTGTGGTCGTCTTCCATGTCACCTTcg         >  3:67298/1‑149 (MQ=255)
                                                                                                                                      gAGGCATAGCCTCGCCTGGTCTCGCGTCATTCCCTTCAGTGGTGCCTGTTTTTCGCAATCTCGATGTCGTTCATCAGGGCTTTTGCGGTGGCCAGCAGGTAATCCGCTGCGCCCAGCATCTTGTGGTCGTCTTCCATGTCACCTTcgcg       <  7:143146/149‑1 (MQ=255)
                                                                                                                                           aTAGCCTCGCCTGGTCTCGCGTCATTCCCTTCAGTGGTGCCTGTTTTTCGCAATCTCGATGTCGTTCATCAGGGCTTTTGCGGTGGCCAGCAGGTAATCCGCTGCGCCCAGCATCTTGTGGTCGTCTTCCATGTCACCTTCGCGTACCa  <  2:82318/149‑1 (MQ=255)
                                                                                                                                             aGCCTCGCCTGGTCTCGCGTCATTCCCTTCAGTGGTGCCTGTTTTTCGCAATCTCGATGTCGTTCATCAGGGCTTTTGCGGTGGCCAGCAGGTAATCCGCTGCGCCCAGCATCTTGTGGTCGTCTTCCATGTCACCt            <  6:99159/137‑1 (MQ=255)
                                                                                                                                             aGCCTCGCCTGGTCTCGCGTCATTCCCTTCAGTGGTGCCTGTTTTTCGCAATCTCGATGTCGTTCATCAGGGCTTTTGCGGTGGCCAGCAGGTAATCCGCTGCGCCCAGCATCTTGTGGTCGTCTTCCATGTCACCt            >  5:99159/1‑137 (MQ=255)
                                                                                                                                             |                                                                                                                                                  
CATGGGATGGATGGCACCGGCGGTGCCGGTGTTCGCGGGCGCGCCCGCTCCCACAGTGTGGGGTGCTTGGCCAATGCAGTGGTTTGCCCTCGCCACGTCTGGCCTGATGAAGCCCATTTCTTCTGACTGCAAGCGAGGCATAGCCTCGCCTGGTCTCGCGTCATTCCCTTCAGTGGTGCCTGTTTTTCGCAATCTCGATGTCGTTCATCAGGGCTTTTGCGGTGGCCAGCAGGTAATCCGCTGCGCCCAGCATCTTGTGGTCGTCTTCCATGTCACCTTCGCGTACCA  >  NC_002947/2037436‑2037723

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