Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,900,710 Δ1 bp 7.5% coding (1261/4026 nt) 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,900,7100A. 7.5% 246.3 / 14.1 80coding (1261/4026 nt)PP_1703assimilatory nitrate reductase/sulfite reductase
Reads supporting (aligned to +/- strand):  ref base A (35/39);  new base . (4/2);  total (39/41)
Fisher's exact test for biased strand distribution p-value = 4.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.57e-01

ACCCGGCAGCCTGCGCAGCGAGGGGGGCGAGGGTACGGTGGCGCTGTTCGGGCAGATGAAAGCCGGTGAAGTGAAAGCCTGCTGGGTCATCTGCAGCAACCCGGTGGCCAGTGTCGCCAACCGTCAGCAAGTGATCGACGGCCTGCGCAAGGCCGAACTGGTCATTACCCAGGACGCCTTCCTCGACACCGAAACCAACCGTTACGCCGACATCCTGCTGCCGGCCGCGCTGTGGGCCGAAGGCGAGGGTGTGATGATCAACAGCGAGCGCAACCTGACCCTGATGCCACGCGCTGT  >  NC_002947/1900562‑1900858
                                                                                                                                                       |                                                                                                                                                    
aCCCGGCAGCCTGCGCAGCGAGGGGGGCGAGGGTACGGTGGCGCTGTTCGGGCAGATGAAAGCCGGTGAAGTGAAAGCCTGCTGGGTCATCTGCAGCAACCCGGTGGCCAGTGTCGCCAACCGTCAGCAAGTGATCGACGGCCTGCGCa                                                                                                                                                      <  1:94528/149‑1 (MQ=255)
   cGGCAGCCTGCGCAGCGAGGGGGGCGAGGGTACGGTGGCGCTGTTCGGGCAGATGAAAGCCGGTGAAGTGAAAGCCTGCTGGGTCATCTGCAGCAACCCGGTGGCCAGTGTCGCCAACCGTCAGCAAGTGATCGACGGCCTGCGCAAgg                                                                                                                                                   <  8:224594/149‑1 (MQ=255)
          cTGCGCAGCGAGGGGGGCGAGGGTACGGTGGCGCTGTTCGGGCAGATGAAAGCCGGTGAAGTGAAAGCCTGCTGGGTCATCTGCAGCAACCCGGTGGCCAGTGTCGCCAACCGTCAGCAAGTGATCGACGGCCTTGCGCAGGCCgt                                                                                                                                               <  5:332076/146‑2 (MQ=255)
          cTGCGCAGCGAGGGGGGCGAGGGTACGGTGGCGCTGTTCGGGCAGATGAAAGCCGGTGAAGTGAAAGCCTGCTGGGTCATCTGCAGCAACCCGGTGGCCAGTGTCGCCAACCGTCAGCAAGTGATCGACGGCCTTGCGCAGGCCgt                                                                                                                                               >  6:332076/1‑145 (MQ=255)
             cgcAGCGAGGGGGGCGAGGGTACGGTGGCGCTGTTCGGGCAGATGAAAGCCGGTGAAGTGAAAGCCTGCTGGGTCATCTGCAGCAACCCGGTGGCCAGTGTCGCCAACCGTCAGCAAGTGATCGACGGCCTTGCGCAGGCCgtcgat                                                                                                                                           <  6:216377/147‑6 (MQ=255)
              gcAGCGAGGGGGGCGAGGGTACGGTGGCGCTGTTCGGGCAGATGAAAGCCGGTGAAGTGAAAGCCTGCTGGGTCATCTGCAGCAACCCGGTGGCCAGTGTCGCCAACCGTCAGCAAGTGATCGACGGCCTGCGCAAGGCCGAACTGGTc                                                                                                                                        <  4:164938/149‑1 (MQ=255)
               cAGCGAGGGGGGCGAGGGTACGGTGGCGCTGTTCGGGCAGATGAAAGCCGGTGAAGTGAAAGCCTGCTGGGTCATCTGCAGCAACCCGGTGGCCAGTGTCGCCAACCGTCAGCAAGTGATCGACGGCCTGCGCAAGGCCGAACTGGTCa                                                                                                                                       <  7:390163/149‑1 (MQ=255)
                    aGGGGGGCGAGGGTACGGTGGCGCTGTTCGGGCAGATGAAAGCCGGTGAAGTGTAAGCCTGCTGGGTCATCTGCAGCAACCCGGTGGCCAGTGTCGCCAACCGTCAGCAAGTGATCGACGGCCTGCGCAAGGCCGAACTGGTCATTAcc                                                                                                                                  <  1:356311/149‑1 (MQ=255)
                     ggggggCGAGGGTACGGTGGCGCTGTTCGGGCAGATGAAAGCCGGTGAAGTGAAAGCCTGCTGGGTCATCTGCAGCAACCCGGTGGCCAGTGTCGCCAACCGTCAGCAAGTGATCGACGGCCTGCGCAAGGCCGAACTGGTCATTAccc                                                                                                                                 <  4:348965/149‑1 (MQ=255)
                     ggggggCGAGGGTACGGTGGCGCTGTTCGGGCAGATGAAAGCCGGTGAAGTGAAAGCCTGCTGGGTCATCTGCAGCAACCCGGTGGCCAGTGTCGCCAACCGTCAGCAAGTGATCGACGGCCTGCGCAAGGCCGAACTGGTCATTAcc                                                                                                                                  >  2:88666/1‑148 (MQ=255)
                       ggggCGAGGGTACGGTGGCGCTGTTCGGGCAGATGAAAGCCGGTGAAGTGAAAGCCTGCTGGGTCATCTGCAGCAACCCGGTGGCCAGTGTCGCCAATCCGTCAGCAAGTGATCGACGGCCTGCGCAAGGCCGAACTGGTCATTACCCa                                                                                                                                <  1:107427/149‑1 (MQ=255)
                          gCGAGGGTACGGTGTCGCTGTTCGGGCAGATGAAAGCCGGTGAAGTGAAAGCCTGCTGGGTCATCTGCAGCAACCCGGTGGCCAGTGTCGCCAACCGTCAGCAAGTGATCGACGGCCTGCGCAAGGCCGAACTGGTCATTACCCAGGAc                                                                                                                            >  6:11940/1‑149 (MQ=255)
                          gCGAGGGTACGGTGGCGCTGTTCGGGCAGATGAAAGCCGGTGAAGTGAAAGCCTGCTGGGTCATCTGCAGCAACCCGGTGGCCAGTGTCGCCAACCGTCAGCAAGTGATCGACGGCCTGCGCAAGGCCGAACTGGTCATTACCCAGGAc                                                                                                                            <  3:36806/149‑1 (MQ=255)
                           cGAGGGTACGGTGGCGCTGTTCGGGCAGATGAAAGCCGGTGAAGTGAAAGCCTGCTGGGTCATCTGCAGCAACCCGGTGGCCAGTGTCGCCAACCGTCAGCAAGTGATCGACGGCCTGCGCAAGGCCGAACTGGTCATTACCCAGGACg                                                                                                                           >  6:4288/1‑149 (MQ=255)
                                gTACGGTGGCGCTGTTCGGGCAGATGAAAGCCGGTGAAGTGAAAGCCTGCTGGGTCATCTGCAGCAACCCGGTGGCCAGTGTCGCCAACCGTCAGCAAGTGATCGACGGCCTGCGCAAGGCCGAACTGGTCATTACCCAGGACGCCTTc                                                                                                                      >  2:255670/1‑149 (MQ=255)
                                 tACGGTGGCGCTGTTCGGGCAGATGAAAGCCGGTGAAGTGAAAGCCTGCTGGGTCATCTGCAGCAACCCGGTGGCCAGTGTCGCCAACCGTCAGCAAGTGATCGACGGCCTGCGCAAGGCCGAACTGGTCATTACCCAGGACGCCTTcc                                                                                                                     >  2:57897/1‑149 (MQ=255)
                                   cGGTGGCGCTGTTCGGGCAGATGAAAGCCGGTGAAGTGAAAGCCTGCTGGGTCATCTGCAGCAACCCGGTGGCCAGTGTCGCCAACCGTCAGCAAGTGATCGACGGCCTGCGCAAGGCCGAACTGGTCATTACCCAGGACGCCTTCCt                                                                                                                    <  2:20609/148‑1 (MQ=255)
                                   cGGGGGCGCTGTTCGGGCAGATGAAAGCCGGTGAAGTGAAAGCCTGCTGGGTCATCTGCAGCAACCCGGTGGCCAGTGTCGCCAACCGTCAGCAAGTGATCGACGGCCTGCGCAAGGCCGAACTGGTCATTACCCAGGACGCCTTCCt                                                                                                                    <  8:107825/148‑1 (MQ=255)
                                     gTGGCGCTGTTCGGGCAGATGAAAGCCGGTGAAGTGAAAGCCTGCTGGGTCATCTGCAGCAACCCGGTGGCCAGTGTCGCCAACCGTCAGCAAGTGATCGACGGCCTGCGCAAGGCCGAACTGGTCATTACCCAGGACGCCTTCCTCGa                                                                                                                 <  5:8877/149‑1 (MQ=255)
                                      tGGCGCTGTTCGGGCAGATGAAAGCCGGTGAAGTGAAAGCCTGCTGGGTCATCTGCAGCAACCCGGTGGCCAGTGTCGCCAACCGTCAGCAAGTGATCGACGGCCTGCGCAAGGCCGAACTGGTCATTACCCAGGACGCCTTCCTCGac                                                                                                                <  4:297414/149‑1 (MQ=255)
                                      tGGCGCTGTTCGGGCAGATGAAAGCCGGTGAAGTGAAAGCCTGCTGGGTCATCTGCAGCAACCCGGTGGCCAGTGTCGCCAACCGTCAGCAAGTGATCGACGGCCTGCGCAAGGCCGAACTGGTCATTACCCAGGACGCCATCCTCGac                                                                                                                <  5:4288/149‑1 (MQ=255)
                                        gcgcTGTTCGGGCAGATGAAAGCCGGTGAAGTGAAAGCCTGCTGGGTCATCTGCAGCAACCCGGTGGCCAGTGTCGCCAACCGTCAGCAAGTGATCGACGGCCTGCGCAAGGCCGAACTGGTCATTACCCAGGACGCCTTCCTCGaca                                                                                                               <  8:347100/148‑1 (MQ=255)
                                        gcgcTGTTCGGGCAGATGAAAGCCGGTGAAGTGAAAGCCTGCTGGGTCATCTGCAGCAACCCGGTGGCCAGTGTCGCCAACCGTCAGCAAGTGATCGACGGCCTGCGCAAGGCCGAACTGGTCATTACCCAGGACGCCTTCCTCGaca                                                                                                               <  7:294300/148‑1 (MQ=255)
                                         cgcTGTTCGGGCAGATGAAAGCCGGTGAAGTGAAAGCCTGCTGGGTCATCTGCAGCAACCCGGTGGCCAGTGTCGCCAACCGTCAGCAAGTGATCGACGGCCTGCGCAAGGCCGAACTGGTCATTACCCATGACGCCTTCCTCGACAcc                                                                                                             <  6:360242/149‑1 (MQ=255)
                                         cgcTGTTCGGGCAGATGAAAGCCGGTGAAGTGAAAGCCTGCTGGGTCATCTGCACCAACCCGGTGGCCAGTGTCGCCAACCGTCAGCAAGTGATCGACGGCCTGCGCAAGGCCGAACTGGTCATTACCCAGGACGCCTTCCTCGACAcc                                                                                                             <  2:336973/149‑1 (MQ=255)
                                          gcTGTTCGGGCAGATGAAAGCCGGTGAAGTGAAAGCCTGCTGGGTCATCTGCAGCAACCCGGTGGCCAGTGTCGCCAACCGTCAGCAAGTGATCGACGGCCTGCGCAAGGCCGAACTGGTCATTACCCAGGACGCCTTCCTCGACAcc                                                                                                             >  6:352100/1‑148 (MQ=255)
                                              ttCGGGCAGATGAAAGCCGGTGAAGTGAAAGCCTGCTGGGTCATCTGCAGCAACCCGGTGGCCAGTGTCGCCAACCGTCAGCAAGTGATCGACAGGCCTGCGCAAGGCCGAACTGGTCATTACCCAGGACGCCTTCCTCGACACCGaaa                                                                                                         <  1:119450/149‑1 (MQ=255)
                                                 gggCAGATGAAAGCCGGTGAAGTGAAAGCCTGCTGGGTCATCTGCAGCAACCCGGTGGCCAGTGTCGCCAACCGTCAGCAAGTGATCGACGGCCTGCGCAAGGCCGAACTGGTCATTACCCAGGACGCCTTCCTCGACACCGAaaccaa                                                                                                     >  8:370517/1‑149 (MQ=255)
                                                 gggCAGATGAAAGCCGGTGAAGTGAAAGCCTGCTGGGTCATCTGCAGCAACCCGGTGGCCAGTGTCGCCAACCGTCAGCAAGTGATCGACGGCCTGCGCAAGGCCGAACTGGTCATTACCCAGGACGCCTTCCTCGACACCGAaaccaa                                                                                                     >  5:143859/1‑149 (MQ=255)
                                                   gCAGATGAAAGCCGGTGAAGTGAAAGCCTGCTGGGTCATCTGCAGCAACCCGGTGGCCAGTGTCGCCAACCGTCAGCAAGTGATCGACGGCCTGCGCAAGGCCGAACTGGTCATTACCCAGGACGCCTTCCTCGaca                                                                                                               >  5:408169/1‑137 (MQ=255)
                                                   gCAGATGAAAGCCGGTGAAGTGAAAGCCTGCTGGGTCATCTGCAGCAACCCGGTGGCCAGTGTCGCCAACCGTCAGCAAGTGATCGACGGCCTGCGCAAGGCCGAACTGGTCATTACCCAGGACGCCTTCCTCGaca                                                                                                               <  6:408169/137‑1 (MQ=255)
                                                      gATGAAAGCCGGTGAAGTGAAAGCCTGCTGGGTCATCTGCAGCAACCCGGTGGCCAGTGTCGCCAACCGTCAGCAAGTGATCGACGGCCTGCGCAAGGCCGAACTGGTCATTACCCAGGACGCCTTCCTCGACACCGAAACCAACCGtt                                                                                                <  5:114483/149‑1 (MQ=255)
                                                          aaaGCCGGTGAAGTGAAAGCCTGCTGGGTCATCTGCAGCAACCCGGTGGCCAGTGTCGCCAACCGTCAGCAAGTGATCGACGGCCTGCGCAAGGCCGAACTGGTCATTACCCAGGACGCCTTCCTCGACACCGAAACCAACCGTTACGc                                                                                            >  6:272853/1‑149 (MQ=255)
                                                           aaGCCGGTGAAGTGAAAGCCTGCTGGGTCATCTGCAGCAACCCGGTGGCCAGTGTCGCCAACCGTCAGCAAGTGATCGACGGCCTGCGCAAGGCCGAACTGGTCATTACCCAGGACGCCTTCCTCGACACCGAAACCAACCGTTACGcc                                                                                           >  5:67299/1‑149 (MQ=255)
                                                               cGGTGAAGTGAAAGCCTGCTGGGTCATCTGCAGCAACCCGGTGGCCAGTGTCGCCAACCGTCAGCAAGTGATCGACGGCCTGCGCAAGGCCGAACTGGTCATTACCCAGGACGCCTTCCTCGACACCGAAACCAACCGTTACGCCGACa                                                                                       >  5:109452/1‑149 (MQ=255)
                                                                 gtgaagtgaaAGCCTGCTGGGTCATCTGCAGCAACCCGGTGGCCAGTGTCGCCAACCGTCAGCAAGTGATCGACCGACTGCGCAAGGCCGAACTGGTCATTACCCAGGACGCCTTCCTCGACACCGAAACCAACCGTAAACCCGACaa                                                                                      >  8:156105/1‑147 (MQ=255)
                                                                 gtgaagtgaaAGCCTGCTGGGTCATCTGAAGCAACCCGTTGTCCAGTTTCGCCAACCGTCAGCAAGTGATCTACGGCCTTCGCAAGGCCGAACTGGTCATTACCCAGTACGGCTTCCTCGACACCGAAACCAACCGTTACGCCGACAt                                                                                      <  7:156105/148‑1 (MQ=255)
                                                                    aagtgaaAGCCTGCTGGGTCATCTGCAGCAACCCGGTGGCCAGTGTCGCCAACCGTCAGCAAGTGATCGACGGCCTGCGCAAGGCCGAACTGGTCATTACCCAGGACGCCTTCCTCGACACCGAAACCAACCGTTACGCCGACATCctg                                                                                  >  2:177049/1‑149 (MQ=255)
                                                                       tgaaAGCCTGCTGGGTCATCTGCAGCAACCCGGTGGCCAGTGTCGCCAACCGTCAGCAAGTGATCGACGGCCTGCGCAAGGCCGAACTGGTCATTACCCAGGACGCCTTCCTCGACACCGAACCGACCCGGTACGCCGCCATCctgctg                                                                               >  6:258907/1‑149 (MQ=255)
                                                                       tgaaAGCCTGCTGGGTCATCTGCAGCAACCCGGTGGCCAGTGTCGCCAACCGTCAGCAAGTGATCGACGGCCTGCGCAAGGCCGAACTGGTCATTACCCAGGACGCCTTCCTCGACACCGAAACCAACCGTTACGCCGACATCctgctg                                                                               <  6:143859/149‑1 (MQ=255)
                                                                                      tCATCTGCAGCAACCCGGTGGCCAGTGTCGCCAACCGTCAGCAAGTGATCGACGGCCTGCGCAAGGCCGAACTGGTCATTACCCAGGACGCCTTCCTCGACACCGAAACCAACCGTTACGCCGACATCCTGCTGCCGGCCGCGCTGTgg                                                                >  2:376784/1‑149 (MQ=255)
                                                                                         tCTGCAGCAACCCGGTGGCCAGTGTCGCTAACCGTCAGCAAGTGATCGACGGCCTGCGCAAGGCCGAACTGGTCATTACCCAGGACGCCTTCCTCGACACCGAAACCAACAGTTACGCCGa                                                                                         >  2:30965/1‑121 (MQ=255)
                                                                                         tCTGCAGCAACCCGGTGGCCAGTGTCGCTAACCGTCAGCAAGTGATCGACGGCCTGCGCAAGGCCGAACTGGTCATTACCCAGGACGCCTTCCTCGACACCGAAACCAACAGTTACGCCGa                                                                                         <  1:30965/121‑1 (MQ=255)
                                                                                             cagcaACCCGGTGGCCAGTGTCGCCAACCGTCAGCAAGTGATCGACGGCCTGCGCAAGGCCGAACTGGTCATTACCCAGGACGCCTTCCTCGACACCGAAACCAACCGTTACGCCGACATCCTGCTGCCGGCCGCGCTGTGGGCCGaa                                                          <  5:272853/148‑1 (MQ=255)
                                                                                              agcaACCCGGTGGCCAGTGTCGCCAACCGTCAGCAAGTGATCGACGGCCTGCGCAAGGCCGAACTGGTCATTACCCAGGACGCCTTCCTCGACACCGAAACCAACCGTTACGCCGACATCCTGCTGCCGGCCGCGCTGTGGGCCGAAgg                                                        <  7:59461/149‑1 (MQ=255)
                                                                                                caACCCGGTGGCCAGTGTCGCCAACCGTCAGCAAGTGATCGACGGCCTGCGCAAGGCCGAACTGGTCATTACCCAGGACGCCTTCCTCGACACCGAAACCAACCGTTACGCCGACATCCTGCTGCCGGCCGCGCTGTGGGCCGAAGGc                                                       >  5:101167/1‑148 (MQ=255)
                                                                                                 aaCCCGGTGGCCAGTGTCGCCAACCGTCAGCAAGTGATCGACGGCCTGCGCAAGGCCGAACTGGTCATTACCCAGGACGCCTTCCTCGACACCGAAACCAACCGTTACGCCGACATCCTGCTGCCGGCCGCGCTGTGGGCCGAAGGCGa                                                     <  1:177049/149‑1 (MQ=255)
                                                                                                   cccGGTGGCCAGTGTCGCCAACCGTCAGCAAGTGATCGACGGCCTGCGCAAGGCCGAACTGGTCATTACCCAGGACGCCTTCCTCGACACCGAAACCAACCGTTACGCCGACATCCTGCTGCCGGCCGCGCTGTGGGCCGAAGGCGAgg                                                   <  3:189757/149‑1 (MQ=255)
                                                                                                     cGGTGGCCAGTGTCGCCAACCGTCAGCAAGTGATCGACGGCCTGCGCAAGGCCGAACTGGTCATTACCCAGGACGCCTTCCTCGACACCGAAACCAACCGTTACGCCGACATCCTGCTGCCGGCCGCGCTGTGGGCCGAAGGCGAGGgt                                                 >  6:301696/1‑149 (MQ=255)
                                                                                                        tGGCCAGTGTCGCCAACCGTCAGCAAGTGATCGACGGCCTGCGCAAGGCCGAACTGGTCATTACCCAGGACGCCTTCCTCGACACCGAAACCAACCGTTACGCCGACATCCTGCTGCCGGCCGCGCTGTGGGCCGAAGGCGAGGGtgtg                                              <  1:255670/149‑1 (MQ=255)
                                                                                                          gCCAGTGTCGCCAACCGTCAGCAAGTGATCGACGGCCTGCGCAAGGCCGAACTGGTCATTACCCAGGACGCCTTCCTCGACACCGAAACCAACCGTTACGCCGACATCCTGCTGCCGGCCGCGCTGTGGGCCGAAGGCGAGGGTGtgat                                            <  1:88666/149‑1 (MQ=255)
                                                                                                            cAGTGTCGCCAACCGTCAGCAAGTGATCGACGGCCTGCGCAAGGCCGAACTGGTCATTACCCAGGACGCCTTCCTCGACACCGAAACCAACCGTTACGCCGACATCCTGCTGCCGGCCGCGCTGTGGGCCGAAGGCGAGGGTGtgatga                                          >  2:134149/1‑149 (MQ=255)
                                                                                                            cAGTGTCGCCAACCGTCAGCAAGTGATCGACGGCCTGCGCAAGGCCGAACTGGTCATTACCCAGGACGCCTTCCTCGACACCGAAACCAACCGTTACGCCGACATCCTGCTGCCGGCCGCGCTGTGGGCCGAAGGCGAGGGTGtgat                                            <  7:40159/147‑1 (MQ=255)
                                                                                                            cAGTGTCGCCAACCGTCAGCAAGTGATCGACGGCCTGCGCAAGGCCGAACTGGTCATTACCCAGGACGCCTTCCTCGACACCGAAACCAACCGTTACGCCGACATCCTGCTGCCGGCCGCGCTGTGGGCCGAAGGCGAGGGTGtgat                                            >  8:40159/1‑147 (MQ=255)
                                                                                                             aGTGTCGCCAACCGTCAGCAAGTGATCGACGGCCTGCGCAAGGCCGAACTGGTCATTACCCAGGACGCCTTCCTCGACACCGAAACCAACCGTTACGCCGACATCCTGCTGCCGGCCGCGCTGTGGGCCGAAGGCGAGGGTGTgatgat                                         <  5:301696/149‑1 (MQ=255)
                                                                                                             aGTGTCGCCAACCGTCAGCAAGTGATCGACGGCCTGCGCAAGGCCGAACTGGTCATTACCCAGGACGCCTTCCTCGACACCGAAACCAACCGTTACGCCGACATCCTGCTGCCGGCCGCGCTGTGGGCCGAAGGCGAGGGTGTgatgat                                         <  5:11940/149‑1 (MQ=255)
                                                                                                               tgtCGCCAACCGTCAGCAAGTGATCGACGGCCTGCGCAAGGCCGAACTGGTCATTACCCAGGACGCCTTCCTCGACACCGAAACCAACCGTTACGCCGACATCCTGCTGCCGGCCGCGCTGTGGGCCGAAGGCGAGGGTGTGATGATc                                        >  7:74334/1‑148 (MQ=255)
                                                                                                                  cGCCAACCGTCAGCAAGTGATCGACGGCCTGCGCAAGGCCGAACTGGTCATTACCCAGGACGCCTTCCTCGACACCGAAACCAACCGTTACGCCGACATCCTGCTGCCGGCCGCGCTGTGGGCCGAAGGCGAGGGTGTGATGATCaa                                      >  7:213354/1‑147 (MQ=255)
                                                                                                                     cAACCGTCAGCAAGTGATCGACGGCCTGCGCAAGGCCGAACTGGTCATTACCCAGGACGCCTTCCTCGACACCGAAACCAACCGTTACGCCGACATCCTGCTGCCGGCCGCGCTGTGGGCCGAAGGCGAGGGTGTGATGATCAACagc                                  >  5:48192/1‑148 (MQ=255)
                                                                                                                      aaCCGTCAGCAAGTGATCGACGGCCTGCGCAAGGCCGAACTGGTCATTACCCAGGACGCCTTCCTCGACACCGAAACCAACCGTTACGCCGACATCCTGCTGCCGGCCGCGCTGTGGGCCGAAGGCGAGGGTGTGATGATCAACagc                                  >  4:294795/1‑147 (MQ=255)
                                                                                                                          cGTCAGCAAGTGATCGACGGCCTGCGCAAGGCCGAACTGGTCATTACCCAGGACGCCTTCCTCGACACCGAAACCAACCGTTACGCCGACATCCTGCTGCCGGCCGCGCTGTGGGCCGAAGGCGAGGGTGTGATGATCAACagcgagcg                             >  7:360561/1‑149 (MQ=255)
                                                                                                                           gTCAGCAAGTGATCGACGGCCTGCGCAAGGCCGAACTGGTCATTACCCAGGACGCCTTCCTCGACACCGAAACCAACCGTTACGCCGACATCCTGCTGCCGGCCGCGCTGTGGGCCGAAGGCGAGGGTGTGATGATCAACagcgagcg                             >  8:294603/1‑148 (MQ=255)
                                                                                                                                cAAGTGATCGACGGCCTGCGCAAGGCCGAACTGGTCATTACCCAGGACGCCTTCCTCGACACCGAAACCAACCGTTACGCCGACATCCTGCTGCCGGCCGCGCTGTGGGCCGAAGGCGAGGGTGTGATGATCAACAGCGAGCGCAACCt                       <  6:109452/149‑1 (MQ=255)
                                                                                                                                  aGTGATCGACGGCCTGCGCAAGGCCGAACTGGTCATTACCCAGGACGCCTTCCTCGACACCGAAACCAACCGTTACGCCGACATCCTGCTGCCGGCCGCGCTGTGGGCCGAAGGCGAGGGTGTGATGATCAACa                                    >  1:235787/1‑134 (MQ=255)
                                                                                                                                  aGTGATCGACGGCCTGCGCAAGGCCGAACTGGTCATTACCCAGGACGCCTTCCTCGACACCGAAACCAACCGTTACGCCGACATCCTGCTGCCGGCCGCGCTGTGGGCCGAAGGCGAGGGTGTGATGATCAACa                                    <  2:235787/134‑1 (MQ=255)
                                                                                                                                  aGTGATCGACGGCCTGCGCAAGGCCGAACTGGTCATTACCCAGGACGCCTTCCTCGACACCGAAACCAACCGTTACGCCGACATCCTGCTGCCGGCCGCGCTGTGGGCCGAAGGCGAGGGTGTGATGATCAACAGCGAGCGCAACCTGa                     <  7:294603/149‑1 (MQ=255)
                                                                                                                                   gTGATCGACGGCCTGCGCAAGGCCGAACTGGTCATTACCCAGGACGCCTTCCTCGACACCGAAACCAACCGTTACGCCGACATCCTGCTGCCGGCCGCGCTGTGGGCCGAAGGCGAGGGTGTGATGATCAACAGCGAGCGCAACCTGAc                    <  5:352100/149‑1 (MQ=255)
                                                                                                                                    tGATCGACGGCCTGCGCAAGGCCGAACTGGTCATTACCCAGGACGCCTTCCTCGACACCGAAACCAACCGTTACGCCGACATCCTGCTGCCGGCCGCGCTGTGGGCCGAAGGCGAGGGTGTGATGATCAACAGCGAGCGCAACCTGAcc                   <  6:101167/149‑1 (MQ=255)
                                                                                                                                    tGATCGACGGCCTGCGCAAGGCCGAACTGGTCATTACCCAGGACGCCTTCCTCGACACCGAAACCAACCGTTACGCCGACATCCTGCTGCCGGCCGCGCTGTGGGCCGAAGGCGAGGGTGTGATGATCAACAGCGAGCGCAACCTGAcc                   <  1:57897/149‑1 (MQ=255)
                                                                                                                                    accagttcGGCCTTGCGCAGGCCGAACTGGTCATTACCCAGGACGCCTTCCTCGACACCGAAACCAACCGTTACGCCGACATCCTGCTGCCGGCCGCGCTGTGGGCCGAAGGCGAGGGTGTGATGATCAACAGCGAGCGCAACCTGAcc                   >  7:126809/8‑149 (MQ=255)
                                                                                                                                        cGACGGCCTGCGCAAGGCCGAACTGGTCATTACCCAGGACGCCTTCCTCGACACCGAAACCAACCGTTACGCCGACATCCTGCTGCCGGCCGCGCTGTGCGCCGAAGGCGAGGGTGTGATGATCAACAGCGAGCGCACCCTGAACCt                 >  4:415611/1‑147 (MQ=255)
                                                                                                                                         ttcGGCCTTGCGCAGGCCGAACTGGTCATTACCCAGGACGCCTTCCTCGACACCGAAACCAACCGTTACGCCGACATCCTGCTGCCGGCCGCGCTGTGGGCCGAAGGCGAGGGTGTGATGATCAACAGCGAGCGCAACCTGACCCTGAt              >  7:11740/3‑149 (MQ=255)
                                                                                                                                         gACGGCCTGCGAAAGGCCGAACTGTTCATTACCCAGGCCGCCTTCCTCGACACCGAAACCAACCGTTACGCCGACCTCCTGCTGCCGGCCGCGCTGTGGGCCGAAGGCGAGGGTGTGATGATCAACAGCGAGCGCAACCTGACCCTGAt              >  3:119640/1‑149 (MQ=255)
                                                                                                                                          aCGGCCTGCGCAAGGCCGAACTGGTCATTACCCAGGACGCCTTCCTCGACACCGAAACCAACCGTTACGCCGACATCCTGCTGCCGGCCGCGCTGTGGGCCGAAGGCGAGGGTGTGATGATCAACAGCGAGCGCAACCTGACCCTGATg             >  8:212926/1‑149 (MQ=255)
                                                                                                                                          aCGGCCTGCGCAAGGCCGAACTGGTCATTACCCAGGACGCCTTCCTCGACACCGAAACCAACCGTTACGCCGACATCCTGCTGCCGGCCGCGCTGTGGGCCGAAGGCGAGGGTGTGATGATCAACAGCGAGCGCAACCTGACCCTGATg             >  8:249798/1‑149 (MQ=255)
                                                                                                                                           cGGCCTTGCGCAGGCCGAACTGGTCATTACCCAGGACGCCTTCCTCGACACCGAAACCAACCGTTACGCCGACATCCTGCTGCCGGCCGCGCTGTGGGCCGAAGGCGAGGGTGTGATGATCAACAGCGAGCGCAACCTGACCCTGATGc            >  8:156316/1‑149 (MQ=255)
                                                                                                                                               ccTGCGCAAGGCCGAACTGGTCATTACCCAGGACGCCTTCCTCGACACCGAAACCAACCGTTACGCCGACATCCTGCTACCGGCCGCGCTGTggg                                                               <  4:341419/95‑1 (MQ=255)
                                                                                                                                               ccTGCGCAAGGCCGAACTGGTCATTACCCAGGACGCCTTCCTCGACACCGAAACCAACCGTTACGCCGACATCCTGCTACCGGCCGCGCTGTggg                                                               >  3:341419/1‑95 (MQ=255)
                                                                                                                                                    cgcAAGGCCGAACTGGTCATTACCCAGGACGCCTTCCTCGACACCGAAACCAACCGTTACGCCGACATCCTGCTGCCGGCCGCGCTGTGGGCCGAAGGCGAGGGTGTGATGATCAACAGCGAGCGCAACCTGACCCTGATGCCACgcgc     >  5:276708/1‑149 (MQ=255)
                                                                                                                                                      cAAGGCCTAACTGGTCATTACCGAGGACGCCTTCCTCGACACCGAAACCAACCGTTACGCCGACATCCTGCTGCCGGCCGCGCTGTGGGCCGAAGGCGAGGGTGTGATGATCAACAGCGAGCTCAACCTGa                     <  2:294346/131‑1 (MQ=255)
                                                                                                                                                      cAAGGCCGAACTGGTCATTACACAGGACGCCTTCCTCGACACCGAAACCAACCGTTACGCCGACATCCTGCTGCCGGCCGCGCTGTGGGCCGAAGGCGAGGGTGTGATGATCAACAGCGAGCTCAACCTGa                     >  1:294346/1‑131 (MQ=255)
                                                                                                                                                       aaGGCCGAACTGGTCATTACCCAGGACGCCTTCCTCGACACCGAAACCAACCGTTACGCCGACATACTGCTGCCGGCCGCGCTGTGGGCCGAAGGCGAGGGTGTGATGATCAACAGCGAGCGCAACCTGACCCTGATGCCACGCGCTGt  >  6:220194/1‑149 (MQ=255)
                                                                                                                                                       |                                                                                                                                                    
ACCCGGCAGCCTGCGCAGCGAGGGGGGCGAGGGTACGGTGGCGCTGTTCGGGCAGATGAAAGCCGGTGAAGTGAAAGCCTGCTGGGTCATCTGCAGCAACCCGGTGGCCAGTGTCGCCAACCGTCAGCAAGTGATCGACGGCCTGCGCAAGGCCGAACTGGTCATTACCCAGGACGCCTTCCTCGACACCGAAACCAACCGTTACGCCGACATCCTGCTGCCGGCCGCGCTGTGGGCCGAAGGCGAGGGTGTGATGATCAACAGCGAGCGCAACCTGACCCTGATGCCACGCGCTGT  >  NC_002947/1900562‑1900858

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: