Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,848,876 C→T 6.5% A106V (GCC→GTC)  cysM → cysteine synthase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,848,8760CT 6.5% 181.3 / 2.0 61A106V (GCC→GTC) cysMcysteine synthase B
Reads supporting (aligned to +/- strand):  ref base C (31/26);  new base T (2/2);  total (33/28)
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.98e-01

GAATTGCGGGGCCAGATCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGAAGGCGCCCGTGACCTGGCCGAGAAGCTGCAGGCCGAAGGCCGTGGCCTGGTGCTCGACCAGTTCGCCAACGGCGACAACCCGATCGCCCACTAC  >  NC_002947/1848728‑1849024
                                                                                                                                                      |                                                                                                                                                    
gAATTGCGGGGCCAGATCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGgccgc                                                                                                                                                      <  4:173182/149‑1 (MQ=255)
  tttGCGGGGGCTTAGCAAGCCAGGCGATACCCTGATGGAAGCCACCTCAGGCAACACGGGCATCGCCTTGGCGATGGAGGCGGCGATCAAGGGCTACATGATGATCCTGATCATGCCCGACCACTCCACCGCCGAGCGCACGgccgccc                                                                                                                                                    <  3:47575/148‑2 (MQ=255)
         ggCCAGATCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGc                                                                                                                                             >  3:154222/1‑149 (MQ=255)
           ccAGATCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGTCATGGCGGCCt                                                                                                                                           >  4:297205/1‑149 (MQ=255)
             aGATCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGTCATGGCGGCCttg                                                                                                                                         <  3:297205/149‑3 (MQ=255)
               aTCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATTGCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATg                                                                                                                                        <  1:279309/149‑1 (MQ=255)
                   aGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGcc                                                                                                                                   >  5:46105/1‑149 (MQ=255)
                      cAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAg                                                                                                                                <  6:293641/149‑1 (MQ=255)
                      cAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAg                                                                                                                                >  4:213675/1‑149 (MQ=255)
                       aGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGc                                                                                                                               <  8:174869/149‑1 (MQ=255)
                       aGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGc                                                                                                                               <  5:105651/149‑1 (MQ=255)
                       aGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGc                                                                                                                               <  4:154222/149‑1 (MQ=255)
                          cGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGa                                                                                                                            <  6:1252/149‑1 (MQ=255)
                                ccTGATCGAAGCCACCTCCGGCAACAGGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTgg                                                                                                                      >  2:51360/1‑149 (MQ=255)
                                ccTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTgg                                                                                                                      >  5:52435/1‑149 (MQ=255)
                                ccTGATCGAAGCCACCTCCGGCAACACGGGCAGCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTgg                                                                                                                      >  5:206789/1‑149 (MQ=255)
                                     tCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGAc                                                                                                                  >  2:217916/1‑148 (MQ=255)
                                            cACCTCCGGCAACACGGGCATCGCCCTGGCGGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGc                                                                                                                               >  1:214899/1‑129 (MQ=255)
                                            cACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGc                                                                                                                               <  2:214899/128‑1 (MQ=255)
                                             aCCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAaggag                                                                                                           >  3:81823/1‑147 (MQ=255)
                                                tCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAggg                                                                                                      >  3:136365/1‑149 (MQ=255)
                                                  cGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAgggg                                                                                                     <  2:40412/148‑1 (MQ=255)
                                                    gCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATg                                                                                                  >  2:129526/1‑149 (MQ=255)
                                                    gCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATg                                                                                                  >  3:169514/1‑149 (MQ=255)
                                                      aaCACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACTAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGcc                                                                                                                                   >  1:252319/1‑114 (MQ=255)
                                                      aaCACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGcc                                                                                                                                   <  2:252319/114‑1 (MQ=255)
                                                       acacGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGaa                                                                                               <  1:217916/149‑1 (MQ=255)
                                                               aTCGCCCTGACGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGAAGGCGCCCg                                                                                       <  1:328717/149‑1 (MQ=255)
                                                                tCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGAAGGCGCCCGt                                                                                      <  4:169514/149‑1 (MQ=255)
                                                                    ccTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGAAGGCGCCCGTGAcc                                                                                  >  6:332715/1‑149 (MQ=255)
                                                                    ccTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGAAGGCGCCCGTGAcc                                                                                  >  2:253077/1‑149 (MQ=255)
                                                                        gCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGAAGGCGCCCGTGACCTg                                                                                <  6:297288/147‑1 (MQ=255)
                                                                            aTGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCAGGTGACCAAGGAGGAGGGGATGGAAGGCGCCCGTGACCTGGCCga                                                                           <  1:253077/149‑1 (MQ=255)
                                                                                 ggcggcgATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCTAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGAAGGCGCCCGTGACCTGGCCGAGAAGc                                                                      >  8:312544/1‑149 (MQ=255)
                                                                                       gATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAaggagg                                                                                                          >  2:292879/1‑107 (MQ=255)
                                                                                              ggCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGAAg                                                                                              >  6:224093/1‑112 (MQ=255)
                                                                                              ggCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGAAg                                                                                              <  5:224093/112‑1 (MQ=255)
                                                                                              ggCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGAAGGCGCCCGTGACCTGGCCGAGAAGCTGCAGGCCGAAgg                                                         <  6:46105/149‑1 (MQ=255)
                                                                                                            tCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGAcc                                                                                                                 <  7:48475/80‑1 (MQ=255)
                                                                                                            tCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGAcc                                                                                                                 >  8:48475/1‑80 (MQ=255)
                                                                                                             ccTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGAAGGCGCCCGTGACCTGGCCGAGAAGCTGCAGGCCGAAGGCCGTGGCCTgg                                              <  3:84244/146‑1 (MQ=255)
                                                                                                             ccTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGAAGGCGCCCGTGACCTGGCCGAGAAGCTGCAGGCCGAAGGCCGTGGCCTgg                                              >  4:84244/1‑146 (MQ=255)
                                                                                                             ccTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGAAGGCGCCCGTGACCTGGCCGAGAAGCTGCAGGCCGAAGGCCGTGGCCTGGTGc                                           >  3:41796/1‑149 (MQ=255)
                                                                                                                  tCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGAAGGCGCCCGTGACCTGGCCGAGAAGCTGCAGGCCGAAGGCCGTGGCCTgg                                              >  1:328496/1‑141 (MQ=255)
                                                                                                                  tCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGAAGGCGCCCGTGACCTGGCCGAGAAGCTGCAGGCCGAAGGCCGTGGCCTgg                                              <  2:328496/141‑1 (MQ=255)
                                                                                                                       cccGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGAAGGCGCCCGTGACCTGGCCGAGAAGCTGCAGGCCGAAGGCCGTGGCCTGGTGCTCGACCAGtt                                 >  1:273262/1‑149 (MQ=255)
                                                                                                                           aCAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGAAGGCGCCCGTGACCTGGCCGAGAAGCTGCAGGCCGAAGGCCGTGGCCTGGTGCTCGACCAGTTCGc                              >  5:143764/1‑148 (MQ=255)
                                                                                                                             caCTCCACCGCCGAGCGCACGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGAAGGCGCCCGTGACCTGGCCGAGAAGCTGCAGGCCGAAGGCCGTGGCCTGGTGCTCGACCAGTTCGCCaa                           <  2:291081/148‑1 (MQ=255)
                                                                                                                             aaCTCCACCGCCGAGCGCTAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGAAGGCGCCCGTGACCTGGCCGAGAAGCTGCAGGCCGAAGGCCGTGGCCTGGTGCTCGACCAGTTCGCCaa                           <  7:312544/149‑1 (MQ=255)
                                                                                                                                tCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATTGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGAAGGCGCCAGTGACCTGGCCGAGAAGATGCAGGCCGAAGGCCGTGGCCTGGTGCTCGACCAGTTCGCCAACgg                        >  2:312161/1‑149 (MQ=255)
                                                                                                                                tCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGAAGGCGCCCGTGACCTGGCCGAGAAGc                                                                      >  1:200863/1‑103 (MQ=255)
                                                                                                                                tCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGAAGGCGCCCGTGACCTGGCCGAGAAGc                                                                      <  2:200863/103‑1 (MQ=255)
                                                                                                                                tCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGAAGGCGCCCGTGACCTGGCCGAGAAGCTGCAGGCCGAAGGCCGTGGCCTGGTGCTCGACCAGTTCGCCAACgg                        >  2:114626/1‑149 (MQ=255)
                                                                                                                                  cACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAgggg                                                                                                     <  5:314483/70‑1 (MQ=255)
                                                                                                                                  cACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAgggg                                                                                                     >  6:314483/1‑70 (MQ=255)
                                                                                                                                   aCCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGAAGGCGCCCGTGACCTGGCCGAGAAGCTGCAGGCCGAAGGCCGTGGCCTGGTGCTCGACCAGTTCGCCAACGGCGa                     >  4:186172/1‑149 (MQ=255)
                                                                                                                                      gccgAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGAAGGCGCCCGTGACCTGGCCGAGAAGCTGCAGGCCGAAGGCCGTGGCCTGGTGCTCGACCAGTTCGCCAACGGCGACaa                  >  2:49317/1‑149 (MQ=255)
                                                                                                                                       ccgAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGAAGGCGCCCGTGACCTGGCCGAGAAGCTGCAGGCCGAAGGCCGTGGCCTGGTGCTCGACCAGTTCGCCAACGGCGACAAc                 <  6:206789/149‑1 (MQ=255)
                                                                                                                                        cgAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGAAGGCGCCCGTGACCTGGCCGAGAAGCTGCAGGCCGAAGGCCGTGGCCTGGTGCTCGACCAGTTCGCCAACGGCGACAAcc                <  1:49317/149‑1 (MQ=255)
                                                                                                                                             gcAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGAAGGCGCCCGTGACCTGGCCGAGAAGCTGCAGGCCGAAGGCCGTGGCCTGGTGCTCGACCAGTTCGCCAACGGCGACAACCCGATc           >  6:320216/1‑149 (MQ=255)
                                                                                                                                             gcAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGAAGGCGCCCGTGACCTGGCCGAGAAGCTGCAGGCCGAAGGCCGTGGCCTGGTGCTCGACCAGTTCGCCAACGGCGACAACCCGATc           <  3:186172/149‑1 (MQ=255)
                                                                                                                                             cataGGCCGTCATGGCGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGAAGGCGCCCGTGACCTGGCCGAGAAGCTGCAGGCCGAAGGCCGTg                                                    >  2:321260/4‑108 (MQ=255)
                                                                                                                                             cataGGCCGTCATGGCGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGAAGGCGCCCGTGACCTGGCCGAGAAGCTGCAGGCCGAAGGCCGTg                                                    <  1:321260/105‑1 (MQ=255)
                                                                                                                                                  gccgccATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGAAGGCGCCCGTGACCTGGCCGAGAAGCTGCAGGCCGAAGGCCGTGGCCTgg                                              >  3:230651/1‑109 (MQ=255)
                                                                                                                                                  gccgccATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGAAGGCGCCCGTGACCTGGCCGAGAAGCTGCAGGCCGAAGGCCGTGGCCTgg                                              <  4:230651/109‑1 (MQ=255)
                                                                                                                                                      ccATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGAAGGCGCCCGTGACCTGGCCGAGAAGCTGCAGGCCGAAGGCCGTGGCCTGGTGCTCGACCAGTTCGCCAACGGCGACAACCCGATCGCCCACTac  <  2:273262/149‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
GAATTGCGGGGCCAGATCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGAAGGCGCCCGTGACCTGGCCGAGAAGCTGCAGGCCGAAGGCCGTGGCCTGGTGCTCGACCAGTTCGCCAACGGCGACAACCCGATCGCCCACTAC  >  NC_002947/1848728‑1849024

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: