Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,655,522 A→T 6.1% L416Q (CTG→CAG) ‡ PP_1449 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,655,5220AT 6.1% 199.0 / 2.6 66L416Q (CTG→CAG) ‡PP_1449hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (29/33);  new base T (2/2);  total (31/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 = 9.09e-01

CTTTATCGACATGTGCCCGGCGCTCAGGGTACTGCCCTTTGCGGTTTCCTTGAGCTGGTCGCTATCGCTGTCACCGCGCCATAAATC‑CTTGCGATGACGCACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTCCAGGCCGGCACCGCTGTCGACCGTCAGCTCGTTACCGGCGTCCACCTTCGCGGCGGTCATGCGCATGTCGTCGCCGGACTGCAACAGCACGCTCTCAGTGCCCAGCAATTGACTGCCCTGTTGTTGCTTGTGGGTGGTGGTGGTCTTGCG  >  NC_002947/1655377‑1655670
                                                                                                                                                  |                                                                                                                                                    
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                                                                                                                          <  5:277381/149‑1 (MQ=255)
                         aGGGTACTGCCCTTTGCGGTTTCCTTGAGCTGGTCGCTATCGCTGTCACCGCGCCATAAATC‑CTTGCGATGACGCACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTCCAGGCCGGCACCGCTGTCGACCGTCAGCTc                                                                                                                          <  3:322799/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                                                                                                                                        >  8:344062/1‑120 (MQ=38)
                                     tttGCGGTTTCCTTGAGCTGGTCGCTATCGCTGTCACCGCGCCATAAATC‑CTTGCGATGACGCACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCCTGGACATCCTCGCCGGCATCgaca                                                                                                                                        <  7:344062/123‑4 (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                                                                                                  >  2:224805/1‑149 (MQ=255)
                                                 ttGAGCTGGTCGCTATCGCTGTCACCGCGCCATAAATC‑CTTGCGATGACGCACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTCCAGGCCGGCACCGCTGTCGACCGTCAGCTCGTTACCGGCGTCCACCTTCgcggc                                                                                                  >  4:293891/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                                                                          <  6:243934/149‑1 (MQ=255)
                                                                         ccGCGCCATAAATC‑CTTGCGATGACGCACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTCCAGGCCGGCACCGCTGTCGACCGTCAGCTCGTTACCGGCGTCCACCTTCGCGGCGGTCATGCGCATGTCGTCGCCGGa                                                                          <  3:46938/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)
                                                                                                                                                 cAGGCCGGCACCGCTGTCGACCGTCAGCTCGTTACCGGCGTCCACCTTCGCGGCGGTCATGCGCATGTCGTCGCCGGACTGCAACAGCACGCTCTCAGTGCCCAGCAATTGACTGCCCtg                                <  3:147830/120‑1 (MQ=255)
                                                                                                                                                 cAGGCCGGCACCGCTGTCGACCGTCAGCTCGTTACCGGCGTCCACCTTCGCGGCGGTCATGCGCATGTCGTCGCCGGACTGCAACAGCACGCTCTCAGTGCCCAGCAATTGACTGCCCtg                                >  4:147830/1‑120 (MQ=255)
                                                                                                                                                 cAGGCCGGCACCGCTGTCGACCGTCAGCTCGTTACCGGCGTCCACCTTCGCGGCGGTCATGCGCATGTCGTCGCCGGACTGCAACAGCACGCTCTCAGTGCCCAGCAATTGACTGCCCTGTTGTTGCTTGTGGGTGGTGGTGGTCTTGc   >  5:340856/1‑149 (MQ=255)
                                                                                                                                                  aGGCCGGCACCGCTGTCGACCGTCAGCTCGTTACCGGCGTCCACCTTCGCGGCGGTCATGCGCATGTCGTCGCCGGACTGCAACAGCACGCTCTCAGTGCCCAGCAATTGACTGCCCTGTTGTTGCTTGTGGGTGGTGGTGGTCTTGCg  >  3:175002/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
CTTTATCGACATGTGCCCGGCGCTCAGGGTACTGCCCTTTGCGGTTTCCTTGAGCTGGTCGCTATCGCTGTCACCGCGCCATAAATC‑CTTGCGATGACGCACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTCCAGGCCGGCACCGCTGTCGACCGTCAGCTCGTTACCGGCGTCCACCTTCGCGGCGGTCATGCGCATGTCGTCGCCGGACTGCAACAGCACGCTCTCAGTGCCCAGCAATTGACTGCCCTGTTGTTGCTTGTGGGTGGTGGTGGTCTTGCG  >  NC_002947/1655377‑1655670

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