Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,655,513 A→T 6.2% L419H (CTC→CAC) ‡ PP_1449 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,655,5130AT 6.2% 198.8 / 2.3 64L419H (CTC→CAC) ‡PP_1449hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (27/33);  new base T (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 = 8.75e-01

GTTGTCGCTTGCCTTTATCGACATGTGCCCGGCGCTCAGGGTACTGCCCTTTGCGGTTTCCTTGAGCTGGTCGCTATCGCTGTCACCGCGCCATAAATC‑CTTGCGATGACGCACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTCCAGGCCGGCACCGCTGTCGACCGTCAGCTCGTTACCGGCGTCCACCTTCGCGGCGGTCATGCGCATGTCGTCGCCGGACTGCAACAGCACGCTCTCAGTGCCCAGCAATTGACTGCCCTGTTGTTGCTTGTGGGTGGTG  >  NC_002947/1655365‑1655659
                                                                                                                                                     |                                                                                                                                                  
gTTGTCGCTTGCCTTTATCGACATGTGCCCGGCGCTCAGGGTACTGCCCTTTGCGGTTTCCTTGAGCTGGTCGCTATCGCTGTCACCGCGCCATAAATC‑CTTGCGATGACGCACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCGa                                                                                                                                                    <  7:171247/149‑1 (MQ=255)
          gCCTTTATCGACATGTGCCCGGCGCTCAGGGTACTGCCCTTTGCGGTTTCCTTGAGCTGGTCGCTATCGCTGTCACCGCGCCATAAATC‑CTTGCGATGACGCACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTcc                                                                                                                                            <  8:316966/147‑1 (MQ=255)
            cTTTATCGACATGTGCCCGGCGCTCAGGGTACTGCCCTTTGCGGTTTCCTTGAGCTGGTCGCTATCGCTGTCACCGCGCCATAAATC‑CTTGCGATGACGCACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTCCAGGc                                                                                                                                        >  6:228850/1‑149 (MQ=255)
            cTTTATCGACATGTGCCCGGCGCTCAGGGTACTGCCCTTTGCGGTTTCCTTGAGCTGGTCGCTATCGCTGTCACCGCGCCATAAATC‑CTTGCGATGACGCACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTCCAGGc                                                                                                                                        >  6:354502/1‑149 (MQ=255)
            cTTTATCGACATGTGCCCGGCGCTCAGGGTACTGCCCTTTGCGGTTTCCTTGAGCTGGTCGCTATCGCTGTCACCGCGCCATAAATC‑CTTGCGATGACGCACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTCCAGGc                                                                                                                                        >  7:137367/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                                                                                                            >  5:245960/1‑140 (MQ=255)
                                                 tttGCGGTTTCCTTGAGCTGGTCGCTATCGCTGTCACCGCGCCATAAATC‑CTTGCGATGACGCACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTCCAGGCCGGCACCGCTGTCGACCGTCAGCTCGtt                                                                                                            <  6:245960/140‑1 (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                                                                 <  8:137367/148‑1 (MQ=255)
                                                                                    aCCGCGCCATAAATC‑CTTGCGATGACGCACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTCCAGGCCGGCACCGCTGTCGACCGTCAGCTCGTTACCGGCGTCCACCTTCGCGGCGGTCATGCGCATGTCGTCGCCg                                                                 >  6:311776/1‑148 (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                                >  6:374428/1‑149 (MQ=255)
                                                                                                                     gctgTTCGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTCCAGGCCGGCACCGCTGTCGACCGTCAGCTCGTTACCGGCGTCCACCTTCGCGGCGGTCATGCGCATGTCGTCGCCGGACTGCAACAGCACGCTCTCAGTGCCCAGCAAt                                >  2:34999/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)
                                                                                                                                                     |                                                                                                                                                  
GTTGTCGCTTGCCTTTATCGACATGTGCCCGGCGCTCAGGGTACTGCCCTTTGCGGTTTCCTTGAGCTGGTCGCTATCGCTGTCACCGCGCCATAAATC‑CTTGCGATGACGCACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTCCAGGCCGGCACCGCTGTCGACCGTCAGCTCGTTACCGGCGTCCACCTTCGCGGCGGTCATGCGCATGTCGTCGCCGGACTGCAACAGCACGCTCTCAGTGCCCAGCAATTGACTGCCCTGTTGTTGCTTGTGGGTGGTG  >  NC_002947/1655365‑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: ‑