Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,615,630 G→C 9.4% intergenic (‑109/+22) hldE ← / ← msbA bifunctional heptose 7‑phosphate kinase/heptose 1‑phosphate adenylyltransferase/lipid A export ATP‑binding/permease MsbA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,615,6300GC 9.4% 187.6 / 7.4 63intergenic (‑109/+22)hldE/msbAbifunctional heptose 7‑phosphate kinase/heptose 1‑phosphate adenylyltransferase/lipid A export ATP‑binding/permease MsbA
Reads supporting (aligned to +/- strand):  ref base G (31/26);  new base C (4/2);  total (35/28)
Fisher's exact test for biased strand distribution p-value = 6.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGGGCTTGATCGAAACGCGGCATGGACAACTTCATGGGCAACCCATATGGAAATAATGAACAGGGGCAGGATATTAGCACAGGGTAGTCAATGGCTTTATGGCCAGATCATCCCTGGAAAATTCCCGTGACAATCGAGGCCATGACGGCCCCGATCAAGGAGGGTTCAGGTGATATCGGCCTTGGCCGGCTCATCCAGCCCCATGGCATGCAAACGGGCATAATGGCCATTAGCCGCAAGCAGCTCGGTGTGAGTACCGCGCTCGACCAGGCGGCCCTGGTCCATGACAA  >  NC_002947/5615486‑5615776
                                                                                                                                                 |                                                                                                                                                  
cGGGGCTTGATCGAAACGCGGCATGGACAACTTCATGGGCAACCCATATGGAAATAATGAACAGGGGCAGGATATTAGCACAGGGTAGTCAATGGCTTTGTGGCCAGATCATCCCTGGAAAATTCCCGTGACAATCGAGACCATGACgg                                                                                                                                                <  8:73544/149‑1 (MQ=255)
   ggCTTGATCGAAACGCGGCATGGACAACTTCATGGGCAACCCATATGGAAATAATGAACAGGGGCAGGATATTAGCACAGGGTAGTCAATGGCTTTATGGCCAGATCATCCCTGGAAAATTCCCGTGACAATCGAGGCCATGACGGccc                                                                                                                                             <  2:353394/149‑1 (MQ=255)
        gATCGAAACGCGGCATGGACAACTTCATGGGCAACCCATATGGAAATAATGAACAGGGGCAGGATATTAGCACAGGGTAGTCAATGGCTTTATGGCCAGATCATCCCTGGAAAATTCCCGTGACAATCGAGGCCATGACGGCCCCGAt                                                                                                                                         >  6:77212/1‑148 (MQ=255)
        gATCGAAACGCGGCATGGACAACTTCATGGGCAACCCATATGGAAATAATGAACAGGGGCAGGATATTAGCACAGGGTAGTCAATGGCTTTATGGCCAGATCATCCCTGGAAAATTCCCGTGACAATCGAGGCCATGACGGCCCCGAt                                                                                                                                         >  6:220639/1‑148 (MQ=255)
          tCGAAACGCGGCATGGACAACTTCATGGGCAACCCATATGGAAATAATGAACAGGGGCAGGATATTAGCACAGGGTAGTCAATGGCTTTATGGCCAGATCATCCCTGGAAAATTCCCGTGACAATCGAGGCCATGACGGCCCCGATc                                                                                                                                        >  1:284615/1‑147 (MQ=255)
          tCGAAACGCGGCATGGACAACTTCATGGGCAACCCATATGGAAATAATGAACAGGGGCAGGATATTAGCACAGGGTAGTCAATGGCTTTATGGCCAGATCATCCCTGGAAAATTCCCGTGACAATCGAGGCCATGACGGCCCCGATc                                                                                                                                        >  2:224496/1‑147 (MQ=255)
             aaaCGCGGCATGGACAACTTCATGGGCAACCCATATGGAAATAATGAACAGGGGCAGGATATTAGCACAGGGTAGTCAATGGCTTTATGGCCAGATCATCCCTGGAAAATTCCCGTGACAATCGAGGCCATGACGGCCCCGATCAagga                                                                                                                                   <  2:243237/149‑1 (MQ=255)
             aaaCGCGGCATGGACAACTTCATGGGCAACCCATATGGAAATAATGAACAGGGGCAGGATATTAGCACAGGGTAGTCAATGGCTTTATGGCCAGATCATCCCTGGAAAATTCCCGTGACAATCGAGGCCATGACGGCCCCGATCAag                                                                                                                                     >  4:129924/1‑147 (MQ=255)
              aaCGCGGCATGGACAACTTCATGGGCAACCCATATGGAAATAATGAACAGGGGCAGGATATTAGCACAGGGTAGTCAATGGCTTTATGGCCAGATCATCCCTGGAAAATTCCCGTGACAATCGAGGCCATGACGGCCCCGATCAaggag                                                                                                                                  <  5:77212/149‑1 (MQ=255)
                       tGGACAACTTCATGGGCAACCCATATGGAAATAATGAACAGGGGCAGGATATTAGCACAGGGTAGTCAATGGCTTTATGGCCAGATCATCCCTGGAAAATTCCCGTGACAATCGAGGCCATGACGGCCCCGATCAAGGAGGGTTCAGGt                                                                                                                         <  1:224496/149‑1 (MQ=255)
                                tCATGGGCAACCCATATGGAAATAATGAACAGGGGCAGGATATTAGCACAGGGTAGTCAATGGCTTTATGGCCAGATCATCCCTGGAAAATTCCCGTGACAATCGAGGCCATGACGGCCCCGATCAAGGAGGGTTCAGGTGATATCGGc                                                                                                                >  5:137063/1‑149 (MQ=255)
                                 cATGGGCAACCCATATGGAAATAATGAACAGGGGCAGGATATTAGCACAGGGTAGTCAATGGCTTTATGGCCAGATCATCCCTGGAAAATTCCCGTGACAATCGAGGCCATGACGGCCCCGATCAAGGAGGGTTCAGGTGATATCGGc                                                                                                                <  5:220639/148‑1 (MQ=255)
                                         aCCCATATGGAAATAATGAACAGGGGCAGGATATTAGCACAGGGTAGTCAATGGCTTTATGGCCAGATCATCCCTGGAAAATTCCCGTGACAATCGAGGCCATGACGGCCCCGATCAAGGAGGGTTCAGGTGATATCGGCCTTGGCCgg                                                                                                       <  4:273143/149‑1 (MQ=255)
                                           ccATATGGAAATAATGAACAGGGGCAGGATTTTAGCACAGGGTAGTCAATGGCTTTATGGCCAGATCATCCCTGGAAAATTCCCGTGACAATCGAGGCCATGACGGCCCCGATCAAGGAGGGTTCAGGTGATATCGGCCTTGGCCGGc                                                                                                      >  5:15696/1‑148 (MQ=255)
                                             atatGGAAATAATGAACAGGGGCAGGATATTAGCACAGGGTAGTCAATGGCTTTATGGCCAGATCATCCCTGGAAAATTCCCGTGACAATCGAGGCCATGACGGCCCCGATCAAGGAGGGTTCAGGTGATATCGGCCTTGGCCGGCTCa                                                                                                   >  3:118431/1‑149 (MQ=255)
                                                tGGAAATAATGAACAGGGGCAGGATATTAGCACAGGGTAGTCAATGGCTTTATGGCCAGATCATCCCTGGAAAATTCCCGTGACAATCGAGGCCATGACGGCCCCGATCAAGGAGGGTTCAGGTGATATCGGCCTTGGCCGGCTCATc                                                                                                 <  3:119157/148‑1 (MQ=255)
                                                   aaaTAATGAACAGGGGCAGGATATTAGCACAGGGTAGTCAATGGCTTTATGGCCAGATCATCCCTGGAAAATTCCCGTGACAATCGAGGCCATGACGGCCCCGATCAagga                                                                                                                                   <  3:375849/111‑1 (MQ=255)
                                                   aaaTAATGAACAGGGGCAGGATATTAGCACAGGGTAGTCAATGGCTTTATGGCCAGATCATCCCTGGAAAATTCCCGTGACAATCGAGGCCATGACGGCCCCGATCAagga                                                                                                                                   >  4:375849/1‑111 (MQ=255)
                                                         tGAACCGGGGCAGGAATTTACCACCGGGTAGTCAATGGCTTTATGGCCATAACACCCCTGTAAAATTCCCGGGACAATCCCGGCCATCCCGGCCCCGAACAAGGCGGGTTCCGGTTATATCGGCCTGGGCCGGCTCATCCAGCCCCagg                                                                                       <  6:137063/149‑3 (MQ=255)
                                                         tGAACAGGGTCAGGATATTAGCACAGGGTAGTCAATGGCTTTATGGCCAGATCATCCCTGGAAAATTCCCGTGACAATCGAGGCCATGACGGCCCCGATCAAGGAGGGTTCAGGTGATATCGGCCTTGGCCGGCTCATCCAGCCCCATg                                                                                       <  5:271814/149‑1 (MQ=255)
                                                         tGAACAGGGGCAGGATATTAGCACAGGGTAGTCAATGGCTTTATGGCCAGATCATCCCTGGAAAATTCCCGTGACAATCGAGGCCATGACGGCCCCGATCAAGGAGGGTTCAGGTGATATCGGCCTTGGCCGGCTCATCCAGCCCCATg                                                                                       <  2:284615/149‑1 (MQ=255)
                                                            aCAGGGCCAGGATATTAGCACAGGGTAGTCAATGGCTTTATGGCCAGATCATCCCTGGAAAATTCCCGTGACAATCGAGGCCATGACGGCCCCGATCAAGGAGGGTTCAGGTGATATCGGCCTTGCCCGGCTCATCCAGCCCCATGGc                                                                                     <  7:240817/148‑1 (MQ=255)
                                                              aGGGGCAGGATATTAGCACAGGGTAGTCAATGGCTTTATGGCCAGATCATCCCTGGAAAATTCCCGTGACAATCGAGGCCATGACGGCCCCGATCAAGGAGGGTTCAGGTGATATCGGCCTTGGCCGGCTCATCCAGCCCCATGGCATg                                                                                  >  3:140370/1‑149 (MQ=255)
                                                              aGGGGCAGGATATTAGCACAGGGTAGTCAATGGCTTTATGGCCAGATCATCCCTGGAAAATTCCCGTGACAATCGAGGCCATGACGGCCCCGATCAAGGAGGGTTCAGGTGATATCGGCCTTGGCCGGCTCATCCAGCCCCATGGCATg                                                                                  <  1:127995/149‑1 (MQ=255)
                                                                gggCAGGATATTAGCACAGGGTAGTCAATGGCTTTATGGCCAGATCAGTCCCTGGAAAATTCCCGTGACAATCGAGGCCATGACGGCCCCGATCAAGGAGGGTTCAGGTGATATCGGCCTTGGCCGGCTCATCCAGCCCCATGGCATg                                                                                  >  1:274650/1‑148 (MQ=255)
                                                                gggCAGGATATTAGCACAGGGTAGTCAATGGCTTTATGGCCAGATCATCCCTGGAAAATTCCCGTGACAATCGAGGCCATGACGGCCCCGATCAAGGAGGGTTCAGGTGATATCGGCCTTGGCCGGCTCATCCAGCCCCATGGCATGc                                                                                 >  7:261354/1‑148 (MQ=255)
                                                                gggCAGGATATTAGCACAGGGTAGTCAATGGCTTTATGGCCAGATCATCCCTGGAAAATTCCCGTGACAATCGAGGCCATGACGGCCCCGATCAAGGAGGGTTCAGGTGATATCGGCCTTGGCCGGCTCATCCAGCCCCATGGCATGc                                                                                 >  6:128400/1‑148 (MQ=255)
                                                                   cAGGATATTAGCACAGGGTAGTCAATGGCTTTATGGCCAGATCATCCCTGGAAAATTCCCGTGACAATCGAGGCCATGACGGCCCCGATCAAGGAGGGTTCAGGTGATATCGGCCTTGGCCGGCTCATCCAGCCCCATGGCATGc                                                                                 >  7:257342/1‑145 (MQ=255)
                                                                   cAGGATATTAGCACAGGGTAGTCAATGGCTTTATGGCCAGATCATCCCTGGAAAATTCCCGTGACAATCGAGGCCATGACGGCCCCGATCAAGGAGGGTTCAGGTGATATCGGCCTTGGCCGGCTCATCCAGCCCCATGGCATGc                                                                                 <  8:257342/145‑1 (MQ=255)
                                                                     ggATATTAGCACAGGGTAGTCAATGGCTTTATGGCCAGATCATCCCTGGAAAATTCCCGTGACAATCGAGGCCATGACGGCCCCGATCAAGGAGGGTTCAGGTGATATCGGCCTTGGCCGGCTCATCCAGCCCCATGGCATGCAAACgg                                                                           <  4:140370/149‑1 (MQ=255)
                                                                      gATTTTAGCACAGGGTAGTCAATGGCTTTATGGCCAGATCATCCCTGGAAAATTCCCGTGACAATCGAGGCCATGACGGCCCCGATCAAGGAGGGTTCAGGTGATATCGGCCTTGGCCGGCTCATCCAGCCCCATGGCATGCAAACggg                                                                          <  6:15696/149‑1 (MQ=255)
                                                                         atTAGCACAGGGTAGTCAATGGCTTTATGGCCAGATCATCCCTGGAAAATTCCCGTGACAATCGAGGCCATGACGGCCCCGATCAAGGAGGGTTCAGGTGATATCGGCCTTGGCCGGCTCATCCAGCCCCATGGCATGCAAACgg                                                                           <  5:280297/145‑1 (MQ=255)
                                                                         atTAGCACAGGGTAGTCAATGGCTTTATGGCCAGATCATCCCTGGAAAATTCCCGTGACAATCGAGGCCATGACGGCCCCGATCAAGGAGGGTTCAGGTGATATCGGCCTTGGCCGGCTCATCCAGCCCCATGGCATGCAAACgg                                                                           >  6:280297/1‑145 (MQ=255)
                                                                         atTAGCACAGGGTAGTCAATGGCTTTATGGCCAGATCATCCCTGGAAAATTCCCGTGACAATCGAGGCCATGACGGCCCCGATCAAGGAGGGTTCAGGTGATATCGGCCTTGGCCGGCTCATCCAGCCCCATGGCATGCAAACGGGCa                                                                        >  5:152352/1‑148 (MQ=255)
                                                                          ttAGCACAGGGTAGTCAATGGCTTTATGGCCAGATCATCCCTGGAAAATTCCCGTGACAATCGAGGCCATGACGGCCCCGATCAAGGAGGGTTCAGGTGATATCGGCCTTGGc                                                                                                          <  8:118690/113‑1 (MQ=255)
                                                                          ttAGCACAGGGTAGTCAATGGCTTTATGGCCAGATCATCCCTGGAAAATTCCCGTGACAATCGAGGCCATGACGGCCCCGATCAAGGAGGGTTCAGGTGATATCGGCCTTGGc                                                                                                          >  7:118690/1‑113 (MQ=255)
                                                                                     tAGTCAATGGCTTTATGGCCAGATCATCCCTGGAAAATTCCCGTGACAATCGAGGCCATGACGGCCCCGATCAAGGAGGGTTCAGGTGATATCGGCCTTGGCCGGCTCATCCAGCCCCATGGCATGCAAACGGGCATAATGGCCATTAg                                                           >  7:360084/1‑149 (MQ=255)
                                                                                          aaTGGCTTTATGGCCAGATCATCCCTGGAAAATTCCCGTGACAATCGAGGCCATGACGGCCCCGATCAAGGAGGGTTCAGGTGATATCGGCCTTGGCCGGCTCATCCAGCCCCATGGCATGCAAACGGGCATAATGGCCATTAGCCGc                                                       >  8:265376/1‑148 (MQ=255)
                                                                                          aaTGGCTTTATGGCCAGATCATCCCTGGAAAATTCCCGTGACAATCGAGGCCATGACGGCCCCGATCAAGGAGGGTTCAGGTGATATCGGCCTTGGCCGGCTCATCCAGCCCCATGGCATGCAAACGGGCATAATGGCCATTAGCCGc                                                       >  5:270322/1‑148 (MQ=255)
                                                                                           aTGGCTTTATGGCCAGATCATCCCTGGAAAATTCCCGTGACAATCGAGGCCATGACGGCCCCGATCAAGGAGGGTTCAGGTGATATCGGCCTTGGCCGGCTCATCCAGCCCCATGGCATGCAAACGGGCATAATGGCCATTAGCCGCaa                                                     <  6:152352/149‑1 (MQ=255)
                                                                                              gCTTTATGGCCAGATCATCCCTGGAAAATTCCCGTGACAATCGAGGCCATGACGGCCCCGATCAAGGAGGGTTCAGGTGATATCGGCCTTGGCCGGCTCATCCAGCCCCATGGCATGCAAACGGGCATAATGGCCATTAGCCGCAag                                                    >  3:377483/1‑147 (MQ=255)
                                                                                                tttATGGCCAGATCATCCCTGGAAAATTCCCGTGACAATCGAGGCCATGACGGCCCCGATCAAGGAGGGTTCAGGTGATATCGGCCTTGGCCGGCTCATCCAGCCCCATGGCATGCAAACGGGCATAATGGCCATTAGCCGCAagcagc                                                >  1:341865/1‑149 (MQ=255)
                                                                                                     ggCCAGATCATCCCTGGAAAATTCCCGTGACAATCGAGGCCATGACGGCCCCGATCAAGGAGGGTTCAGGTGATATCGGCCTTGGCCGGCTCATCCAGCCCCATGGCATGCAAACGGGCATAATGGCAATTAGCCGCAAGCAGCTCGgt                                           >  3:178846/1‑149 (MQ=255)
                                                                                                      gCCAGATCATCCCTGGAAAATTCCCGTGACAATCGAGGCCATGACGGCCCCGATCAAGGAGGGTTCAGGTGATATCGGCCTTGGCCGGCTCATCCAGCCCCATGGCATGCAAACGGGCATAATGGCCATTAGCCGCAAGCAGCTCGgtg                                          >  6:259309/1‑149 (MQ=255)
                                                                                                       ccAGATCATCCCTGGAAAATTCCCGTGACAATCGAGGCCATGACGGCCCCGATCAAGGAGGGTTCAGGTGATATCGGCCTTGGCCGGCTCATCCAGCCCCATGGCATGCAAACGGGCATAATGGCCATTAGCCGCAAGCAGCTCGgtgt                                         >  4:274393/1‑149 (MQ=255)
                                                                                                         aGATCATCCCTGGAAAATTCCCGTGACAATCGAGGCCATGACGGCCCCGATCAAGGAGGGTTCAGGTGATATCGGCCTTGGCCGGCTCATCCAGCCCCATGGCATGCAAACGGGCATAATGGCCATTAGCCGCAAGCAGCTCGGTGTGa                                       <  7:265376/149‑1 (MQ=255)
                                                                                                         aGATCATCCCTGGAAAATTCCCGTGACAATCGAGGCCATGACGGCCCCGATCAAGGAGGGTTCAGGTGATATCGGCCTTGGCCGGCTCATCCAGCCCCATGGCATGCAAACGGGCATAATGGCCATTAGCCGCAAGCAGCTCGGTGTGa                                       <  5:128400/149‑1 (MQ=255)
                                                                                                                   cTGGAAAATTCCCGTGACAATCGAGGCCATGACGGCCCCGATCAAGGAGGGTTCAGGTGATATCGGCCTTGGCCGGCTCATCCAGCCCCATGGCATGCAAACGGGCATAATGGCCATTAGCCGCAAGCAGCTCGGTGTGAGTACCgcgc                              >  4:239739/1‑149 (MQ=255)
                                                                                                                              ccGTGACAATCGAGGCCATGACGGCCCCGATCAAGGAGGGTTCAGGTGATATCGGCCTTGGCCGGCTCATCCAGCCCCATGGCATGCAAACGGGCATAATGGCCATTAGCCGCAAGCAGCTCGGTGTGAGTACCGCGCTCGACCAggcg                   >  6:50754/1‑149 (MQ=255)
                                                                                                                              ccGTGACAATCGAGGCCATGACGGCCCCGATCAAGGAGGGTTCAGGTGATATCGGCCTTGGCCGGCTCATCCAGCCCCATGGCATGCAAACGGGCATAATGGCCATTAGCCGCAAGCAGCTCGGTGTGAGTACCGCGCTCGACCAggcg                   >  6:85955/1‑149 (MQ=255)
                                                                                                                               cGTGACAATCGAGGCCATGACGGCCCCGATCAAGGAGGGTTCAGGTGATATCGGCCTTGGCCGGCTCATCCAGCCCCATGGCATGCAAACGGGCATAATGGCCATTAGCCGCAAGCAGTTCGGTGTGAGTACCGCGCTCGACCAggcgg                  >  7:337493/1‑149 (MQ=255)
                                                                                                                               cGTGACAATCGAGGCCATGACGGCCCCGATCAAGGAGGGTTCAGGTGATATCGGCCTTGGCCGGCTCATCCAGCCCCATGGCATGCAAACGGGCATAATGGCCATTAGCCGCAAGCAGCTCGGTGTGAGTACCGCGCTCGACCAggcgg                  <  3:129924/149‑1 (MQ=255)
                                                                                                                                tccttgaTCGGGGCCGTCATGGCCTCGATTAAGGAGGGTTCAGGTGATATCGGCCTTGGCCGGCTCATCCAGCCCCATGGCATGCAAACGGGCATAATGGCCATGGGCCGCAAGCAGCTCGGTGTGAGTACCGCGCTCGACCAggcggc                 >  1:317099/7‑149 (MQ=255)
                                                                                                                                   ttgaTCGGGGCCGTCATGGCCTCGATCAAGGAGGGTTCAGGTGATATCGGCCTTGGACGGCTCATCCAGCCCCATGGCATGCAAACGGGCATAATGGCCATTAGCCGCAAGCAGCTCGGTGTGAGTACAGAGCTCGACCAGGCGGCGCt              >  1:210919/4‑149 (MQ=255)
                                                                                                                                      aTCGGGGCCGTCATGGCCTCGATCAAGGAGGGTTCAGGTGATATCGGCCTTGGCCGGCTCATCCAGCCCCATGGCATGCAAACGGGAATAATGGCCATTAGACGCAAGCAGCTCGGTGTGAGTACCGCGATCGACCCGGCGGCCCTGGt           >  4:20280/1‑149 (MQ=255)
                                                                                                                                      aTCGAGGCCATGACGGCCCCGATCAAGGAGGGTTCAGGTGATATCGGCCTTGGCCGGCTCATCCAGCCCCATGGCATGCAAACGGGCATAATGGCCATTAGCCGCAAGCAGCTCGGTGTGAGTACCGCGCTCGACCAGGCGGCCCTGGt           <  3:239739/149‑1 (MQ=255)
                                                                                                                                       tCGGGGCCGTCATGGCCTCGATCAAGGAGGGTTCAGGTGATATCGGCCTTGGCCGGCTCATCCAGCCCCATGGCATGCAAACGGGCATAATGGCCATTAGCCGCAAGCAGCTCGGTGTGAGTACCGCGCTc                            >  1:255911/1‑131 (MQ=255)
                                                                                                                                       tCGGGGCCGTCATGGCCTCGATCAAGGAGGGTTCAGGTGATATCGGCCTTGGCCGGCTCATCCAGCCCCATGGCATGCAAACGGGCATAATGGCCATTAGCCGCAAGCAGCTCGGTGTGAGTACCGCGCTc                            <  2:255911/131‑1 (MQ=255)
                                                                                                                                             ccATGACGGCCCCGCTCAAGGAGGGTTCAGGTGATATCGGCCTTGGCCGGCTCATCCAGCCCCATGGCATGCAAACGGGCATAATGGCCATTAGCCGCAAGCAGCGAGGGGTGAGTACCGCGCTCGACCAGCCCGCCCAGGGCCATGa     >  2:308573/1‑148 (MQ=255)
                                                                                                                                             ccATGACGGCCCCGATCAAGGAGGGTTCAGGTGATATCGGCCTTGGCCGGCTCATCCAGCCCCATGGCATGCAAACGGGCATAATGGCCATTAGCCGCAAGCAGCTCGGTGTGAGTACCGCGCTCGACCAGGCGGCCCTGGTCCATGAc    <  5:50754/149‑1 (MQ=255)
                                                                                                                                             ccATGACGGCCCCGATCAAGGAGGGTTCAGGTGATATCGGCCTTGGCCGGCTCATCCAGCCCCAGGGCATGCAAACGGGCATAATGGCCATTAGCCCCACGCAGCGCGGTGTGAGTACCGCGCTCGCCCAGGCGGCACTGGACCATgc     >  2:312674/1‑147 (MQ=255)
                                                                                                                                               aTGACGGCCCCGATCAAGGAGGGTTCAGGTGATATCGGCCTTGGCCGGCTCATCCAGCCCCATGGCATGCAAACGGGCATAATGGCCATTAGCCGCAAGCAGCTCGGTGTGAGTACCGCGCTCGACCAGGCGGCCCTGGTCCATGACaa  <  8:261354/149‑1 (MQ=255)
                                                                                                                                               aTGACGGCCCCGATCAAGGAGGGTTCAGGTGATATCGGCCTTGGCCGGCTCATCCAGCCCCATGGCATGCAAACGGGCATAATGGCCATTAGCCGCAAGCAGCTCGGTGTGAGTACCGCGCTCGACCAGGCGGCCCTGGTCCATGACaa  <  8:337493/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
CGGGGCTTGATCGAAACGCGGCATGGACAACTTCATGGGCAACCCATATGGAAATAATGAACAGGGGCAGGATATTAGCACAGGGTAGTCAATGGCTTTATGGCCAGATCATCCCTGGAAAATTCCCGTGACAATCGAGGCCATGACGGCCCCGATCAAGGAGGGTTCAGGTGATATCGGCCTTGGCCGGCTCATCCAGCCCCATGGCATGCAAACGGGCATAATGGCCATTAGCCGCAAGCAGCTCGGTGTGAGTACCGCGCTCGACCAGGCGGCCCTGGTCCATGACAA  >  NC_002947/5615486‑5615776

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