Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,327,297 A→T 9.6% N74I (AAC→ATC) ‡ PP_3797 → membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,327,2970AT 9.6% 184.7 / 10.0 63N74I (AAC→ATC) ‡PP_3797membrane protein
Reads supporting (aligned to +/- strand):  ref base A (25/31);  new base T (3/3);  total (28/35)
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.02e-01

TTCTTTGTCCTGATCGTCTGTTTCACCGGCATGCTCGCCGTGGTCAGTCAGGAAATCGTCTGGCTCGCCGATCCGGCCGCACGCGCCAACAAGCCAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACATGGTGGTGGAGCGCCTGAGCCAGCCGGATGGTTCGCACTTCGCGGTAAAAGCCAATGTCACCCTGCCTGACGGCACCAGCCCGACGCTGTATGTAAACCCCTACACCGGAACCATTCAGGGCAGGAC  >  NC_002947/4327152‑4327441
                                                                                                                                                 |                                                                                                                                                
ttCTTTGTCCTGATCGTCTGTTTCACCGGCATGCTCGCCGTGGTCAGTCAGGAAATCGTCTGGCTCGCCGATCCGGCCGCACGCGCCAACAAGCCAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAg                                                                                                                                              >  1:296145/1‑150 (MQ=255)
 tCTTTGTCCTGATCGTCTGTTTGGCCGGCATGCTCGCCGTGGTCAGTCAGGAAATCGTCTGGCTCGCCGATCCGGCCGCACGCGCTAACAAGCCAGACGGGGACACCGAGCGCATGAGCTTTCAGCAAGTACTGGAGGGCTTGAACaatg                                                                                                                                             <  2:1235399/150‑3 (MQ=255)
 tCTTTGTCCTGATCGTCTGTTTCACCGGCATGCTCGCCGTGGTCAGTCAGGAAATCGTCTGGCTCGCCGATCCGGCCGCACGCGCCAACAAGCCAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGATGCCTTGAACAAgg                                                                                                                                             <  1:806970/150‑1 (MQ=255)
     tGTCCTGATCGTCTGTTTCACCGGCATGCTCGCCGTGGTCAGTCAGGAAATCGTCTGGCTCGCCGATCCGGCCGCACGCGCCAACAAGCCAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGTTCAAGGcctc                                                                                                                                         <  2:1001080/150‑3 (MQ=255)
          tGATCGTCTGTTTCACCGGCATGCTCGCCGTGGTCAGTCAGGAAATCGTCTGGCTCGCCGATCCGGCCGCACGCGCCAACAAGCCAGACGCGGACACCGAGCGCATGAGCTTTCAGCCGGTACAGGGGGCCTTGAACAAGGCCGAGCCgg                                                                                                                                    >  2:180191/1‑150 (MQ=255)
          tGATCGTCTGTTTCACCGGCATGCTCGCCGTGGTCAGTCAGGAAATCGTCTGGCTCGCCGATCCGGCCGCACGCGCCAACAAGCCAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCgg                                                                                                                                    >  1:379908/1‑150 (MQ=255)
           gATCGTCTGTTTCACCGGCATGCTCGCCGTGGTCAGTCAGGAAATCGTCTGGCTCGCCGATCCGGCCGCACGCGCCAACAAGCCAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGa                                                                                                                                   >  1:1312106/1‑150 (MQ=255)
                tCTGTTTCACCGGCATGCTCGCCGTGGTCAGTCAGGAAATCGTGTGGCTCGCCGATCCGGCCGCACGCGCCAACAAGCCAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACAtgg                                                                                                                              <  1:866342/150‑1 (MQ=255)
                    tttCACCGGCATGCTCGCCGTGGTCAGTCAGGAAATCGTCTGGCTCGCCGATCCGGCCGCACGCGCCAACAAGCCAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACAtggtggt                                                                                                                          >  1:1417109/1‑150 (MQ=255)
                     ttcagtggCATGGTCGTCCGGGGCAGTCAGGAAGTTGTCTAGCGCGCGACGCCGTCCCCACGCGCCAACACGTCAGACGCGGACACCGAGCGCATGAGCTGTCAGCAGGTACTGGAGGCCTTGACCAAGGCAAAGCCGGACAtggtgatg                                                                                                                         <  2:277780/144‑1 (MQ=255)
                       cACCGGCATGCTCGCCGTGGTCAGTCAGGAAATCGTCTGGCTCGCCGATCCGGCCGCACGCGCCAACAAGCCAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACATGGTGGTGGa                                                                                                                       <  1:915076/150‑1 (MQ=255)
                           ggCATGCTCGCCGTGGTCAGTCAGGAAATCGTCTGGCTCGCCGATCCGGCCGCACGCGCCAACAAGCCAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACATGGTGGTGGAgcgc                                                                                                                   <  1:940448/150‑1 (MQ=255)
                               tGCTCGCCGTGGTCAGTCAGGCAATCGTATGGCTGTCCTAGCCGTTCGCACGCGCCAACAAGCCAGACGCGGACACCGGGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACATGGTGGTGGAGCGCCTGa                                                                                                               <  2:12998/150‑1 (MQ=255)
                               tGCTCGCCGTGGTCAGTCAGGAAATCGGCTGGCTCGCCGATCCGGCCGCACGCGCCAACAAGCCAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACATGGTGGTGGAGCGCCTGa                                                                                                               <  1:229817/150‑1 (MQ=255)
                                    gCCGTGGTCAGTCAGGAAATCGTCTGGCTCGCCGATCCGGCCGCACGCGCCAACAAGCCAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACATGGTGGTGGAGCGCCTGagccag                                                                                                          >  1:241006/1‑150 (MQ=255)
                                           tcagtcagGAAATCGTCTGGCTCGCCGATCCGGCCGCACGCGCCAACAAGCCAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACATGGTGGTGGAGCGCCTGAGCCAGCCGGATg                                                                                                   <  2:241006/150‑1 (MQ=255)
                                              gtcagGAAATCGTCTGGCTCGCCGATCCGGCCGCACGCGCCAACAAGCCAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGTTCAAGGc                                                                                                                                            <  1:124699/106‑1 (MQ=255)
                                              gtcagGAAATCGTCTGGCTCGCCGATCCGGCCGCACGCGCCAACAAGCCAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGTTCAAGGc                                                                                                                                            >  2:124699/1‑106 (MQ=255)
                                               tcagGCAATCGTCTGGCTCGCCGATCCGGCCGCACGCGCCAACAAGCCAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGcc                                                                                                                                           <  1:1155835/106‑1 (MQ=255)
                                               tcagGAAATCGTCTGGCTCGCCGATCCGGCCGCACGCGCCAACAAGCCAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGcc                                                                                                                                           >  2:1155835/1‑106 (MQ=255)
                                                  ggAAATCGTCTGGCTCGCCGATCCGGCCGCACGCGCCAACAAGCCAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACATGGTGGTGGAGCGCCTGAGCCAGCCGGATGGTTCGCa                                                                                            >  2:1116617/1‑150 (MQ=255)
                                                          tCTGGCTCGCCGATCCGGCCGCACGCGCCAACAAGCCAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACATGGTGGTGGAGCGCCTGAGCCAGCCGGATGGTTc                                                                                               <  1:516022/139‑1 (MQ=255)
                                                          tCTGGCTCGCCGATCCGGCCGCACGCGCCAACAAGCCAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACATGGTGGTGGAGCGCCTGAGCCAGCCGGATGGTTc                                                                                               >  2:516022/1‑139 (MQ=255)
                                                              gCTCGCCGATCCGGCCGCACGCGCCAACAAGCCAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACATGGTGGTGGAGCGCCTGAGCCAGCCGGATg                                                                                                   <  1:61492/131‑1 (MQ=255)
                                                              gCTCGCCGATCCGGCCGCACGCGCCAACAAGCCAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACATGGTGGTGGAGCGCCTGAGCCAGCCGGATg                                                                                                   >  2:61492/1‑131 (MQ=255)
                                                               cTCGCCGATCCGGCCGCACGCGCCAACAAGCCAGACGCGGACACCGCGCGCATGAGGGTGCAGCACGTCCTGTCTGCCGTGTACAAGGCCGAGCCGGACATGGTGGTGGAGCGCCTGAGCCAGCCGGATGGTTCGCACTTCGCGGTaaaa                                                                               <  2:337858/150‑1 (MQ=255)
                                                                tCGCCGATCCGGCCGCACGCGCCAACAAGCCAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACATGGTGGTGGAGCGCCTGAACAAGGCGGATGGTACGGACTACGCGGAaaaac                                                                              >  2:1014497/1‑149 (MQ=255)
                                                                      aTCCGGCCGCACGCGCCAACAAGCCAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGTTCAAGGcctcc                                                                                                                                        >  1:1217004/1‑83 (MQ=255)
                                                                      aTCCGGCCGCACGCGCCAACAAGCCAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGTTCAAGGcctcc                                                                                                                                        <  2:1217004/86‑4 (MQ=255)
                                                                          ggCCGCACGCGCCAACAAGCCAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACATGGTGGTGGAGCGCCTGAGCCAGCCGGATGGTTCGCACTTCGCGGTAAAAGCCAATGTCAc                                                                    >  1:984341/1‑150 (MQ=255)
                                                                           gCCGCACGCGCCAACAAGCCAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACATGGTGGTGGAGCGCCTGAGCCAGCCGGATGGTTCGCACTTCGCGGTAAAAGCCAATGTCAcc                                                                   <  2:1312106/150‑1 (MQ=255)
                                                                                 cgcgCCAACAAGCCAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACATGGTGGTGGAGCGCCTGAGCCAGCCGGATGGTTCGCACTTCGCGGTAAAAGCCAATGTCACcctgcct                                                             <  1:180191/150‑1 (MQ=255)
                                                                                  gcgcCAACAAGCCAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCg                                                                                                                                     >  1:1393790/1‑77 (MQ=255)
                                                                                  gcgcCAACAAGCCAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCg                                                                                                                                     <  2:1393790/77‑1 (MQ=255)
                                                                                    gcCAACAAGCCAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACAtggtg                                                                                                                            <  2:180936/84‑1 (MQ=255)
                                                                                    gcCAACAAGCCAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACAtggtg                                                                                                                            >  1:180936/1‑84 (MQ=255)
                                                                                     ccAACAAGCCAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACATGGTGGTGGAGCGCCTGAGCCAGCCGGATGGTTCGCACTTCGCGGTAAAAGCCAATGTCACCCTGCCTGACg                                                         <  1:1331538/150‑1 (MQ=255)
                                                                                      caacaaGCCAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACATGGTGGTGGAGCGCCTGAGCCAGCCGGATGGTTCGCACTTCGCgg                                                                                    >  1:119217/1‑122 (MQ=255)
                                                                                      caacaaGCCAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACATGGTGGTGGAGCGCCTGAGCCAGCCGGATGGTTCGCACTTCGCgg                                                                                    <  2:119217/122‑1 (MQ=255)
                                                                                         caaGCCAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACATGGTGGTGGAGCGCCTGAGCCAGCCGGATGGTTCGCACTTCGCGGTAAAAGCCAATGTCACCCTGCCTGACGGCAc                                                     <  1:1014497/150‑1 (MQ=255)
                                                                                            gCCAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGACGGACATGGTGGTGGAGCGCCTGAGCCAGCCGGATGGTTCGCACTTCGCGGTAAAAGCCAATGTCACCCTGCCTGACGGCACCAg                                                  <  2:379908/150‑1 (MQ=255)
                                                                                             ccAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACATGGTGGTGGAGCGCCTGAGCCAGCCGGATGGTTCGCACTTCGCGGTAAAAGCCAATGTCACCCTGCCTGACGGCACCAGc                                                 >  2:994952/1‑150 (MQ=255)
                                                                                              cAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACATGGTGGTGGAGCGCCTGAGCCAGCCGGATGGTTCGCACTTCGCGGTAAAAGCCAATGTCACCCTGCCTGACGGCACCAGcc                                                <  1:1116617/150‑1 (MQ=255)
                                                                                                      gACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACATGGTGGTGGAGCGCCTGAGCCAGCCGGATGGTTCGCACTTCGCGGTAAAAGCCAATGTCACCCTGCCTGACGGCACCAGCCCGACGCtg                                        >  1:1213875/1‑150 (MQ=255)
                                                                                                         acCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACATGGTGGTGGAGCGCCTGACACGGAGAAATGTTTCGCACTTCGCGGTAAAAGCCAATGTCACCCTCCCTGACGGACCCAGGCCGCCGCtttat                                     >  1:88045/1‑150 (MQ=255)
                                                                                                           cGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACATGGTGGTGGAGCGCCTGAGCCAGCCGGATGGTTCGCACTTCGCGGTAAAAGCCAAt                                                                         <  2:218463/112‑1 (MQ=255)
                                                                                                           cGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACATGGTGGTGGAGCGCCTGAGCCAGCCGGATGGTTCGCACTTCGCGGTAAAAGCCAAt                                                                         >  1:218463/1‑112 (MQ=255)
                                                                                                               cgcATGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACATGGTGGTGGAGCGCCTGAGCCAGCCGGATGGTTCGCACTTCGCGGTAAAAGCCAATGTCACCCTGCCTGACGGCACCAGCCCGACGCTGTATGTAAAc                               <  1:994952/150‑1 (MQ=255)
                                                                                                                  aTGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACATGGTGGTGGAGCGCCTGAGCCAGCCGGATGGTTCGCACTTCGCGGTAAAAGCCAATGTCACCCTGCCTGACGGCACCa                                                   <  2:64699/127‑1 (MQ=255)
                                                                                                                  aTGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACATGGTGGTGGAGCGCCTGAGCCAGCCGGATGGTTCGCACTTCGCGGTAAAAGCCAATGTCACCCTGCCTGACGGCACCa                                                   >  1:64699/1‑127 (MQ=255)
                                                                                                                  aTGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACATGGTGGTGGAGCGCCTGAGCCAGCCGGATGGTTCGCACTTCGCGGTAAAAGCCAATGTCACCCTGCCTGACGGCACCACACCTAGGTTGGAGGTAAAcccc                            <  2:1417109/150‑1 (MQ=255)
                                                                                                                        tttCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACATGGTGGTGGAGCGCCTGAGCCAGCCGGATGGTTCGCACTTCGCGGTAAAAGCCAATGTCACCCTGCCTGACGGCAc                                                     <  1:774204/119‑1 (MQ=255)
                                                                                                                        tttCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACATGGTGGTGGAGCGCCTGAGCCAGCCGGATGGTTCGCACTTCGCGGTAAAAGCCAATGTCACCCTGCCTGACGGCAc                                                     >  2:774204/1‑119 (MQ=255)
                                                                                                                        tttCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACATGGTGGTGGAGCGCCTGAGCCAGCCGGATGGTTCGCACTTCGCGGTAAAAGCCAATGTCACCCTGCCTGACGGCACCAGCCCGACGCtgtat                                     <  1:70748/135‑1 (MQ=255)
                                                                                                                        tttCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACATGGTGGTGGAGCGCCTGAGCCAGCCGGATGGTTCGCACTTCGCGGTAAAAGCCAATGTCACCCTGCCTGACGGCACCAGCCCGACGCtgtat                                     >  2:70748/1‑135 (MQ=255)
                                                                                                                            agcagGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACATGGTGGTGGAGCGCCTGAGCCAGCCGGATGGTTCGCACTTCGCGGTAAACCCCAATGTCACCATGATTGATGGGTCCATCCCGACGCTGTATGTAAACCCCTACACCCGaa                  <  2:49692/150‑1 (MQ=255)
                                                                                                                                 gTACTGGAGGCCTTGAACAAGGCCGAGCCGGACATGGGGGTGGAGCGCCTGAGCCAGCCGGATGGTTCGCACTTCGCGGTAAAAGCCAATGTCACCCTGCCTGACGGCACCAGCCCGACGCTGTATGTAAACCCCTACACCGGAACCAtt             >  1:609041/1‑150 (MQ=255)
                                                                                                                                       tcggCCTTGTTCAAGGCCGAGCCGGACATGGTGGTGGAGCGCCTGAGCCAGCCGGATGGTTCGCACTTCGCGGTAAAAGCCAATGTCACCCTGCCTGACGGCACCAGCCCGACGCTGTATGTAAACCCCTACACCGGAACCATTCAGGGc       >  1:1105012/3‑150 (MQ=255)
                                                                                                                                        aGGCCTTGAACAAGGCCGAGCCGGACATGGTGGTGGAGCGCCTGAGCCAGCCGGATGGTTCGCACTTCGCGGTAAAAGCCAATGTCACcctgcct                                                             >  1:668140/1‑95 (MQ=255)
                                                                                                                                        aGGCCTTGAACAAGGCCGAGCCGGACATGGTGGTGGAGCGCCTGAGCCAGCCGGATGGTTCGCACTTCGCGGTAAAAGCCAATGTCACcctgcct                                                             <  2:668140/95‑1 (MQ=255)
                                                                                                                                        aGGCCTTGAACAAGGCCGAGCCGGACATGGTGGTGGAGCGCCTGAGCCAGCCGGATGGTTCGCACTTCGCGGTAAAAGCCAATGTCACCCTGCCTGACGGCACCAGCCCGACGCTGTATGTAAACCCCTACACCGGAACCATTCAGGGCa      <  2:296145/150‑1 (MQ=255)
                                                                                                                                        aGGCCTTGAACAAGGCCGAGCCGGACATGGTGGTGGAGCGCCTGAGCCAGCCGGATGGTTCGCACTTCGCGGTAAAAGCCAATGTCACCCTGCCTGACGGCACCAGCCCGACGCTGTATGTAAACCCCTACACCGGAACCATTCAGGGCa      <  1:1130815/150‑1 (MQ=255)
                                                                                                                                            cTTGAACAAGGCCGAGCCGGACATGGTGGTGGAGCGCCTGAGCCAGCCGGATGGTTCGCACTTCGCGGTAAAAGCCAATGTCACCCTGCCTGACGGCACCAGCCCGACGCTGTATGTAAACCCCTACACCGGAACCATTCAGGGCAGGAc  >  2:25047/1‑150 (MQ=255)
                                                                                                                                                aaCAAGGCCGAGCCGGACATGGTGGTGGAGCGCCTGAGCCAGCCGGATGGTTCGCACTTCGCGGTAAAAGCCAATGTCACcctgcc                                                              >  2:889799/1‑86 (MQ=255)
                                                                                                                                                aaCAAGGCCGAGCCGGACATGGTGGTGGAGCGCCTGAGCCAGCCGGATGGTTCGCACTTCGCGGTAAAAGCCAATGTCACcctgcc                                                              <  1:889799/86‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
TTCTTTGTCCTGATCGTCTGTTTCACCGGCATGCTCGCCGTGGTCAGTCAGGAAATCGTCTGGCTCGCCGATCCGGCCGCACGCGCCAACAAGCCAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACATGGTGGTGGAGCGCCTGAGCCAGCCGGATGGTTCGCACTTCGCGGTAAAAGCCAATGTCACCCTGCCTGACGGCACCAGCCCGACGCTGTATGTAAACCCCTACACCGGAACCATTCAGGGCAGGAC  >  NC_002947/4327152‑4327441

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