Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 817,397 G→A 11.6% A257A (GCC→GCT PP_0703 ← transmembrane sensor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947817,3970GA11.6% 154.6 / 9.3 60A257A (GCC→GCTPP_0703transmembrane sensor
Reads supporting (aligned to +/- strand):  ref base G (24/29);  new base A (3/4);  total (27/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.05e-01

CTGTACTCGACCCGCACGGGCAGGCTGCGGGTGAGCAGTTGCAGGATTTGCCCGGTGTCAGCCAACTGATAGGTGCCTGAAACCTTCAAGCCCGCCACTTCTGGAGCGCAGCGCAACAGGCCCGGTCGGTAGCGGCCGAGTTCGCTGGCAAACTCGGCAAGCGGCATCTGCTGCACGCTGAGTACACCTTCGGTCCAACCCCAGGGGTCGTTGTGCAGTGTTCTGGCCTGAATGCTGCCGTCGGCTGTGACTTGCAAGGTTTCACCGGGCAGCACTTGCCTCGCTGAG  >  NC_002947/817251‑817538
                                                                                                                                                  |                                                                                                                                              
cTGTACTCGACCCGCACGGGCAGGCTGCGGGTGAGCAGTTGCAGGATTTGCCCGGTGTCAGCCAACTGATAGGTGCCTGAAACCTTCAAGCCCGCCACTTCTGGAGCGCAGCGCAACAGGCCCGGTCGGTAGCGGCCGAGTTCGCTgg                                                                                                                                               <  3:89223/148‑1 (MQ=255)
cTGTACTCGACCCGCACGGGCAGGCTGCGGGTGAGCAGTTGCAGGATTTGCCCGGTGTCAGCCAACTGATAGGTGCCCGAAACCTTCAAGCCCGCCACTTCTGGAGCGCAGCGCAACAGGCCCGGTCGGTAGCGGCCGAGTTCGCTGGc                                                                                                                                              <  2:45558/149‑1 (MQ=255)
      tCGACCCGCACGGGCAGGCTGCGGGTGAGCAGTTGCAGGATTTGCCCGGTGTCACCCAACTGATAGGTGCCTGAAACCTTCAAGCCCGCCAGTTCTGGAGCGCAGCGCAACAGGCCCGGTCGGTAGCGGCCGAGTTCGCTGGCAAACTc                                                                                                                                        >  1:242125/1‑149 (MQ=255)
        gACCCGCACGGGCAGGCTGCGGGTGAGCAGTTGCAGGATTTGCCCGGTGTCAGCCAACTGATAGGTGCCTGAAACCTTCAAGCCCGCCACTTCTGGAGCGCAGCGCAACAGGCCCGGTCGGTAGCGGCCGAGTTTGCCAGCGAACTc                                                                                                                                        <  1:24158/147‑1 (MQ=255)
        gACCCGCACGGGCAGGCTGCGGGTGAGCAGTTGCAGGATTTGCCCGGTGTCAGCCAACTGATAGGTGCCTGAAACCTTCAAGCCCGCCACTTCTGGAGCGCAGCGCAACAGGCCCGGTCGGTAGCGGCCGAGTTTGCCAGCGAACTc                                                                                                                                        >  2:24158/1‑147 (MQ=255)
        gACCCGCACGGGCAGGCTGCGGGTGAGCAGTTGCAGGATTTGCCCGGTGTCAGCCAACTGATAGGTGCCTGAAACCTTCAAGCCCGCCACTTCTGGAGCGCAGCGCAACAGGCCCGGTCGGTAGCGGCCGAGTTTGCCAGCGAACTCgg                                                                                                                                     <  1:153433/149‑1 (MQ=255)
        gACCCGCACGGGCAGGCTGCGGGTGAGCAGTTGCAGGATTTGCCCGGTATCAGCCAACTGATAGGTGCCTGAAACCTTCAAGCCCGCCACTTCTGGAGCGCAGCGCAACAGGCCCGGTCGGTAGCGGCCGAGTTTGCCAGCGAACTCgg                                                                                                                                     <  2:71463/149‑1 (MQ=255)
        gACCCGCACGGGCAGGCTGCGGGTGAGCAGTTGCAGGATTTGCCCGGTATCAGCCAACTGATAGGTGCCTGAAACCTTCAAGCCCGCCACTTCTGGAGCGCAGCGCAACAGGCCCGGTCGGTAGCGGCCGAGTTTGCCAGCGAACTCgg                                                                                                                                     >  1:71463/1‑149 (MQ=255)
           ccGCACGGGCAGGCTGCGGGTGAGCAGTTGCAGGATTTGCCCGGTGTCAGCCAACTGATAGGTGCCTGAAACCTTCAAGCCCGCCACTTCTGGAGCGCAGCGCAACAGGCCCGGTCGGTAGCGGCCGAGTTCGCTGGCAAACTCGGc                                                                                                                                    >  1:138166/1‑147 (MQ=255)
             gCACGGGCAGGCTGCGGGTGAGCAGTTGCAGGATTTGCCCGGTGTCAGCCAACTGATAGGTGCCTGAAACCTTCAAGCCCGCCACTTCTGGAGCGCAGCGCAACAGGCCCGGTCGGTAGCGGCCGAGTTCGCTGGCAAACTCGGCAAGc                                                                                                                                >  3:125577/1‑149 (MQ=255)
             gCACGGGCAGGCTGCGGGTGAGCAGTTGCAGGATTTGCCCGGTGTCAGCCAACTGATAGGTGCCTGAAACCTTCAAGCCCGCCACTTCTGGAGCGCAGCGCAACAGGCCCGGTCGGTAGCGGCCGAGTTCGCTGGCAAACTCGGCAAGc                                                                                                                                <  5:118590/149‑1 (MQ=255)
                        cTGCGGGTGAGCAGTTGCAGGATTTGCCCGGTGTCAGCCAACTGATAGGTGCCTGAAACCTTCAAGCCCGCCACTTCTGGAGCGCAGCGCAACAGGCCCGGTCGGTAGCGGCCGAGTTCGCTGGCAAACTCGGCAAGCGGCATctgctg                                                                                                                     <  4:29886/149‑1 (MQ=255)
                            gggTGAGCAGTTGCAGGATTTGCCCGGTGTCAGCCAACTGATAGGTGCCTGACACCTTCAAGCCCGCCACTTCTGGAGCGAAGCGCAACACGCCAGGTCGGTAGCGGCCGAGTTCGCTGGCAAACTCGGCAAGCGGCATCTGCTGCaag                                                                                                                 >  5:99910/1‑147 (MQ=255)
                            gggTGAGCAGTTGCAGGATTTGCCCGGTGTCAGCCAACTGATAGGTGCCTGAAACCTTCAAGCCCGCCACTTCTGGAGCGCAGCGCAACAGGCCCGGTCGGTAGCGGCCGAGTTCGCTGGCAAACTCGGCAAGCGGCATCTGCTGCACg                                                                                                                 >  7:232554/1‑149 (MQ=255)
                               tGAGCAGTTGCAGGATTTGCCCGGTGTCAGCCAACTGATAGGTGCCTGAAACCTTCAAGCCCGCCACTTCTGGAGCGCAGCGCAACAGGCCCGGTCGGTAGCGGCCGAGTTCGCTGGCAAACTCGGCAAGCGGCATCTGCTGCACGCTg                                                                                                              <  1:114887/149‑1 (MQ=255)
                                      ttGCAGGATTTGCCCGGTGTCAGCCAACTGATAGGTGCCTGAAACCTTCAAGCCCGCCACTTCTGGAGCGCAGCGCAACAGGCCCGGTCGGTAGCGGCCGAGTTTGCCAGCGAACTCGGc                                                                                                                                    <  2:32168/120‑1 (MQ=255)
                                      ttGCAGGATTTGCCCGGTGTCAGCCAACTGATAGGTGCCTGAAACCTTCAAGCCCGCCACTTCTGGAGCGCAGCGCAACAGGCCCGGTCGGTAGCGGCCGAGTTTGCCAGCGAACTCGGc                                                                                                                                    >  1:32168/1‑120 (MQ=255)
                                        gCAGGATTTGCCCGGTGTCAGCCAACTGATAGGTGCCTGAAACCTTCAAGCCCGCCACTTCTGGAGCGCAGCGCAACAGGCCCGGTCGGTAGCGGCCGAGTTCGCTGGCAAACTCGGCAAGCGGCATCTGCTGCACGCTGAGTACACCt                                                                                                     >  3:163075/1‑149 (MQ=255)
                                         cAGGATTTGCCCGGTGTCAGCCAACTGATAGGTGCCTGAAACCTTCAAGCCCGCCACTTCTGGAGCGCAGCGCAACAGGCCCGGTCGGTAGCGGCCGAGTTCGCTGGCAAACTCgg                                                                                                                                     <  5:68388/116‑1 (MQ=255)
                                         cAGGATTTGCCCGGTGTCAGCCAACTGATAGGTGCCTGAAACCTTCAAGCCCGCCACTTCTGGAGCGCAGCGCAACAGGCCCGGTCGGTAGCGGCCGAGTTCGCTGGCAAACTCgg                                                                                                                                     >  6:68388/1‑116 (MQ=255)
                                         cAGGACTTGCCCGGTGGCAGCCAACTGATAGGGGCCTGACACCTTCCAGCCCGCCACTTCTGGAGCGCAGCGCAACAGGCCCGGTCGGGAGCGGCCGAGTTCGCTGGCAAACTCGGCAAGCGGCATCTGCTGCACGCTGAGTACACCtt                                                                                                    <  6:102129/149‑1 (MQ=255)
                                              tttGCCCGGTGTCAGCCAACTGATAGGTGCCTGAAACCTTCAAGCCCGCCACTTCTGGAGCGCAGCGCAACAGGCCCGGTCGGTAGCGGCCGAGTTCGCTGGCAAACTCGGCAAGCGGCATCTGCTGCACGCTGAGTACACCTTCGGTc                                                                                               <  7:278472/149‑1 (MQ=255)
                                              tttGCCCGGTGTCAGCCAACTGATAGGTGCCTGAAACCTACAAGACCGCCACTTCTGGAGCGCAGCGCAACAGGCCCGGTCGGTAGCGGCCGAGTTCGCTGGCAAACTCGGCAAGCGGCATCTGCTGCACGCTGAGTACACCTTCGGTc                                                                                               >  7:216463/1‑149 (MQ=255)
                                                    cGGTGTCAGCCAACTGATAGGTGCCTGAAACCTTCAAGCCCGCCACTTCTGGAGCGCAGCGCAACAGGCCCGGTCGGTAGCGGCCGAGTTCGCTGGCAAACTCGGCAAGCGGCATCTGCTGCACGCTGAGTacac                                                                                                       >  3:159596/1‑135 (MQ=255)
                                                    cGGTGTCAGCCAACTGATAGGTGCCTGAAACCTTCAAGCCCGCCACTTCTGGAGCGCAGCGCAACAGGCCCGGTCGGTAGCGGCCGAGTTCGCTGGCAAACTCGGCAAGCGGCATCTGCTGCACGCTGAGTacac                                                                                                       <  4:159596/135‑1 (MQ=255)
                                                    cGGTGTCAGCCAACTGATAGGTGCCTGAAACCTTCAAGCCCGCCACTTCTGGAGCGCAGCGCAACAGGCCCGGTCGGTAGCGGCCGAGTTCGCTGGCAAACTCGGCAAGCGGCATCTGCTGCACGCTGAGTACACCTTCGGTCCAAccc                                                                                         >  7:252939/1‑149 (MQ=255)
                                                         tCAGCCATCTGTTAGGTGCCTGAAACCTTCAAGCCCGCCACTTCTGGAGCGCAGCGCATCAGGCCCGGTCGGTAGCGGCCGAGTTCGCTGGCAAACTCGGCAAGCGGCATCTGCTGCACGCTGAGTACACCTTCGGTCCAACCCCAgg                                                                                     <  8:216463/148‑1 (MQ=255)
                                                               aaCTGATAGGTGCCTGAAACCTTCAAGCCCGCCACTTCTGGAGCGCAGCGCAACAGGCCCGGTCGGTAGCGGCCGAGTTCGCTGGCAAACTCGGCAAGCGGCATCTGCTGCACGCTGAGTACACCTTCGGTCCAACCCCAGGGGTCGtt                                                                              <  4:125577/149‑1 (MQ=255)
                                                               aaCTGATAGGTGCCTGAAACCTTCAAGCCCGCCACTTCTGGAGCGCAGCGCAACAGGCCCGGTCGGTAGCGGCCGAGTTCGCTGGCAAACTCGGCAAGCGGCATCTGCTGCACGCTGAGTACACCTTCGGTCCAACCCCAGGGGTCGtt                                                                              <  2:138166/149‑1 (MQ=255)
                                                                  tGATAGGTGCCTGAAACCTTCAAGCCCGCCACTTCTGGAGCGCAGCGCAACAGGCCCGGTCGGTAGCGGCCGAGTTCGCTGGCAAACTCGGCAAGCGGCATCTGCTGGACGCTGAGTACACCTTCGGTCCAACCCCAGGGGTCGTtgtg                                                                           >  5:275847/1‑149 (MQ=255)
                                                                   gATAGGTGCCTGAAACCTTCAAGCCCGCCACTTCTGGAGCGCAGCGCAACAGGCCCGGTCGGTAGCGGCCGAGTTCGCTGGCAAACTCGGCAAGCGGCATCTGCTGCACGCTGAGTACACCTTCGGTCCAACCCCAGGGGTCGTtgtg                                                                           >  3:97260/1‑148 (MQ=255)
                                                                    aTAGGTGCCTGAAACCTTCAAGCCCGCCACTTCTGGAGCGCAGCGCAACAGGCCCGGTCGGTAGCGGCCGAGTTCGCTGGCAAACTCGGCAAGCGGCATCTGCTGCACGCTGAGTACACCTTCGGTCCAACCCCAGGGGTCGTTGTTCa                                                                         <  2:242125/149‑1 (MQ=255)
                                                                    aTAGGTGCCTGAAACCTTCAAGCCCGCCACTTCTGGAGCGCAGCGCAACAGGCCCGGTCGGTAGCGGCCGAGTTCGCTGGCAAACTCGGCAAGCGGCATCTGCTGCACGCTGAGTACACCTTCGGTCCAACCCCAGGGGTCGTTGTGCa                                                                         <  3:293171/149‑1 (MQ=255)
                                                                           ccTGAAACCTTCAAGCCCGCCACTTCTGGAGCGCAGCGCAACAGGCCCGGTCGGTAGCGGCCGAGTTCGCTGGCAAACTCGGCAAGCGGCATCTGCTGCACGCTGAGTACACCTTCGGTCCAACCCCAGGGGTCGTtgtg                                                                           <  6:176663/140‑1 (MQ=255)
                                                                           ccTGAAACCTTCAAGCCCGCCACTTCTGGAGCGCAGCGCAACAGGCCCGGTCGGTAGCGGCCGAGTTCGCTGGCAAACTCGGCAAGCGGCATCTGCTGCACGCTGAGTACACCTTCGGTCCAACCCCAGGGGTCGTtgtg                                                                           >  5:176663/1‑140 (MQ=255)
                                                                             tGAAACCTTCAAGCCCGCCACTTCTGGAGCGCAGCGCAACAGGCCCGGTCGGTAGCGGCCGAGTTCGCTGGCAAACTCGGCAAGCGGCATCTGCTGCACGCTGAg                                                                                                            >  6:16397/1‑105 (MQ=255)
                                                                             tGAAACCTTCAAGCCCGCCACTTCTGGAGCGCAGCGCAACAGGCCCGGTCGGTAGCGGCCGAGTTCGCTGGCAAACTCGGCAAGCGGCATCTGCTGCACGCTGAg                                                                                                            <  5:16397/105‑1 (MQ=255)
                                                                             tGAAACCTTCAAGCCCGCCACTTCTGGAGCGCAGCGCAACAGGCCCGGTCGGTAGCGGCCGAGTTCGCTGGCAAACTCGGCAAGCGGCATCTGCTGCACGCTGAGTACACCTTCGGTCCAACCCCAg                                                                                      <  1:195896/127‑1 (MQ=255)
                                                                             tGAAACCTTCAAGCCCGCCACTTCTGGAGCGCAGCGCAACAGGCCCGGTCGGTAGCGGCCGAGTTCGCTGGCAAACTCGGCAAGCGGCATCTGCTGCACGCTGAGTACACCTTCGGTCCAACCCCAg                                                                                      >  2:195896/1‑127 (MQ=255)
                                                                                   cTTCAAGCCCGCCACTTCTGGAGCGCAGCGCAACAGGCCCGGTCGGTAGCGGCCGAGTTCGCTGGCAAACTCGGCAAGCGGCATCTGCTGCACGCTGAGTACACCTTCGGTCCAACCCCAGGGGTCGTTGTGCAGTGTTCTGGCCTGaa                                                          <  8:232554/149‑1 (MQ=255)
                                                                                    ttCAAGCCCGCCACTTCTGGAGCGCAGCGCAACAGGCCCGGTCGGTAGCGGCCGAGTTCGCTGGCAAACTCGGCAAGCGGCATCTGCTGCACGCTGAGTACACCTTCGGTCCAACCCCAGGGGTCGTTGTGCAGTGTTCTGGCCTGAAt                                                         >  8:31972/1‑149 (MQ=255)
                                                                                       aaGCCCGCCACTTCTGGAGCGCAGCGCAACAGGCCCGGTCGGTAGCGGCCGAGTTCGCTGGCAAACTCGGCAAGCGGCATCTGCTGCACGCTGAGTACACCTTCGGTCCAACCCCAGGGGTCGTTGTGCAGTGTTCTGGCCTGAAtgc                                                       <  4:97260/148‑1 (MQ=255)
                                                                                       aaGCCCGCCACTTCTGGAGCGCAGCGCAACAGACCCGGTCGGTAGCGGCCGAGTTCGCTGGCAAACTCGGCAAGCGGCATCTGCTGCACGCTGAGTACACCTTCGGTCCAACCCCAGGGGTCGTTGTGCAGTGTTCTGGCCTGAAtgc                                                       <  6:253772/148‑1 (MQ=255)
                                                                                        aGCCCGCCACTTCTGGAGCGCAGCGCAACAGGCCCGGTCGGTAGCGGCCGAGTTCGCTGGCAAACTCGGCAAGCGGCATCTGCTGCACGCTGAGTACACCTTCGGTCCAACCCCAGGGGTCGTTGTGCAGTGTTCTGGCCTGAAtgctg                                                     <  6:275847/149‑1 (MQ=255)
                                                                                        aGCCCGCCACTTCTGGAGCGCAGCGCAACAGGCCCGGTCGGTAGCGGCCGAGTTCGCTGGCAAACTCGGCAAGCGGCATCTGCTGCACGCTGAGTACACCTTCGGTCCAACCCCAGGGGTCGTTGTGCAGTGTTCTGGCCTGAAtgctg                                                     >  2:231337/1‑149 (MQ=255)
                                                                                          cccGCCACTTCTGGAGCGCAGCGCAACAGGCCCGGTCGGTAGCGGCCGAGTTCGCTGGCAAACTCGGCAAGCGGCATCTGCTGCACGCTGAGTACACCTTCGGTCCAACCCCAGGGGTCGTTGTGCAGTGTTCTGGCCTGAATGCTGcc                                                   >  2:45911/1‑149 (MQ=255)
                                                                                           ccGCCACTTCTGGAGCGCAGCGCAACAGGCCCGGTCGGTAGCGGCCGAGTTCGCTGGCAAACTCGGCAAGCGGCATCTGCTGCACGCTGAGTACACCTTCGGTCCaa                                                                                            >  5:27769/1‑107 (MQ=255)
                                                                                           ccGCCACTTCTGGAGCGCAGCGCAACAGGCCCGGTCGGTAGCGGCCGAGTTCGCTGGCAAACTCGGCAAGCGGCATCTGCTGCACGCTGAGTACACCTTCGGTCCaa                                                                                            <  6:27769/107‑1 (MQ=255)
                                                                                                 cTTCTGGAGCGCAGCGCAACAGGCCCGGTCGGTAGCGGCCGAGTTCGCTGGCAAACTCGGCAAGCGGCATCTGCTGCACGc                                                                                                                >  6:156399/1‑81 (MQ=255)
                                                                                                 cTTCTGGAGCGCAGCGCAACAGGCCCGGTCGGTAGCGGCCGAGTTCGCTGGCAAACTCGGCAAGCGGCATCTGCTGCACGc                                                                                                                <  5:156399/81‑1 (MQ=255)
                                                                                                 cTTCTGGAGCGCAGCGCAACAGGCCCGGTCGGTAGCGGCCGAGTTCGCTGGCAAACTCGGCAAGCGGCATCTGCTGCACGCTGAGTACACCTTCGGTCCAACCCCAGGGGTCGTTGTGCAGTGTTCTGGCCTGAATGCTGCCGTCGGCt                                            <  4:163075/149‑1 (MQ=255)
                                                                                                       gAGCGCAGCGCAACAGGCCCGGTCGGTAGCGGCCGAGTTCGCTGGCAAACTCGGCAAGCGGCATCTGCTGCACGCTGAGTACACCTTCGGTCCAACCCCAGGGGTCGTTGTGCAGTGTTCTGGCCTGAATGCTGCCGTCGGCTGTGACt                                      >  5:47485/1‑149 (MQ=255)
                                                                                                         gcgcagcgcaACAGGCCCGGTCGGTAGCGGCCGAGTTCGCTGGCAAACTCGGCAAGCGGCATCTGCTGCACGCTGAGTACACCTTCGGTCCAACCCCAGGGGTCGTTGTGCAGTGTTCTGGCCTGAATGCTGCCGTCGGCTGTGACTTg                                    <  8:252939/149‑1 (MQ=255)
                                                                                                             agcgcaACAGGCCCGGTCGGTAGCGGCCGAGTTCGCTGGCAAACTCGAGCAAGCGGCATCTGCTGCACGCTGAGTACACCtt                                                                                                    <  1:42338/82‑1 (MQ=255)
                                                                                                                         ccGGTCGGTAGCGGCCGAGTTCGCTGGCAAACTCGGCAAGCGGCATCTGCTGCACGCTGAGTACACCTTCGGTCCAACCCCAGGGGTCGTTGTGCAGTGTTCTGGCCTGAATGCTGCCGTCgg                                              <  6:135673/123‑1 (MQ=255)
                                                                                                                         ccGGTCGGTAGCGGCCGAGTTCGCTGGCAAACTCGGCAAGCGGCATCTGCTGCACGCTGAGTACACCTTCGGTCCAACCCCAGGGGTCGTTGTGCAGTGTTCTGGCCTGAATGCTGCCGTCgg                                              >  5:135673/1‑123 (MQ=255)
                                                                                                                         ccGGTCGGTAGCGGCCGAGTTCGCTGGCAAACTCGGCAAGCGGCATCTGCTGCACGCTGAGTACACCTTCGGTCCAACCCCAGGGGTCGTTGTGCAGTGTTCTGGCCTGAATGCTGCCGTCGGCTGTGACTTGCAAGGTTTCACCGGgc                    <  6:47485/149‑1 (MQ=255)
                                                                                                                          cggtcggtAGCGGCCGAGTTCGCTGGCAAACTCGGCAAGCGGCATCTGCTGCACTCTGAGTACACCTTCGGTCCAACCCCAGGGGTCGTTGTGCAGTGTTCTGGCCTGAATGCTGCCGTCGGCTGTGACTTGCAAGGTTTCACCGGgca                   >  5:257798/1‑149 (MQ=255)
                                                                                                                                gtAGCGGCCGAGTTCGCTGGCAAACTCGGCAAGCGGCATCTGCTGCACGCTGAGTACACCTTCGGTCCAACCCCAGGGGTCGTTGTGCAGTGTTCTGGCCTGAATGCTGCCGTCGGCTGTGACTTGCAAGGTTTCACCGGGCAGCACtt             <  1:231337/149‑1 (MQ=255)
                                                                                                                                       ccGAGTTCGCTGGCAAACTCGGCAACCGGCATCTGCTCCACGCTGAGTACACCTTCGGTCCAACCCCAGGGGTCGTTGTGCAGTGTTCTGGCCTGAATGCTGCCGTCGGCTGTGACTTGCAAGGTTTCACCGGGCAGCACTTGCCTCGc      >  1:255586/1‑149 (MQ=255)
                                                                                                                                           gTTCGCTGGCAAACTCGGCAAGCGGCATCTGCTGCACGCTGAGTACACCTTCGGTCCAACCCCAGGGGTCGTTGTGCAGTGTTCTGGCCTGAATGCTGCCGTCGGCTGTGACTTGCAAGGTTTCACCGGGCAGCACTTGCCTCGCTGAg  <  4:39701/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                              
CTGTACTCGACCCGCACGGGCAGGCTGCGGGTGAGCAGTTGCAGGATTTGCCCGGTGTCAGCCAACTGATAGGTGCCTGAAACCTTCAAGCCCGCCACTTCTGGAGCGCAGCGCAACAGGCCCGGTCGGTAGCGGCCGAGTTCGCTGGCAAACTCGGCAAGCGGCATCTGCTGCACGCTGAGTACACCTTCGGTCCAACCCCAGGGGTCGTTGTGCAGTGTTCTGGCCTGAATGCTGCCGTCGGCTGTGACTTGCAAGGTTTCACCGGGCAGCACTTGCCTCGCTGAG  >  NC_002947/817251‑817538

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: