Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 58,747 T→A 14.6% intergenic (‑658/+30) 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,7470TA14.6% 184.9 / 16.3 74intergenic (‑658/+30)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 = 1.00e+00

GGTGAGTAGGGATTTCTCGGCTTTTCTGTGAAAAGCACCCCTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTCACATACAGCCAACGATGATCACGAGGTGGCGCATACAGGTGGCGCGCTTTC  >  NC_002947/58600‑58894
                                                                                                                                                    |                                                                                                                                                    
ggTGAGTAGGGATTTCTCGGCTTTTCTGTGAAAAGCACCCCTTCATTTCACAAATAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCgctg                                                                                                                                                     >  1:246046/1‑149 (MQ=255)
    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)
                                                                                                                                                    |                                                                                                                                                    
GGTGAGTAGGGATTTCTCGGCTTTTCTGTGAAAAGCACCCCTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTCACATACAGCCAACGATGATCACGAGGTGGCGCATACAGGTGGCGCGCTTTC  >  NC_002947/58600‑58894

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: