Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 58,750 T→A 14.8% intergenic (‑661/+27) PP_5424 ← / ← PP_0050 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00294758,7500TA14.8% 182.8 / 21.5 74intergenic (‑661/+27)PP_5424/PP_0050hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (37/26);  new base A (4/7);  total (41/33)
Fisher's exact test for biased strand distribution p-value = 2.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.65e-01

AGTAGGGATTTCTCGGCTTTTCTGTGAAAAGCACCCCTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTCACATACAGCCAACGATGATCACGAGGTGGCGCATACAGGTGGCGCGCTTTCCA  >  NC_002947/58604‑58896
                                                                                                                                                   |                                                                                                                                                   
aGTAGGGATTTCTCGGCTTTTCTGTGAAAAGCACCCCTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCagcagc                                                                                                                                                   >  7:139783/1‑149 (MQ=255)
   aGGGATTTCTCGGCTTTTCTGTGAAAAGCACCCCTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCgg                                                                                                                                                 <  7:199165/148‑1 (MQ=255)
      gATTTCTCGGCTTTTCTGTGAAAAGCACCCCTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCAGCAGCGGGCtt                                                                                                                                             >  2:62674/1‑149 (MQ=255)
       aTTTCTCGGCTTTTCTGTGAAAAGCACCCCTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTGTTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCttt                                                                                                                                            >  6:326495/1‑149 (MQ=255)
        aaggggtgCTTTTCACAGAAAAGCACCCCTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAAGGTAATATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCtttt                                                                                                                                           >  1:269170/8‑149 (MQ=255)
           ctgggCTTTTCTGTGAAAAGCAGCCCTTCATTTCACAAAAAAATTGGCTAAGATTCCAATTTCACGGGTTCTGGCCTATCACCTATTTTGTGAACCTAGAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCAGCAGCGGGCTTTTTa                                                                                                                                         <  4:46620/145‑1 (MQ=255)
           ctcGTCTTTTCTGTGAAAAGCACCCCTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGGCACTAAAAAGCCCGCAGCAGCGGGCTTTTTa                                                                                                                                         <  8:237034/148‑1 (MQ=255)
           ctcGGCTTTTCTGTGAAAAGCACCCCTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCAGCAGCGGGCTTTTTa                                                                                                                                         <  8:139783/148‑1 (MQ=255)
           atcGGCTTTTCTGTTAAAAGCACCCCTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGTGTCTGGCCTATCACCTATTTTTTGAACCTACAAGGTTTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTa                                                                                                                                         <  8:45849/147‑1 (MQ=255)
            tcGGCTTTTCTGTGAAAAGCACCCCTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGt                                                                                                                                       >  2:95352/1‑149 (MQ=255)
             gtgCTTTTCACAGAAAAGCACCCCTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTg                                                                                                                                      >  2:172025/3‑149 (MQ=255)
              ggCTTTTCTGTGAAAAGCACCCCTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAAGGTGCTAATTCGCCAGACACTAAAAAGCCCGCAGCAGCGGGCTTTTTa                                                                                                                                         <  1:62674/146‑1 (MQ=255)
                 ttttCTGTGAAAAGCACCCCTTCATTTCACAAAAAAAATGGCTAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAAGGTGTAATTAGCCAGACACTAAAAAGCCCGCTGCTGCGGGGTTTTTAGTGACTg                                                                                                                                  >  4:174631/1‑149 (MQ=255)
                      tgtgAAAAGCACCCCTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCaaa                                                                                                                              >  2:209100/1‑148 (MQ=255)
                        tgAAAAGCACCCCTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGGCACTAAAAAGCCCGCAGCAGCGGGCTTTTTAGTGTCTggcgaa                                                                                                                             <  4:199866/147‑6 (MQ=38)
                          aaaaGCACCCCTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAAt                                                                                                                          >  2:71036/1‑148 (MQ=255)
                          aaaaGCACCCCTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATc                                                                                                                         <  3:251160/149‑1 (MQ=255)
                           aaaGCACCCCTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCAGCAGCGGGCttt                                                                                                                                            <  6:248883/129‑1 (MQ=255)
                           aaaGCACCCCTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCAGCAGCGGGCttt                                                                                                                                            >  5:248883/1‑129 (MQ=255)
                               cACCCCTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAg                                                                                                                   <  2:69051/149‑1 (MQ=255)
                                  cccTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTa                                                                                                                                         >  5:256708/1‑125 (MQ=255)
                                  cccTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTa                                                                                                                                         <  6:256708/125‑1 (MQ=255)
                                  cccTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCg                                                                                                                >  2:82575/1‑149 (MQ=255)
                                         tttCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGaa                                                                                                         >  6:81187/1‑149 (MQ=255)
                                                           gAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGc                                                                                       >  6:583/1‑149 (MQ=255)
                                                           gAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGACCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCAGGTGAGAATGTTAACTATcaca                                                                                           >  6:587/1‑145 (MQ=255)
                                                             aGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTg                                                                                     <  5:81187/149‑1 (MQ=255)
                                                                   cAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGc                                                                                                                                     <  6:277970/96‑1 (MQ=255)
                                                                   cAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGc                                                                                                                                     >  5:277970/1‑96 (MQ=255)
                                                                   cAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATg                                                                               <  1:95352/149‑1 (MQ=255)
                                                                            ggggTCTGGCCTATCACCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAg                                                                       >  4:223709/1‑148 (MQ=255)
                                                                              ggTCTGGCCTATCACCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTACGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGaa                                                                     <  1:71036/149‑1 (MQ=255)
                                                                                    gCCTATCACCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGt                                                               <  2:246046/148‑1 (MQ=255)
                                                                                    gCCTATCACCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTc                                                              >  7:3999/1‑149 (MQ=255)
                                                                                      gtATCAGCTATTTTGTGAACCTACAAGGTGTAATTCGCCCGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCACACAATGAGTAGCCGGTGAGAATGTTCACGAGCACACCGCTGGCATTGGCCACCAGACCTTAGTGGc                                                            >  4:40314/2‑149 (MQ=255)
                                                                                           aCCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCAtt                                                        >  5:59409/1‑148 (MQ=255)
                                                                                           aCCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTc                                                       <  1:172025/149‑1 (MQ=255)
                                                                                           aCCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTc                                                       <  2:269170/149‑1 (MQ=255)
                                                                                            ccTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTca                                                      >  8:243980/1‑149 (MQ=255)
                                                                                            ccTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTca                                                      >  1:316116/1‑149 (MQ=255)
                                                                                             ctgttttGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTcac                                                     <  5:326495/146‑1 (MQ=255)
                                                                                             cTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTca                                                      >  1:13413/1‑148 (MQ=255)
                                                                                             cTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGACCACCAGAACATAGTCGCCATTca                                                      >  7:160448/1‑148 (MQ=255)
                                                                                             cTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGACCACCAGAACATAGTCGCCATTca                                                      >  3:8912/1‑148 (MQ=255)
                                                                                                  ttGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTCACATAc                                                 <  6:59409/148‑1 (MQ=255)
                                                                                                   tgtgAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTCACATACa                                                <  7:243980/148‑1 (MQ=255)
                                                                                                   tgtgAAACTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCAATTCACATACa                                                <  6:241817/148‑1 (MQ=255)
                                                                                                     tgAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGACCACCAGAACATAGTCGCCATTCACATACAGcc                                             <  4:8912/149‑1 (MQ=255)
                                                                                                      gAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGACCACCAGAACATAGTCGCCATTCACATACAGcc                                             <  8:160448/148‑1 (MQ=255)
                                                                                                          cTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGGAATGGCCACCAGAACATAGTCGCCATTCACATACAGCCAACga                                        >  2:221846/1‑149 (MQ=255)
                                                                                                            aCAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTc                                                       >  2:134307/1‑132 (MQ=255)
                                                                                                            aCAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTCACATACAGCCAACgat                                       >  5:168775/1‑148 (MQ=255)
                                                                                                                ggTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTCACATACAGCCAACGATGATCa                                  >  8:104363/1‑149 (MQ=255)
                                                                                                                ggTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTCACATACAGCCAACGATGATCa                                  >  8:174072/1‑149 (MQ=255)
                                                                                                                ggTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCATGCAAACAATAAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGGAGACATTCAAAAACAGCAAAAgaatag                                    >  4:241725/1‑143 (MQ=255)
                                                                                                                  tgtAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAATCAATCAGTAGCCGGTGAGAATGTTGACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTCACATACAGCCAACGATGATCACg                                >  7:129068/1‑149 (MQ=255)
                                                                                                                  tgtAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGATTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTCACATACAGCCAACGATGATCACg                                <  6:168775/149‑1 (MQ=255)
                                                                                                                   gtAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTCACATACAGCCAACGATGATCAc                                 >  3:181300/1‑147 (MQ=255)
                                                                                                                    ctAAATCGCCAGACACTAAAAACCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCAtt                                                        <  4:46938/123‑1 (MQ=255)
                                                                                                                     tAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTCACATACAGCCAACGATGATCACGa                               <  2:316116/148‑1 (MQ=255)
                                                                                                                      aaTTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTcaca                                                    <  1:256750/126‑1 (MQ=255)
                                                                                                                      aaTTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTcaca                                                    >  2:256750/1‑126 (MQ=255)
                                                                                                                       aTTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTCACATACAGCCAACGATGATCACGAGGt                            <  8:3999/149‑1 (MQ=255)
                                                                                                                        ttCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTCACATACa                                                >  4:92899/1‑128 (MQ=255)
                                                                                                                        ttCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTCACATACa                                                <  3:92899/128‑1 (MQ=255)
                                                                                                                           gCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTCACATACAGCCACCGATGATCACGAGGTGgcg                        >  3:38407/1‑149 (MQ=255)
                                                                                                                               cgatcTAAAAAGCCGGCAGCAGCGGGCTTTTTAGTGCCTTCAAACTATCAGTTGGCGGTGAGTATGTTCACGATCACACCGATGTCAATGGCCACCAGAAAATAGTCGACATTCACATACAGGCAACGATGATCACTATTTGGCGCATa                    <  8:178046/145‑1 (MQ=255)
                                                                                                                                    tAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTCACATACCGCCAACGATGATCACGAGGTg                           >  7:261304/1‑137 (MQ=255)
                                                                                                                                    tAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTCACATACCGCCAACGATGATCACGAGGTg                           <  8:261304/137‑1 (MQ=255)
                                                                                                                                    tAAAAAGCCCGCAGCAGCGGGCTTTTTAGTGTCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTCACATACAGCCAACGATGATCACGAGGTGGCGCATACAGGt               >  2:310708/1‑149 (MQ=255)
                                                                                                                                      aaaaGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTCACATACAGCCAACGATGATCACGAGGTGGCGCATACAg                 <  2:247549/145‑1 (MQ=255)
                                                                                                                                      aaaaGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTCACATACAGCCAACGATGATCACGAGGTGGCGCATACAg                 >  1:247549/1‑145 (MQ=255)
                                                                                                                                      aaaaGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTCACATACAGCCAACGATGATCACGAGGTGGCGCATACAGGTgg             >  6:327107/1‑149 (MQ=255)
                                                                                                                                      aaaaGCCCGCAGCAGCGGGCTTTTTAGTGCCTGCAAAAAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGATGGCAATGGCCACCAGAACATAGTCGCCATTCACATACAGCCAACGATGATCACGAGGTGGCGCATACAgagcg             >  7:178046/1‑145 (MQ=255)
                                                                                                                                             cGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTCACATACAGCCAACGATGATCACGAGGTGGCGCATACAGGTGGCGCGCtt      >  6:158473/1‑149 (MQ=255)
                                                                                                                                               ctgctgCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTCACATACAGCCAACGATGATCACGAGGTGGCGCATACAGGTGGCGCGCTTTc    <  5:158473/149‑1 (MQ=255)
                                                                                                                                                  ctgcGGTCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGATCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTCACATACAGCCAACGATTATCACGAGGTGGCGCATACAGGTGGCGCGCTTTCca  <  5:587/148‑1 (MQ=255)
                                                                                                                                                  ctgcGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCCTTCACATACAGCCAACGATGATCACGAGGTGGCGCATACAGGTGGCGCGCTTTCca  <  5:583/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
AGTAGGGATTTCTCGGCTTTTCTGTGAAAAGCACCCCTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTCACATACAGCCAACGATGATCACGAGGTGGCGCATACAGGTGGCGCGCTTTCCA  >  NC_002947/58604‑58896

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: