Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 58,750 T→A 13.1% 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,7500TA13.1% 186.3 / 17.1 75intergenic (‑661/+27)PP_5424/PP_0050hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (32/33);  new base A (6/4);  total (38/37)
Fisher's exact test for biased strand distribution p-value = 7.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGTAGGGATTTCTCGGCTTTTCTGTGAAAAGCACCCCTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTCACATACAGCCAACGATGATCACGAGGTGGCGCATACAGGTGGCGCGCTTTCC  >  NC_002947/58604‑58895
                                                                                                                                                  |                                                                                                                                                 
aGTAGGGATTTCTCGGCTTTTCTGTGAAAAGCACCCCTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCtgctgc                                                                                                                                                 <  2:38916/149‑1 (MQ=255)
 gTAGGGATTTCTCGGCTTTTCTGTGAAAAGCACCCCTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGTGGTCTGGCCTATCACCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCtgctgc                                                                                                                                                 <  4:166314/148‑1 (MQ=255)
 gTAGGGATTTCTCGGCTTTTCTGTGAAAAGCACCCCTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCg                                                                                                                                                <  1:207357/149‑1 (MQ=255)
 gTAGGGATTTCTCGGCTTTTCTGTGAAAAGCACCCCTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCAGCAGCg                                                                                                                                                >  5:186745/1‑149 (MQ=255)
  tAGGGATTTCTCGGCTTTTCTGTGAAAAGCACCCCTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCAGCAGCgg                                                                                                                                               >  7:205298/1‑149 (MQ=255)
        tttCTCGGCTTTTCTGTGAAAAGCACCCCTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCtttt                                                                                                                                         <  6:198085/149‑1 (MQ=255)
         ttCTCGGCTTTTCTGTGAAAAGCACCCCTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCAGCAGCGGGCtttt                                                                                                                                         <  8:205298/148‑1 (MQ=255)
          tctcAGCGTTTCTGTGAAAAGCACCCCTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCttttt                                                                                                                                        <  2:11877/148‑1 (MQ=255)
           ctcGGCTTTTCTGTGAAAAGCACCCCTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCtttt                                                                                                                                         >  2:200966/1‑146 (MQ=255)
           ctcGGCTTTTCTGTGAAAAGCACCCCTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTa                                                                                                                                       <  1:200966/148‑1 (MQ=255)
           ctcGGCTTTTCTGTGAAAAGCACCCCTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCAGCAGCGGGCTTTTTa                                                                                                                                       <  6:186745/148‑1 (MQ=255)
            acGGCTATTCTGTGAAAAGCACCCCTTCATTTCACAAAAAAAATGGCCAAGATTCCAATTTCACCGGGTCTGGCCTATCACCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGt                                                                                                                                     <  8:140675/148‑1 (MQ=255)
                 ttttCTGTGAAAAGCACCCCTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTg                                                                                                                                >  7:163199/1‑149 (MQ=255)
                      tgtgAAAAGCACCCCTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTa                                                                                                                                       >  1:72594/1‑137 (MQ=255)
                      tgtgAAAAGCACCCCTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTa                                                                                                                                       <  2:72594/137‑1 (MQ=255)
                           aaaGCACCCCTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAAt                                                                                                                        >  5:68798/1‑147 (MQ=255)
                                  cccTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCg                                                                                                               >  4:17543/1‑149 (MQ=255)
                                    cTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGt                                                                                                             >  4:110303/1‑149 (MQ=255)
                                       cATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTga                                                                                                           <  6:68798/148‑1 (MQ=255)
                                         tttCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGaa                                                                                                        >  7:28146/1‑149 (MQ=255)
                                             acaAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGtt                                                                                                    <  2:89487/149‑1 (MQ=255)
                                             acaAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGtt                                                                                                    >  3:10515/1‑149 (MQ=255)
                                             acaAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCAGCAGCGGGCTTTTTa                                                                                                                                       <  4:62719/114‑1 (MQ=255)
                                             acaAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCAGCAGCGGGCTTTTTa                                                                                                                                       >  3:62719/1‑114 (MQ=255)
                                                   aaaaTGGCGAAGATTCCTATTTCACGGGGTCTTGCCTATCACCTATTTTGTGAACCTACAAGTTGTTATTCTCGAGACACTAAAAAGGCAGCTGCTGCTGGCTTTTTAGTGCCTGCATACAATCAGTAGCCGGTGAGAATGTTCACGAt                                                                                              <  5:44456/149‑1 (MQ=255)
                                                   aaaaTGGCGAAGATTCCAATTTCACTGGGTCTGTCCTATCACCTATTTTGTGAACCTACAAGTTGTAATTCGCCAGACTCTAAAAATCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGAt                                                                                              <  5:40395/149‑1 (MQ=255)
                                                     aaTGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATca                                                                                            <  3:123995/149‑1 (MQ=255)
                                                            aaGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCt                                                                                     <  3:17543/149‑1 (MQ=255)
                                                            aaGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCt                                                                                     <  1:160238/149‑1 (MQ=255)
                                                             aGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTg                                                                                    <  8:28146/149‑1 (MQ=255)
                                                              gATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTgg                                                                                   <  3:110303/149‑1 (MQ=255)
                                                                 tCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCaa                                                                                >  2:14810/1‑149 (MQ=255)
                                                                  ccAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTg                                                                                                            <  2:243582/120‑1 (MQ=255)
                                                                  ccAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTg                                                                                                            >  1:243582/1‑120 (MQ=255)
                                                                   cAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACAcc                                                                                        <  7:191758/139‑1 (MQ=255)
                                                                   cAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACAcc                                                                                        >  8:191758/1‑139 (MQ=255)
                                                                    aaTTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATgg                                                                             >  7:94586/1‑149 (MQ=255)
                                                                     aTTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGc                                                                            >  8:254135/1‑149 (MQ=255)
                                                                     aTTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGc                                                                            >  1:48361/1‑149 (MQ=255)
                                                                                    gCCTATCAACTATTTTGTGATCCTACAAGGTTTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTTAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTc                                                             <  6:256333/149‑1 (MQ=255)
                                                                                     ccTATCACCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCg                                                            <  1:14810/149‑1 (MQ=255)
                                                                                           aCCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCAtt                                                       >  6:245715/1‑148 (MQ=255)
                                                                                           aCCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCAtt                                                       >  3:157994/1‑148 (MQ=255)
                                                                                           aCCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTc                                                      <  2:257047/149‑1 (MQ=255)
                                                                                           aCCTATTTTGTGAACCTACAAGGTGTAAATCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTATTAGTGCCTGCAAACAATCAGTTGCCGGTGAGAATGTTCACTATCACACCGCTGGAAATGGCAACAAGAAGATAGACGCCAtt                                                       >  8:197959/1‑148 (MQ=255)
                                                                                            ccTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTca                                                     >  5:29333/1‑149 (MQ=255)
                                                                                             cTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTcac                                                    <  2:48361/149‑1 (MQ=255)
                                                                                             cTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTca                                                     >  6:47969/1‑148 (MQ=255)
                                                                                              tATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTcaca                                                   >  6:11526/1‑149 (MQ=255)
                                                                                               aTTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTcaca                                                   >  1:154970/1‑148 (MQ=255)
                                                                                                   tgtgAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTCACATACa                                               <  8:163199/148‑1 (MQ=255)
                                                                                                   tgtgAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTCACATACa                                               <  4:10515/148‑1 (MQ=255)
                                                                                                   tgtgAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCACCATTCACATACAg                                              >  2:114078/1‑149 (MQ=255)
                                                                                                       aaCCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATg                                                                              >  7:162787/1‑113 (MQ=255)
                                                                                                       aaCCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATg                                                                              <  8:162787/113‑1 (MQ=255)
                                                                                                             cAAGGTGTAATTCGCCAGACACTAAAAAGCACGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTCACATACAGCCAACgat                                      <  2:154970/147‑1 (MQ=255)
                                                                                                              aaGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTCACATACAGCCAACgatgat                                   <  4:157994/149‑1 (MQ=255)
                                                                                                               aGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTGACAAACAGCCAACgatgat                                   >  7:229704/1‑148 (MQ=255)
                                                                                                                ggTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTCACATACAGCCAACGATGATCa                                 >  8:156876/1‑149 (MQ=255)
                                                                                                                    tAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTCACATACAGCCAACGATGATCACGa                              <  8:94586/148‑1 (MQ=255)
                                                                                                                      aTTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTCACATACAGc                                             <  8:216854/131‑1 (MQ=255)
                                                                                                                      aTTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACAGCTGGCAATGGCCACCAGAACATAGTCGCCATTCACATACAGc                                             >  7:216854/1‑131 (MQ=255)
                                                                                                                          gCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTCACATACAGCCAACGATGATCACGAGGTGgc                        <  5:47969/148‑1 (MQ=255)
                                                                                                                           ccAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTCACATACAGCCAACGATGATCACGAgg                            >  6:178727/1‑143 (MQ=255)
                                                                                                                           ccAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTCACATACAGCCAACGATGATCACGAgg                            <  5:178727/143‑1 (MQ=255)
                                                                                                                              ggcacTAAAAAGCCCGCAGCAGCGGGCTTTTTAGTGTCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAATCGCCATTCACATACAGCCAACGATGATCACGAGGTGGCGCAt                    >  3:206404/3‑148 (MQ=255)
                                                                                                                                 acTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTCACATACAGCCAACGATGATCACGAGGTGGCGCATACAg                <  7:156876/149‑1 (MQ=255)
                                                                                                                                  cttaaaaGCGCTCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTCACATACAGGCAACGATGATCACGAGGTGTCGCTTACAgg               <  6:275331/146‑1 (MQ=255)
                                                                                                                                   tAAAAAGCCCGCAGCAGCGGGCTTTTTAGTGTCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTCACATACAGCCAACGATGATCACGAGGTGGCGCATACAGGt              >  7:103616/1‑149 (MQ=255)
                                                                                                                                     aaaaGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTc                                                                                                   >  5:50830/1‑62 (MQ=255)
                                                                                                                                     aaaaGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTc                                                                                                   <  6:50830/62‑1 (MQ=255)
                                                                                                                                     aaaaGCCCGCAGCAGCGGGCTTTTTAGTGTCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTCACATACAGCCAACgatg                                     >  8:246614/1‑124 (MQ=255)
                                                                                                                                     aaaaGCCCGCAGCAGCGGGCTTTTTAGTGTCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTCACATACAGCCAACgatg                                     <  7:246614/124‑1 (MQ=255)
                                                                                                                                       aaGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTCACATACAGCCAACGATGATCACGAGGTGGCGCATACAGGTGgcg          >  8:53869/1‑149 (MQ=255)
                                                                                                                                          cccGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTCACATACAGCCAACGATGATCACGAGGTGGCGCATACAGGTGgcgcgc       >  8:185781/1‑149 (MQ=255)
                                                                                                                                          cccGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTCACATACAGCCAACGATGATCACGAGGTGGCGCATACAGGTGgcgcgc       >  3:224803/1‑149 (MQ=255)
                                                                                                                                            cGCTGCTGCGGGCTTTTAAGAGCCTGCAAACAATCAGAAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATACTCGCCATTCACATACAGCCAACGATGATCACGAGGTGGCGCATACAGCTGGCGCGCtt     >  8:97401/1‑149 (MQ=255)
                                                                                                                                                 ctgcGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTCACATACAGCCAACGATGATCACGAGGTGGCGCATACAGGTGGCGCGCTTTcc  <  7:254135/147‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
AGTAGGGATTTCTCGGCTTTTCTGTGAAAAGCACCCCTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTCACATACAGCCAACGATGATCACGAGGTGGCGCATACAGGTGGCGCGCTTTCC  >  NC_002947/58604‑58895

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: