Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,073,609 G→T 5.5% L135I (CTC→ATC)  rutA ← pyrimidine oxygenase, FMN‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,073,6090GT 5.5% 251.3 / 3.1 72L135I (CTC→ATC) rutApyrimidine oxygenase, FMN‑dependent
Reads supporting (aligned to +/- strand):  ref base G (31/37);  new base T (2/2);  total (33/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.59e-01

TCGCCTTTAAAATCGCTTTTCCCCGTGCCCCACAGGTCGCGCAGCACCTGAACATATTCGGTGAGATAGTCGTAACGACGGGAGAAATAGTCATCGCCAGGCCAGATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGTCACGAGGTTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGCCAGCCCCGCCATCAAGGTGAA  >  NC_000913/1073461‑1073753
                                                                                                                                                    |                                                                                                                                                
tCGCCTTTAAAATCGCTTTTCCCCGTGCCCCACAGGTCGCGCAGCACCTGAACATATTCGGTGAGATAGTCGTAACGACGGGAGAAATAGTCATCGCCAGGCCAGATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGTCACGAg                                                                                                                                                  >  8:234846/1‑149 (MQ=255)
  gCCTTTAAAATCGCTTTTCCCCGTGCCCCACAGGTCGCGCAGCACCTGAACATATTCGGTGAGATAGTCGTAACGACGGGAGAAATAGTCATCGCCAGGCCAGATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGTCACGAGGt                                                                                                                                                >  5:281503/1‑149 (MQ=255)
      ttAAAATCGCTTTTCCCCGTGCCCCACAGGTCGCGCAGCACCTGAACATATTCGGTGAGATAGTCGTAACGACGGGAGAAATAGTCATCGCCAGGCCAGATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGTCACGAGGt                                                                                                                                                >  8:226373/1‑145 (MQ=255)
      ttAAAATCGCTTTTCCCCGTGCCCCACAGGTCGCGCAGCACCTGAACATATTCGGTGAGATAGTCGTAACGACGGGAGAAATAGTCATCGCCAGGCCAGATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGTCACGAGGt                                                                                                                                                <  7:226373/145‑1 (MQ=255)
      ttAAAATCGCTTTTCCCCGTGCCCCACAGGTCGCGCAGCACCTGAACATATTCGGTGAGATAGTCGTAACGACGGGAGAAATAGTCATCGCCAGGCCAGATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGTCACGAGGTTGAc                                                                                                                                            >  8:1614/1‑149 (MQ=255)
      ttAAAATCGCTTTTCCCCGTGCCCCACAGGTCGCGCAGCACCCGAACATATTCGGTGAGATAGTCGTAACGACGGGAGAAATAGTCATCGCCAGGCCAGATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGTCACGAGGTTGAc                                                                                                                                            >  2:332157/1‑149 (MQ=255)
         aaaTCGCTTTTCCCCGTGCCCCACAGGTCGCGCAGCACCTGAACATATTCGGTGAGATAGTCGTAACGACGGGAGAAATAGTCATCGCCAGGCCAGATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGTCACGAGGTTGACGc                                                                                                                                          <  5:183014/148‑1 (MQ=255)
              gCTTTTCCCCGTGCCCCACAGGTCGCGCAGCACCCGAACATATTCGGTGAGATAGTCGTAACGACGGGAGAAATAGTCATCGCCAGGCCAGATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGTCACGAGGTTGACGCCAAAAc                                                                                                                                    <  1:332157/149‑1 (MQ=255)
                   tCCCCGTGCCCCACAGGTCGCGCAGCACCTGAACATATTCGGTGAGATAGTCGTAACGACGGGAGAAATAGTCATCGCCAGGCCAGATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGTCACGAGGTTGACGCCAAAAcgccc                                                                                                                                <  8:380528/148‑1 (MQ=255)
                   tCCCCGTGCCCCACAGGTCGCGCAGCACCTGAACATATTCGGTGAGATAGTCGTAACGACGGGAGAAATAGTCATCGCCAGGCCAGATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGTCACGAGGTTGACGCCAAAAcgccc                                                                                                                                <  4:203780/148‑1 (MQ=255)
                      ccGTGCCCCACAGGTCGCGCAGCACCTGAACATATTCGGTGAGCTAGTCGTAACGACGGGAGAAATAGTCATCGCCGGGCCAGATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGTCACGAGGTTGACGCCAAAAcgcccgc                                                                                                                              >  5:90907/1‑147 (MQ=255)
                      ccGTGCCCCACAGGTCGCGCAGCACCTGAACATATTCGGTGAGATAGTCGTAACGACGGGAGAAATAGTCATCGCCAGGCCAGATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGTCACGAGGTTGACGCCAAAACGCCCGCCa                                                                                                                            <  4:296144/149‑1 (MQ=255)
                        gTGCCCCACAGGTCGCGCAGCACCTGAACATATTCGGTGAGATAGTCGTAACGACGGGAGAAATAGTCATCGCCAGGCCAGATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGTCACGAGGTTGACGCCAAAACGCCCGCCaga                                                                                                                          <  5:55666/149‑1 (MQ=255)
                        gTGCCCCACAGGTCGCGCAGCACCTGAACATATTCGGTGAGATAGTCGTAACGACGGGAGAAATAGTCATCGCCAGGCCAGATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGTCACGAGGTTGACGCCAAAACGCCCGCCa                                                                                                                            >  6:133584/1‑147 (MQ=255)
                           ccccACAGGTCGCGCAGCACCTGAACATATTCGGTGAGATAGTCGTAACGACGGGAGAAATAGTCATCGCCAGGCCAGATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGTCACGAGGTTGACGCCAAAACGCCCGCCAGAGAt                                                                                                                       <  3:29920/149‑1 (MQ=255)
                           ccccACAGGTCGCGCAGCACCTGAACATATTCGGTGAGATAGTCGTAACGACGGGAGAAATAGTCATCGCCAGGCCAGATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGTCACGAGGTTGACGCCAAAACGCCCGCCAGAGAt                                                                                                                       <  3:308041/149‑1 (MQ=255)
                           ccccACAGGTCGCGCAGCACCTGAACATATTCGGTGAGATAGTCGTAACGACGGGAGAAATAGTCATCGCCAGGCCAGATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGTCACGAGGTTGACGCAAAAACGCCCGCCAGAGAt                                                                                                                       <  4:153621/149‑1 (MQ=255)
                            cccACAGGTCGCGCAGCACCTGAACATATTCGGTGAGATAGTCGTAACGACGGGAGAAATAGTCATCGCCAGGCCAGATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGTCACGAGGTTGACGCCAAAACGCCCGCCAGAGATg                                                                                                                      <  8:105267/149‑1 (MQ=255)
                            cccACAGGTCGCGCAGCACCTGAACATATTCGGTGAGATAGTCGTAACGACGGGAGAAATAGTCATCGCCAGGCCAGATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGTCACGAGGTTGACGCCAAAACGCCCGCCAGAGATg                                                                                                                      <  5:143168/149‑1 (MQ=255)
                             ccACAGGTCGCGCAGCACCTGAACATATTCGGTGAGATAGTCGTAACGACGGGAGAAATAGTCATCGCCAGGCCAGATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGTCACGAGGtt                                                                                                                                               >  8:61193/1‑123 (MQ=255)
                             ccACAGGTCGCGCAGCACCTGAACATATTCGGTGAGATAGTCGTAACGACGGGAGAAATAGTCATCGCCAGGCCAGATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGTCACGAGGtt                                                                                                                                               <  7:61193/123‑1 (MQ=255)
                             ccACAGGTCGCGCAGCACCTGAACATATTCGGTGAGATAGCCGTAACGACGGGAGAAATAGTCATCGCCAGGCCAGATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGTCACGAGGTTGACGCCAAAACGCCCGCCAGAGATg                                                                                                                      >  1:309549/1‑148 (MQ=255)
                              cacaGGTCGCGCAGCACCTGAACATATTCGGTGAGATAGTCGTAACGACGGGAGAAATAGTCATCGCCAGGCCAGATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGTCACGAGGTTGACGCCAAAACGCCCGCCAGAGATg                                                                                                                      >  6:174842/1‑147 (MQ=255)
                               acaGGTCGCGCAGCACCTGAACATATTCGGTGAGATAGTCGTAACGACGGGAGAAATAGTCATCGCCAGGCCAGATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGTCACGAGGTTGACGCCAAAACGCCCGCCAGAGATGGAg                                                                                                                   >  7:365559/1‑149 (MQ=255)
                                      gcgcAGCACCTGAACATATTCGGTGAGATAGTCGTAACGACGGGAGAAATAGTCATCGCCAGGCCAGATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGTCACGAGGTTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATg                                                                                                             <  6:281503/148‑1 (MQ=255)
                                         cagcaCCTGAACATATTCGGTGAGATAGTCGTAACGCCGGGAGAAATAGTCATCGCCAGGCCAGATACCCCTCTGCTCATACTCGGGCTTTTGCCAGCCAGTCACGAGGTTGACGCCAAAACGACCGCCAGAGATGGAGTCGATGGTTg                                                                                                         >  2:68016/1‑149 (MQ=255)
                                            caCCTGAACATATTCGGTGAGATAGTCGTAACGACGGGAGAAATAGTCATCGCCAGGCCAGATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGTCACGAGGTTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCgg                                                                                                      <  7:1614/149‑1 (MQ=255)
                                              ccTGAAAATATTCGGTGAGAAAGTCGTAACGACGGGAGAAATAGTCATCGCCAGGCCAGATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGTCACGAGGTTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGcc                                                                                                    <  6:296025/149‑1 (MQ=255)
                                               cTGAACATATTCGGTGAGATAGTCGGAACGACGGGAGAAATAGTCATCGCCAGGCCAGATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGTCACGAGGTGGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCa                                                                                                   <  3:11860/149‑1 (MQ=255)
                                                 gAACATATTCGGTGAGATAGTCGTAACGACGGGAGAAATAGTCATCGCCAGGCCAGATACCCACCTGCTCATACTCGGGCTTTTGCCAGCCAGTCACGAGGTTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATa                                                                                                 >  2:42616/1‑149 (MQ=255)
                                                     atatTCGGTGAGATAGTCGTAACGACGGGAGAAATAGTCATCGCCAGGCCAGATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGTCACGAGGTTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACggg                                                                                             <  1:134013/149‑1 (MQ=255)
                                                       atTCGGTGAGATAGTCGTAACGACGGGAGAAATAGTCATCGCCAGGCCAGATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGTCACGATGTTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGc                                                                                            >  5:300892/1‑148 (MQ=255)
                                                       atTCGGTGAGATAGTCGTAACGACGGGAGAAATAGTCATCGCCAGGCCAGATACCCATCTGCTCATACTCGGGCTTTTGACAGCCAGTCACGATGTTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGc                                                                                            >  5:110706/1‑148 (MQ=255)
                                                          cGGTGAGATAGTCGTAACGACGGGAGAAATAGTCATCGCCAGGCCAGATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGTCACGAGGTTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGcgac                                                                                         <  8:365559/148‑1 (MQ=255)
                                                               agaTAGTCGTAACGACGGGAGAAATAGTCATCGCCAGGCCAGATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGTCACGAGGTTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGc                                                                                   <  6:331270/149‑1 (MQ=255)
                                                                gaTAGTCGTAACGACGGGAGAAATAGTCATCGCCAGGCCAGATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGTCACGAGGTTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATCGTTGCGGCCATACGGGCGACGATTTcc                                                                                  >  2:299018/1‑149 (MQ=255)
                                                                   aGTCGTAACGACGGGAGAAATAGTCATCGCCAGGCCAGATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGTCACGAGGTTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCg                                                                                 >  7:206945/1‑147 (MQ=255)
                                                                     tCGTAACGACGGGAGAAATAGTCATCGCCAGGCCAGATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGTCACGAGGTTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAgg                                                                             >  2:332413/1‑149 (MQ=255)
                                                                     tCGTAACGACGGGAGAAATAGTCATCGCCAGGCCAGATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGTCACGAGGTTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAgg                                                                             >  3:106383/1‑149 (MQ=255)
                                                                       gTAACGACGGGAGAAATAGTCATCGCCAGGCCAGATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGTCACGAGGTTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGt                                                                            >  8:174410/1‑148 (MQ=255)
                                                                               gggAGAAATAGTCATCGCCGGGCCAGATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGTCACGAGGTTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGGTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAg                                                                   <  6:90907/149‑1 (MQ=255)
                                                                               gggAGAAATAGTCATCGCCAGGCCAGATACCCATCTGCTCATACTCGGGCGTTTGCCAGCCACTCACGAGGATGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAg                                                                   <  6:48922/149‑1 (MQ=255)
                                                                                ggAGAAATAGTCATCGCCAGGCCAGATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGTCACGAGGTTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGtt                                                                      >  4:42079/1‑145 (MQ=255)
                                                                                ggAGAAATAGTCATCGCCAGGCCAGATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGTCACGAGGTTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGtt                                                                      <  3:42079/145‑1 (MQ=255)
                                                                                 gagaAATAGTCATCGCCAGGCCAGATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGTCACGAGGTTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAg                                                                   >  8:248048/1‑147 (MQ=255)
                                                                                 gagaAATAGTCATCGCCAGGCCAGATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGTCACGAGGTTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAg                                                                   <  7:248048/147‑1 (MQ=255)
                                                                                   gaAATAGTCATCGCCAGGCCAGATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGTCACGAGGTTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTgg                                                               <  2:92433/149‑1 (MQ=255)
                                                                                    aaaTAGTCATCGCCAGGCCAGATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGTCACGAGGTTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGgc                                                              <  1:133280/149‑1 (MQ=255)
                                                                                         gTCATCGCCAGGCCAGATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGTCACGAGGTTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCTGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGcagcag                                                         <  3:20629/149‑1 (MQ=255)
                                                                                          tCATCGCCAGGCCAGATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGTCACGAGGTTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGt                                                        >  5:319045/1‑149 (MQ=255)
                                                                                            aTCGCCAGGCCAGATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGTCACGAGGTTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTg                                                       <  4:149964/148‑1 (MQ=255)
                                                                                               gccaggccagATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGTCACGAGGTTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGt                                                   <  2:145970/149‑1 (MQ=255)
                                                                                                 caggccagATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGTCACGATGTTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTaa                                                 <  6:110706/149‑1 (MQ=255)
                                                                                                 caggccagATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGTCACGATGTTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTaa                                                 <  6:300892/149‑1 (MQ=255)
                                                                                                   ggccagATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGTCACGAGGTTGACGCCAAACCGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAAt                                               <  1:332413/149‑1 (MQ=255)
                                                                                                     ccagATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGTCACGAGGTTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCt                                             >  1:226339/1‑149 (MQ=255)
                                                                                                       agctACCCATCTGCTCATACTCGGGCTGTTGCCTGCCAGTCACGAGGTTGACGCCAATTCGGCCGCCAGTGAGGGAGTCGATGGTTGCGGCCATACGGGCGACTTTTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGa                                           <  5:136914/146‑1 (MQ=255)
                                                                                                           aCCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGTCACGAGGTTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATgc                                       >  3:4253/1‑149 (MQ=255)
                                                                                                                    cTCATACTCGGGCTTTTGCCAGCCAGTCACGAGGTTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAgg                                 <  4:106383/146‑1 (MQ=255)
                                                                                                                     tCATACTCGGGCTTTTGCCAGCCAGTCACGAGGTTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCAc                             >  8:296024/1‑149 (MQ=255)
                                                                                                                        tACTCGGGCTTTTGCCAGCCAGTCACGAGGTTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGgc                          >  6:229976/1‑149 (MQ=255)
                                                                                                                          cTCGGGCTTTTGCCAGCCAGTCACGAGGTTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGgccg                        >  1:368717/1‑149 (MQ=255)
                                                                                                                           tCGGGCTTTTGCCAGCCAGTCACGAGGTTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGa                                           >  3:138315/1‑129 (MQ=255)
                                                                                                                           tCGGGCTTTTGCCAGCCAGTCACGAGGTTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGa                                           <  4:138315/129‑1 (MQ=255)
                                                                                                                                 ttttGCCAGCCAGTCACGAGGTTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCgcccgccc                 >  5:31925/1‑149 (MQ=255)
                                                                                                                                   ttGCCAGCCAGTCACGAGGTTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGCCAGcccc                <  5:133584/148‑1 (MQ=255)
                                                                                                                                       cagccagTCACGAGGTTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGCCAGCCCCGCCa            >  2:32024/1‑148 (MQ=255)
                                                                                                                                         gccagTCACGAGGTTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGCCAGCCCCGCCATc          >  6:200856/1‑148 (MQ=255)
                                                                                                                                           cagTCACGAGGTTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGCCAGCCCCGCCATCAAg       <  2:368717/149‑1 (MQ=255)
                                                                                                                                           cagTCACGAGGTTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGCCAGCCCCGCCATCAAg       <  8:337834/149‑1 (MQ=255)
                                                                                                                                             gTCACGAGGTTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGCCAGCCCCGCCATCAAGGt     <  5:174842/149‑1 (MQ=255)
                                                                                                                                                aCGAGGTTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGCCAGCCCCGCCATCAAGGTGaa  >  7:158634/1‑149 (MQ=255)
                                                                                                                                                aCGAGGTTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGCCAGCCCCGCCATCAAGGTGaa  >  8:8238/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
TCGCCTTTAAAATCGCTTTTCCCCGTGCCCCACAGGTCGCGCAGCACCTGAACATATTCGGTGAGATAGTCGTAACGACGGGAGAAATAGTCATCGCCAGGCCAGATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGTCACGAGGTTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGCCAGCCCCGCCATCAAGGTGAA  >  NC_000913/1073461‑1073753

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