Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,011,632 G→A 5.2% A58T (GCC→ACC)  yqeA → putative amino acid kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,011,6320GA 5.2% 246.7 / 2.1 77A58T (GCC→ACC) yqeAputative amino acid kinase
Reads supporting (aligned to +/- strand):  ref base G (34/39);  new base A (2/2);  total (36/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.82e-01

CCTGGGCGGGAATGCGCTGGGCGATGACCTGGCCGGGCAAATGAAAGCGGTAAAAATTACTTCTCAGGCAATTGTTGATTTAATTGCCCAGGGACATGAAGTTATCGTCACTCATGGTAACGGGCCGCAGGTGGGCATGATTAACCAGGCCTTTGAAGCTGCCGCAAAAACTGAAGCGCACTCGCCAATGCTGCCGATGTCTGTCTGTGTTGCATTAAGCCAGGGTTATATTGGTTACGATCTACAAAACGCTTTAAGAGAAGAACTGCTTTCTCGGGGCATTAATAAACCTGTAG  >  NC_000913/3011484‑3011779
                                                                                                                                                    |                                                                                                                                                   
ccTGGGCGGGAATGCGCTGGGCGATGACCTGGCCGGGCAAATGAAAGCGGTAAAAATTACTTCTCAGGCAATTGTTGATTTAATTGCCCAGCGACATGAAGTTATCGTCACTCATCGTAACGGGCCGCAGGTGGGCATGATTAACCAgg                                                                                                                                                     >  6:168108/1‑149 (MQ=255)
     gCGGGAATGCGCTGGGCGATGACCTGGCCGGGCAAATGAAAGCGGTAAAAATTACTTCTCAGGCAATTGTTGATTTAATTGCCCAGGGACATGAAGTTATCGTCACTCATGGTAACGGGCCGCAGGTGGGCATGATTAACCAGGCCttt                                                                                                                                                >  5:215751/1‑149 (MQ=255)
       gggAATGCGCTGGGCGATGACCTGGCCGGGCAAATGAAAGCGGTAAAAATTACTTCTCAGGCAATTGTTGATTTAATTGCCCAGGGACATGAAGTTATCGTCACTCATGGTAACGGGCCGCAGGTGGGCATGATTAACCAGGCCTTTg                                                                                                                                               >  1:152865/1‑148 (MQ=255)
         gAATGCGCTGGGCGATGACCTGGCCGGGCAAATGAAAGCGGTAAAAATTACTTCTCAGGCAATTGTTGATTTAATTGCCCAGGGACATGAAGTTATCGTCACTCATGGTAACGGGCCGCAGGTGGGCATGATTAACCAGGCCTTTGAAg                                                                                                                                            >  5:295215/1‑149 (MQ=255)
              cgcTGGGCGATGACCTGGCCGGGCAAATGAAAGCGGTAAAAATTACTTCTCAGGCAATTGTTGATTTAATTGCCCAGGGACATGAAGTTATCGTCACTCATGGTAACGGGCCGCAGGTGGGCATGATTAACCAGGCCTTTGa                                                                                                                                              >  8:352932/1‑142 (MQ=255)
              cgcTGGGCGATGACCTGGCCGGGCAAATGAAAGCGGTAAAAATTACTTCTCAGGCAATTGTTGATTTAATTGCCCAGGGACATGAAGTTATCGTCACTCATGGTAACGGGCCGCAGGTGGGCATGATTAACCAGGCCTTTGa                                                                                                                                              <  7:352932/142‑1 (MQ=255)
               gcTGGGCGATGACCTGGCCGGGCAAATGAAAGCGGTTAAAATTACTTCTCAGGCAATTGTTGATTTAATTGCCCAGGGACATGAAGTTATCGTCACTCATGGTAACGGGCCGCAGGTGGGCATGATTAACCAGGCCTTTGAAGCTGCCg                                                                                                                                      >  3:239700/1‑149 (MQ=255)
                cTGGGCGATGACCTGGCCGGGCAAATGAAAGCGGTAAAAATTACTTCTCAGGCAATTGTTGATTTAATTGCCCAGGGACATGAAGTTATCGTCACTCATGGTAACGGGCCGCAGGTGGGCATGATTAACCAGGCCTTTGAAGCTGCCGc                                                                                                                                     <  6:295215/149‑1 (MQ=255)
                 tGGGCGATGACCTGGCCGGGCAAATGAAAGCGGTAAAAATTACTTCTCAGGCAATTGTTGATTTAATTGCCCAGGGACATGAAGTTATCGTCACTCATGGTAACGGGCCGCAGGTGGGCATGATTAACCAGGCCTTTGAAGCTGCCGc                                                                                                                                     >  8:220246/1‑148 (MQ=255)
                   ggCGATGACCTGGCCGGGCAAATGAAAGCGGTAAAAATTACTTCTCAGGCAATTGTTGATTTAATTGCCCAGGGACATGAAGTTATCGTCACTCATGGTAACGGGCCGCAGGTGGGCATGATTAACCAGGCCTTTGAAGCTGCCGCaaa                                                                                                                                  <  8:26681/149‑1 (MQ=255)
                     cGATGACCTGGCCGGGCAAATGAAAGCGGTAAAAATTACTTCTCAGGCAATTGTTGATTTAATTGCCCAGGGACATGAAGTTATCGTCACTCATGGTAACGGGCCGCAGGTGGGCATGATTAACCAGGCCTTTGAAGCTGCCGCaaaaa                                                                                                                                <  7:127489/149‑1 (MQ=255)
                      gATGACCTGGCCGGGCAAATGAAAGCGGTAAAAATTACTTCTCAGGCAATTGTTGATTTAATTGCCCAGGGACATGAAGTTATCGTCACTCATGGTAACGGGCCGCAGGTGGGCATGATTAACCAGGCCTTTGAAGCTGCCGCAAAAAc                                                                                                                               >  1:210125/1‑149 (MQ=255)
                            cTGGCCGGGCAAATGAAAGCGGTAAAAATTACTTCTCAGGCAATTGTTGATTTAATTGCCCAGGGACATGAAGTTATCGTCACTCATGGTAACGGGCCGCAGGTGGGCATGATTAACCAGGCCTTTGAAGCTGCCGCAAAAACTGAAgc                                                                                                                         >  6:17746/1‑149 (MQ=255)
                            cTGGCCGGGCAAATGAAAGCGGTAAAAATTACTTCTCAGGCAATTGTTGATTTAATTGCCCAGGGACATGAAGTTATCGTCACTCATGGTAACGGGCCGCAGGTGGGCATGATTAACCAGGCCTTTGAAGCTGCCGCAAAAACTGAAgc                                                                                                                         >  8:303674/1‑149 (MQ=255)
                            cTGGCCGGGCAAATGAAAGCGGTAAAAATTACCTCTCAGGCAATTGTTGATTTAATTGCCCAGGGACATGAAGTTATCGTCACTCATGGTAACGGGCCGCAGGTGGGCATGATTAACCAGGCCTTTGAAGCTGCCGCAAAAACTGAAgc                                                                                                                         <  2:312034/149‑1 (MQ=255)
                                      aaaTGAAAGCGGTAAAAATTACTTCTCAGGCAATTGTTGATTTAATTGCCCAGGGACATGAAGTTATCGTCACTCATGGTAACGGGCCGCAGGTGGGCATGATTAACCAGGCCTTTGAAGCTGCCGCAAAAACTGAAg                                                                                                                          <  6:176897/138‑1 (MQ=255)
                                      aaaTGAAAGCGGTAAAAATTACTTCTCAGGCAATTGTTGATTTAATTGCCCAGGGACATGAAGTTATCGTCACTCATGGTAACGGGCCGCAGGTGGGCATGATTAACCAGGCCTTTGAAGCTGCCGCAAAAACTGAAg                                                                                                                          >  5:176897/1‑138 (MQ=255)
                                       aaTGAAAGCGGTAAAAATTACATCTCAGGCAATTGTTGATTTAATTGCCCAGGGAAATGAAGTTATCGTCACTCATGGTAAGGGGCCGCAGGTGGGCATGATTAACAAGGCCTTTGAAGCTGCCGCAAAAACTGAAGCGCACACGCCaa                                                                                                              <  5:168108/149‑1 (MQ=255)
                                             aGCGGTAAAAATTACTTCTCAGGCAATTGTTGATTTAATTGCCCAGGGACATGAAGTTATCGTCACTCATGGTAACGGGCCGCAGGTGGGCATGATTAACCAGGCCTTTGAAGCTGCCGCAAAAACTGAAGCGCACTCGCCAAtgctgc                                                                                                        <  2:363651/149‑1 (MQ=255)
                                                    aaaaTTACTTCTCAGGCAATTGTTGATTTAATTGCCCAGGGACATGAAGTTATCGTCACTCATGGTAACGGGCCGCAGGTGGGCATGATTAACCAGGCCTTTGAAGCTGCCGCAAAAACTGAAGCGCACTCGCCAATGCTGCCGAtgt                                                                                                  <  8:287769/148‑1 (MQ=255)
                                                       aTTACTTCTCAGGCAATTGTTGATTTAATTGCCCAGGGACATGAAGTTATCGTCACTCATGGTAACGGGCCGCAGGTGGGCATGATTAACCAGGCCTTTGAAGCTGCCGCAAAAACTGAAGCGCACTCGCCAATGCTGCCGAtgtctgt                                                                                              >  4:106541/1‑149 (MQ=255)
                                                            ttCTCAGGCAATTGTTGATTTAATTGCCCAGGGACATGAAGTTATCGTCACTCATGGTAACGGGCCGCAGGTGGGCATGATTAACCAGGCCTTTGAAGCTGCCGCAAAAACTGAAGCGCACTCGCCAATGCTGCCGATGTCTGTCtgtg                                                                                         >  5:233129/1‑149 (MQ=255)
                                                             tctcAGGCAATTGTTGATTTAATTGCCCAGGGACATGAAGTTATCGTCACTCATGGTAACGGGCCGCAGGTGGGCATGATTAACCAGACCTTTGAAGCTGCCGCAAAAACTGAAGCGCACTCGCCAATGCTGCCGATGTCTGTCTgtgt                                                                                        >  3:35338/1‑149 (MQ=255)
                                                              ctcAGGCAATTGTTGATTTAATTGCCCAGGGACATGAAGTTATCGTCACTCATGGTAACGGGCCGCAGGTGGGCATGATTAACCAGGCCTTTGAAGCTGCCGCAAAAACTg                                                                                                                             >  6:269949/1‑111 (MQ=255)
                                                              ctcAGGCAATTGTTGATTTAATTGCCCAGGGACATGAAGTTATCGTCACTCATGGTAACGGGCCGCAGGTGGGCATGATTAACCAGGCCTTTGAAGCTGCCGCAAAAACTg                                                                                                                             <  5:269949/111‑1 (MQ=255)
                                                              ctcAGGCAATTGTTGATTTAATTGCCCAGGGACATGAAGTTATCGTCACTCATGGTAACGGGCCGCAGGTGGGCATGATTAACCAGGCCTTTGAAGCTGCCGCAAAAACTGAAGCGCACTCGCCAATGCTGCCGATGTCTGTCTgtgt                                                                                        <  6:233129/148‑1 (MQ=255)
                                                              ctcAGGCAATTGTTGATTTAATTGCCCAGGGACATGAAGTTATCGTCACTCATGGTAACGGGCCGCAGGTCGGCATGATTAACCAGGCCTTTGAAGCTGCCGCAAAAACTGAAGCGCACTCGCCAATGCTGCCGATGTCTGTCTgtgta                                                                                       >  2:256432/1‑148 (MQ=255)
                                                                   gCAATTGTTGATTTAATTGCCCAGGGACATGAAGTTATCGTCACTCATGGTAACGGGCCGCAGGTGGGCATGATTAACCAGGCCTTTGAAGCTGCCGCAAAAACTGAAGCGCACTCGCCAATGCTGCCGATGTCTGTCTGTGTTGCAt                                                                                   >  5:113728/1‑148 (MQ=255)
                                                                   gCAATTGTTGATTTAATTGCCCAGGGACATGAAGTTATCGTCACTCATGGTAACGGGCCGCAGGTGGGCATGATTAACCAGGCCTTTGAAGCTGCCGCAAAAACTGAAGCGCACTCGCCAATGCTGCCGATGTCTGTCTGTGTTGCAt                                                                                   >  3:228900/1‑148 (MQ=255)
                                                                     aaTTGTTGATTTAATTGCCCAGGGACATGAAGTTATCGTCACTCATGGTAACGGGCCGCAGGTGGGCATGATTAACCAGGCCTTTGAAGCTGCCGCAAAAACTGAAGCGCACTCGCCAATGCTGCCGATGTCTGTCTGTGTTGCATTaa                                                                                <  1:40328/149‑1 (MQ=255)
                                                                       ttgttgATTTAATTGCCCAGGGACATGAAGTTATCGTCACTCATGGTAACGGGCCGCAGGTGGGCATGATTAACCAGGCCTTTGAAGCTGCCGCAAAAACTGAAGCGCACTCGCCAATGCTGCCGATGTCTGTCTGTGTTGCATTAAGc                                                                              <  2:52826/149‑1 (MQ=255)
                                                                       ttgttgATTTAATTGCCCAGGGACATGAAGTTATCGTCACTCATGGTAACGGGCCGCAGGTGGGCATGATTAACCAGGCCTTTGAAGCTGCCGCAAAAACTGAAGCGCACTCGCCAATGCTGCCGATGTCTGTCTGTGTTGCATTAAGc                                                                              <  3:106541/149‑1 (MQ=255)
                                                                          ttgATTTAATTGCCCAGGGACATGAAGTTATCGTCACTCATGGTAACGGGCCGCAGGTGGGCATGATTAACCAGGCCTTTGAAGCTGCCGCAAAAACTGAAGCGCACTCGCCAATGCTGCCGATGTCTGTCTGTGTTGCATTAAGCCAg                                                                           <  7:225783/149‑1 (MQ=255)
                                                                          ttgATTTAATTGCCCAGGGACATGAAGTTATCGTCACTCATGGTAACGGGCCGCAGGTGGGCATGATTAACCAGACCTTTGAAGCTGCCGCAAAAACTGAAGCGCACTCGCCAATGCTGCCGATGTCTGTCTGTGTTGCATTAAGCCa                                                                            <  4:35338/148‑1 (MQ=255)
                                                                             aTTTAATTGCCCAGGGACATGAAGTTATCGTCACTCATGGTAACGGGCCGCAGGTGGGCATGATTAACCAGACCTTTGAAGCTGCCGCAAAAACTGAAGCGCACTCGCCAATGCTGCCGATGTCTGTCTGTGTTGCATTAAGCCAGGGt                                                                        >  5:358082/1‑149 (MQ=255)
                                                                              tttAATTGCCCAGGGACATGAAGTTATCGTCACTCATGGTAACGGGCCGCAGGTGGGCATGATTAACCAGGCCTTTGAAGCTGCCGCAAAAACTGAAGCGCACTCGCCAATGCTGCCGATGTCTGTCTGTGTTGCATTAAGCCAGGGtt                                                                       <  8:354062/149‑1 (MQ=255)
                                                                               ttAATTGCCCAGGGACATGAAGTTATCGTCACTCATGGTAACGGGCCGCAGGTGGGCATGATTAACCAGGCCTTTGAAGCTGCCGCAAAAACTGAAGCGCACTCGCCAATGCTGCCGATGTCTGTCTGTGTTGCATTAAGCCAGGGTta                                                                      >  5:322230/1‑149 (MQ=255)
                                                                                       ccAGGGACATGAAGTTATCGTCACTCATGGTAACGGGCCGCAGGTGGGCATGATTAACCAGGCCTTTGAAGCTGCCGCAAAAACTGAAGCGCACTCGCCAATGCTGCCGATGTCTGTCTGTGTTGCATTAAGCCAGGGTTATATTGGt                                                               >  4:146342/1‑148 (MQ=255)
                                                                                        cAGGGACATGAAGTTATCGGATCTCATGGTTACGGGCCGCAGGTGGGCATGATTAACCAGGCCTTTGAAGCTGCCGCAAATACTGAAGCGCACTCGCCAATGCTGCCGATGTCTGTCTTTGTTGCATTATGACAGGGTTATATTGGTTa                                                             <  2:224582/149‑1 (MQ=255)
                                                                                         aGGGACATGAAGTTATCGTCACTCATGGTAACGGGCCGCAGGTGGGCATGATTAACCAGGCCTTTGAAGCTGCCGCAAAAACTGAAGCGCACTCGCCAATGCTGCCGATGTCTGTCTGTGTTGCATTAAGCCAGGGTTATATTGGTTAc                                                            >  1:101784/1‑149 (MQ=255)
                                                                                         aGGGACATGAAGTTATCGTCACTCATGGTAACGGGCCGCAGGTGGGCATGATTAACCAGGCCTTTGAAGCTGCCGCAAAAACTGAAGCGCACTCGCCAATGCTGCCGATGTCTGTCTGTGTTGCATTAAGCCAGGGTTATATTGGTTAc                                                            >  7:13905/1‑149 (MQ=255)
                                                                                               aTGAAGTTATCGTCACTCATGGTAACGGGCCGCAGGTGGGCATGATTAACCAGGCCTTTGAAGCTGCCGCAAAAACTGAAGCGCACTCGCCAATGCTGCCGATGTCTGTCTGTGTTGCATTAAGCCAGGGTTATATTGGTTACGATCTa                                                      <  7:176999/149‑1 (MQ=255)
                                                                                                tGAAGTTATCGTCACTCATGGTAACGGGCCGCAGGTGGGCATGATTAACCAGGCCTTTGAAGCTGCCGCAAAAACTGAAGCGCACTCGCCAATGCTGCCGATGTCTGTCTGTGTTGCATTAAGCCAGGGTTATATTGGTTACGATCTa                                                      >  5:229166/1‑148 (MQ=255)
                                                                                                tGAAGTTATCGTCACTCATGGTAACGGGCCGCAGGTGGGCATGATTAACCAGGCCTTTGAAGCTGCCGCAAAAACTGAAGCGCACTCGCCAATGCTGCCGATGTCTGTCTGTGTTGCATTAAGCCAGGGTTATATTGGTTACGATCTAc                                                     <  5:63636/149‑1 (MQ=255)
                                                                                                 gAAGTTATCGTCACTCATGGTAACGGGCCGCAGGTGGGCATGATTAACCAGGCCTTTGAAGCTGCCGCAAAAACTGAAGCGCACTCGCCAATGCTGCCGATGTCTGTCTGTGTTGCATTAAGCCAGGGTTATATTGGTTACGATCTAc                                                     >  2:269167/1‑148 (MQ=255)
                                                                                                    gTTATCGTCACTCATGGTAACGGGCCGCAGGTGGGCATGATTAACCAGGCCTTTGAAGCTGCCGCAAAAACTGAAGCGCACTCGCCAATGCTGCCGATGTCTGTCTGTGTTGCATTAAGCCAGGGTTATATTGGTTACGATCTACaaaa                                                 <  7:220246/149‑1 (MQ=255)
                                                                                                       aTCGTCACTCATGGTAACGGGCCGCAGGTGGGCATGATTAACCAGGCCTTTGAAGCTGCCGCAAAAACTGAAGCGCACTCGCCAATGCTGCCGATGTCTGTCTGTGTTGCATTAAGCCAGGGTTATATTGGTTACGATCTACAAAACGc                                              >  2:88823/1‑149 (MQ=255)
                                                                                                        tCGTCACTCATGGTAACGGGCCGCAGGTGGGCATGATTAACCAGGCCGTTGAAGCTGCCGCAAAAATTGAAGCGCACTCGCCAATGCTGCCGATGTCTGTCTGTGTTGCATTAAGCCAGGGTTATATTGGTTACGATCTACAAAACGCt                                             <  6:215751/149‑1 (MQ=255)
                                                                                                             aCTCATGGTAACGGGCCGCAGGTGGGCATGATTAACCAGGCCTTTGAAGCTGCCGCAAAAACTGAAGCGCACTCGCCAATGCTGCCGATGTCTGTCTGTGTTGCATTAAGCCAGGGTTATATTGGTTACGATCTACAAAACGCTTTaa                                         <  2:101784/148‑1 (MQ=255)
                                                                                                             aCTCATGGTAACGGGCCGCAGGTGGGCATGATTAACCAGGCCTTTGAAGCTGCCGCAAAAACTGAAGCGCACTCGCCAATGCTGCCGATGGCTGTCTGTGTTGCATTAAGCCAGGGTTATATTGGTTACGATCTACAAAACGCTTTaa                                         <  6:113728/148‑1 (MQ=255)
                                                                                                               tCATGGTAACGGGCCGCAGGTGGGCATGATTAACCAGGCCTTTGAAGCTGCCGCAAAAACTGAAGCGCACTCGCCAATGCTGCCGATGTCTGTCTGTGTTGCATTAAGCCAGGGTTATATTGGTTACGATCTACAAAACGCTTTaa                                         >  2:118040/1‑146 (MQ=255)
                                                                                                               tCATGGTAACGGGCCGCAGGTGGGCATGATTAACCAGGCCTTTGAAGCTGCCGCAAAAACTGAAGCGCACTCGCCAATGCTGCCGATGTCTGTCTGTGTTGCATTAAGCCAGGGTTATATTGGTTACGATCTACAAAACGCTTTaa                                         >  7:12789/1‑146 (MQ=255)
                                                                                                                cATGGTAACGGGCCGCAGGTGGGCATGATTAACCAGGCCTTTGAAGCTGCCGCAAAAACTGAAGCGCACTCGCCAATGCTGCCGATGTCTGTCTGTGTTGCATTAAGCCAGGGTTATATTGGTTACGATCTACAAAACGCTTTAagag                                      <  8:12789/148‑1 (MQ=255)
                                                                                                                 agggTAACGGGCCGCAGGTGGGCATGCGTAACCAGGCCTTTGAAGCTGCCGCAAAAACTGCAGCGCACTCGCCAATGCTGCCGATGTCTGTCTGTGTTGCATTAAGCCTGGGTTATATTGGTTACGATCTACAAAACGCTTTAAgaga                                     <  4:239700/146‑1 (MQ=255)
                                                                                                                 aTGGTAACGGGCCGCAGGTGGGCATGATTAACCAGGCCTTTGAAGCTGCCGCAAAAACTGAAGCGCACTCGCCAATGCTGCCGATTTCTGTCTGTGTTGCATTAAGCCAGGGTTATATTGGTTACGATCTACAAAACGCTTTAAGagaa                                    <  6:144684/149‑1 (MQ=255)
                                                                                                                 aTGGTAACGGGCCGCAGGTGGGCATGATTAACCAGGCCTTTGAAGCTGCCGCAAAAACTGAAGCGCACTCGCCAATGCTGCCGATGTCTGTCTGTGTTGCATTAAGCCAGGGTTATATTGGTTACGATCTACAAAACGCTTTAagag                                      >  2:239200/1‑147 (MQ=255)
                                                                                                                 aTGGTAACGGGCCGCAGGTGGGCATGATTAACCAGGCCTTTGAAGCTGCCGCAAAAACTGAAGCGCACTCGCCAATGCTGCCGATGTCTGTCTGTGTTGCATTAAGCCAGGGTTATATTGGTTACGATCTACAAAACGCTTTAAGagaa                                    <  6:322230/149‑1 (MQ=255)
                                                                                                                  tGGTAACGGGCCGCAGGTGGGCATGATTAACCAGGCCTTTGAAGCTGCCGCAAAAACTGAAGCGCACTCGCCAATGCTGCCGATGTCTGTCTGTGTTGCATTAAGCCAGGGTTATATTGGTTACGATCTACAAAACGCTTTAAGagaag                                   <  4:228900/149‑1 (MQ=255)
                                                                                                                   ggTAACGGGCCGCAGGTGGGCATGATTAACCAGGCCTTTGAAGCTGCCGCAAAAACTGAAGCGCACTCGCCAATGCTGCCGATGTCTGTCTGTGTTGCATTAAGCCAGGGTTATATTGGATACGATCTACAAAACGCTTTAAGagaag                                   >  1:347158/1‑148 (MQ=255)
                                                                                                                     tAACGGGCCGCAGGTGGGCATGATTAACCAGGCCTTTGAAGCTGCCGCAAAAACTGAAGCGCACTCGCCAATGCTGCCGATGTCTGTCTGTGTTGCATTAAGCCAGGGTTATATTGGTTACGATCTACAAAACGCTTTAAGAGAAGAAc                                >  7:96480/1‑149 (MQ=255)
                                                                                                                     tAACGGGCCGCAGGTGGGCATGATTAACCAGGCCTTTGAAGCTGCCGCAAAAACTGAAGCGCACTCGCCAATGCTGCCGATGTCTGTCTGTGTTGCATTAAGCCAGGGTTATATTGGTTACGATCTACAAAACGCTTTAAGAGAAGAAc                                <  1:88823/149‑1 (MQ=255)
                                                                                                                         gggCCGCAGGTGGGCATGATTAACCAGGCCTTTGAAGCTGCCGCAAAAACTGAATCGCACTCGCCAATGCTGCCGATGTCTGTCTGTGTTGCATTAAGCCAGGGTTATATTGGTTACGATCTACAAAACGCTTTAAGAGAAGAACTGCt                            <  2:210125/149‑1 (MQ=255)
                                                                                                                            ccGCAGGTGGGCATGATTAACCAGGCCTTTGAAGCTGCCGCAAAAACTGAAGCGCACTCGCCAATGCTGCCGATGTCTGTCTGTGTTGCATTAAGCCAGGGTTATATTGGTTACGATCTACAAAACGCTTTAAGAGAAGAACTGCTTtc                         <  7:303674/149‑1 (MQ=255)
                                                                                                                             cGCAGGTGGGCATGATTAACCAGGCCTTTGAAGCTGCCGCAAAAACTGAAGCGCACTCGCCAATGCTGCCGATGTCTGTCTGTGTTGCATTAAGCCAGGGTTATATTGGTTACGATCTACAAAACGCTTTAAGAGAAGAACTGCTTtct                        >  7:143924/1‑149 (MQ=255)
                                                                                                                                 ggTGGGCATGATTAACCAGACCTTTGAAGCTGCCGCAAAAACTGAAGCGCACTCGCCAATGCTGCCGATGTCTGTCTGTGTTGCATTAAGCCAGGGTTATATTGGTTACGATCTACAAAACGCTTTAAGAGAAGAACTGCTTTCTCggg                    <  6:358082/149‑1 (MQ=255)
                                                                                                                                  gTGGGCATGATTAACCAGGCCTTTGAAGCTGCCGCAAAAACTGAAGCGCACTCGCCAATGCTGCCGATGTCTGTCTGTGTTGCATTAAGCCAGGGTTATATTGGTTACGATCTACAAAACGCTTTAAGAGAAGAACTGCTTTCTCgggg                   >  2:12975/1‑149 (MQ=255)
                                                                                                                                     ggCATGATTAACCAGGCCTTTGAAGCTGCCGCAAAAACTGAAGCGCACTCGCCAATGCTGCCGATGTCTGTCTGTGTTGCATTAAGCCAGGGTTATATTGGTTACGATCTACAAAACGCTTTAAGAGAAGAACTGCTTTCTCGGGGCAt                <  1:269167/149‑1 (MQ=255)
                                                                                                                                       cATGATTAACCAGGCCTTTGAAGCTGCCGCAAAAACTGAAGCGCACTCGCCAATGCTGCCGATGTCTGTCTGTGTTGCATTAAGCCAGGGTTATATTGGTTACGATCTACAAAACGCTTTAAGAGAAGAACTGCTTTCTCGGGGCATta              <  5:17746/149‑1 (MQ=255)
                                                                                                                                        aTGATTAACCAGGCCTTTGAAGCTGCCGCAAAAACTGAAGCGCACTCGCCAATGCTGCCGATGTCTGTCTGTGTTGCATTAAGCCAGGGTTATATTGGATACGATCTACAAAACGCTTTAAGAGAAGAACTGCTTTCTCGGGGCATtaa             <  2:347158/149‑1 (MQ=255)
                                                                                                                                         tGATTAACCAGGCCTTTGAAGCTGCCGCAAAAACTGAAGCGCACTCGCCAATGCTGCCGATGTCTGTCTGTGTTGCATTAAGCCAGGGTTATATTGGTTACGATCTACAACACGCTTTaa                                         >  1:278629/1‑120 (MQ=255)
                                                                                                                                         tGATTAACCAGGCCTTTGAAGCTGCCGCAAAAACTGAAGCGCACTCGCCAATGCTGCCGATGTCTGTCTGTGTTGCATTAAGCCAGGGTTATATTGGTTACGATCTACAAAACGCTTTaa                                         <  2:278629/120‑1 (MQ=255)
                                                                                                                                         tGATTAACCAGGCCTTTGAAGCTGCCGCAAAAACTGAAGCGCACTCGCCAATGCTGCCGATGTCTGTCTGTGTTGCATTAAGCCAGGGTTATATTGGTTACGATCTACAAAACGCTTTAAGAGAAGAACTGCTTTCTCGGGGCATtaa             >  3:147140/1‑148 (MQ=255)
                                                                                                                                           aTTAACCAGGCCTTTGAAGCTGCCGCAAAAACTGAAGCGCACTCGCCAATGCTGCCGATGTCTGTCTGTGTTGCATTAAGCCAGGGTTATATTGGTTACGATCTACAAAACGCTTTAAGAGAAGAACTGCTTTCTCGGGGCATtaataa          <  8:143924/149‑1 (MQ=255)
                                                                                                                                              aaCCAGGCCTTTGAAGCTGCCGCAAAAACTGAAGCGCACTCGCCAATGCTGCCGATGGCTGTCTGTGTTGCATTAAGCCAGGGTTATATTGGTTACGATCTACAAAACGCTTTAAGAGACGAACTGCTTTCTCGGGGCATTAATAAAc        <  2:152865/148‑1 (MQ=255)
                                                                                                                                               aCCAGGCCTTTGAAGCTGCCGCAAAAACTGAAGCGCACTCGCCAATGCTGCCGATGTCTGTCTGTGTTGCATTAAGCCAGGGTTATATTGGTTACGATCTACAAAACGCTTTAAGAGAAGAACTGCTTTCTCGGGGCATTAATAAACCt      >  2:265726/1‑149 (MQ=255)
                                                                                                                                                  aGGCCTTTGAAGCTGCCGCAAAAACTGAAGCGCACTCGCCAATGCTGCCGATGTGTGTCTGTGTTGCATTAAGCCAGGGTTATATTGGTTACGATCTACAAAACGCTTTAAGAGAAGAACTGCTTTCTCGGGGCATTAATAAACCTg     >  6:347979/1‑147 (MQ=255)
                                                                                                                                                  aGGCCTTTGAAGCTGCCGCAAAAACTGAAGCGCACTCGCCAATGCTGCCGATGTCTGTCTGTGTTGCATTAAGCCAGGGTTATATTGGTTACGATCTACAAAACGCTTTAAGAGAAGAACTGCTTTCTCGGGGCATTAATAAACCTGTa   <  1:239200/149‑1 (MQ=255)
                                                                                                                                                  aGGCCTTTGAAGCTGCCGCAAAAACTGAAGCGCACTCGCCAATGCTGCCGATGTCTGTCTGTGTTGCATTAAGCCAGGGTTATATTGGTTACGATCTACAAAACGCTTTAAGAGAAGAACTGCTTTCTCGGGGCATTAATAAACCTGTa   <  1:118040/149‑1 (MQ=255)
                                                                                                                                                   ggCCTTTGAAGCTGCCGCAAAAACTGAAGCGCACTCGCCAATGCTGCCGATGTCTGTCTGTGTTGCATTAAGCCAGGGTTATATTGGTTAAGATCTACAAAACGCTTTAAGAGAAGAACTGCTTTCTCGGGGCATTAATAAACCTGTAg  >  8:181762/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CCTGGGCGGGAATGCGCTGGGCGATGACCTGGCCGGGCAAATGAAAGCGGTAAAAATTACTTCTCAGGCAATTGTTGATTTAATTGCCCAGGGACATGAAGTTATCGTCACTCATGGTAACGGGCCGCAGGTGGGCATGATTAACCAGGCCTTTGAAGCTGCCGCAAAAACTGAAGCGCACTCGCCAATGCTGCCGATGTCTGTCTGTGTTGCATTAAGCCAGGGTTATATTGGTTACGATCTACAAAACGCTTTAAGAGAAGAACTGCTTTCTCGGGGCATTAATAAACCTGTAG  >  NC_000913/3011484‑3011779

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: