Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,083,011 C→A 6.7% D419Y (GAC→TAC)  alaS ← alanine‑‑tRNA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,083,0110CA 6.7% 218.2 / 4.2 72D419Y (GAC→TAC) alaSalanine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base C (38/29);  new base A (2/3);  total (40/32)
Fisher's exact test for biased strand distribution p-value = 6.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.88e-01

CTTGCCCTGGCCCTCGGTGGCGTCGTAACCGAGGAACTCGGTGGCGCTGTCGACCTTGACCAGGCTGTTGTAATCCATGCCGAAGGCGCTGGCGGAGCGGGCGCGCTCACGCTGGGCGTCCATTTCACGCTCGAAGCCGGCTTCGTCGATGGTCAGCTCACGCTCACGAGCGATGTCGGCAGTCAGGTCCATCGGGAAGCCATAGGTGTCGTACAGCTTGAACACCACGTCGCCCGGTACCACGTCGCCTTTCAGCTGGGCCAGGTCTTGCTCGAGGATGCGCAGGCCCTGCTC  >  NC_002947/5082865‑5083158
                                                                                                                                                  |                                                                                                                                                   
cTTGCCCTGGCCCTCGGTGGCGTCGTAACCGAGGAACTCGGTGGCGCTGTCGACCTTGACCAGGCTGTTGTAATCCATGCCGAAGGCGCTGGCGGAGCGGGCGCGCTCACGCTGGGCGTCCATTTCACGCTCGAAGCCGGCTTCGTCGa                                                                                                                                                   >  3:29062/1‑149 (MQ=255)
  tGCCCTGGCCCTCGGTGGCGTCGTAACCGAGGAACTCGGTGGCGCTGTCGACCTTGACCAGGCTGTTGTAAGCCATGCCGAAGGCGCTGGGGGAGCGGGCGCGCTCCCGCGGGGCGGCCATTTCACGCGCGAAGCCGGCTTCGTCGATg                                                                                                                                                 >  3:30087/1‑149 (MQ=255)
           ccTCGGTGGCGTCGTAACCGAGGAACTCGGTGGCGCTGTCGACCTTGACCAGGCTGTTGTAATCCATGCCGAAGGCGCTGGCGGAGCGGGCGCGCTCACGCTGGGCGTCCATTTCACGCTCGAAGCCGGCTTCGTCGATGGTCAgctca                                                                                                                                        >  5:364847/1‑149 (MQ=255)
           ccTCGGTGGCGTCGTAACCGAGGAACTCGGTGGCGCTGTCGACCTTGACCAGGCTGTTGTAATCCATGCCGAAGGCGCTGGCGGAGCGGGCGCGCTCACGCTGGGCGTCCATTTCACGCTCGAAGCCGGCTTCGTCGATGGTCAgctca                                                                                                                                        >  7:68593/1‑149 (MQ=255)
            cTCGGTGGCGTCGTAACCGAGGAACTCGGTGGCGCTGTCGACCTTGACCAGGCTGTTGTAATCCATGCCGAAGGCGCTGGCGGAGCGGGCGCGCTCACGCTGGGCGTCCATTTCACGCTCGAAGCCGGCTTCGTCGATGGTCAgctcac                                                                                                                                       >  2:225810/1‑149 (MQ=255)
            cTCGGTGGCGTCGTAACCGAGGAACTCGGTGGCGCTGTCGACCTTGACCAGGCTGTTGTAATCCATGCCGAAGGCGCTGGCGGAGCGGGCGCGCTCACGCTGGGCGTCCATTTCACGCTCGAAGCCGGCTTCGTAGATGGTCAgctcac                                                                                                                                       >  3:288809/1‑149 (MQ=255)
            cTCGGTGGCGTCGTAACCGAGGAACTCGGTGGCGCTGTCGACCTTGACCAGGCTGTTGTAATCCATGCCGAAGGCGCTGGCGGAGCGGGCGCGCTCACGCTGGGCGTCCATTTCACGCTCGAAGCCGGCTTCGTAGATGGTCAgctcac                                                                                                                                       >  7:296374/1‑149 (MQ=255)
                 tGGCGTCGTAACCGAGGAACTCGGTGGCGCTGTCGACCTTGACCAGGCTGTTGTAATCCATGCCGAAGGCGCTGGCGGAGCGGGCGCGCTCACGCTGGGCGTCCATTTCACGCTCGAAGCCGGCTTCGTAGATGGTCAgctcacgctca                                                                                                                                  <  8:296374/149‑1 (MQ=255)
                 tGGCGTCGTAACCGAGGAACTCGGTGGCGCTGTCGACCTTGACCAGGCTGTTGTAATCCATGCCGAAGGCGCTGGCGGAGCGGGCGCGCTCACGCTGGGCGTCCATTTCACGCTCGAAGCCGGCTTCGTAGATGGTCAgctcacgctca                                                                                                                                  <  4:288809/149‑1 (MQ=255)
                      tcgtAACCGAGGAACTCGGTGGCGCTGTCGACCTTGACCAGACTGTTGTAATCCATGCCGAAGGCGCTGGCGGAGCGGGCGCGCTCACGCTGGGCGTCCATTTCACGCTCGAAGCCGGCTTCGTCGATGGTCAGCTCACGCTCACGAGc                                                                                                                             <  8:71555/149‑1 (MQ=255)
                          aaCCGAGGAACTCGGTGGCGCTGTCGACCTTGACCAGGCTGTTGTAATCCATGCCGAAGGCGCTGGCGGAGCGGGCGCGCTCACGCTGGGCGTCCATTTCACGCTCGAAGCCGGCTTCGTCGATGGTCAGCTCACGCTCACGAGCGATg                                                                                                                         >  4:310442/1‑149 (MQ=255)
                          aaCCGAGGAACTCGGTGGCGCTGACGACCTTGACCAGGCTGTTGTAATCCATGCCGAAGGCGCTGGCGGAGCGGGCGCGCTCACGCTGGGCGTCCATTTCACGCTCGAAGCCGGCTTCGTCGATGGTCAGCTCACGCTCACGAGCGATg                                                                                                                         >  6:19873/1‑149 (MQ=255)
                             cGAGGAACTCGGTGGCGCTGTCGACCTTGACCAGGCTGTTGTAATCCATGCCGAAGGCGCTGGCGGAGCGGGCGCGCTCACGCTGGGCGTCCATTTCACGCTCGAAGCCGGCTTCGTCGATGGTCAGCTCACGCTCACGAGCGATGTCg                                                                                                                      >  3:296781/1‑149 (MQ=255)
                              gAGGAACTCGGTGGCGCTGTCGACCTTGACCAGGCTGTTGTAATCCATGCCTAAGGCGCTGGCGGAGCGGGCGCGCTCACGCTGGGCGTCCATTTCACGCTCGAAGCCGGCTTCGTCGATGGTCAGCTCACGCTCACGAGCGATGTCgg                                                                                                                     <  4:355195/149‑1 (MQ=255)
                              gAGGAACTCGGTGGCGCTGTCGACCTTGACCAGGCTGTTGTAATCCATGCCGAAGGCGCTGGCGGAGCGGGCGCGCTCACGCTGGGCGTCCATTTCACGCTCGAAGCCGGCTTCGTCGATGGTCAGCTCACGCTCACGAGCGATGTCgg                                                                                                                     >  1:84799/1‑149 (MQ=255)
                              gAGGAACTCGGTGGCGCTGTCGACCTTGACCAGGCTGTTGTAATCCATCCCGAAGGCGCTGGCGGAGCGGGCGCGCTCACGCTGGGCGTCCATTTCACGCACGAAGCCGGCTTCGTCGATGGTCAGCTCACGCTCACGAGCGATGTCgg                                                                                                                     >  7:357446/1‑149 (MQ=255)
                                ggAACTCGGTGGCGCTGTCGACCTTGACCAGGCTGTTGTAATCCATGCCGAAGGCGCTGGCGGAGCGGGCGCGCTCACGCTGGGCGTCCATTTCACGCTCGAAGCCGGCTTCGTCGATGGTCAGCTCACGCTCACGAGCGATGTCgg                                                                                                                     >  1:385005/1‑147 (MQ=255)
                                ggAACTCGGTGGCGCTGTCGACCTTGACCAGGCTGTTGTAATCCATGCCGAAGGCGCTGGCGGAGCGGGCGCGCTCACGCTGGGCGTCCATTTCACGCTCGAAGCCGGCTTCGTCGATGGTCAGCTCACGCTCACGAGCGATGTCGGCa                                                                                                                   <  7:77817/149‑1 (MQ=255)
                                  aaCTCGGTGGCGCTGTCGACCTTGACCAGGCTGTTGTAATCCATGCCGAAGGCGCTGGCGGAGCGGGCGCGCTCACGCTGGGCGTCCATTTCACGCTCGAAGCCGGCTTCGTCGATGGTCAGCTCACGCTCACGAGCGATGTCGGCAGt                                                                                                                 <  8:5256/149‑1 (MQ=255)
                                      cGGTGGCGCTGTCGACCTTGACCAGGCTGTTGTAATCCATGCCGAAGGCGCTGGCGGAGCGGGCGCGCTCACGCTGGGCGTCCATTTCACGCTCGAAGCCGGCTTCGTCGATGGTCAGCTCACGCTCACGAGCGATGTCGGCAGTCAgg                                                                                                             >  2:365287/1‑149 (MQ=255)
                                       ggTGGCGCTGTCGACCTTGACCAGGCTGTTGTAATCCATGCCGAAGGCGCTGGCGGAGCGGGCGCGCTCACGCTGGGCGTCCATTTCACGCTCGAAGCCGGCTTCGTCGATGGTCAGCTCACGCTCACGAGCGATGTCGGCAGTCAgg                                                                                                             <  8:19166/148‑1 (MQ=255)
                                         tGGCGCTGTCGACCTTGACCAGGCTGTTGTAATCCATGCCGAAGGCGCTGGCGGAGCGGGCGCGCTCACGCTGGGCGTCCATTTCACGCTCGAAGCCGGCTTCGTCGATGGTCAGCTCACGCTCACGAGCGATGTCGGCAGTCAGGTc                                                                                                           >  3:163504/1‑148 (MQ=255)
                                            cgcTGTCGACCTTGACCAGGCTGTTGTAATCCATGCCGAAGGCGCTGGCGGAGCGGGCGCGCTCACGCTGGGCGTCCATTTCACGCTCGAAGCCGGCTTCGTCGATGGTCAGCTCACGCTCACGAGCGATGTCGGCAGTCAGGTCCATc                                                                                                       >  2:72272/1‑149 (MQ=255)
                                             gcTGTCGACCTTGACCAGGCTGTTGTAATCCATGCCGAAGGCGCTGGCGGAGCGGGCGCGCTCACGCTGGGCGTCCATTTCACGCTCGAAGCCGGCTTCGTCGATGGTCAGCTCACGCTCACGAGCGATGTCGGCAGTCAGGTCCATCg                                                                                                      >  4:126487/1‑149 (MQ=255)
                                              cTGTCGACCTTGACTATGTTGTTGTAATCCATGCCTAATGCGCTGTCGTATCTTGCGCGATCACGCTTGGCGTCCATTTCACGCTCGAAGCTGGGTTCGTAGATTGTCATCTCAGTCTCACGAGCGATGTCGTCAGTCATGTACATCgg                                                                                                     <  4:296781/149‑1 (MQ=255)
                                               tGTCGACCTTGTCCAGGGTGTTGTAATCCTTGCCGAAGGCGCTGGCGGAGCGGGCGCGCTCACGCTGGGCGTCCATTTCACGCTCGTAGCCGGCTTCGTCGATGGTCAGCTCACGCTCACGAGCGATGTCGGCAGTCAGGTCGATCggg                                                                                                    <  8:357446/149‑1 (MQ=255)
                                                 tCGACCTTGACCAGGCTGTTGTAATCCATGCCGAAGGCGCTGGCGGAGCGGGCGCGCTCACGCTGGGCGTCCATTTCACGCTCGAAGCCGGCTTCGTCGATGGTCAGCTCACGCTCACGAGCGATGTCGGCAGTCAGGTCCATCGGGa                                                                                                   <  4:29062/148‑1 (MQ=255)
                                                 gcGACCTTGACCAGGCTGTTGTACTCCATGCCGAAGGCGCTGGCGGAGCGGGCGCGCGCACGCTGGGCGTCCATTTCACGCGCGAAGCCGGCTGCGTCGATGGTCAGCTCACGCTCACGAGCGATGTCGGCAGTCAGGTCCATCggg                                                                                                    <  4:30087/146‑1 (MQ=255)
                                                   gACCTTGACCAGGCTGTTGTAATCCATGCCGAAGGCGCTGGCGGAGCGGGCGCGCTCACGCTGGGCGTCCATTTCACGCTCGAAGCCGGCTTCGTCGATGGTCAGCTCACGCTCACGAGCGATGTCg                                                                                                                      >  6:339971/1‑127 (MQ=255)
                                                         gACCAGGCTGTTGTAATCCATGCCGAAGGCGCTGGCGGAGCGGGCGCGCTCACGCTGGGCGTCCATTTCACGCTCGAAGCCGGCTTCGTCGATGGTCAGCTCACGCTCACGAGCGATGTCGGCAGTCAGGTCCATCGGGAAGCCATAgg                                                                                          >  2:171199/1‑149 (MQ=255)
                                                                    tgtAATCCATGCCGAAGGCGCTGGCGGAGCGGGCGCGCTCACGCTGGGCGTCCATTTCACGCTCGAAGCCGGCTTCGTCGATGGTCAGCTCACGCTCACGAGCGATGTCGGCAGTCAGGTCCATCGGGAAGCCATAGGTGTCGTACAg                                                                                <  8:240770/148‑1 (MQ=255)
                                                                    tgtAATCCATGCCGAAGGCGCTGGCGGAGCGGGCGCGCTCACGCTGGGCGTCCATTTCACGCTCGAAGCCGGCTTCGTCGATGGTCAGCTCACGCTCACGAGCGATGTCGGCAGTCAGGTCCATCGGGAAGCCATAGGTGTCGTACAGc                                                                               >  3:388282/1‑149 (MQ=255)
                                                                    tgtAATCCATGCCGAAGGCGCTGGCGGAGCGGGCGCGCTCACGCTGGGCGTCCATTTCACGCTCGAAGCCGGCTTCGTCGATGGCCAGCTCACGCTCACGAGCGATGTCGGCAGTCAGGTCCATCGGGAAGCCATAGGTGTCGTACAGc                                                                               >  6:305367/1‑149 (MQ=255)
                                                                     gtAATCCATGCCGAAGGCGCTGGCGGAGCGGGCGCGCTCACGCTGGGCGTCCATTTCACGCTCGAAGCCGGCTTCGTCGATGGTCAGCTCACGCTCACGAGCGATGTCGGCAGTCAGGTCCATCGGGAAGCCATAGGTGTCGTACAGc                                                                               <  6:364847/148‑1 (MQ=255)
                                                                       aaTCCATGCCGAAGGCGCTGGCGGAGCGGGCGCGCTCACGCTGGGCGTCCATTTCACGCTCGAAGCCGGCTTCGTCGATGGCCAGCTCACGCTCACGAGCGATGTCGGCAGTCAGGTCCATCGGGAAGCCATAGGTGTCGTACAGCtt                                                                             <  5:305367/148‑1 (MQ=255)
                                                                        aTCCATGCCGAAGGCGCTGGCGGAGCGGGCGCGCTCACGCTGGGCGTCCATTTCACGCTCGAAGCCGGCTTCGTCGATGGTCAGCTCACGCTCACGAGCGATGTCGGCAGTCAGGTCCATCGGGAAGCCATAGGTGTCGTACAGCTTg                                                                            <  5:19873/148‑1 (MQ=255)
                                                                         tCCATGCCGAAGGCGCTGGCGGAGCGGGCGCGCTCACGCTGGGCGTCCATTTCACGCTCGAAGCCGGCTTCGTCGATGGTCAGCTCACGCTCACGAGCGATGTCGGCAGTCAGGTCCATCGGGAAGCCATAGGTGTCGTACAGCTTGaa                                                                          <  1:365287/149‑1 (MQ=255)
                                                                         tCCATGCCGAAGGCGCTGGCGGAGCGGGCGCGCTCACGCTGGGCGTCCATTTCACGCTCGAAGCCGGCTTCGTCGATGGTCAGCTCACGCTCACGAGCGATGTCGGCAGTCAGGTCCATCGGGAAGCCATAGGTGTCGTACAGCTTGaa                                                                          >  5:72974/1‑149 (MQ=255)
                                                                          ccATGCCGAAGGCGCTGGCGGAGCGGGCGCGCTCACGCTGGGCGTCCATTTCACGCTCGAAGCCGGCTTCGTCGATGGTCAGCTCACGCTCACGAGCGATGTCGGCAGTCAGGTCCATCGGGAAGCCATAGGTGTCGTACAGCTTGaa                                                                          >  5:132290/1‑148 (MQ=255)
                                                                             tGCCGAAGGCGCTGGCGGAGCGGGCGCTCTCACGCTGGGCGTCCATTTCACGCTCGAAGCCGGCTTCGTCGATGGTCAGCTCACGCTCACGAGCGATGTCGGCAGTCAGGTCCATCGGGAAGCCATAGGTGTCGTACAGCTTGAAcacc                                                                      <  6:72974/149‑1 (MQ=255)
                                                                             tGCCGAAGGCGCTGGCGGAGCGGGCGCGCTCACGCTGGGCGTCCATTTCACGCTCGAAGCCGGCTTCGTCGATGGTCAGCTCACGCTCACGAGCGATGTCGGCAGTCAGGTCCATCGGGAAGCCATAGGTGTCGTACAGCTTGAacac                                                                       >  4:231957/1‑148 (MQ=255)
                                                                               ccGAAGGCGCTGGCGGAGCGGGCGCGCTCACGCTGGGCGTCCATTTCACGCTCGAAGCCGGCTTCGTCGATGGTCAGCTCACGCTCACGAGCGATGTCGGCAGTCAGGTCCATCGGGAAGCCATAGGTGTCGTACAGCTTGAAcaccac                                                                    >  6:168410/1‑149 (MQ=255)
                                                                               ccGAAGGCGCTGGCGGAGCGGGCGCGCTCACGCTGGGCGTCCATTTCACGCTCGAAGCCGGCTTCGTCGATGGTCAGCTCACGCTCACGAGCCATGTCGGCAGTCAGGTCCAACGGGAAGCCAAAGGTGTCGTACAGCTTGAAcaccac                                                                    >  1:245273/1‑149 (MQ=255)
                                                                                cGAAGGCGCTGGCGGAGCGGGCGCGCTCACGCTGGGCGTCCATTTCACGCTCGAAGCCGGCTTCGTCGATGGTCAGCTCACGCTCACGAGCGATGTCGGCAGTCAGGTCCATCGGGAAGCCATAGGTGTCGTACAGCTTGAACACCACg                                                                   >  4:68710/1‑149 (MQ=255)
                                                                                 gAAGGCGCTGGCGGAGCGGGCGCGCTCACGCTGGGCGTCCATTTCACGCTCGAAGCCGGCTTCGTCGATGGTCAACTCCCGCTCACGAGCGATGTCGGCAGTCAGGTCCATCGGGAAGCCATAGGTGTCGTACAGCTTGAACACCACGt                                                                  >  2:12736/1‑149 (MQ=255)
                                                                                   aGGCGCTGGCGGAGCGGGCGCGCTCACGCTGGGCGTCCATTTCACGCTCGAAGCCGGCTTCGTCGATGGTCAGCTCACGCTCACGAGCGATGTCGGCAGTCAGGTCCATCGGGAAGCCATAGGTGTc                                                                                      <  5:375270/127‑1 (MQ=255)
                                                                                   aGGCGCTGGCGGAGCGGGCGCGCTCACGCTGGGCGTCCATTTCACGCTCGAAGCCGGCTTCGTCGATGGTCAGCTCACGCTCACGAGCGATGTCGGCAGTCAGGTCCATCGGGAAGCCATAGGTGTc                                                                                      >  6:375270/1‑127 (MQ=255)
                                                                                           gCGGAGCGGGCGCGCTCACGCTGGGCGTCCATTTCACGCTCGAAGCCGGCTTCGTCGATGGTCAGCTCACGCTCACGAGCGATGTCGGCAGTCAGGTCCATCGGGAAGCCATAGGTGTCGTACAGCTTGAACACCACGTCGCCCGGt                                                          <  8:68593/147‑1 (MQ=255)
                                                                                              gAGCGGGCGCGCTCACGCTGGGCGTCCATTTCACGCTCGAAGCCGGCTTCGTCGATGGTCAGCTCACGCTCACGAGCGATGTCGGCAGTCAGGTCCATCGGGAAGCCATa                                                                                            >  2:155256/1‑110 (MQ=255)
                                                                                              gAGCGGGCGCGCTCACGCTGGGCGTCCATTTCACGCTCGAAGCCGGCTTCGTCGATGGTCAGCTCACGCTCACGAGCGATGTCGGCAGTCAGGTCCATCGGGAAGCCATa                                                                                            <  1:155256/110‑1 (MQ=255)
                                                                                               aGCGGGCGCGCTCACGCTGGGCGTCCATTTCACGCTCGAAGCCGGCTTCGTCGATGGTCAGCTCACGCTCACGAGCGATGTCGGCAGTCAGGTCCATCggg                                                                                                    >  8:232354/1‑101 (MQ=255)
                                                                                               aGCGGGCGCGCTCACGCTGGGCGTCCATTTCACGCTCGAAGCCGGCTTCGTCGATGGTCAGCTCACGCTCACGAGCGATGTCGGCAGTCAGGTCCATCggg                                                                                                    <  7:232354/101‑1 (MQ=255)
                                                                                                        gcTCACGCTGGGCGTCCAATTCACGCTCGAAGCCGGCTTCGTCGATGGTCAGCTCACGCTCACGAGCGATGTCGGCAGTCAGGTCCATCCGGAAGCCATAGGTGTCGTACAGCTTGCACACCACGTCCCCCGCTACCACGTCGCCATTc                                           >  8:105697/1‑149 (MQ=255)
                                                                                                         cTCACGCTGGGCGTCCATTTCACGCTCGAAGCCGGCTTCGTCGATGGTCAGCTCACGCTCACGAGCGATGTCGGCAGTCAGGTCCATCGGGAAGCCATAGGTGTCGTACAGCTTGAACACCACGTCGCCCGGTACCACGTCGCCTTTCa                                          <  2:84799/149‑1 (MQ=255)
                                                                                                            aCGCTGGGCGTCCATTTCACGCTCGAAGCCGGCTTCGTCGATGGTCAGCTCACGCTCACGAGCGATGTCGGCAGTCAGGTCCATCGGGAAGCCATAGGTGTCGTACAGCTTGAACACCACGTCGCCCGGTACCACGTCGCCTTTCAg                                         <  1:297615/147‑1 (MQ=255)
                                                                                                            aCGCTGGGCGTCCATTTCACGCTCGAAGCCGGCTTCGTCGATGGTCAGCTCACGCTCACGAGCGATGTCGGCAGTCAGGTCCATCGGGAAGCCATAGGTGTCGTACAGCTTGAACACCACGTCGCCCGGTACCACGTCGCCTTTCAGCt                                       >  3:230619/1‑149 (MQ=255)
                                                                                                              gCTGGGCGTCCATTTCACGCTCGAAGCCGGCTTCGTCGATGGTCAGCTCACGCTCACGAGCGATGTCGGCAGTCAGGTCCTTCGGGAAGCCATAGGTGTCGTACAGCTTGAACACCACGTCGCCCGGTACCACGTCGCCTTTCAGCTgg                                     >  5:188718/1‑149 (MQ=255)
                                                                                                              gCTGGGCGTCCATTTCACGCTCGAAGCCGGCTTCGTCGATGGTCAACTCACGCTCACGAGCGATGTCGGCAGTCAGGTCCATCGGGAAGCCATAGGTGTCGTACAGCTTGAACACCACGTCGCCCGGTACCACGTCGCCTTTCAGCTgg                                     <  1:12736/149‑1 (MQ=255)
                                                                                                                  ggCGTCCATTTCACGCTCGAAGCCGGCTTCGTCGATGGTCAGCTCACGCTCACGAGCGATGTCGGCAGTCAGGTCCATCGGGAAGCCATAGGTGTCGTACAGCTTGAACACCACGTCGCCCGGTACCACGTCGCCTTTCAGCTGGGc                                   >  2:198618/1‑147 (MQ=255)
                                                                                                                     gTCCATTTCACGCTCGAAGCCGGCTTCGTCGATGGTCAGCTCACGCTCACGAGCGATGTCGGCAGTCAGGTCCATCGGGAAGCCATAGGTGTCGTACAGCTTGAACACCACGTCGCCCGGTACCACGTCGCCTTTCAGCTGGGCCAGGt                              >  5:391865/1‑149 (MQ=255)
                                                                                                                        cATTTCACGCTCGAAGCCGGCTTCGTCGATGGTCAGCTCACGCTCACGAGCGATGTCGGCAGTCAGGTCCATCGGGAAGCCATAGGTGTCGTACAGCTTGAACACCACGTCGCCCGGTACCACGTCGCCTTTCAGCTGGGCCAGGTCtt                           <  3:310442/149‑1 (MQ=255)
                                                                                                                        cATTTCACGCTCGAAGCCGGCTTCGTCGATGGTCAGCTCACGCTCACGAGCGATGTCGGCAGTCAGGTCCATCGGGAAGCCATAGGTGTCGTACAGCTTGAACACCACGTCGCCCGGTACCACGTCGCCTTTCAGCTGGGCCAGGTCtt                           <  2:385005/149‑1 (MQ=255)
                                                                                                                              aCGCTCGAAGCCGGCTTCGTCGATGGTCAGCTCACGCTCACGAGCGATGTCGGCAGTCAGGTCCTTCGGGAAGCCATAGGTGTCGTACAGCTTGAACACCACGTCGCCCGGTACCACGTCGCCTTTCAGCTGGGCCAGGTCTTGCTCGa                     <  6:188718/149‑1 (MQ=255)
                                                                                                                              aCGCTCGAAGCCGGCTTCGTCGATGGTCAGCTCACGCTCACGAGCGATGTCGGCAGTCAGGTCCATCGGGAAGCCATAGGTGTCGTACAGCTTGAACACCACGTCGCCCGGTACCACGTCGCCTTTCAGCTGGGCCAGGTCTTGCTc                       >  4:248468/1‑147 (MQ=255)
                                                                                                                               cGCTCGAAGCCGGCTTCGGCGTTGGTCAGCTCACGCTCACGAGCGATGTCGGCAGTCAGGTCCATCGGGAAGCCATAGGTGTCGTACAGCTTGAACACCACGTCGCCCGGTACCACGTCGCGTTTCAGCTGGGCCAGGTCTTGCTCGAg                    <  2:245273/149‑1 (MQ=255)
                                                                                                                                gCTCGAAGCCGGCTTCGTCGATGGTCAGCTCACGCTCACGAGCGATGTCGGCAGTCAGGTCCATCGGGAAGCCATAGGTGTCGTACAGCTTGAACACCACGTCGCCCGGTACCACGTCGCCTTTCAGCTGGGCCAGGTCTTGCTCGAg                    >  4:387008/1‑148 (MQ=255)
                                                                                                                                  tCGAAGCCGGCTTCGTCGATGGTCAGCTCACGCTCACGAGCGATGTCGGCAGTCAGGTCCATCGGGAAGCAATAGGTGTCGTACAGCTTGAACACCACGTCGCCCGGTACCACGTCGCCTTTCAGCTGGGCCAGGTCTTCCTCGAGGAt                 >  3:242907/1‑149 (MQ=255)
                                                                                                                                    gAAGCCGGCTTCGTCGATGGTCAGCTCACGCTCACGAGCGATGTCGGCAGTCAGGTCCATCGGGAAGCCATAGGTGTCGTACAGCTTGAACACCACGTCGCCCGGTACCACGTCGCCTTTCAGCTGGGCCAGGTCTTGCTCGAGGAt                 <  1:171199/147‑1 (MQ=255)
                                                                                                                                     aaGCCGGCTTCGTCGATGGTCAGCTCACGCTCACGAGCGATGTCGGCAGTCAGGTCCATCGGGAAGCCATAGGTGGCGTACAGCTTGAACACCACGTCGCCCGGTACCACGTCGCCTTTCAGCTGGGCCAGGTCTTGCTCGAGGATgcg              <  3:348001/149‑1 (MQ=255)
                                                                                                                                      aGCCGGCTTCGTCGATGGTCAGCTCACGCTCACGAGCGATGTCGGCAGTCAGGTCCATCGGGAAGCCATAGGTGTCGTACAGCTTGAACACCACGTCGCCCGGTACCACGTCGCCTTTCAGCTGGGCCAGGTCTTGCTCGAGGATgcgc             <  1:225810/149‑1 (MQ=255)
                                                                                                                                          ggCTTCGTCGATGGTCAGCTCACGCTCACGAGCGATGTCGGCAGTCAGGTCCATCGGGAAGCCATAGGTGTCGTACAGCTTGAACACCACGTCGCCCGGTACCACGTCGCCTTTCAGCTGGGCCAGGTCTTGCTCGAGGATGCGCAGGc         <  1:198618/149‑1 (MQ=255)
                                                                                                                                          ggCTTCGTCGATGGTCAGCTCACGCTCACGAGCGATGTCGGCAGTCAGGTCCATCGGGAAGCCATAGGTGTCGTACAGCTTGAACACCACGTCGCCCGGTACCACGTCGCCTTTCAGCTGGGCCAGGTCTTGCTCGAGGATGCGCAGGc         >  2:312563/1‑149 (MQ=255)
                                                                                                                                               cgtcgATGGTCAGCTCACGCTCACGAGCGATGTCGGCAGTCAGGTCCATCGGGAAGCCATAGGTGTCGTACAGCTTGAACACCACGTCGCCCGGTACCACGTCGCCTTTCAGCTGGGCCAGGTCTTGCTCGAGGATGCGCAGGCCCTg     <  3:231957/148‑1 (MQ=255)
                                                                                                                                                 tcgATGGTCAGCTCACGCTCACGAGCGATGTCGGCAGTCAGGTCCATCGGGAAGCCATAGGTGTCGTACAGCTTGAACACCACGTCGCCCGGTACCACGTCGCCTTTCAGCTGGGCCAGGTCTTGCTCGAGGATGCGCAGGCCCTGCTc  >  5:189834/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
CTTGCCCTGGCCCTCGGTGGCGTCGTAACCGAGGAACTCGGTGGCGCTGTCGACCTTGACCAGGCTGTTGTAATCCATGCCGAAGGCGCTGGCGGAGCGGGCGCGCTCACGCTGGGCGTCCATTTCACGCTCGAAGCCGGCTTCGTCGATGGTCAGCTCACGCTCACGAGCGATGTCGGCAGTCAGGTCCATCGGGAAGCCATAGGTGTCGTACAGCTTGAACACCACGTCGCCCGGTACCACGTCGCCTTTCAGCTGGGCCAGGTCTTGCTCGAGGATGCGCAGGCCCTGCTC  >  NC_002947/5082865‑5083158

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