Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,899,807 A→T 7.9% K120* (AAG→TAG) ‡ PP_1703 → assimilatory nitrate reductase/sulfite reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,899,8070AT 7.9% 189.7 / 4.6 63K120* (AAG→TAG) ‡PP_1703assimilatory nitrate reductase/sulfite reductase
Reads supporting (aligned to +/- strand):  ref base A (28/30);  new base T (2/3);  total (30/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

AGCCAGCCGCCAGTGCGCACCAGCATGGACCAGGCCATCGCTGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTACCAGGACTTCGAGCACGCCGACGTATTTCT  >  NC_002947/1899660‑1899955
                                                                                                                                                   |                                                                                                                                                    
agccagccGCCAGTGCGCACCAGCATGGACCAGGCCATCGCTGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCcaacaa                                                                                                                                                     <  4:293378/149‑1 (MQ=255)
 gccagccGCCAGTGCGCACCAGCATGGACCAGGCCATCGCTGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAg                                                                                                                                                    >  5:212276/1‑149 (MQ=255)
   cagccGCCAGTGCGCACCAGCATGGACCAGGCCATCGCTGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGtt                                                                                                                                                  >  3:107018/1‑149 (MQ=255)
    agccGCCAGTGCGCACCAGCATGGACCAGGCCATCGCTGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTg                                                                                                                                                 <  3:175694/149‑1 (MQ=255)
    agccGCCAGTGCGCACCAGCATGGACCAGGCCATCGCTGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTg                                                                                                                                                 <  7:240236/149‑1 (MQ=255)
         ccAGTGCGCACCAGCATGGACCAGGCCATCGCTGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACTTGTTGGCCag                                                                                                                                            <  1:217635/149‑2 (MQ=255)
           aGTGCGCACCAGCATGGACCAGGCCATCGCTGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAgg                                                                                                                                          <  2:69743/149‑1 (MQ=255)
              gcgcACCAGCATGGACCAGGCCATCGCTGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTg                                                                                                                                                 <  8:30264/139‑1 (MQ=255)
              gcgcACCAGCATGGACCAGGCCATCGCTGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTg                                                                                                                                                 >  7:30264/1‑139 (MQ=255)
               cgcACCAGCATGGACCAGGCCATCGCTGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCt                                                                                                                                       <  5:298543/148‑1 (MQ=255)
                    cAGCATGGACCAGGCCATTGCTGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATcc                                                                                                                                 <  7:330022/149‑1 (MQ=255)
                    cAGCATGGACCAGGCCATCGCTGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAgg                                                                                                                                          <  1:168570/140‑1 (MQ=255)
                    cAGCATGGACCAGGCCATCGCTGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAgg                                                                                                                                          >  2:168570/1‑140 (MQ=255)
                    cAGCATGGACCAGGCCATCGCTGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATcc                                                                                                                                 >  7:308196/1‑149 (MQ=255)
                      gCATGGACCAGGCCATCGCTGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTCGCCAACAAGTTGGCCAAGGGCTTTATCCGt                                                                                                                               <  1:284535/149‑1 (MQ=255)
                      gCATGGACCAGGCCATCGCTGCCACCGCCGAGCCTCTGCTCGGCATCATCGACCAGCACGTCCCGGATGCCGTGTCGCTGTACGTGTAGGGGAAAATGTCGCTTGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGt                                                                                                                               <  5:341869/149‑1 (MQ=255)
                       cATGGACCAGGCCATCGCTGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTa                                                                                                                              <  8:270812/149‑1 (MQ=255)
                           gACCAGGCCATCGCTGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATTCGTAcgc                                                                                                                           <  8:309131/148‑1 (MQ=255)
                            aCCAGGCCATCGCTGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCt                                                                                                                                       <  6:748/135‑1 (MQ=255)
                            aCCAGGCCATCGCTGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCt                                                                                                                                       >  5:748/1‑135 (MQ=255)
                                    aTCGCTGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGa                                                                                                                 <  6:212276/149‑1 (MQ=255)
                                     tCGCTGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAg                                                                                                                <  4:57062/149‑1 (MQ=255)
                                       gCTGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGt                                                                                                               <  4:107018/148‑1 (MQ=255)
                                        cTGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTcc                                                                                                             <  4:2557/149‑1 (MQ=255)
                                           ccACCGCCGAGTGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCAAAGGGATTTATACGTACGCGCCACATCGAGTCCAAc                                                                                                          >  6:51201/1‑149 (MQ=255)
                                               cgccgAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTcgc                                                                                                      <  8:301186/149‑1 (MQ=255)
                                                             gcgGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAAc                                                                                                          >  1:21682/1‑131 (MQ=255)
                                                             gcgGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAAc                                                                                                          <  2:21682/131‑1 (MQ=255)
                                                               ggCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGgccag                                                                                      >  7:189371/1‑149 (MQ=255)
                                                                  atcatcGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGgccagcgc                                                                                   <  4:119095/149‑1 (MQ=255)
                                                                     atcGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATgg                                                                                          <  7:284914/139‑1 (MQ=255)
                                                                     atcGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATgg                                                                                          >  8:284914/1‑139 (MQ=255)
                                                                     atcGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGgccagcgccag                                                                                >  4:51290/1‑149 (MQ=255)
                                                                     atcGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGgccagcgccag                                                                                >  5:323688/1‑149 (MQ=255)
                                                                      tcGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGgccagcgccag                                                                                >  2:119195/1‑148 (MQ=255)
                                                                        gACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCcagcag                                                                             >  6:79708/1‑149 (MQ=255)
                                                                                 ggCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACaa                                                                    <  6:323688/149‑1 (MQ=255)
                                                                                              tGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCcac                                                                                                                      <  8:94910/86‑1 (MQ=255)
                                                                                              tGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCcac                                                                                                                      >  7:94910/1‑86 (MQ=255)
                                                                                              tGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGgc                                                       >  4:23389/1‑149 (MQ=255)
                                                                                              tGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGgc                                                       <  2:33960/149‑1 (MQ=255)
                                                                                                      tACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATgg                                               >  2:131522/1‑149 (MQ=255)
                                                                                                       aCGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGc                                              >  8:179686/1‑149 (MQ=255)
                                                                                                        cGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTg                                                         <  6:190814/137‑1 (MQ=255)
                                                                                                        cGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTg                                                         >  5:190814/1‑137 (MQ=255)
                                                                                                           gtCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCAc                                          >  7:324004/1‑149 (MQ=255)
                                                                                                              cGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTg                                       >  6:287178/1‑149 (MQ=255)
                                                                                                                     aTGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTa                                <  5:51201/149‑1 (MQ=255)
                                                                                                                        tCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTgg                                                                          >  2:342309/1‑104 (MQ=255)
                                                                                                                        tCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTgg                                                                          <  1:342309/104‑1 (MQ=255)
                                                                                                                        tCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTAcc                              <  1:131522/148‑1 (MQ=255)
                                                                                                                                ggCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTACCAGGACTTc                      >  7:2109/1‑148 (MQ=255)
                                                                                                                                   ccAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTACCAGGACTTCGAg                   >  4:21353/1‑148 (MQ=255)
                                                                                                                                       cccttGGCCAACTTGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTACCAGGa                          <  2:346037/133‑1 (MQ=255)
                                                                                                                                       cccttGGCCAACTTGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTACCAGGa                          >  1:346037/5‑137 (MQ=255)
                                                                                                                                         cttGGCCAACTTGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTACCAGGACTTCGa                    <  2:226598/139‑1 (MQ=255)
                                                                                                                                         cttGGCCAACTTGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTACCAGGACTTCGa                    >  1:226598/3‑141 (MQ=255)
                                                                                                                                           tGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGACGCTTGGCGCCGATGGCCCACCTGGCAGTTACCAGGACTTCGAGCACGCCGAc          >  7:13705/1‑149 (MQ=255)
                                                                                                                                            ggCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTACCAGGACTTCGAGCACGCCGACg         >  5:368204/1‑149 (MQ=255)
                                                                                                                                             gCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTACCAGGACTTCGAg                   <  8:61077/138‑1 (MQ=255)
                                                                                                                                             gCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTACCAGGACTTCGAGCACGCCGAc          <  3:23389/147‑1 (MQ=255)
                                                                                                                                             gCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAACCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTACCAGGACTTCGAg                   >  7:61077/1‑138 (MQ=255)
                                                                                                                                              ccAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTACCAGGACTTCGAGCACGCCGACGTa       >  3:301155/1‑149 (MQ=255)
                                                                                                                                              ccAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTACCAGGACTTCGAGCACGCCGACGTa       <  8:324004/149‑1 (MQ=255)
                                                                                                                                               caacaaGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATg                                                <  3:130469/107‑1 (MQ=255)
                                                                                                                                               caacaaGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATg                                                >  4:130469/1‑107 (MQ=255)
                                                                                                                                                   aaGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTACCAGGACTTCGAGCACGCCGACGTATTTCt  >  1:126693/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
AGCCAGCCGCCAGTGCGCACCAGCATGGACCAGGCCATCGCTGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTACCAGGACTTCGAGCACGCCGACGTATTTCT  >  NC_002947/1899660‑1899955

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