Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,212,472 C→T 7.2% intergenic (+36/+58) PP_4593 → / ← PP_4594 hypothetical protein/cystathionine gamma‑synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,212,4720CT 7.2% 207.2 / 4.0 69intergenic (+36/+58)PP_4593/PP_4594hypothetical protein/cystathionine gamma‑synthase
Reads supporting (aligned to +/- strand):  ref base C (31/33);  new base T (2/3);  total (33/36)
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.80e-01

TAGCCGCCGAAGGCAACGGCACCAACGACGAAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGACAACCGCAGCAGG  >  NC_002947/5212327‑5212617
                                                                                                                                                  |                                                                                                                                                 
tAGCCGCCGAAGGCAACGGCACCAACGACGAAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAg                                                                                                                                                >  5:27523/1‑149 (MQ=255)
 aGCCGCCGAAGGCAACGGCACCAACGACGAAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAga                                                                                                                                               >  6:217331/1‑148 (MQ=255)
  gccgccGAAGGCAACGGCACCAACGACGAAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGcc                                                                                                                                              >  5:144870/1‑149 (MQ=255)
   ccgccgAAGGCAACGGCACCAACGACGAAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGGAAAAAGGCCAAGATAGAACGGATCATTAATCCGCCCAAATAATGACAGGCATCTATTGAATCCACCGGCCCTTTCGCGGGTAAAcc                                                                                                                                              >  1:323384/1‑149 (MQ=255)
    cgccgAAGGCAACGGCACCAACGACGAAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCg                                                                                                                                            >  3:219447/1‑149 (MQ=255)
       cgAAGGCAACGGCACCAACGACGAAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTc                                                                                                                                         <  2:314802/149‑1 (MQ=255)
        gAAGGCAACGGCACCAACGACGAAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTcc                                                                                                                                        >  7:245923/1‑149 (MQ=255)
            gCAACGGCACCAACGACGAAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTaca                                                                                                                                    <  2:23086/149‑1 (MQ=255)
             cAACGGCACCAACGACGAAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGTAAACCCGCTCCTacga                                                                                                                                   <  3:234439/149‑3 (MQ=255)
                       aaCGACGAAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCCATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGa                                                                                                                         <  1:313963/149‑1 (MQ=255)
                         cgacgaAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCCTTTCGCGGGTAAACCCGCTCCTACGACGGTGTAGGAGc                                                                                                                       <  2:323384/149‑1 (MQ=38)
                           acgaAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCgg                                                                                                                     >  6:319579/1‑149 (MQ=255)
                            cgaAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCggg                                                                                                                    <  7:133114/149‑1 (MQ=255)
                             gaAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGTAAACCCGCTCCTACGACGGTGTAGGAGCgggc                                                                                                                   >  1:60210/1‑148 (MQ=38)
                               aGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGttt                                                                                                                 <  5:217331/149‑1 (MQ=255)
                                  cGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGTAAACCCGCTCCTACGACGGTGTAGGAGCGGGCTTGcc                                                                                                              <  2:60210/149‑1 (MQ=37)
                                  cGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTAcc                                                                                                              <  6:144870/149‑1 (MQ=255)
                                    tCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTa                                                                                                                                      <  2:326129/123‑1 (MQ=255)
                                    tCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTa                                                                                                                                      >  1:326129/1‑123 (MQ=255)
                                      aaCCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCcgc                                                                                                           <  7:171275/148‑1 (MQ=255)
                                        ccGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTa                                                                                                                            <  8:150070/129‑1 (MQ=255)
                                        ccGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTa                                                                                                                            >  7:150070/1‑129 (MQ=255)
                                          gCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAg                                                                                                      >  8:206917/1‑149 (MQ=255)
                                          gCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAg                                                                                                      >  4:211918/1‑149 (MQ=255)
                                               cgcAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCt                                                                                                 <  2:405671/149‑1 (MQ=255)
                                                  aaGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCg                                                                                              >  5:22783/1‑149 (MQ=255)
                                                    ggAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACcgtcg                                                                                                                              >  1:32047/1‑115 (MQ=255)
                                                    ggAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACcgtcg                                                                                                                              <  2:32047/115‑1 (MQ=255)
                                                          gTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTcac                                                                                      <  3:211918/149‑1 (MQ=255)
                                                            tGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTa                                                                                                                >  7:47188/1‑121 (MQ=255)
                                                            tGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTa                                                                                                                <  8:47188/121‑1 (MQ=255)
                                                            tGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTcaca                                                                                     >  6:118813/1‑148 (MQ=255)
                                                                     aaGGACAAGATCGATCGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGAc                                                                              >  3:69498/1‑147 (MQ=255)
                                                                      aGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCt                                                                                                 <  2:206430/127‑1 (MQ=255)
                                                                      aGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCt                                                                                                 >  1:206430/1‑127 (MQ=255)
                                                                        gACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCcg                                                                                                            >  1:119453/1‑114 (MQ=255)
                                                                        gACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCcg                                                                                                            <  2:119453/114‑1 (MQ=255)
                                                                          cAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCt                                                                       <  6:22783/149‑1 (MQ=255)
                                                                          cAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCt                                                                       >  4:389766/1‑149 (MQ=255)
                                                                            aGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACcgtc                                                                                                                               >  8:357099/1‑91 (MQ=255)
                                                                            aGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACcgtc                                                                                                                               <  7:357099/91‑1 (MQ=255)
                                                                                  aaCGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCt                                                                       >  3:233303/1‑141 (MQ=255)
                                                                                  aaCGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCt                                                                       <  4:233303/141‑1 (MQ=255)
                                                                                               tCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGcc                                                              >  1:307531/1‑137 (MQ=255)
                                                                                               tCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGcc                                                              <  2:307531/137‑1 (MQ=255)
                                                                                               tCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTcc                                                  <  4:69498/149‑1 (MQ=255)
                                                                                                  gCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTc                                                >  6:110311/1‑148 (MQ=255)
                                                                                                   cccAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGa                                              >  5:264042/1‑149 (MQ=255)
                                                                                                      aaaTCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCa                                           <  5:319579/149‑1 (MQ=255)
                                                                                                       aaTCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAg                                          <  8:245923/149‑1 (MQ=255)
                                                                                                         tCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGt                                                                                                                   >  8:373842/1‑74 (MQ=255)
                                                                                                         tCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGt                                                                                                                   <  7:373842/74‑1 (MQ=255)
                                                                                                             gACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACg                                    >  3:193180/1‑149 (MQ=255)
                                                                                                              aCAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACgg                                   <  4:193180/149‑1 (MQ=255)
                                                                                                                  gCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAAt                               >  8:95732/1‑149 (MQ=255)
                                                                                                                   cATCTATTGAATCCACCGGCCTCTTCGCCGGCAACCCCGCTCCTCCACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGAAATc                              <  8:267679/149‑1 (MQ=255)
                                                                                                                     tCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAgcg                                                                         >  3:364748/1‑104 (MQ=255)
                                                                                                                     tCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAgcg                                                                         <  4:364748/104‑1 (MQ=255)
                                                                                                                         ttGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTc                           <  6:27523/146‑1 (MQ=255)
                                                                                                                         ttGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGAt                        <  5:118813/149‑1 (MQ=255)
                                                                                                                         ttGAATCCACCGGCCTCCTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCg                                               >  6:210165/1‑126 (MQ=255)
                                                                                                                                 aCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCa                                      >  3:19119/1‑127 (MQ=255)
                                                                                                                                 aCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCa                                      <  4:19119/127‑1 (MQ=255)
                                                                                                                                 aCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTAcc                 >  2:195785/1‑148 (MQ=255)
                                                                                                                                  ccGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGa               <  5:110311/149‑1 (MQ=255)
                                                                                                                                              cgGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGACAACCGcagcag   <  7:95732/149‑1 (MQ=255)
                                                                                                                                               gggCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATc                              <  3:84019/121‑1 (MQ=255)
                                                                                                                                               gggCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATc                              >  4:84019/1‑121 (MQ=255)
                                                                                                                                               gggCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGACAACCGCAGCAgg  <  3:389766/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
TAGCCGCCGAAGGCAACGGCACCAACGACGAAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGACAACCGCAGCAGG  >  NC_002947/5212327‑5212617

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: