Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 34,258 T→A 5.1% L232Q (CTG→CAG)  czrSA → metal‑sensor histidine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00294734,2580TA 5.1% 244.8 / 2.1 78L232Q (CTG→CAG) czrSAmetal‑sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base T (43/31);  new base A (2/2);  total (45/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 = 8.46e-01

CTGTGCGCAGAGGCCTGCTACCACTGACCAAGTTCCTCAAGATCGCCTCGAAAATCTCAACCGAGAGCCTCGACCATCGCCTTCCAACTGATGGGATGCCGGCAGAACTCAAAATGCTCGCCGACGGCATCAACATCATGCTGGCTCGCCTGGATGATGGGGTTCAGCAGCTCTCGCAGTTTTCTGACGACCTTGCTCACGAACTTCGAGCACCGCTTTCCAACCTGATGGGCAAGGCCCAGGTTGCCTTGACCAGAGAACGATCGCTTTCCGAGTACCGGGAGGTT  >  NC_002947/34117‑34403
                                                                                                                                             |                                                                                                                                                 
cTGTGCGCAGAGGCCTGCTACCACTGACCAAGTTCCTCAAGATCGCCTCGAAAATCTCAACCGAGAGCCTCGACCATCGCCTTCCAACTGATGGGATGCCGGCAGAACTCAAAATGCTCGCCGACGGCATCAACATCATGCTGGCTCGc                                                                                                                                            >  7:245821/1‑149 (MQ=255)
     cgcAGAGGCCTGCTACCACTGACCAAGTTCCTCAAGATCGCCTCGAAAATCTCAACCGAGAGCCTCGACCATCGCCTTCCAACTGATGGGATGCCGGCAGAACTCAAAATGCTCGCCGACGGCATCAACATCATGCTGGCTCGCCTGga                                                                                                                                       <  5:82798/149‑1 (MQ=255)
      gcAGAGGCCTGCTACCACTGACCAAGTTCCTCAAGATCGCCTCGAAAATCTCAACCGAGAGCCTCGACCATCGCCTTCCAACTGATGGGATGCCGGCAGAACTCAAAATGCTCGCCGACGGCATCAACATCATGCTGGCTCGCCTGgat                                                                                                                                      >  5:209144/1‑149 (MQ=255)
       cAGAGGCCTGCTACCACTGACCAAGTTCCTCAAGATCGCCTCGAAAATCTCAACCGAGAGCCTCGACCATCGCCTTCCAACTGATGGGATGCCGGCAGAACTCAAAATGCTCGCCGACGGCATCAACATCATGCTGGCTCGCCTGgat                                                                                                                                      >  2:183316/1‑148 (MQ=255)
       cAGAGGCCTGCTACCACTGACCAAGTTCCTCAAGATCGCCTCGAAAATCTCAACCGAGAGCCTCGACCATCGCCTACCAACTGATGGGATGCCGGCAGAACTCAAAATGCTCGCCGACGGCATCAACATCATGCTGGCTCGCCTGgat                                                                                                                                      >  2:135966/1‑148 (MQ=255)
            gCCTGCTACCACTGACCAAGTTCCTCAAGATCGCCTCGAAAATCTCAACCGAGAGCCTCGACCATCGCCTTCCAACTGATGGGATGCCGGCAGAACTCAAAATGCTCGCCGACGGCATCAACATCATGCTGGCTCGCCTGGATGATggg                                                                                                                                <  1:183316/149‑1 (MQ=255)
            gCCTGCTACCACTGACCAAGTTCCTCAAGATCGCCTCGAAAATCTCAACCGAGAGCCTCGACCATCGCCTTCCAACTGATGGGATGCCGGCAGAACTCAAAATGCTCGCCGACGGCATCAACATCATGCTGGCTCGCCTGGATGATggg                                                                                                                                <  8:262346/149‑1 (MQ=255)
                 cTACCACTGACCAAGTTCCTCAAGATCGCCTCGAAAATCTCAACCGAGAGCCTCGACCATCGCCTTCCAACTGATGGGATGCCGGCAGAACTCAAAATGCTCGCCGACGGCATCAACATCATCCAGGCGAGCCAGCatgatgttgat                                                                                                                             <  5:35774/147‑6 (MQ=255)
                  tACCACTGACCAAGTTCCTCAAGATCGCCTCGAAAATCTCAACCGAGAGCCTCGACCATCGCCTTCCAACTGATGGGATGCCGGCAGAACTCAAAATGCTCGCCGACGGCATCAACATCATGCTGGCTCGCCTGGATGATGGGGTTcag                                                                                                                          >  1:29238/1‑149 (MQ=255)
                     cACTGACCAAGTTCCTCAAGATCGCCTCGAAAATCTCAACCGAGAGCCTCGACCATCGCCTTCCAACTGATGGGATGCCGGCAGAACTCAAAATGCTCGCCGACGGCATCAACATCATGCTGGCTCGCCTGGATGATGGGGTTcagcag                                                                                                                       >  6:108076/1‑149 (MQ=255)
                     cACTGACCAAGTTCCTCAAGATCGCCTCGAAAATCTCAACCGAGAGCCTCGACCATCGCCTTCCAACTGATGGGATGCCGGCAGAACTCAAAATGCTCGCCGACGGCATCAACATCATCCAGGCGAGCCAGCatgatgttgat                                                                                                                             >  8:284703/1‑138 (MQ=255)
                     cACTGACCAAGTTCCTCAAGATCGCCTCGAAAATCTCAACCGAGAGCCTCGACCATCGCCTTCCAACTGATGGGATGCCGGCAGAACTCAAAATGCTCGCCGACGGCATCAACATCATCCAGGCGAGCCAGCatgatgttgat                                                                                                                             <  7:284703/143‑6 (MQ=255)
                          aCCAAGTTCCTCAAGATCGCCTCGAAAATCTCAACCGAGAGCCTCGACCATCGCCTTCCAACTGATGGGATGCCGGCAGAACTCAAAATGCTCGCCGACGGCATCAACATCATGCTGGCTCGCCTGGATGATGGGGTTCAGCAGctc                                                                                                                    <  6:209144/147‑1 (MQ=255)
                          aCCAAGTTCCTCAAGATCGCCTCGAAAATCTCAACCGAGAGCCTCGACCATCGCCTTCCAACTGATGGGATGCCGGCAGAACTCAAAATGCTCGCCGACGGCATCAACATCATGCTGGCTCGCCTGGATGATGGGGTTCAGCAGCtctc                                                                                                                  >  3:98540/1‑149 (MQ=255)
                          aCCAAGTTCCTCAAGATCGCCTCGAAAATCTCAACCGAGAGCCTCGACCATCGCCTTCCAACTGATGGGATGCCGGCAGAACTCAAAATGCTCGCCGACGGCATCAACATCATGCTGGCTCGCCTGGATGATGGGGTTCAGCAGCtctc                                                                                                                  >  3:80037/1‑149 (MQ=255)
                           ccAAGTTCCTCAAGATCGCCTCGAAAATCTCAACCGAGAGCCTCGACCATCGCCTTCCAACTGATGGGATGCCGGCAGAACTCAAAATGCTCGCCGACGGCATCAACATCATGCTGGCTCGCCTGGATGATGGGGTTCAGCAGCTCTCg                                                                                                                 <  4:80037/149‑1 (MQ=255)
                            cAAGTTCCTCAAGATCGCCTCGAAAATCTCAACCGAGAGCCTCGACCATCGCCTTCCAACTGATGGGATGCCGGCAGAACTCAAAATGCTCGCCGACGGCATCAACATCATGCTGGCTCGCCTGGATGATGGGGTTCAGCAGCTCTCg                                                                                                                 >  8:52388/1‑148 (MQ=255)
                             aaGTTCCTCAAGATCGCCTCGAAAATCTCAACCGAGAGCCTCGCCCATCGCCTTCCAACTGATGGGATGCCGGCAGAACTCAAAATGCTCGCCGACGGCATCAACATCATGCTGGCTCGCCTGGATGATGGGGTTCAGCAGCTCTCGCa                                                                                                               >  5:378587/1‑149 (MQ=255)
                                   cTCAAGATCGCCTCGAAAATCTCAACCGAGAGCCTCGACCATCGCCTTCCAACTGATGGGATGCCGGCAGAACTCAAAATGCTCGCCGACGGCATCAACATCATGCTGGCTCGCCTGGATGATGGGGTTCAGCAGCTCTCGCAGtttt                                                                                                          <  2:29238/148‑1 (MQ=255)
                                    tCAAGATCGCCTCGAAAATCTCAACCGAGAGCCTCGACCATCGCCTTCCAACTGATGGGATGCCGGCAGAACTCAAAATGCTCGCCGACGGCATCAACATCATGCTGGCTCGCCTGGATGATGGGGTTCAGCAGCTCTCGCAGTTTTCt                                                                                                        >  1:125251/1‑149 (MQ=255)
                                      aaGATCGCCTCGAAAATCTCCACCGCGAGCCTCGACCATCGCCTTCCAACTGATGGGATGCCGGCAGCACTCAAAATGCTCGCCGACGGCATCAACATCATGCTGGCTCGCCTGGATGATGGGGTTCAGCAGCTCTCGCAGTTTTCTgc                                                                                                      >  5:72736/1‑148 (MQ=255)
                                        gATCGCCTCGAAAATCTCAACCGAGAGCCTCGACCATCGCCTTCCAACTGATGGGATGCCGGCAGAACTCAAAATGCTCGCCGACGGCATCAACATCATGCTGGCTCGCCTGGATGATGGGGTTCAGCAGCTCTCGCAGTTTTCTgac                                                                                                     >  5:189657/1‑148 (MQ=255)
                                               tCGAAAATCTCAACCGAGAGCCTCGACCATCGCCTTCCAACTGATGGGATGCCGGCAGAACTCAAAATGCTCGCCGACGGCATCAACATCATGCTGGCTCGCCTGGATGATGGGGTTCAGCAGCTCTCGCa                                                                                                               <  4:177061/131‑1 (MQ=255)
                                               tCGAAAATCTCAACCGAGAGCCTCGACCATCGCCTTCCAACTGATGGGATGCCGGCAGAACTCAAAATGCTCGCCGACGGCATCAACATCATGCTGGCTCGCCTGGATGATGGGGTTCAGCAGCTCTCGCa                                                                                                               >  3:177061/1‑131 (MQ=255)
                                               tCGAAAATCTCAACCGAGAGCCTCGACCATCGCCTTCCAACTGATGGGATGCCGGCAGAACTCAAAATGCTCGCCGACGGCATCAACATCATGCTGGCTCGCCTGGATGATGGGGTTCAGCAGCTCTCGCAGTTTTCTGACGACCTTGc                                                                                             >  6:131550/1‑149 (MQ=255)
                                                         cAACCGAGAGCCTCGACCATCGCCTTCCAACTGATGGGATGCCGGCAGAACTCAAAATGCTCGCCGACTGCATCAACATCATGCTGGCTCGCCTGGATGATGGGGTTCAGCAGCTCTCGCAGTTTTCTGACCACCTTGCTCACGAACtt                                                                                   >  7:267088/1‑149 (MQ=255)
                                                          aaCCGAGAGCCTCGACCATCGCCTTCCAACTGATGGGATGCCGGGAGAACTCAAAATGCTCGCCGACGGCATCAACATCATGCTGGCTCGCCTGGATGATGGGGTTCAGCAGCTCTCGCAGTTTTCTGACGACCTTGCTCACGAACTTc                                                                                  >  1:79992/1‑149 (MQ=255)
                                                          aaCCGAGAGCCTCGACCATCGCCTTCCAACTGATGGGATGCCGGCATAAATCAATATGCTCGCCGACGGCATCAACATCATGCTGGCTCGCCTGGATGATGGGGTTCAGCAGCTCTCGCAGTTTTCTGACGACCTTGCTCACGTACTTc                                                                                  <  7:143836/149‑1 (MQ=255)
                                                          aaCCGAGAGCCTCGACCATCGCCTTCCAACTGATGGGATGCCGGCAGAACTCAAAATGCTCGCCGACGGCATCAACATCATGCTGGCTCGCCTGGATGATGGGGTTCAGCAGCTCTCGCAGTTTTCTGACGACCTTGCTCACGAACTTc                                                                                  <  3:135351/149‑1 (MQ=255)
                                                          aaCCGAGAGCCTCGACCATCGCCTTCCAACTGATGGGATGCCGGCAGAACTCAAAATGCTCGCCGACGGCATCAACATCATGCTGGCTCGCCTGGATGATGGGGTTCAGCAGCTCTCGCAGTTTTCTGACGACCTTGCTCACGAACTTc                                                                                  >  4:135351/1‑149 (MQ=255)
                                                            ccGAGAGCCTCGACCATCGCCTTCCAACTGATGGGATGCCGGCAGAACTCAAAATGCTCGCCGACGGCATCAACATCATGCTGGCTCGCCTGGATGATGGGGTTCAGCAGCTCTCGCAGTTTTCTGACGACCTTGCTCACGAACTTCGa                                                                                <  1:369491/149‑1 (MQ=255)
                                                             cGAGAGCCTCGACCATCGCCTTCCAACTGATGGGATGCCGGCAGAACTCAAAATGCTCGCCGACGGCATCAACATCATGCTGGCTCGCCTGGATGATGGGGTTCAGCAGCTCTCGCAGTTTTCTGACGACCTTGCTCACGAACTTCGAg                                                                               <  6:189657/149‑1 (MQ=255)
                                                             cGAGAGCCTCGACCATCGCCTTCCAACTGATGGGATGCCGGCAGAACTCAAAATGCTCGCCGACGGCATCAACATCATGCTGGCTCGCCTGGATGATGGGGTTCAGCAGCTCTCGCAGTTTTCTGACGACCTTGCTCACGAACTTCGAg                                                                               >  2:130778/1‑149 (MQ=255)
                                                                  gCCTCGACCATCGCCTTCCAACTGATGGGATGCTGGCAGAACTCAAAATGCTCGCCGACGGCATCAACATCATGCTGGCTCGCCTGGATGATGGGGTTCAGCAGCTCTCGCAGTTTTCTGACGACCTTGCTCACGAACTTCGAGCAcc                                                                           <  1:346182/148‑1 (MQ=255)
                                                                       gACCATCGCCTTCCAACTGATGGGATGCCGGCAGAACTCAAAATGCTCGCCGACGGCATCAACATCATGCTGGCTCGCCTGGATGATGGGGTTCAGCAGCTCTCGCAGTTTTCTGACGACCTTGCTCACGAACTTCGAGCACCGCTTTc                                                                     <  2:125251/149‑1 (MQ=255)
                                                                          cATCGCCTTCCAACTGATGGGATGCCGGCAGAACTCAAAATGCTCGCCGAAGGCATCAACATCATGCTGGCTCGCCTGGATGATGGGGTTCAGCAGCTCTCGCAGTTTTCTGACGACCTTGCTCACGAACTTCGAGCACCGCTTTCCaa                                                                  >  6:220824/1‑149 (MQ=255)
                                                                           aTCGCCTTCCAACTGATGGGATGCCGGCAGAACTCAAAATGCTCGCCGACGGCATCAACATCATGCTGGCTCGCCTGGATGATGGGGTTCAGCAGCTCTCGCAGTTTTCTGACGACCTTGCTCACGAACTTCGAGCACCGCTTTCCAAc                                                                 <  1:307494/149‑1 (MQ=255)
                                                                            tCGCCTTCCAACTGATGGGATGCCGGCAGAACTCAAAATGCTCGCCGACGGCATCAACATCATGCTGGCTCGCCTGGATGATGGGGTTCAGCAGCTCTCGCAGTTTTCTGACGACCTTGCTCACGAACTTCGAGCACCGCTTTCCAAcc                                                                >  1:334563/1‑149 (MQ=255)
                                                                             cGCCTTCCAACTGATGGGATGCCGGCAGAACTCAAAATGCTCGCCGACGGCATCAACATCATGCTGGCTCGCCTGGATGATGGGGTTCAGCAGCTCTCGCAGTTTTCTGACGACCTTGCTCACGAACTTCGAGCACCGCTTTCCAACCt                                                               >  6:89248/1‑149 (MQ=255)
                                                                               ccTTCCAACTGATGGGATGCCTGCAGAACTCAAAATGCTCGCCGACGGCATCAACATCATGCTGTCTCGCCTTGATGATGGGGTTCAGCAGCTCTCGCAGTTTTCTGACGACCTTGCTCACGAACTTCGAGCACCGCTTTTCAACCTGa                                                             <  7:153004/149‑1 (MQ=255)
                                                                                 ttCCAACTGATGGGATGCCGGCAGAACTCAAAATGCTCGCCGACGGCATCAACATCATGCTGGCTCGCCTGGATGATGGGGTTCAGCAGCTCTCGCAGTTTTCTGACGACCTTGCTCACGAACTTCGAGCACCGCTTTCCAACCTGATg                                                           >  5:342286/1‑149 (MQ=255)
                                                                                   ccAACTGATGGGATGCCGGCAGAACTCAAAATGCTCGCCGACGGCATCAACATCATGCTGGCTCGCCTGGATGATGGGGTTCAGCAGCTCTCGCAGTTTTCTAACGACCTTGCTCACGAACTTCGAGCACCGCTTTCCAACCTgggcc                                                          <  8:245821/148‑5 (MQ=255)
                                                                                    cAACTGATGGGATGCCGGCAGAACTCAAAATGCTCGCCGACGGCATCAACATCATGCTGGCTCGCCTGGATGATGGGGTTCAGCAGCTCTCGCAGTTTTCTGACGACCTTGCTCACGAACTTCGAGCACCGCTTTCCAACCTGATggg                                                         >  8:39443/1‑148 (MQ=255)
                                                                                     aaCTGATGGGATGCCGGCAGAACTCAAAATGCTCGCCGACGGCATCAACATCATGCTGGCTCGCCTGGATGATGGGGTTCAGCAGCTCTCGCAGTTTTCTGACGACCTTGCTCACGAACTTCGAGCACCGCTTTCCAACCTGATGGGc                                                        >  2:31211/1‑148 (MQ=255)
                                                                                     aaCTGATGGGATGCCGGCAGAACTCAAAATGCTCGCCGACGGCATCAACATCATGCTGGCTCGCCTGGATGATGGGGTTCAGCAGCTCTCGCAGTTTTCTGACGACCTTGCTCACGAACTTCGAGCACCGCTTTCCAACCTGATGGGc                                                        >  1:267528/1‑148 (MQ=255)
                                                                                     aaCTGATGGGATGCCGGCAGAACTCAAAATGCTCGCCGACGGCATCAACATCATGCTGGCTCGCCTGGATGATGGGGTTCAGCAGCTCTCGCAGTTTTCTGACGACCTTGCTCACAAACTTCGAGCACCGCTTTCCAACCTGATGGGc                                                        >  4:383536/1‑148 (MQ=255)
                                                                                      aCTGATGGGATGCCGGCAGAACTCAAAATGCTCGCCGACGGCATCAACATCATGCTGGCTCGCCTGGATGATGGGGTTCAGCAGCTCTCGCAGTTTTCTGACGACCTTGCTCACGAACTTCGAGCACCGCTTTCCAACCTGATGGGCaa                                                      <  3:121626/149‑1 (MQ=255)
                                                                                        tGATGGGATGCCGGCAGAACTCAAAATGCTCGCCGACGGCATCAACATCATGCTGGCTCGCCTGGATGATGGGGTTCAGCAGCTCTCGCAGTTTTCTGACGACCTTGCTCACGAACTTCGAGCACCGCTTTCCAACCTGATGGGCAAgg                                                    >  5:157104/1‑149 (MQ=255)
                                                                                          aTGGGATGCCTGCAGATCTAAAAATGCTCGCCGACGGCATCAACATCATGCTGGCTCGCCTGGATGATGGGGTTCAGCAGCTCTCGCAGTTTTCTGACGACGTTGCTCTCGAACTTCGAGCACCGCTTTCCAACCTGATGGGCAAgg                                                    <  2:375838/147‑1 (MQ=255)
                                                                                            gggATGCCGGCAGAACTCAAAATGCTCGCCGACGGCATCAACATCATGCTGGCTCGCCTGGATGATGGGGTTCAGCAGCTCTCGCAGTTTTCTGACGACCTTGCTCACGAACTTCGAGCACCGCTTTCCAACCTGATGGGCAAGGcc                                                  >  2:257407/1‑147 (MQ=255)
                                                                                             ggATGCCGGCAGAACTCAAAATGCTCGCCGACGGCATCAACATCATGCTGGCTCGCCTGGATGATGGGGTTCAGCAGCTCTCGCAGTTTTCTGACGACCTTGCTCACGAACTTCGAGCACCGCTTTCCAACCTGATGGGCAAGGCCCAg                                               <  4:98540/149‑1 (MQ=255)
                                                                                             ggATGCCGGCAGAACTCAAAATGCTCGCCGACGGCATCAACATCATGCTGGCTCGCCTGGATGATGGGGTTCAGCAGCTCTCGCAGTTTTCTGACGACCTTGCTCACAAACTTCGAGCACCGCTTTCCAACCTGATGGGCAAGGCCCa                                                <  3:383536/148‑1 (MQ=255)
                                                                                              gATGCCGGCAGAACTCAAAATGCTCGCCGACGGCATCAACATCATGCTGGCTCGCCTGGATGATGGGGTTCAGCAGCTCTCGCAGTTTTCTGACGACCTTGCTCACGAACTTCGAGCACCGCTTTCCAACCTGATGGTCAAGGcaacct                                              >  8:238024/1‑144 (MQ=255)
                                                                                              gATGCCGGCAGAACTCAAAATGCTCGCCGACGGCATCAACATCATGCTGGCTCGCCTGGATGATGGGGTTCAGCAGCTCTCGCAGTTTTCTGACGACCTTGCTCACGAACTTCGAGCACCGCTTTCCAACCTGATGGGCAAGGCCCAgg                                              >  1:211757/1‑149 (MQ=255)
                                                                                               aTGCCGGCTGAACTCAAAATGCTCGCCGACGGCATCAACATCATGCTGGCTCGCCTGGATGATGGGGTTCAGCAGCACTCGCAGATTTCTGACGACCTTGCTCACGAACTTCGAGCTCCGGTTTCGAACATGATGGGCAAGGCCAAGGt                                             >  7:304399/1‑149 (MQ=255)
                                                                                               aTGCCGGCAGAACTCAAAATGCTCGCCGACGGCATCAACATCATGCTGGCTCGCCTGGATGATGGGGTTCAGCAGCTCTCGCAGTTTTCTGACGACCTTGCTCACGAACTTCGAGCACCGCTTTCCAACCTGATGGGCAAGGCCCAGGt                                             >  4:267090/1‑149 (MQ=255)
                                                                                                    ggCAGAACTCAAAATGCTCGCCGACGGCATCAACATCATGCTGGCTCGCCTGGATGATGGGGTTCAGCAGCTCTCGCAGTTTTCTGACGCCCTTGCTCACGAACTTCGAGCACCGCTTTCCAACCTGATGGGCAAGGCCCAGGTTGCCt                                        >  5:280093/1‑149 (MQ=255)
                                                                                                    ggCAGAACTCAAAATGCTCGCCGACGGCATCAACATCATGCTGGCTCGCCTGGATGATGGGGTTCAGCAGCTCTCGCAGTTTTCTGACGACCTTGCTCACGAACTTCGAGCACCGCTTTCCAACCTGATGGGCAAGGCCCAGGTTGCCt                                        >  4:194812/1‑149 (MQ=255)
                                                                                                        gAACTCAAAATGCTCGCCGACGGCATCAACATCATGCTGGCTCGCCTGGATGATGGGGTTCAGCAGCTCTCGCAGTTTTCTGACGACCTTGCTCACGAACTTCGAGCACCGCTTTCCAACCTGATGGGCAAGGCCCAGGTTGCCTTGAc                                    <  2:79992/149‑1 (MQ=255)
                                                                                                        gAACTCAAAATGCTCGCCGACGGCATCAACATCATGCTGGCTCGCCTGGATGATGGGGTTCAGCAGCTCTCGCAGTTTTCTGACGACCTTGCTCACGAACTTCGAGCACCGCTTTCCAACCTGATGGGCAAGGCCCAGGTTGCCTTGAc                                    >  6:293056/1‑149 (MQ=255)
                                                                                                         aaCTCAAAATGCTCGCCGACGGCATCAACATCATGCTGGCTCGCCTGGATGATGGGGTTCAGCAGCTCTCGCAGTTTTCTGACGCCCTTGCTCACGAACTTCGAGCACCGCTTTCCAACCTGATGGGCAAGGCCCAGGTTGCCTTGAc                                    >  7:174870/1‑148 (MQ=255)
                                                                                                           cTCAAAATGCTCGCCGACGGCATCAACATCATGCTGGCTCGCCTGGATGATGGGGTTCAGCAGCTCTCGCAGTTTTCTGACGACCTTGCTCACGAACTTCGAGCACCGCTTTCCAACCTGATGGGCAAGGCCCAGGTTGCCTTGAc                                    >  7:240020/1‑146 (MQ=255)
                                                                                                              aaaaTGCTCGCCGACGGCATCAACATCATGCTGGCTCGCCTGGATGATGGGGTTCAGCAGCTCTCGCAGTTTTCTGACGACCTTGCTCACGAACTTCGAGCACCGCTTTCCAACCTGATGGGCAAGGCCCAGGTTGCCTTGACCagag                               >  4:83455/1‑148 (MQ=255)
                                                                                                              aaaaTGCTCGCCGACGGCATCAACATCATGCTGGCTCGCCTGGATGATGGGGTTCAGCAGCTCTCGCAGTTTTCTGACGACCTTGCTCACGAACTTCGAGCACCGCTTTCCAACCTGATGGGCAAGGCCCAGGTTGCCTTGACCaga                                <  3:194812/147‑1 (MQ=255)
                                                                                                                    atCGCCCACGGCAGCAACATGATGCTTGCTCGGCTGGATGATGGGGTGCAGCAGCTCTCGCAGTTTTGTGACGACCTTGATCTCGAACTTCGAGCACCGCTTTCCACCCTGATGGGCAAGGCGCAGGTTGCCTTGACGAGAgatcgatc                        <  8:304399/148‑1 (MQ=255)
                                                                                                                       gCCGACGGCATCAACATCATGCTGGCTCGCCTGGATGATGGGGTTCAGCAGCTCTCGCAGTTTTCTGACGACCTTGCTCACGAACTTCGAGCACCGCTTTCCAACCTGATGGGCAAGGCCCAGGTTGCCTTGACCAGAGAACGATCg                       <  1:130778/147‑1 (MQ=255)
                                                                                                                       gCCGACGGCATCAACATCATGCTGGCTCGCCTGGATGATGGGGTTCAGCAGCTCTCGCAGTTTTCTGACGACCTTGCTCACGAACTTCGAGCACCGCTTTCCAACCTGATGGGCAAGGCCCAGGTTGCCTTGACCAGAGAACGATCGc                      <  2:211757/148‑1 (MQ=255)
                                                                                                                        ccGACGGCATCAACATCATGCTGGCTCGCGTGGATGATGGGGTTCAGCAGCTCTCGCAGTTTTCTGACGACCTTGCTCACGAACTTCGAGCACCGCTTTCCAACCTGATGGGCAAGGCCCAGGTTGCCTTGACCAGAGAACGATCGCtt                    <  5:293056/149‑1 (MQ=255)
                                                                                                                        ccGAAGGCATCTACTTCATGCTGGCTCGCCTGGATGATGGGGTTCAGCAGCTCTGGCAGTTTTCTGACGACCTTTCTCACTAACTTCGAGCACCGCTTTCGAACCTGATGGGCAAGGCCCAGGTTGCCTTGACCAGAGAACGATCGCtt                    <  5:400222/149‑1 (MQ=255)
                                                                                                                           aCGGCATCAACATCATGCTGGCTCGCCTGGATGATGGGGTTCAGCAGCTCTCGCAGTTTTCTGACGACCTTGCTCACGAACTTCGAGCACCGCTTTCCACCCTGATGGGCAAGGCCCAGGTTGCCTTGACCAGAGAACGATCGCTTTcc                 <  6:72736/149‑1 (MQ=255)
                                                                                                                             ggCATCAACATCATGCTGGCTCGCCTGGATGATGGGGTTCAGCAGCTCTCGCAGTTTTCTGACGACCTTGCTCACGAACTTCGAGCACCGCTTTCCAACCTGATGGGCAAGGCCCAGGTTGCCTTGACCAGAGAACGATCGCTTTCCGa               <  6:378587/149‑1 (MQ=255)
                                                                                                                             ggCATCAACATCATGCTGGCTCGCCTGGATGATGGGGTTCAGCAGCTCTCGCAGTTTTCTGACGACCTTGCTCACGAACTTCGAGCACCGCTTTCCAACCTGATGGGCAAGGCCCAGGTTGCCTTGACCAGAGAACGATCGCTTTCCGa               <  7:52388/149‑1 (MQ=255)
                                                                                                                                  ccccatcatCCAGGCGAGCCAGCATGATGGGGTTCAGCAGCTCTCGCAGTTTTCTGACGACCTTGCTCACGAACTTCGAGCACCGCTTTCCAACCTGATGGGCAAGGCCCAGGTTGCCTTGACCAGAGAACGATCGCTTTCCGAGTAcc          >  1:23521/4‑149 (MQ=255)
                                                                                                                                   aaCATCATGCTGGCTCGCCTGGATGATGGGGTTCAGCAGCTCTCGCAGTTTTCTGACGACCTTGCTCACGAACTTCGAGCACCGCTTTCCAACCTGATGGGCAAGGCCCAGGTTGCCTTGACCAGAGAACGATCGCTTTCCGAGTACCg         >  3:194600/1‑149 (MQ=255)
                                                                                                                                    aCATCATGCTGGCTCGCCTGGATGATGGGGTTCAGCAGCTCTCGCAGTTTTCTGACGACCTTGCTCACGAACTTCGAGCACCGCTTTCCAACCTGATGGGCAAGGCCCAGGTTGCCTTGACCAGAGAACGATCGCTTTCCGAGTACCgg        >  1:20587/1‑149 (MQ=255)
                                                                                                                                      atcatGCTGGCTCGCCTGGATGATGGGGTTCAGCAGCTCTCGCAGTTTTCTGACGACCTTGCTCACGAACTTCGAGCACCGCTTTCCAACCTGATTGGCAAGGCCCAGGTTGCCTTGACCAGAGAACGATCGCTTTCCGAGTACCgg        >  4:274351/1‑147 (MQ=255)
                                                                                                                                       tcatGCTGGCTCGCCTGGATGATGGGGTTCAGCAGCTCTCGCAGTTTTCTGACGACCTTGCTCACGAACTTCGAGCACCGCTTTCCAACCTGATGGGCAAGGCCCAGGTTGCCTTGACCAGAGAACGATCGCTTTCCGAGTACCGGGAg     <  5:131550/149‑1 (MQ=255)
                                                                                                                                          tGCTGGCTCGCCTGGATGATGGGGTTCAGCAGCTCTCGCAGTTTTCTGACGACCTTGCTCACGAACTTCGAGCACCGCTTTCCAACCTGATGGGCAAGGCCCAGGTTGCCTTGACCAGAGAACGATCGCTTTCCGAGTACCGGGAGGtt  >  2:127823/1‑149 (MQ=255)
                                                                                                                                             |                                                                                                                                                 
CTGTGCGCAGAGGCCTGCTACCACTGACCAAGTTCCTCAAGATCGCCTCGAAAATCTCAACCGAGAGCCTCGACCATCGCCTTCCAACTGATGGGATGCCGGCAGAACTCAAAATGCTCGCCGACGGCATCAACATCATGCTGGCTCGCCTGGATGATGGGGTTCAGCAGCTCTCGCAGTTTTCTGACGACCTTGCTCACGAACTTCGAGCACCGCTTTCCAACCTGATGGGCAAGGCCCAGGTTGCCTTGACCAGAGAACGATCGCTTTCCGAGTACCGGGAGGTT  >  NC_002947/34117‑34403

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: