Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,655,517 T→C 6.2% I418V (ATC→GTC)  PP_1449 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,655,5170TC 6.2% 186.4 / 2.3 64I418V (ATC→GTC) PP_1449hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (27/33);  new base C (2/2);  total (29/35)
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.77e-01

GCCTTTATCGACATGTGCCCGGCGCTCAGGGTACTGCCCTTTGCGGTTTCCTTGAGCTGGTCGCTATCGCTGTCACCGCGCCATAAATC‑CTTGCGATGACGCACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTCCAGGCCGGCACCGCTGTCGACCGTCAGCTCGTTACCGGCGTCCACCTTCGCGGCGGTCATGCGCATGTCGTCGCCGGACTGCAACAGCACGCTCTCAGTGCCCAGCAATTGACTGCCCTGTTGTTGCTTGTGGGTGGTG  >  NC_002947/1655375‑1655659
                                                                                                                                               |                                                                                                                                              
gCCTTTATCGACATGTGCCCGGCGCTCAGGGTACTGCCCTTTGCGGTTTCCTTGAGCTGGTCGCTATCGCTGTCACCGCGCCATAAATC‑CTTGCGATGACGCACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTcc                                                                                                                                            <  8:316966/147‑1 (MQ=255)
  cTTTATCGACATGTGCCCGGCGCTCAGGGTACTGCCCTTTGCGGTTTCCTTGAGCTGGTCGCTATCGCTGTCACCGCGCCATAAATC‑CTTGCGATGACGCACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTCCAGGc                                                                                                                                        >  7:137367/1‑149 (MQ=255)
  cTTTATCGACATGTGCCCGGCGCTCAGGGTACTGCCCTTTGCGGTTTCCTTGAGCTGGTCGCTATCGCTGTCACCGCGCCATAAATC‑CTTGCGATGACGCACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTCCAGGc                                                                                                                                        >  6:354502/1‑149 (MQ=255)
  cTTTATCGACATGTGCCCGGCGCTCAGGGTACTGCCCTTTGCGGTTTCCTTGAGCTGGTCGCTATCGCTGTCACCGCGCCATAAATC‑CTTGCGATGACGCACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTCCAGGc                                                                                                                                        >  6:228850/1‑149 (MQ=255)
   tttATCGACATGTGCCCGGCGCTCAGGGTACTGCCCTTTGCGGTTTCCTTGAGCTGGTCGCTATCGCTGTCACCGCGCCATAAATC‑CTTGCGATGACGCACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTCCAGGcc                                                                                                                                       >  7:180302/1‑149 (MQ=255)
   tttATCGACATGTGCCCGGCGCTCAGGGTACTGCCCTTTGCGGTTTCCTTGAGCTGGTCGCTATCGCTGTCACCGCGCCATAAATC‑CTTGCGATGACGCACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTCCAGGcc                                                                                                                                       >  7:349325/1‑149 (MQ=255)
          aCATGTGCCCGGCGCTCAGGGTACTGCCCTTTGCGGTTTCCTTGAGCTGGTCGCTATCGTTGTCACCGCGCCATAAATC‑CTTGCGATGACGCACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTCCAGGCCGGCACCg                                                                                                                                <  7:59472/149‑1 (MQ=255)
            aTGTGCCCGGCGCTCAGGGTACTGCCCTTTGCGGTTTCCTTGAGCTGGTCGCTATCGCTGTCACCGCGCCATAAATC‑CTTGCGATGACGCACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCCTGGACATCCTCGCCGGc                                                                                                                                    <  1:227171/143‑1 (MQ=255)
            aTGTGCCCGGCGCTCAGGGTACTGCCCTTTGCGGTTTCCTTGAGCTGGTCGCTATCGCTGTCACCGCGCCATAAATC‑CTTGCGATGACGCACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCCTGGACATCCTCGCCGGc                                                                                                                                    >  2:227171/1‑143 (MQ=255)
                gCCCGGCGCTCAGGGTACTGCCCTTTGCGGTTTCCTTGAGCTGGTCGCTATCGCTGTCACCGCGCCATAAATC‑CTTGCGATGACGCACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTCCAGGCCGGCACCGCTGTc                                                                                                                           >  4:322799/1‑148 (MQ=255)
                 cccGGCGCTCAGGGTACTGCCCTTTGCGGTTTCCTTGAGCTGGTCGCTATCGCTGTCACCGCGCCATAAATC‑CTTGCGATGACGCACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTCCAGGCCGGCACCGCTGTCGa                                                                                                                         <  3:117242/149‑1 (MQ=255)
                  ccGGCGCTCAGGGTACTGCCCTTTGCGGTTTCCTTGAGCTGGTCGCTATCGCTGTCACCGCGCCATAAATC‑CTTGCGATGACGCACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTCCAGGCCGGCACCGCTGTCGAc                                                                                                                        >  4:354048/1‑149 (MQ=255)
                   cGGCGCTCAGGGTACTGCCCTTTGCGGTTTCCTTGAGCTGGTCGCTATCGCTGTCACCGCGCCATAAATC‑CTTGCGATGACGCACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTCCAGGCCGGCACCGCTGTCGAcc                                                                                                                       <  5:228850/149‑1 (MQ=255)
                      cgcTCAGGGTACTGCCCTTTGCGGTTTCCTTGAGCTGGGCGCTATCGCTGTCACCGCGCCATACATC‑CTTGCGATGACGCACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTCCAGGCCGGCACCGCTGTCGACCGTc                                                                                                                    <  4:324106/149‑1 (MQ=255)
                        cTCAGGGTACTGCCCTTTGCGGTTTCCTTGAGCTGGTCGCTATCGCTGTCACCGCGCCATAAATC‑CTTGCGATGACGCACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTCCAGGCCGGCACCGCTGTCGACCGTCAg                                                                                                                  >  5:243934/1‑149 (MQ=255)
                           aGGGTACTGCCCTTTGCGGTTTCCTTGAGCTGGTCGCTATCGCTGTCACCGCGCCATAAATC‑CTTGCGATGACGCACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTCCAGGCCGGCACCGCTGTCGACCGTCAGCt                                                                                                                <  8:180302/148‑1 (MQ=255)
                           aGGGTACTGCCCTTTGCGGTTTCCTTGAGCTGGTCGCTATCGCTGTCACCGCGCCATAAATC‑CTTGCGATGACGCACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTCCAGGCCGGCACCGCTGTCGACCGTCAGCTc                                                                                                               <  3:322799/149‑1 (MQ=255)
                           aGGGTACTGCCCTTTGCGGTTTCCTTGAGCTGGTCGCTATCGCTGTCACCGCGCCATAAATC‑CTTGCGATGACGCACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTCCAGGCCGGCACCGCTGTCGACCGTCAGCTc                                                                                                               <  5:277381/149‑1 (MQ=255)
                           aGGGTACTGCCCTTTGCGGTTTCCTTGAGCTGGTCGCTATCGCTGTCACCGCGCCATAAATC‑CTTGCGATGACGCACCTGCGGTTCGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTCCAGCCCGGCACCGCTGTCGACCGTCAGCTc                                                                                                               <  4:291008/149‑1 (MQ=255)
                                 cTGCCCTTTGCGGTTTCCTTGAGCTGGTCGCTATCGCTGTCACCGCGCCATAAATC‑CTTGCGATGACGCACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTCCAGGCCGGCACCGCTGTCGACCGTCAGCTCGTTAc                                                                                                          <  5:354502/148‑1 (MQ=255)
                                       tttGCGGTTTCCTTGAGCTGGTCGCTATCGCTGTCACCGCGCCATAAATC‑CTTGCGATGACGCACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTCCAGGCCGGCACCGCTGTCGACCGTCAGCTCGtt                                                                                                            <  6:245960/140‑1 (MQ=255)
                                       tttGCGGTTTCCTTGAGCTGGTCGCTATCGCTGTCACCGCGCCATAAATC‑CTTGCGATGACGCACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTCCAGGCCGGCACCGCTGTCGACCGTCAGCTCGtt                                                                                                            >  5:245960/1‑140 (MQ=255)
                                       tttGCGGTTTCCTTGAGCTGGTCGCTATCGCTGTCACCGCGCCATAAATC‑CTTGCGATGACGCACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCCTGGACATCCTCGCCGGCATCgaca                                                                                                                             <  7:344062/123‑4 (MQ=38)
                                       tttGCGGTTTCCTTGAGCTGGTCGCTATCGCTGTCACCGCGCCATAAATC‑CTTGCGATGACGCACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCCTGGACATCCTCGCCGGCATCgaca                                                                                                                             >  8:344062/1‑120 (MQ=38)
                                            ggTTTCCTTGAGCTGGTCGCTATCGCTGTCACCGCGCCATAAATC‑CTTGCGATGACGCACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTCCAGGCCGGCACCGCTGTCGACCGTCAGCTCGTTACCGGCGTCCACCt                                                                                              >  2:40278/1‑149 (MQ=255)
                                              tttCCTTGAGCTGGTCGCTATCGCTGTCACCGCGCCATAAATC‑CTTGCGATGACGCACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTCCAGGCCGGCACCGCTGTCGACCGTCAGCTCGTTACCGGCGTCCACCTTc                                                                                            >  3:169326/1‑149 (MQ=255)
                                              tttCCTTGAGCTGGTCGCTATCGCTGTCACCGCGCCATAAATC‑CTTGCGATGACGCACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGACGAGGATGTCCAGGCCGGCACCGCTGTCGACCGTCAGCTCGTTACCGGCGTCCACCTTc                                                                                            >  8:137849/1‑149 (MQ=255)
                                                  attGAGCTGGTCGCTATCGCTGTCACCGCGCCATAAATC‑CTTGCGATGACGCACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTCCAGGCCGGCACCGCTGTCGACCGTCAGCTCGTTACCGGCGTCCCCCTTCgcgg                                                                                        >  4:252368/2‑149 (MQ=255)
                                                   ttGAGCTGGTCGCTATCGCTGTCACCGCGCCATAAATC‑CTTGCGATGACGCACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTCCAGGCCGGCACCGCTGTCGACCGTCAGCTCGTTACCGGCGTCCACCTTCgcggc                                                                                       >  4:293891/1‑149 (MQ=255)
                                                   ttGAGCTGGTCGCTATCGCTGTCACCGCGCCATAAATC‑CTTGCGATGACGCACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTCCAGGCCGGCACCGCTGTCGACCGTCAGCTCGTTACCGGCGTCCACCTTCgcggc                                                                                       >  2:224805/1‑149 (MQ=255)
                                                           gTCGCTATCGCTGTCACCGCGCCATAAATC‑CTTGCGATGACGCACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTCCAGGCCGGCACCGCTGTCGACCGTCAGCTCGTTACCGGCGTCCACCTTCGCGGCGGTCATgc                                                                               <  1:224805/149‑1 (MQ=255)
                                                              gCTATCGCTGGCACCGCGCCATAATTC‑CTTGCGATGACGCTCGTGCTGTTTGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTCCAGGCCGGCACCGCTGTCGACCGTCAGCTCGTTACCGGCGTCCACCTTCGCGGCGGTCATgcgc                                                                             <  5:152533/148‑1 (MQ=255)
                                                               cTATCGCTGTCACCGCGCCATAAATCGCTTGCGATGACGCACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTCCAGGCCGGCACCGCTGTCGACCGTCAGCTCGTTACCGGCGTCCACCTTCGCGGCGGTCATGCGCa                                                                            >  1:192264/1‑149 (MQ=255)
                                                               cTATCGCTGTCACCGCGCCATAAATC‑CTTGCGATGACGCACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTCCAGGCCGGCACCGCTGTCGACCGTCAGCTCGTTACCGGCGTCCACCTTCGCGGCGGTCATGCGCAt                                                                           <  5:48062/149‑1 (MQ=255)
                                                                  tCGCTGTCACCGCGCCATAAATC‑CTTGCGATGACGCACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTCCAGGCCGGCACCGCTGTCGACCGTCAGCTCGTTACCGGCGTCCACCTTCGCGGCGGTCATGCGCATg                                                                          <  3:293891/147‑1 (MQ=255)
                                                                     cTGTCACCGCGCCATAAATC‑CTTGCGATGACGCACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTCCAGGCCGGCACCGCTGTCGACCGTCAGCTCGTTACCGGCGTCCACCTTCGCGGCGGTCATGCGCATgtcgtc                                                                     <  5:153512/149‑1 (MQ=255)
                                                                     cTGTCACCGCGCCATAAATC‑CTTGCGATGACGCACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTCCAGGCCGGCACCGCTGTCGACCGTCAGCTCGTTACCGGCGTCCACCTTCGCGGCGGTCATGCGCATgtcgtc                                                                     >  4:161201/1‑149 (MQ=255)
                                                                      tGTCACCGCGCCATAAATC‑CTTGCGATGACGCACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTCCAGGCCGGCACCGCTGTCGACCGTCAGCTCGTTACCGGCGTCCACCTTCGCGGCGGTCATGCGCATGtcgtcg                                                                    >  6:220458/1‑149 (MQ=255)
                                                                          aCCGCGCCATAAATC‑CTTGCGATGACGCACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTCCAGGCCGGCACCGCTGTCGACCGTCAGCTCGTTACCGGCGTCCACCTTCGCGGCGGTCATGCGCATGTCGTCGCCg                                                                 >  6:311776/1‑148 (MQ=255)
                                                                          aCCGCGCCATAAATC‑CTTGCGATGACGCACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTCCAGGCCGGCACCGCTGTCGACCGTCAGCTCGTTACCGGCGTCCACCTTCGCGGCGGTCATGCGCATGTCGTCGCCg                                                                 <  8:137367/148‑1 (MQ=255)
                                                                           ccGCGCCATAAATC‑CTTGCGATGACGCACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTCCAGGCCGGCACCGCTGTCGACCGTCAGCTCGTTACCGGCGTCCACCTTCGCGGCGGTCATGCGCATGTCGTCGCCGGa                                                               <  3:46938/149‑1 (MQ=255)
                                                                           ccGCGCCATAAATC‑CTTGCGATGACGCACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTCCAGGCCGGCACCGCTGTCGACCGTCAGCTCGTTACCGGCGTCCACCTTCGCGGCGGTCATGCGCATGTCGTCGCCGGa                                                               <  6:243934/149‑1 (MQ=255)
                                                                            cgcgCCATAAATC‑CTTGCGATGACGCACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTCCAGGCCGGCACCGCTGTCGACCGTCAGCTCGTTACCGGCGTCCACCTTCGCGGCGGTCATGCGCATGTCGTCGCCGGAc                                                              >  6:258667/1‑149 (MQ=255)
                                                                                 cATAAATC‑CTTGCGATGACGCACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTCCAGGCCGGCACCGCTGTCGACCGTCAGCTCGTTACCGGCGTCCACCTTCGCGGCGGTCATGCGCATGTCGTCGCCGGACTGCaa                                                         <  3:354048/149‑1 (MQ=255)
                                                                                   tAAATC‑CTTGCGATGACGCACCTGCTGTGCGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTCCAGGCCGGCACCGCTGTCGACCGTCAGCTCGTTACCGGTGTCCACCTTCGCGGCGGTCATGCGCATGGCGTCGCCGGACTGCAACa                                                       <  1:267468/149‑1 (MQ=255)
                                                                                     aaTC‑CTTGCGATGACGCACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTCCAGGCCGGCACCGCTGTCGACCGTCAGCTCGTTACCGGCGTCCACCTTCGCGGCGGTCATGCGCATGTCGTCGCCGGACTGCAACAGc                                                     <  3:161201/149‑1 (MQ=255)
                                                                                              cGATGACGCACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTCCAGGCCGGCACCGCTGTCGACCGTCAGCTCGTTACCGGCGTCCACCTTCGCGGCGGTCATGCGCATGTCGTCGCCGGACTGCAACAGCACGctct                                              <  5:311776/148‑1 (MQ=255)
                                                                                              cGATGACGCACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTCCAGGCCGGCACCGCTGTCGACCGTCAGCTCGTTACCGGCGTCCACCTTCGCGGCGGTCATGCGCATGTCGTCGCCGGACTGCAACAGCACGCtctc                                             <  8:349325/149‑1 (MQ=255)
                                                                                               gATGACGCACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTCCAGGCCGGCACCGCTGTCGACCGTCAGCTCGTTACCGGCGTCCACCTTCGCGGCGGTCATGCGCATGTCGTCGCCGGACTGCAACAGCACGCTCTCa                                            <  7:204204/149‑1 (MQ=255)
                                                                                                      cACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTCCAGGCCGGCACCGCTGTCGACCGTCAGCTCGTTACCGGCGTCCACCTTCGCGGCGGTCATGCGCATGTCGTCGCCGGACTGCAACAGCACGCTCTCAGTGCCCa                                     <  2:81645/149‑1 (MQ=255)
                                                                                                        ccTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTCCAGGCCGGCACCGCTGTCGACCGTCAGCTCGTTACCGGCGTCCACCTTCGCGGCGGTCATGCGCATGTCGTCGCCGGACTGCAACAGCACGCTCTCAGTGCCCAg                                    >  3:242865/1‑148 (MQ=255)
                                                                                                           gctgTTCGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTCCAGGCCGGCACCGCTGTCGACCGTCAGCTCGTTACCGGCGTCCACCTTCGCGGCGGTCATGCGCATGTCGTCGCCGGACTGCAACAGCACGCTCTCAGTGCCCAGCAAt                                >  2:34999/1‑149 (MQ=255)
                                                                                                           gctgTTCGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTCCAGGCCGGCACCGCTGTCGACCGTCAGCTCGTTACCGGCGTCCACCTTCGCGGCGGTCATGCGCATGTCGTCGCCGGACTGCAACAGCACGCTCTCAGTGCCCAGCAAt                                >  6:374428/1‑149 (MQ=255)
                                                                                                                      cGTTTGCTGTCGATGCCGACGAGGATGTCCAGGCCGGCACCGCTGTCGACCGTCAGCTCGTTACCGGCGTCCACCTTCGCGGCGGTCATGCGCATGTCGTCGCCGGACTGCAACAGCACGCTCTCAGTGCCCAGCAATTGACTGCCCtg                     <  7:137849/149‑1 (MQ=255)
                                                                                                                            cTGTCGATGCCGGCGAGGATGTCCAGGCCGGCACCGCTGTCGACCGTCAGCTCGTTACCGGCGTCCACCTTCGCGGCTGTCATGCGCATGTCGTCGCCGGACTGCAACAGCACGCTCTCAGTGCCCAGCAATTGACTGCCCTGttgttg               <  5:258667/149‑1 (MQ=255)
                                                                                                                            cTGTCGATGCCGGCGAGGATGTCCAGGCCGGCACCGCTGTCGACCGTCAGCTCGTTACCGGCGTCCACCTTCGCGGCGGTCATGCGCATGTCGTCGCCGGACTGCAACAGCACGCTCTCAGTGCCCAGCAATTGACTGCCCTGttgttg               <  1:182360/149‑1 (MQ=255)
                                                                                                                             tGTCGATGCCGGCGAGGATGTCCAGGCCGGCACCGCTGTCGACCGTCAGCTCGTTACCGGCGTCCACCTTCGCGGCGGTCATGCGCATGTCGTCGCCGGACTGCAACAGCACGCTCTCAGTGCCCAGCAATTGACTGCCCTGTTGTTGc              >  4:331939/1‑149 (MQ=255)
                                                                                                                                 gATGCCGGCGAGGATGTCCAGGCCGGCACCGCTGTCGACCGTCAGCTCGTTACCGGCGTCCACCTTCGCGGCGGTCATGCGCATGTCGTCGCCGGACTGCAACAGCACGCTCTCAGTGCCCAGCAATTGACTGCCCTGTTGTTGCTtgt          <  4:169326/149‑1 (MQ=255)
                                                                                                                                  aTGCCGGCGAGGATGTCCAGGCCGGCACCGCTGTCGACCGTCAGCTCGTTACCGGCGTCCACCTTCGCGGCGGTCATGCGCATGTCGTCGCCGGACTGCAACAGCACGCTCTCAGTGCCCAGCAATTGACTGCCCTGTTGTTGCTtgtg         <  4:242865/149‑1 (MQ=255)
                                                                                                                                  aTGCCGGCGAGGATGTCCAGGCCGGCACCGCTGTCGACCGTCAGCTCGTTACCGGCGTCCACCTTCGCGGCGGTCATGCGCATGTCGTCGCCGGACTGCAACAGCACGCTCTCAGTGCCCAGCAATTGACTGCCCTGTTGTTGCTtgt          <  7:57594/148‑1 (MQ=255)
                                                                                                                                  aTGCCGGCGAGGATGTCCAGGCCGGCACCGCTGTCGACCGTCAGCTCGTTACCGGCGTCCACCTTCGCGGCGGTCATGCGCATGTCGTCGCCGGACTGCAACAGCACGCTCTCAGTGCCCAGCAATTGACTGCCCTGTTGTTGCTtgt          >  8:57594/1‑148 (MQ=255)
                                                                                                                                      cGGCGAGGATGTCCAGGCCGGCACCGCTGTCGACCGTCAGCTCGTTACCGGCGTCCACCTTCGCGGCGGTCATGCGCATGTCGTCGCCGGACTGCAACAGCACGCTCTCAGTGCCCAGCAATTGACTGCCCTGTTGTTGCTTgtgggtg     >  6:351647/1‑149 (MQ=255)
                                                                                                                                       ggCGAGGATGTCCAGGCCGGCACCGCTGTCGACCGTCAGCTCGTTACCGGCGTCCACCTTCGCGGCGGTCATGCGCATGTCGTCGCCGGACTGCAACAGCACGCTCTCAGTGCCCAGCAATTGACTGCCCTGTTGTTGCTTgtgggtgg    >  2:315235/1‑149 (MQ=255)
                                                                                                                                         cGAGGATGTCCAGGCCGGCACCGCTGTCGACCGTCAGCTCGTTACCGGCGTCCACCTTCGCGGCGGTCATGCGCATGTCGTCGCCGGACTGCAACAGCACGCTCTCAGTGCCCAGCAATTGACTGCCCTGTTGTTGCTTGTGggtggtg  <  5:374428/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                              
GCCTTTATCGACATGTGCCCGGCGCTCAGGGTACTGCCCTTTGCGGTTTCCTTGAGCTGGTCGCTATCGCTGTCACCGCGCCATAAATC‑CTTGCGATGACGCACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTCCAGGCCGGCACCGCTGTCGACCGTCAGCTCGTTACCGGCGTCCACCTTCGCGGCGGTCATGCGCATGTCGTCGCCGGACTGCAACAGCACGCTCTCAGTGCCCAGCAATTGACTGCCCTGTTGTTGCTTGTGGGTGGTG  >  NC_002947/1655375‑1655659

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: ‑