Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,543,595 C→G 9.0% A10A (GCC→GCG rnz → ribonuclease Z

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,543,5950CG 9.0% 231.9 / 11.7 78A10A (GCC→GCGrnzribonuclease Z
Reads supporting (aligned to +/- strand):  ref base C (36/35);  new base G (3/4);  total (39/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 = 4.88e-01

CAGCAAGCTCATCTGGCGACAGGATAACTAGCGGGCGGTCGACGGCGGGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGCCTGCTGCACACCCCGCTTTCGATCCGCGACCTGCGCGCGATCTT  >  NC_002947/4543448‑4543741
                                                                                                                                                   |                                                                                                                                                  
cAGCAAGCTCATCTGGCGACAGGATAACTAGCGGGCGGTCGACGGCGGGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTccgccg                                                                                                                                                   >  8:165810/1‑149 (MQ=255)
   cAAGCTCATCTGGCGACAGGATAACTAGCGGGCGGTCGACGGCGGGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCgggg                                                                                                                                                >  4:149490/1‑149 (MQ=255)
    aaGCTCATCTGGCGACAGGATAACTAGCGGGCGGTCGACGGCGGGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGt                                                                                                                                               >  1:114094/1‑149 (MQ=255)
       cTCATCTGGCGACAGGATAACTAGCGGGCGGTCGACGGCGGGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCCCGGCGGAGGTGcc                                                                                                                                            <  3:31614/149‑1 (MQ=255)
          aTCTGGCGACAGGATAACTAGCGGGCGGTCGACGGCGGGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCccac                                                                                                                                         >  6:262597/1‑149 (MQ=255)
          aTCTGGCGACAGGATAACTAGCGGGCGGTCGACGGCGGGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCCCGGCGGAGGTGccaa                                                                                                                                          <  5:234494/148‑3 (MQ=255)
              ggCGACAGGATAACTAGCGGGCGGTCGAGGGCGGGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAg                                                                                                                                     >  1:212113/1‑149 (MQ=255)
               gCGACAGGATAACTAGCGGGCGGTCGACGGCGGGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAgg                                                                                                                                    <  7:165810/149‑1 (MQ=255)
                cGACAGGATAACTAGCGGGCGGTCGACGGCGGGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGc                                                                                                                                   <  2:212113/149‑1 (MQ=255)
                   cAGGATAACTAGCGGGCGGTCGACGGCGGGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACg                                                                                                                                >  7:260799/1‑149 (MQ=255)
                      gATAACTAGCGGGCGGTCGACGGCGGGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCaa                                                                                                                             >  5:182514/1‑149 (MQ=255)
                         aaCTAGCGGGCGGTCGACGGCGGGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTCGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGt                                                                                                                          <  8:310213/149‑1 (MQ=255)
                          aCTAGCGGGCGGTCGACGGCGGGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTc                                                                                                                         <  4:141780/149‑1 (MQ=255)
                            tAGCGGGCGGTCGACGGCGGGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAg                                                                                                                       >  4:315144/1‑149 (MQ=255)
                            tAGCGGGCGGTCGACGGCGGGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCAGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAg                                                                                                                       >  3:227054/1‑149 (MQ=255)
                                 ggcggTCGACGGCGGGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCca                                                                                                                  >  8:99024/1‑149 (MQ=255)
                                    ggTCGACGGCGGGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGTCCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAg                                                                                                               >  1:189844/1‑149 (MQ=255)
                                    ggTCGACGGCGGGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAg                                                                                                               <  7:59425/149‑1 (MQ=255)
                                    ggTCGACGGCGGGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAg                                                                                                               >  1:42798/1‑149 (MQ=255)
                                        gACGGCGGGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGt                                                                                                           <  7:99024/149‑1 (MQ=255)
                                               ggTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGaa                                                                                                    <  2:114094/149‑1 (MQ=255)
                                               ggTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGaa                                                                                                    <  2:42798/149‑1 (MQ=255)
                                                     aaTCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGc                                                                                              >  5:42235/1‑149 (MQ=255)
                                                      aTCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCg                                                                                             >  4:76777/1‑149 (MQ=255)
                                                         cGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGaa                                                                                                    >  6:86391/1‑139 (MQ=255)
                                                         cGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGaa                                                                                                    <  5:86391/139‑1 (MQ=255)
                                                          ggCAAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAg                                                                                         >  3:33315/1‑149 (MQ=255)
                                                             tAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAgc                                                                                                                      >  8:23075/1‑117 (MQ=255)
                                                             tAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAgc                                                                                                                      <  7:23075/117‑1 (MQ=255)
                                                              aaGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAg                                                                                                                                     <  4:80793/101‑1 (MQ=255)
                                                              aaGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAg                                                                                                                                     >  3:80793/1‑101 (MQ=255)
                                                                ggTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTg                                                                                   <  3:76777/149‑1 (MQ=255)
                                                                    tgCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTAt                                                                               <  3:149490/149‑1 (MQ=255)
                                                                     gcgcCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATc                                                                              <  4:33315/149‑1 (MQ=255)
                                                                      cgcCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCt                                                                             <  1:514/149‑1 (MQ=255)
                                                                      cgcCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCt                                                                             <  4:108946/149‑1 (MQ=255)
                                                                          ccGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTc                                                                         <  6:42235/149‑1 (MQ=255)
                                                                             tATCACCTTGGCCGCTTCAGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGAt                                                                      <  4:227054/149‑1 (MQ=255)
                                                                               tCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGAtt                                                                     <  8:260799/148‑1 (MQ=255)
                                                                                 aCCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCCCCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTgcg                                                                  >  3:15272/1‑149 (MQ=255)
                                                                                  ccTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCTGCCACTGGTATCTGGTCGATTgcgg                                                                 >  2:120071/1‑149 (MQ=255)
                                                                                       gCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAg                                                            <  3:315144/149‑1 (MQ=255)
                                                                                        ccGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAgg                                                           >  1:358085/1‑149 (MQ=255)
                                                                                           cTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCAc                                                        <  3:124417/149‑1 (MQ=255)
                                                                                                ggCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTa                                                                                <  4:123659/120‑1 (MQ=255)
                                                                                                ggCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTa                                                                                >  3:123659/1‑120 (MQ=255)
                                                                                                 gCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCa                                                  <  5:262597/149‑1 (MQ=255)
                                                                                                   cTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCAcc                                                       >  1:36076/1‑142 (MQ=255)
                                                                                                   cTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCAcc                                                       <  2:36076/142‑1 (MQ=255)
                                                                                                    tGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAAt                                                                                                                            >  6:119592/1‑72 (MQ=255)
                                                                                                    tGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTAt                                                                               >  1:318186/1‑117 (MQ=255)
                                                                                                    tGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTAt                                                                               <  2:318186/117‑1 (MQ=255)
                                                                                                    tGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCg                                               >  5:248225/1‑149 (MQ=255)
                                                                                                     gAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGc                                              >  5:244464/1‑149 (MQ=255)
                                                                                                        aCAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGATGCCAGCGGCAGCCACTGGt                                                                                 <  4:338814/111‑1 (MQ=255)
                                                                                                        aCAGGAGCCTGACCATCGACCTGCACTTCCTTGGCACCTCCCCCTGCCTGCCCACCAACGCACGCAATGACAGCGCCACAGCCGTGATCCAAGCCAGCGGCAGCCACtcga                                                                                 >  3:338814/1‑108 (MQ=255)
                                                                                                          aGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCCCCCAGCACCGCctgct                                         >  6:268039/1‑149 (MQ=255)
                                                                                                                     cATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGCCTGCTGCACACCCCGc                              >  7:347427/1‑149 (MQ=255)
                                                                                                                         gACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGCCTGCTGCACACCCCGCTTTc                          <  6:244464/149‑1 (MQ=255)
                                                                                                                           ccTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGCCTGCTGCACACCCCGCTTTCGa                        >  3:340316/1‑149 (MQ=255)
                                                                                                                            ctgctgTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGATGCCAGCGGCTGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGCCTGCTGCACACCCCGCTTTCGAt                       <  1:120071/149‑1 (MQ=255)
                                                                                                                            ctgctgTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGcc                                             >  2:236945/1‑127 (MQ=255)
                                                                                                                            ctgctgTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGcc                                             <  1:236945/127‑1 (MQ=255)
                                                                                                                              gctgTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTg                                                                            >  3:364275/1‑94 (MQ=255)
                                                                                                                              gctgTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTg                                                                            <  4:364275/94‑1 (MQ=255)
                                                                                                                               ctgTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGGCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGCCTGCTGCACACCCCGCTTTCGATc                      <  6:182514/147‑1 (MQ=255)
                                                                                                                                  ttggtgggCACCCCGGCGGAGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGCctgct                                         >  2:111125/7‑125 (MQ=255)
                                                                                                                                  ttggtgggCACCCCGGCGGAGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGCctgct                                         <  1:111125/119‑1 (MQ=255)
                                                                                                                                    ggtgggCACCCCGGCGGAGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGCCTGCTGCACACCCCGCTTTCGATCCGCGAcc               >  3:107661/5‑149 (MQ=255)
                                                                                                                                     gtgggCACCCCGGCGGAGGTGCCAACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGCCTGCTGCACACCCCGCTTTCGATCcgc                   <  1:194399/141‑1 (MQ=255)
                                                                                                                                     gtgggCACCCCGGCGGAGGTGCCAACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGCCTGCTGCACACCCCGCTTTCGATCcgc                   >  2:194399/4‑144 (MQ=255)
                                                                                                                                      ttGGCACCTCCGCCGGGGTGCCCACCAAGTCACGCTATGTCTGCGCCACAGCCGTGATCGTAGAGATCGGCAGCAACTGGTATGTGGTCGATTGCGGCTaa                                                             <  2:15891/101‑1 (MQ=255)
                                                                                                                                      ttGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGaa                                                             >  1:15891/1‑101 (MQ=255)
                                                                                                                                        ggCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTAGATTGCGGCGAAGGCACCCATCACCGCCTGCTGCACACCCCGCTTTCGATCCGCGACCTgc            <  2:358085/148‑1 (MQ=255)
                                                                                                                                             cTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGCCTGCTGCACACCCCGCTTTCGATCCGCGACCTgcgcg         >  7:294696/1‑146 (MQ=255)
                                                                                                                                             cTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGCCTGCTGCACACCCCGCTTTCGATCCGCGACCTgcgcg         <  8:294696/146‑1 (MQ=255)
                                                                                                                                             cTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGCCTGCTGCACACCCCGCTTTCGATCCGCGACCTGCGCGCGa      >  7:12384/1‑149 (MQ=255)
                                                                                                                                              tCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGCCTGCTGCACACCCCGCTTTCGATcc                     >  6:49682/1‑133 (MQ=255)
                                                                                                                                              tCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGCCTGCTGCACACCCCGCTTTCGATcc                     <  5:49682/133‑1 (MQ=255)
                                                                                                                                               ccgccgGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGCctgc                                          <  4:269918/111‑1 (MQ=255)
                                                                                                                                               ccgccgGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGCctgc                                          >  3:269918/1‑111 (MQ=255)
                                                                                                                                                 gccgGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGCCTGCTGCACACCCCGCTTTCGATCCGCGACCTGCGCGCGATCtt  <  8:347427/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CAGCAAGCTCATCTGGCGACAGGATAACTAGCGGGCGGTCGACGGCGGGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGCCTGCTGCACACCCCGCTTTCGATCCGCGACCTGCGCGCGATCTT  >  NC_002947/4543448‑4543741

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: