Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,900,704:1 +T 7.8% coding (1255/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,7041.T 7.8% 297.8 / 19.9 77coding (1255/4026 nt)PP_1703assimilatory nitrate reductase/sulfite reductase
Reads supporting (aligned to +/- strand):  ref base . (33/38);  new base T (4/2);  total (37/40)
Fisher's exact test for biased strand distribution p-value = 4.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.50e-01

CCACCCGGCAGCCTGCGCAGCGAGGGGGGCGAGGGTACGGTGGCGCTGTTCGGGCAGATGAAAGCCGGTGAAGTGAAAGCCTGCTGGGTCATCTGCAGCAACCCGGTGGCCAGTGTCGCCAACCGTCAGCAAGTGATCGACGGCCTGCGCAAGGCCGAACTGGTCATTACCCAGGACGCCTTCCTCGACACCGAAACCAACCGTTACGCCGACATCCTGCTGCCGGCCGCGCTGTGGGCCGAAGGCGAGGGTGTGATGATCAACAGCGAGCGCAACCTGACCCTGATGC  >  NC_002947/1900560‑1900848
                                                                                                                                                   |                                                                                                                                                
ccACCCGGCAGCCTGCGCAGCGAGGGGGGCGAGGGTACGGTGGCGCTGTTCGGGCAGATGAAAGCCGGTGAAGTGAAAGCCTGCTGGGTCATCTGCAGCAACCCGGTGGCCAGTGTCGCCAACCGTCAGCAAGTGATCGACGGCCTgcg                                                                                                                                              >  2:356311/1‑149 (MQ=255)
  aCCCGGCAGCCTGCGCAGCGAGGGGGGCGAGGGTACGGTGGCGCTGTTCGGGCAGATGAAAGCCGGTGAAGTGAAAGCCTGCTGGGTCATCTGCAGCAACCCGGTGGCCAGTGTCGCCAACCGTCAGCAAGTGATCGACGGCCTGCGCa                                                                                                                                            <  1:94528/149‑1 (MQ=255)
     cGGCAGCCTGCGCAGCGAGGGGGGCGAGGGTACGGTGGCGCTGTTCGGGCAGATGAAAGCCGGTGAAGTGAAAGCCTGCTGGGTCATCTGCAGCAACCCGGTGGCCAGTGTCGCCAACCGTCAGCAAGTGATCGACGGCCTGCGCAAgg                                                                                                                                         <  8:224594/149‑1 (MQ=255)
            cTGCGCAGCGAGGGGGGCGAGGGTACGGTGGCGCTGTTCGGGCAGATGAAAGCCGGTGAAGTGAAAGCCTGCTGGGTCATCTGCAGCAACCCGGTGGCCAGTGTCGCCAACCGTCAGCAAGTGATCGACGGCCTTGCGCAGGCCgt                                                                                                                                     >  6:332076/1‑145 (MQ=255)
            cTGCGCAGCGAGGGGGGCGAGGGTACGGTGGCGCTGTTCGGGCAGATGAAAGCCGGTGAAGTGAAAGCCTGCTGGGTCATCTGCAGCAACCCGGTGGCCAGTGTCGCCAACCGTCAGCAAGTGATCGACGGCCTTGCGCAGGCCgt                                                                                                                                     <  5:332076/146‑2 (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                                                                                                     <  7:294300/148‑1 (MQ=255)
                                          gcgcTGTTCGGGCAGATGAAAGCCGGTGAAGTGAAAGCCTGCTGGGTCATCTGCAGCAACCCGGTGGCCAGTGTCGCCAACCGTCAGCAAGTGATCGACGGCCTGCGCAAGGCCGAACTGGTCATTACCCAGGACGCCTTCCTCGaca                                                                                                     <  8:347100/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                                                                                                     <  6:408169/137‑1 (MQ=255)
                                                     gCAGATGAAAGCCGGTGAAGTGAAAGCCTGCTGGGTCATCTGCAGCAACCCGGTGGCCAGTGTCGCCAACCGTCAGCAAGTGATCGACGGCCTGCGCAAGGCCGAACTGGTCATTACCCAGGACGCCTTCCTCGaca                                                                                                     >  5:408169/1‑137 (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                                                     >  3:341419/1‑95 (MQ=255)
                                                                                                                                                 ccTGCGCAAGGCCGAACTGGTCATTACCCAGGACGCCTTCCTCGACACCGAAACCAACCGTTACGCCGACATCCTGCTACCGGCCGCGCTGTggg                                                     <  4:341419/95‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
CCACCCGGCAGCCTGCGCAGCGAGGGGGGCGAGGGTACGGTGGCGCTGTTCGGGCAGATGAAAGCCGGTGAAGTGAAAGCCTGCTGGGTCATCTGCAGCAACCCGGTGGCCAGTGTCGCCAACCGTCAGCAAGTGATCGACGGCCTGCGCAAGGCCGAACTGGTCATTACCCAGGACGCCTTCCTCGACACCGAAACCAACCGTTACGCCGACATCCTGCTGCCGGCCGCGCTGTGGGCCGAAGGCGAGGGTGTGATGATCAACAGCGAGCGCAACCTGACCCTGATGC  >  NC_002947/1900560‑1900848

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: