Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,938,839 C→T 8.7% E73K (GAG→AAG)  PP_1738 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,938,8390CT 8.7% 222.1 / 8.0 81E73K (GAG→AAG) PP_1738hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (39/33);  new base T (3/4);  total (43/38)
Fisher's exact test for biased strand distribution p-value = 7.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.64e-01

ACAGCTTTTTCAGGTATTCATTGATGCTGCTGCTGTAGGGCGCGCCGGTGCCGGCGATGATGAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACGGGTTGGTAAGAGGGTAGCCATAGCTGGCAGCGTCGATATC  >  NC_002947/1938692‑1938986
                                                                                                                                                   |                                                                                                                                                   
aCAGCTTTTTCAGGTATTCATTGATGCTGCTGCTGTAGGGCGCGCCGGTGCCGGCGATGATGAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTcg                                                                                                                                                    >  5:378142/1‑149 (MQ=255)
aCAGCTTTTTCAGGTATTCATTGATGCTGCTGCTGTAGGGCGCGCCGGTGCCGGCGATGATGAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTc                                                                                                                                                     <  1:379525/148‑1 (MQ=255)
aCAGCTTTTTCAGGTATTCATTGATGCTGCTGCTGTAGGGCGCGCCGGTGCCGGCGATGATGAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTc                                                                                                                                                     <  6:238982/148‑1 (MQ=255)
 cAGCTTTTTCAGGTATTCATTGATGCTGCTGCTGTAGGGCGCGCCGGTGCCGGCGATGATGAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTcgc                                                                                                                                                   <  4:26483/149‑1 (MQ=255)
   gCTTTTTCAGGTATTCATTGATGCTGCTGCTGTAGGGCGCGCCGGTGCCGGCGATGATGAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCgcgc                                                                                                                                                 <  2:13717/149‑1 (MQ=255)
   gCTTTTTCAGGTATTCATTGATGCTGCTGCTGTAGGGCGCGCCGGTGCCGGCGATGATGAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCgcgc                                                                                                                                                 <  2:378129/149‑1 (MQ=255)
   gCTTTTTCAGGTATTCATTGATGCTGCTGCTGTAGGGCGCGCCGGTGCCGGCGATGATGAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCgcgc                                                                                                                                                 >  6:326995/1‑149 (MQ=255)
     tttttCAGGTATTCATTGATGCTGCTGCTGTAGGGCGCGCCGGTGCCGGCGATGATGAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTc                                                                                                                                               >  7:365017/1‑149 (MQ=255)
         tCAGGTATTCATTGATGCTGCTGCTGTAGGGCGCGCCGGTGCCGGCGATGATGAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCgg                                                                                                                                            <  5:90421/148‑1 (MQ=255)
            ggTATTCATTGATGCTGCTGCTGTAGGGCGCGCCGGTGCCGGCGATGATGAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGgcg                                                                                                                                         >  5:321022/1‑148 (MQ=255)
                ttCATTGATGCTGCTGCTGTAGGGCGCGCCGGTGCCGGCGATGATGAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGt                                                                                                                                    >  6:95786/1‑149 (MQ=255)
                  cATTGATGCTGCTGCTGTAGGGCGCGCCGGTGCCGGCGATGATGAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGtt                                                                                                                                   >  6:169364/1‑148 (MQ=255)
                   aTTGATGCTGCTGCTGTAGGGCGCGCCGGTGCCGGCGATGATGAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGtt                                                                                                                                   >  3:25269/1‑147 (MQ=255)
                   aTTGATGCTGCTGCTGTAGGGCGCGCCGGTGCCGGCGATGATGAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCa                                                                                                                                 <  1:313704/149‑1 (MQ=255)
                   aTTGATGCTGCTGCTGTAGGGCGCGCCGGTGCCGGCGATGATGAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCa                                                                                                                                 <  3:174264/149‑1 (MQ=255)
                   aTTGATGCTGCTGCTGTAGGGCGCGCCGGTGCCGGCGATGATGAAGATAAGCGGCGCTTCCTGATCCTCGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCATGAAGTTGTAGGGCAACGAGAACTCGCGATACGGGCGCAAGTACa                                                                                                                                 >  6:279855/1‑149 (MQ=255)
                   aTTGATGCTGCTGCTGTAGGGCGCGCCGGTGCCGGCGAGGATGAAGAGCAGCGGCGCTTCGTGATCGTGGCGGGCCAGGCGGTACTTGAGCTGCTTCACTGGCCAGAAGTGGTCGGGCGACGTGGGGTCGGGCTCGGGgcggcgggg                                                                                                                                   >  6:202405/1‑141 (MQ=255)
                      gATGCTGCTGCTGTAGGGCGCGCCGGTGCCGGCGATTATGAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAgg                                                                                                                               <  7:364503/148‑1 (MQ=255)
                        tgctgctgcTGTAGGGCGCGCCGGTGCCGGCGATGATGAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTTGCGCCCGGAGCGCGAGTTCacgttg                                                                                                                            <  6:174005/149‑6 (MQ=255)
                        tgctgctgcTGTAGGGCGCGCCGGTGCCGGCGATGATGAAGATCAGCGGCGCGTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCGTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTg                                                                                                                            <  8:155912/149‑1 (MQ=255)
                         gctgctgctGTAGGGCGCGCCGGTGCCGGCGATGATGAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTcc                                                                                                                                              >  6:51913/1‑130 (MQ=255)
                         gctgctgctGTAGGGCGCGCCGGTGCCGGCGATGATGAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTcc                                                                                                                                              <  5:51913/130‑1 (MQ=255)
                         gctgctgctGTAGGGCGCGCCGGTGCCGGCGATGATGAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTg                                                                                                                            >  6:253257/1‑148 (MQ=255)
                         gctgctgctGTAGGGCGCGCCGGTGCCGGCGATGATGAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTg                                                                                                                            >  8:364749/1‑148 (MQ=255)
                                  gTAGGGCGCGCCGGTGCCGGCGATGATGAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACt                                                                                                                  >  5:64101/1‑149 (MQ=255)
                                   tAGGGCGCGCCGGTGCCGGCGATGATGAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGAc                                                                                                                   <  1:370884/147‑1 (MQ=255)
                                   tAGGGCGCGCCGGTGCCGGCGATGATGAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACt                                                                                                                  <  6:184728/148‑1 (MQ=255)
                                    catcatcgcCGGCACCGGCGATGATGAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTACACTGCCCAGAAGTTTTCGGGCAACGTGAACTCGCGCTCCCGGCGCAAGATCAGGCTGTAGTCGGACt                                                                                                                  >  4:270200/7‑147 (MQ=255)
                                     gggCGCGCCGGTGCCGGCGATGATGAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAg                                                                                                               <  5:317977/149‑1 (MQ=255)
                                     gggCGCGCCGGTGCCGGCGATGATGAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAg                                                                                                               <  5:326995/149‑1 (MQ=255)
                                               gTGCCGGCGATGATGAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCAtc                                                                                                     >  5:15207/1‑149 (MQ=255)
                                                   cGGCGATGATGAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCAtcgtcg                                                                                                 <  5:95786/149‑1 (MQ=255)
                                                      cGATGATGAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTg                                                                                              <  7:364749/149‑1 (MQ=255)
                                                           attaaGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCa                                                                                          <  5:305927/145‑1 (MQ=255)
                                                                 aTCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAgg                                                                                   <  8:365017/149‑1 (MQ=255)
                                                                   cAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAggg                                                                                  <  3:34443/148‑1 (MQ=255)
                                                                    aGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGgc                                                                                 <  6:378142/148‑1 (MQ=255)
                                                                      cggcgCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAATTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCt                                                                              >  8:175252/1‑149 (MQ=255)
                                                                             tCGTGATCCTGGCGTGCCAGGCGGTACTTGAGCTTCTTCTCTTCGCAGAAGTTGTCGGTCAACTTGATCTCGCGCTCTGTTGGCAAGTTCAGGCTGTAGTGTTTCTGAGTGTTTTCATCGTCGGTGGGCTGCTCAGGGGGCtttttcg                                                                        <  8:8051/148‑6 (MQ=255)
                                                                               gTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCTGGCGCTAGTTCAGGCTGTAGTCGGACTGAGTCATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCcggc                                                                      <  5:279855/148‑1 (MQ=255)
                                                                               gTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCcgtcg                                                                     >  5:16420/1‑149 (MQ=255)
                                                                               gTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCcggcg                                                                     <  6:64101/149‑1 (MQ=255)
                                                                               gTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCcggcg                                                                     >  1:230412/1‑149 (MQ=255)
                                                                                tGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCcggcgg                                                                    >  2:92687/1‑149 (MQ=255)
                                                                                 gATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCggcggc                                                                   >  6:164914/1‑149 (MQ=255)
                                                                                   tCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCgt                                                                 >  6:304016/1‑149 (MQ=255)
                                                                                     cTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCgtgg                                                               <  2:230412/149‑1 (MQ=255)
                                                                                     cTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCgtg                                                                <  5:169364/148‑1 (MQ=255)
                                                                                      tGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGggtggt                                                              >  8:67494/1‑149 (MQ=255)
                                                                                       ggCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCcggc                                                                      <  1:60068/140‑1 (MQ=255)
                                                                                       ggCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCcggc                                                                      >  2:60068/1‑140 (MQ=255)
                                                                                                ggCGGTACTTGAGCTTCTTCACTGCCCAGAAGTGGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGAtggtg                                                    <  5:304016/149‑1 (MQ=255)
                                                                                                          gAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACg                                          >  7:364623/1‑149 (MQ=255)
                                                                                                          gAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACg                                          >  6:292246/1‑149 (MQ=255)
                                                                                                                 ttcACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACGGGTTgg                                    >  2:238831/1‑148 (MQ=255)
                                                                                                                  tcACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACGGGTTGGt                                   >  7:216696/1‑148 (MQ=255)
                                                                                                                   cACTGCCCAGAAGTTGTCGGGCAACGTGAACTCCCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACGGGTTGGt                                   <  4:25269/147‑1 (MQ=255)
                                                                                                                    aCTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACGGGTTGGt                                   >  1:45993/1‑146 (MQ=255)
                                                                                                                     cTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACGGGTTGGTaa                                 >  4:376498/1‑147 (MQ=255)
                                                                                                                     cTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACGGGTTGGTAaga                               <  6:15207/149‑1 (MQ=255)
                                                                                                                     cTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACGGGTTGGTAaga                               <  8:216696/149‑1 (MQ=255)
                                                                                                                      tGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGaa                                            <  5:335222/135‑1 (MQ=255)
                                                                                                                      tGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGaa                                            >  6:335222/1‑135 (MQ=255)
                                                                                                                         ccAGAAGTTGTCGGGCAACGTGAACTAGCGCTCCGGGCGCAAGTTCAGGCTCTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGCCGCTTGTCCGGCGGCGTTGATGAGATGGAGTCCTCGAACGGATTGCTAAGAggt                            >  6:271324/1‑147 (MQ=255)
                                                                                                                           aGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTTGCGATGGTGGCCTCGAACGGGTTGGTAAGAGGGTAg                         <  3:270200/149‑1 (MQ=255)
                                                                                                                           aGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACGGGTTGGTAAGAGGGTAg                         >  7:46875/1‑149 (MQ=255)
                                                                                                                                ttGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCAACGTCGCTGGGTAGCTCAGGTCGCTTTTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACGGGTTGGTAAGAGGGTAGCCAt                     <  5:271324/148‑1 (MQ=255)
                                                                                                                                   tCGGGCAACGTGAATTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACGGGTTGGTAAGAGGGTAGCCATAGCt                 <  7:175252/149‑1 (MQ=255)
                                                                                                                                   tCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACGGGTTGGTAAGAGGCTACCCATAGCt                 >  7:187053/1‑149 (MQ=255)
                                                                                                                                     gggCCACGTGATCTAGCGCTCCGGGCGCAAGTTCAGGCTGTAGGCGGACTGAGTGTTTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACGGGTTGGTAAGAGGGTAGCCATAGCTgg               <  8:187053/149‑1 (MQ=255)
                                                                                                                                      ggCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACGGGTTGGTAAGAGGGTAGCCATAGCTgg               >  8:238409/1‑148 (MQ=255)
                                                                                                                                       gCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACGGGTTGGTAAGAGGGTAGCCATAGCTgg               >  3:154261/1‑147 (MQ=255)
                                                                                                                                        cAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACGGGTTGGTAAGAGGGTAGCCATAGCTGGCAg            >  5:381906/1‑149 (MQ=255)
                                                                                                                                        cAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACGGGTTGGTAAGAGGGTAGCCATAGCTGGCAg            >  5:387072/1‑149 (MQ=255)
                                                                                                                                         aaCGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACGGGTTGGTAAGAGGGTAGCCATAGCTGGCAGc           <  5:292246/149‑1 (MQ=255)
                                                                                                                                           cctGAACTTGCGCCCGGAGCGCGAGTTCACGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCcggcg                                                                     >  3:238923/3‑89 (MQ=25)
                                                                                                                                           cctGAACTTGCGCCCGGAGCGCGAGTTCACGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCcggcg                                                                     <  4:238923/87‑1 (MQ=25)
                                                                                                                                            gTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACGGGTTGGTAAGAGGGTAGCCATAGCTGGCAGCGTc        >  4:43714/1‑149 (MQ=255)
                                                                                                                                            gTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACGGGTTGGTAAGAGGGTAGCCATAGCTGGCAGCGTc        >  2:270179/1‑149 (MQ=255)
                                                                                                                                             tGAACTTGCGCCCGGAGCGCGAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGACAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACGGGTTg                                     <  7:5429/119‑1 (MQ=255)
                                                                                                                                             tGAACTTGCGCCCGGAGCGCGAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGACAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACGGGTTg                                     >  8:5429/1‑119 (MQ=255)
                                                                                                                                              gAACTTGCGCCCGGAGCGCGAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACGGGt                                       <  4:178728/116‑1 (MQ=255)
                                                                                                                                              gAACTTGCGCCCGGAGCGCGAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACGGGt                                       >  3:178728/1‑116 (MQ=255)
                                                                                                                                                 cTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGCCGTGGTGGCGATGGTGGCCTCGAACGGGTTGGTAAGAGGGTAGCCATAGCTGGCAGCGTCGatct   >  7:161803/1‑147 (MQ=255)
                                                                                                                                                  tCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACGGGTTGGTAAGAGGGTAGCCATAGCTGGCAGCGTCGATATc  >  5:351976/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
ACAGCTTTTTCAGGTATTCATTGATGCTGCTGCTGTAGGGCGCGCCGGTGCCGGCGATGATGAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACGGGTTGGTAAGAGGGTAGCCATAGCTGGCAGCGTCGATATC  >  NC_002947/1938692‑1938986

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