Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,425,141 G→C 6.8% intergenic (‑199/‑601) PP_5647 ← / → PP_3917 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,425,1410GC 6.8% 183.7 / 4.6 59intergenic (‑199/‑601)PP_5647/PP_3917hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (29/26);  new base C (2/2);  total (31/28)
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.82e-01

ATTTCGCTTAAATGCGACCGGTGGAAAATGATCGGCCCTAATCCTTGCTGATAGATATCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCATTGGATGCCCTCGATGGGTCCGGTGGGCCTTGTTTCTGCGCGCCTCGAAA  >  NC_002947/4424998‑4425289
                                                                                                                                               |                                                                                                                                                    
aTTTCGCTTAAATGCGACCGGTGGAAAATGATCGGCCCTAATCCTTGCTGATAGATATCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGgggct                                                                                                                                                 >  5:292252/1‑149 (MQ=255)
       ttAAATGCGACCGGTGGAAAATGATCGGCCCTAATCCTTGCTGATAGATATCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCgg                                                                                                                                                     >  1:330232/1‑138 (MQ=255)
       ttAAATGCGACCGGTGGAAAATGATCGGCCCTAATCCTTGCTGATAGATATCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCgg                                                                                                                                                     <  2:330232/138‑1 (MQ=255)
       ttAAATGCGACCGGTGGAAAATGATCGGCCCTAATCCTTGCTGATAGATATCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTc                                                                                                                                          >  1:318615/1‑149 (MQ=255)
       ttAAATGCGACCGGTGGAAAATGATCGGCCCTAATCCTTGCTGATAGATATCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTc                                                                                                                                          >  2:149778/1‑149 (MQ=255)
         aaaTGCGACCGGTGGAAAATGATCGGCCCTAATCCTTGCTGATAGATATCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCt                                                                                                                                         <  3:5273/148‑1 (MQ=255)
           aTGCGACCGGTGGAAAATGATCGGCCCTAATCCTTGCTGATAGATATCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTc                                                                                                                                      <  3:89124/149‑1 (MQ=255)
                 ccGGTGGAAAATGATCGGCCCTATTCCTTGCTGATAGATATCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATAACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCATAAAGCCCAGCGGGGGGCTGGGCTTCTGTcagga                                                                                                                                  >  8:66296/1‑145 (MQ=255)
                            tGATCGGCCCTAATCCTTGCTGATAGATACCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGa                                                                                                                     <  5:5834/149‑1 (MQ=255)
                                  gCCCTAATCCTTGCTGATAGATATCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGc                                                                                                               <  5:93724/149‑1 (MQ=255)
                                  gCCCTAATCCTTGCTGATAGATATCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGc                                                                                                               >  6:93724/1‑149 (MQ=255)
                                    ccTAATCCTTGCTGATAGATATCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTcagca                                                                                                                                  <  6:211402/128‑1 (MQ=255)
                                    ccTAATCCTTGCTGATAGATATCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTcagca                                                                                                                                  >  5:211402/1‑128 (MQ=255)
                                        aTCCTTGCTGATAGATATCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCaa                                                                                                         <  1:70461/149‑1 (MQ=255)
                                           cTTGCTGATAGATATCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAAtt                                                                                                       <  1:314480/148‑1 (MQ=255)
                                                 gatagataTCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTg                                                                                                >  4:122950/1‑149 (MQ=255)
                                                  atagataTCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTgg                                                                                               >  4:154796/1‑149 (MQ=255)
                                                  atagataTCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTgg                                                                                               >  8:179510/1‑149 (MQ=255)
                                                      atatCCGCGGGAGAGCGCCAGTAGAGAACCGAGCCGGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAATAGCCCAGCGCGTGGCTGGGCTTCAGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGc                                                                                           <  8:90360/149‑1 (MQ=255)
                                                         tCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTa                                                                                         >  4:170398/1‑148 (MQ=255)
                                                         tCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTAc                                                                                        <  6:292252/149‑1 (MQ=255)
                                                               ggACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAAt                                                                                  <  1:221929/149‑1 (MQ=255)
                                                                 aCAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAAt                                                                                  >  7:274935/1‑147 (MQ=255)
                                                                 aCAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAAt                                                                                  <  8:274935/147‑1 (MQ=255)
                                                                     cgcCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGgcgc                                                                            <  3:122950/149‑1 (MQ=255)
                                                                      gcCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCAGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCt                                                                           >  8:272798/1‑149 (MQ=255)
                                                                         aGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTtgt                                                                        <  7:39558/149‑1 (MQ=255)
                                                                           tAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGa                                                                      >  2:54380/1‑149 (MQ=255)
                                                                           tAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGa                                                                      <  1:149778/149‑1 (MQ=255)
                                                                              agaACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAg                                                                   >  5:81704/1‑149 (MQ=255)
                                                                                 aCCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCta                                                                >  3:165250/1‑149 (MQ=255)
                                                                                   cGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATa                                                                    <  2:54965/143‑1 (MQ=255)
                                                                                   cGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATa                                                                    >  1:54965/1‑143 (MQ=255)
                                                                                        ccGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCt                                                         <  1:54380/149‑1 (MQ=255)
                                                                                          gATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCt                                                                                                                                         >  3:53111/1‑67 (MQ=255)
                                                                                          gATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCt                                                                                                                                         <  4:53111/67‑1 (MQ=255)
                                                                                          gATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCg                                                       >  2:87061/1‑149 (MQ=255)
                                                                                             cACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGt                                                                                                             >  5:302131/1‑92 (MQ=255)
                                                                                             cACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGt                                                                                                             <  6:302131/92‑1 (MQ=255)
                                                                                               cTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTTTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCAtt                                                  <  7:272798/149‑1 (MQ=255)
                                                                                                   gCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGCCTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCATTGGAt                                              >  7:312307/1‑149 (MQ=255)
                                                                                                    cccGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATg                                                                                 >  1:148700/1‑113 (MQ=255)
                                                                                                    cccGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATg                                                                                 <  2:148700/113‑1 (MQ=255)
                                                                                                     ccGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCATTGGATg                                             <  3:170398/148‑1 (MQ=255)
                                                                                                               gacagaAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCATTGGATGCCCTCGATggg                                  >  6:329454/1‑149 (MQ=255)
                                                                                                                     aGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACCCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCATTGGATGCCCTCGATGGGTCCgg                             <  7:179510/148‑1 (MQ=255)
                                                                                                                          aaGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCATTGGATGCCCTCGATGGGTCCGGTgg                          <  8:227631/146‑1 (MQ=255)
                                                                                                                          aaGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCATTGGATGCCCTCGATGGGTCCGGTgg                          >  7:227631/1‑146 (MQ=255)
                                                                                                                           aGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAAt                                                                                  >  6:184343/1‑89 (MQ=255)
                                                                                                                           aGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAAt                                                                                  <  5:184343/89‑1 (MQ=255)
                                                                                                                              tgacagaaGCCCAGCCCCGCGCTGGGCTTTTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCATTGGATGCCCTCGATGGGTCCGGTGGGCCTTGt                   >  5:264881/7‑149 (MQ=255)
                                                                                                                                acagaaGCCCAGCCCCGCGCTGGGCTTTTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCATTGGATGCCCTCGATGGGTCCGGTGGGCCTTGttt                 <  6:264881/145‑1 (MQ=255)
                                                                                                                                 aaaaaGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCATTGGATGCCCTCGATGGGTCCGGTGGGCCTTGTTTc                >  2:279873/1‑149 (MQ=255)
                                                                                                                                 aaaaaGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCATTGGATGCCCTCGATGGGTCCGGTGGGCCTTGTTTc                >  6:236797/1‑149 (MQ=255)
                                                                                                                                    aaGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCATTGGATGCCCTCGATGGGTCCGGTGGGCCTTGTTTCTgc             <  3:154796/149‑1 (MQ=255)
                                                                                                                                    aaGCCCAGCCCCGCGCTGGGCTTTTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAg                                                                   <  4:257850/95‑1 (MQ=39)
                                                                                                                                    aaGCCCAGCCCCGCGCTGGGCTTTTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAg                                                                   >  3:257850/1‑95 (MQ=39)
                                                                                                                                      gCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTgg                                                                                               >  5:68/1‑65 (MQ=255)
                                                                                                                                      gCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTgg                                                                                               <  6:68/65‑1 (MQ=255)
                                                                                                                                       cccAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCATTGGATGCCCTCGATGGGTCCGGTGGGCCTTGTTTCTgcgcg          >  4:120805/1‑149 (MQ=255)
                                                                                                                                           gcgcGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCATTGGATGCCCTCGATGGGTCCGGTGGGCCTTGTTTCTGCGCGCCTc      >  7:65929/1‑149 (MQ=255)
                                                                                                                                               ggggCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCATTGGATGCCCTCGATGGGTCCGGTGGGCCTTGTTTCTGCGCGCCTCGaaa  <  8:65929/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
ATTTCGCTTAAATGCGACCGGTGGAAAATGATCGGCCCTAATCCTTGCTGATAGATATCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCATTGGATGCCCTCGATGGGTCCGGTGGGCCTTGTTTCTGCGCGCCTCGAAA  >  NC_002947/4424998‑4425289

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