Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,479,813 A→G 9.8% P388P (CCA→CCG ybdR → Zn‑dependent oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,479,8130AG 9.8% 188.4 / 8.8 72P388P (CCA→CCGybdRZn‑dependent oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (30/34);  new base G (3/4);  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 = 9.52e-01

CTTCCCGAACTGCTGGAGCACATCGAAGCAGGGCGCCTGCAACCGGAACTGATCGTTACTCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTGGGCCATCGTCAGCGTGTTCCGGAGC  >  NC_002947/4479667‑4479955
                                                                                                                                                  |                                                                                                                                               
cTTCCCGAACTGCTGGAGCACATCGAAGCAGGGCGCCTGCAACCGGAACTGATCGTTACTCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAgg                                                                                                                                               >  6:82820/1‑149 (MQ=255)
cTTCCCGAACTGCTGGAGCACATCGAAGCAGGGCGCCTGCAACCGGAACTGATCGTTACTCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAgg                                                                                                                                               <  8:100454/149‑1 (MQ=255)
 ttCCCGAACTGCTGGAGCACATCGAAGCAGGGCGCCTGCAACCGGAACTGATCGTTACTCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCGGGt                                                                                                                                              >  1:416490/1‑149 (MQ=255)
            ctgGAGCACATCGAAGCAGGGCGCCTGCAACCGGAACTGATCGTTACTCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCgg                                                                                                                                   >  4:338984/1‑149 (MQ=255)
               gAGCACATCGAAGCAGGGCGCCTGCAACCGGAACTGATCGTTACTCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCGGGTGCAGCACCTGGCa                                                                                                                                 <  2:416490/148‑1 (MQ=255)
                   acaTCGAAGCAGGGCGCCTGCAACCGGAACTGATCGTTACTCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTg                                                                                                                            >  6:289265/1‑149 (MQ=255)
                   acaTCGAAGCAGGGCGCCTGCAACCGGAACTGATCGTTACTCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTg                                                                                                                            <  2:406350/149‑1 (MQ=255)
                    caTCGAAGCAGGGCGCCTGCAACCGGAACTGATCGTTACTCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCGCCCGGCACCTTgg                                                                                                                           >  2:413229/1‑149 (MQ=255)
                    caTCGAAGCAGGGCGCCTGCAACCGGAACTGATCGTTACTCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTgg                                                                                                                           <  1:419414/149‑1 (MQ=255)
                     aTCGAAGCAGGGCGCCTGCAACCGGAACTGTTCGTTACTCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCCGCACCTTggg                                                                                                                          >  1:259253/1‑149 (MQ=255)
                     aTCGAAGCAGGGCGCCTGCAACCGGAACTGATCGTTACTCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCGGGTGCAGCACCTGGCAcgagg                                                                                                                            >  6:234528/1‑143 (MQ=255)
                     aTCGAAGCAGGGCGCCTGCAACCGGAACTGATCGTTACTCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCGGGTGCAGCACCTGGCAcgagg                                                                                                                            <  5:234528/147‑5 (MQ=255)
                        gAAGCAGGGCGCCTGCAACCGGAACTGATCGTTACTCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGccc                                                                                                                       >  1:194875/1‑149 (MQ=255)
                         aaGCAGGGCGCCTGCAACCGGAACTGATCGTTACTCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCt                                                                                                                      <  5:149792/149‑1 (MQ=255)
                         aaGCAGGGCGCCTGCAACCGGAACTGATCGTTACTCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCt                                                                                                                      >  3:167078/1‑149 (MQ=255)
                          aGCAGGGCGCCTGCAACCGGAACTGATCGTTACTCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTg                                                                                                                     >  3:21106/1‑149 (MQ=255)
                               ggCGCCTGCAACCGGAACTGATCGTTACTCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCCGCATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGaccac                                                                                                                <  1:260423/149‑1 (MQ=255)
                               ggCGCCTGCAACCGGAACTGATCGTTACTCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGcccac                                                                                                                >  4:292907/1‑149 (MQ=255)
                                    cTGCAACCGGAACTGATCGTTACTCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTggg                                                                                                                          >  3:266839/1‑134 (MQ=255)
                                    cTGCAACCGGAACTGATCGTTACTCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTggg                                                                                                                          <  4:266839/134‑1 (MQ=255)
                                         aCCGGAACTGATCGTTACTCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTagag                                                                                                        >  6:222941/1‑147 (MQ=255)
                                         aCCGGAACTGATCGTTACTCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAgagga                                                                                                      <  6:177766/149‑1 (MQ=255)
                                          ccGGAACTGATCGTTACTCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAgaggag                                                                                                     >  7:176567/1‑149 (MQ=255)
                                          ccGGAACTGATCGTTACTCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAgaggag                                                                                                     <  5:149461/149‑1 (MQ=255)
                                           cGGAACTGATCGTTACTCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGaggagg                                                                                                    <  5:222941/149‑1 (MQ=255)
                                                    tCGTTACTCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGaa                                                                                                 >  2:433850/1‑142 (MQ=255)
                                                    tCGTTACTCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGaa                                                                                                 <  1:433850/142‑1 (MQ=255)
                                                      gTTACTCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGa                                                                                        <  3:338984/149‑1 (MQ=255)
                                                       ttACTCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGaa                                                                                       >  5:9417/1‑149 (MQ=255)
                                                       ttACTCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGaa                                                                                       <  1:30258/149‑1 (MQ=255)
                                                         aCTCACCGCTTAGCCCTAGCAGACGCGGCCAGTGGCTACAAAATGTTCGCTCAAAACCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCGTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAAc                                                                                      <  8:317550/148‑1 (MQ=255)
                                                         aCTCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAAc                                                                                      <  5:177721/148‑1 (MQ=255)
                                                          cTCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACg                                                                                     <  3:13551/148‑1 (MQ=255)
                                                           tCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACgg                                                                                    >  4:103706/1‑148 (MQ=255)
                                                               cGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGaa                                                                                       <  3:343676/141‑1 (MQ=255)
                                                               cGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGaa                                                                                       >  4:343676/1‑141 (MQ=255)
                                                                gCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGc                                                                              <  3:103706/149‑1 (MQ=255)
                                                                 cTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGc                                                                              >  4:352332/1‑148 (MQ=255)
                                                                  ttAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGACAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAAAGGTCTCGc                                                                              >  2:193144/1‑147 (MQ=255)
                                                                        cTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGaccac                                                                                                                >  4:5469/1‑108 (MQ=255)
                                                                        cTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGaccac                                                                                                                <  3:5469/108‑1 (MQ=255)
                                                                           gaagaGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCGCCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCt                                                                   <  1:413229/149‑1 (MQ=255)
                                                                           gaagaGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCt                                                                   <  2:259253/149‑1 (MQ=255)
                                                                           gaagaGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCt                                                                   <  2:194875/149‑1 (MQ=255)
                                                                           gaagaGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCt                                                                   <  3:382267/149‑1 (MQ=255)
                                                                            aagaGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAAGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCt                                                                   >  5:407429/1‑148 (MQ=255)
                                                                            aagaGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAAGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCt                                                                   <  6:407429/148‑1 (MQ=255)
                                                                             agagGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTgg                                                                 <  6:9417/149‑1 (MQ=255)
                                                                              gagGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGa                                                                <  8:176567/149‑1 (MQ=255)
                                                                                 gcggcCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAAATGATCATGAACGGTCTCGCCAGTTACTTCTGGAtt                                                              <  1:296311/149‑1 (MQ=255)
                                                                                 gcggcCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTg                                                             >  2:124001/1‑149 (MQ=255)
                                                                                 gcggcCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTg                                                             <  5:82820/149‑1 (MQ=255)
                                                                                         gggCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCg                                                     >  2:305333/1‑149 (MQ=255)
                                                                                          ggCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCgg                                                    <  5:289265/149‑1 (MQ=255)
                                                                                           gCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCgg                                                    >  5:171653/1‑148 (MQ=255)
                                                                                              aCAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAAt                                                >  3:409695/1‑149 (MQ=255)
                                                                                               cAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAAt                                                >  5:405092/1‑148 (MQ=255)
                                                                                               cAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGCTACGGCAAtc                                               <  5:260888/149‑2 (MQ=255)
                                                                                                  aaTGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAAttctt                                            <  4:167078/149‑1 (MQ=255)
                                                                                                      ttCGATAAAAAGGAGGATAACTGCCGCAAGGGCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTg                                         <  6:301314/148‑1 (MQ=255)
                                                                                                          aTCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCAt                                    >  3:315356/1‑149 (MQ=255)
                                                                                                           tCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCAGTAGCGGCAATTCTTCTGCTCATc                                   >  2:417284/1‑149 (MQ=255)
                                                                                                              aaaGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGc                                        <  2:400263/141‑1 (MQ=255)
                                                                                                              aaaGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGc                                        >  1:400263/1‑141 (MQ=255)
                                                                                                                       tAACTGCCGCCACGTCCTCCTCGTGCCCCGAGCCGCCCCCGGCACCTTGGGCCCTGCCCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGACAATTCTTCTGCTCATCGACCTGTGGGcc                       <  6:238367/149‑1 (MQ=255)
                                                                                                                              cGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTGGGCCATCGTCa                <  1:417284/149‑1 (MQ=255)
                                                                                                                                       cccaaggTGCCGGGTGCAGCACCTGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTGGGCCATCGTCAGCGTGTTcc       <  1:175758/143‑1 (MQ=255)
                                                                                                                                       aTCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTGGGCCATCGTCAGCGTGTTc        <  1:124001/148‑1 (MQ=255)
                                                                                                                                         ccTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTGGGCCATCGTCCGCGTGTTCCgg     >  6:403253/1‑149 (MQ=255)
                                                                                                                                           tCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTGGGCCATCGTCAGCGTGTTCCGGAg   >  3:103556/1‑149 (MQ=255)
                                                                                                                                              tGCCGGGTGCAGCACCTGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTGGGCCATCGTCAGCGTGTTCCGGAGc  >  3:20684/1‑147 (MQ=255)
                                                                                                                                              tGCCGGGTGCAGCACCTGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTGGGCCATCGTCAGCGTGTTCCGGAGc  <  4:20684/147‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
CTTCCCGAACTGCTGGAGCACATCGAAGCAGGGCGCCTGCAACCGGAACTGATCGTTACTCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTGGGCCATCGTCAGCGTGTTCCGGAGC  >  NC_002947/4479667‑4479955

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: