Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,049,303 G→T 7.3% intergenic (‑32/+28) PP_1822 ← / ← folEA‑I hypothetical protein/GTP cyclohydrolase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,049,3030GT 7.3% 173.0 / 2.1 57intergenic (‑32/+28)PP_1822/folEA‑Ihypothetical protein/GTP cyclohydrolase I
Reads supporting (aligned to +/- strand):  ref base G (28/24);  new base T (3/2);  total (31/26)
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.95e-01

CTTGCTCGATGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGG‑‑‑CTCGGGGCCAGCAACCAATCACTTGATCA‑‑GGCTGAGGAA‑‑‑CTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGACACCGCGCATCATCATGCACATG  >  NC_002947/2049157‑2049451
                                                                                                                                                  |                                                                                                                                                             
cTTGCTCGATGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGaa                                                                                                                                                             <  3:19169/149‑1 (MQ=255)
      cGATGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCg                                                                                                                                                       >  3:364125/1‑149 (MQ=255)
      cGATGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCg                                                                                                                                                       <  6:352513/149‑1 (MQ=255)
          gCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGG‑‑‑CTc                                                                                                                                               >  2:2040/1‑149 (MQ=255)
              gCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGG‑‑‑CTCgggg                                                                                                                                           >  4:295711/1‑149 (MQ=255)
               ccTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGTAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGG‑‑‑CTCGGGGc                                                                                                                                          >  2:102628/1‑149 (MQ=255)
                   ggcggcggcGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGG‑‑‑CTCGGGGCCa                                                                                                                                        <  3:295711/147‑1 (MQ=255)
                      ggcggcggcGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGG‑‑‑CTCGGGGCCAGCAAc                                                                                                                                   <  4:118330/149‑1 (MQ=255)
                      ggcggcgTCGCGTGGGCGCTTGGCCGGGCGGCTGATCCAGTCCCCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGAGCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGG‑‑‑CTCGGGGCCAGCAAc                                                                                                                                   <  6:390865/149‑1 (MQ=255)
                       gcggcgTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGG‑‑‑CTCGGGGCCAGCAAcc                                                                                                                                  <  3:267077/149‑1 (MQ=255)
                             tCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGG‑‑‑CTCGGGGCCAGCAACCAATCAc                                                                                                                            >  6:206173/1‑149 (MQ=255)
                               gcgTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGG‑‑‑CTCGGGGCCAGCAACCAATCACt                                                                                                                           <  1:75784/148‑1 (MQ=255)
                                 gTTTGCGCTGGGCCGCGCGCCTGACCCCGCCACCCACCCCGACCCGCGGGCCGCGGCCACCCCCCCGGGCTTTCCCGACCAGCGGCACCCCCTCAACCCCAAACCCCCACCCGGACGGGGCGG‑‑‑CTCGGGGCCAGCAACCAATCACTTGa                                                                                                                        <  6:60270/149‑1 (MQ=255)
                                  tttGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGG‑‑‑CTCGGGGCCAGCAACCAATCACTTGAt                                                                                                                       >  2:315386/1‑149 (MQ=255)
                                       gcTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGG‑‑‑CTCGGGGCCAGCAACCAATCACTTGATCA‑‑GGc                                                                                                                >  6:383822/1‑149 (MQ=255)
                                           ggCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGG‑‑‑CTCGGGGCCAGCAACCAATCACTTGATCA‑‑GGCTGa                                                                                                             <  5:218085/148‑1 (MQ=255)
                                               gggCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGG‑‑‑CTCGGGGCCAGCAACCAATCACTTGATCA‑‑GGCTGAGGAA‑‑‑c                                                                                                     <  8:13107/149‑1 (MQ=255)
                                                   ggcTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGG‑‑‑CTCGGGGCCAGCAAc                                                                                                                                   >  2:352153/1‑120 (MQ=255)
                                                   ggcTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGG‑‑‑CTCGGGGCCAGCAAc                                                                                                                                   <  1:352153/120‑1 (MQ=255)
                                                    gcTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGG‑‑‑CTCGGGGCCAGCAACCAATCACTTGATCA‑‑GGCTGAGGAA‑‑‑CTCGCt                                                                                                <  1:2040/149‑1 (MQ=255)
                                                       gATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGG‑‑‑CTCGGGGCCAGCAACCAATCACTTGATCA‑‑GGCTGAGGAA‑‑‑CTCGCTGCg                                                                                             >  2:138818/1‑149 (MQ=255)
                                                                aCCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGG‑‑‑CTCGGGGCCAGCAACCAATCACTTGATCA‑‑GGCTGAGGAA‑‑‑CTCGCTGCGGGTggcggc                                                                                    >  4:231663/1‑149 (MQ=255)
                                                                aCCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGG‑‑‑CTCGGGGCCAGCAACCAATCACTTGATCA‑‑GGCTGAGGAA‑‑‑CTCGCTGCGGGTggcggc                                                                                    <  5:206173/149‑1 (MQ=255)
                                                                   gAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGG‑‑‑CTCGGGGCCAGCAACCAATCACTTGATCA‑‑GGCTGAGGAA‑‑‑CTCGCTGCgg                                                                                            <  8:218816/138‑1 (MQ=255)
                                                                   gAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGG‑‑‑CTCGGGGCCAGCAACCAATCACTTGATCA‑‑GGCTGAGGAA‑‑‑CTCGCTGCgg                                                                                            >  7:218816/1‑138 (MQ=255)
                                                                         gATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGG‑‑‑CTCGGGGCCAGCAACCAATCACTTGATAA‑‑GGCTGAGGAA‑‑‑CTCGCTGCGGGTGGCGGCGTTGTcgcg                                                                           >  3:205269/1‑149 (MQ=255)
                                                                                  gCCGAGGCCAACCCGCAGGGTTTTGTCGTTCATGGGCAACTCCTTAAATGCAAAACCTGACCCTTAAGGGTCGCGCCTCTCGGGGCCAGCAACCAATCACTTGATCA‑‑GGCTGAGGAA‑‑‑CTCGCTGCGGGTGGCGGCGTTGTCGCGGaa                                                                        <  1:315473/147‑1 (MQ=255)
                                                                                  gCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGG‑‑‑CTCGGGGCCAGCAACCAATCACTTGATCA‑‑GGCTGa                                                                                                             <  3:73962/109‑1 (MQ=255)
                                                                                  gCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGG‑‑‑CTCGGGGCCAGCAACCAATCACTTGATCA‑‑GGCTGa                                                                                                             >  4:73962/1‑109 (MQ=255)
                                                                                          cAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGG‑‑‑CTCGGGGCCAGCAACCAATCACTTGATCA‑‑GGCTGAGGAA‑‑‑CTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCa                                                          >  1:57288/1‑149 (MQ=255)
                                                                                               cGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGG‑‑‑CTCGTGGCCAGCAACCAATCACTTGATCA‑‑GGCTGAGGAA‑‑‑CTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAg                                                     <  4:364125/149‑1 (MQ=255)
                                                                                                gCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGG‑‑‑CTCGGGGCCAGCAACCAATCACTTGATCA‑‑GGCTGAGGAA‑‑‑CTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAgg                                                    <  4:355035/149‑1 (MQ=255)
                                                                                                     gCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGG‑‑‑CTCGGGGCCAGCAACCAATCACTTGATCA‑‑GGCTGAGGAA‑‑‑CTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATc                                               <  3:231663/149‑1 (MQ=255)
                                                                                                           aCGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGG‑‑‑CTCGGGGCCAGCAACCAATCACTTGATCA‑‑GGCTGAGGAA‑‑‑CTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTc                                         >  3:333231/1‑149 (MQ=255)
                                                                                                                     gCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGG‑‑‑CTCGGGGCCAGCAACCAATCACTTGATCA‑‑GGCTGAGGAA‑‑‑CTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCt                               >  5:184406/1‑149 (MQ=255)
                                                                                                                      cAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGG‑‑‑CTCGGGGCCAGCAACCAATCACTTGATCA‑‑GGCTGAGGAA‑‑‑CTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCt                               >  6:112280/1‑148 (MQ=255)
                                                                                                                       aaCTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGG‑‑‑CTCGGGGCCAGCAACCAATCACTTGATCA‑‑GGCTGAGGAA‑‑‑CTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTtc                             >  7:366241/1‑149 (MQ=255)
                                                                                                                         cTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGG‑‑‑CTCGGGGCCAGCAACCAATCACTTGATCA‑‑GGCTGAGGAA‑‑‑CTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTtctc                           <  2:57288/149‑1 (MQ=255)
                                                                                                                         cTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGG‑‑‑CTCGGGGCCAGCAACCAATCACTTGATCA‑‑GGCTGAGGAA‑‑‑CTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTtctc                           <  1:138818/149‑1 (MQ=255)
                                                                                                                          tCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGG‑‑‑CTCGGGGCCAGCAACCAATCACTTGATCA‑‑GGCTGAGGAA‑‑‑CTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCt                               <  4:280949/144‑1 (MQ=255)
                                                                                                                          tCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGG‑‑‑CTCGGGGCCAGCAACCAATCACTTGATCA‑‑GGCTGAGGAA‑‑‑CTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCt                               >  3:280949/1‑144 (MQ=255)
                                                                                                                            cTGAAATGCAAAAGCCGACCCTGAAGGGTCGG‑‑‑CTCGGGGCCAGCAACCAATCACTTGATCA‑‑GGCTGAGGAA‑‑‑CTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGac                        >  6:28515/1‑149 (MQ=255)
                                                                                                                             tGAAATGCAAAAGCCGACCCTGAAGGGTCGG‑‑‑CTCGGGGCCAGCAACCAATCACTTGATc                                                                                                                      >  5:320822/1‑59 (MQ=255)
                                                                                                                             tGAAATGCAAAAGCCGACCCTGAAGGGTCGG‑‑‑CTCGGGGCCAGCAACCAATCACTTGATc                                                                                                                      <  6:320822/59‑1 (MQ=255)
                                                                                                                             tGAAATGCAAAAGCCGACCCTGAAGGGTCGG‑‑‑CTCGGGGCCAGCAACCAATCACTTGATCA‑‑GGCTGAGGAA‑‑‑CTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGaca                       >  2:131481/1‑149 (MQ=255)
                                                                                                                                 aatcAAAATCCGCCCCTTTAGGGTGGG‑‑‑GTCGGGGGCATCAACCAATAACTTGATTATCCGCTGAGGAACTCCTCCCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGa                         <  1:315386/145‑1 (MQ=255)
                                                                                                                                 aTGCAAAAGCCGACCCTGAAGGGTCGG‑‑‑CTCGGGGCCAGCAACCAATCACTTGATCA‑‑GGCTGAGGAA‑‑‑CTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGACACcgc                   >  5:25219/1‑149 (MQ=255)
                                                                                                                                     aaaaGCCGACCCTGAAGGGTCGG‑‑‑CTCGGGGCCAGCAACCAATCACTTGATCA‑‑GGCTGAGGAA‑‑‑CTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGACACCGCGcat               >  6:19960/1‑149 (MQ=255)
                                                                                                                                     aaaaGCCGACCCTGAAGGGTCGG‑‑‑CTCGGGGCCAGCAACCAATCACTTGATCA‑‑GGCTGAGGAA‑‑‑CTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGACACCGCGcat               >  4:303592/1‑149 (MQ=255)
                                                                                                                                       gaGCCGACCCTTCAGGGTCGG‑‑‑CTCGGGGCCAGCAACCAATCACTTGATCA‑‑GGCTGAGGAA‑‑‑CTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGACACCGCGcatca             >  3:117069/2‑149 (MQ=255)
                                                                                                                                         gCCGACCCTTCAGGGTCGG‑‑‑CTCGGGGCCAGCAACCAATCACTTGATCA‑‑GGCTGAGGAA‑‑‑CTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGATTTCTGCTTCTCGACACCGCGcatcat            >  1:203498/1‑148 (MQ=255)
                                                                                                                                         gCCGACCCTGAAGGGTCGG‑‑‑CTCGGGGCCAGCAACCAATCACTTGATCA‑‑GGCTGAGGAA‑‑‑CTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGACACCGCGcatcatc           <  5:383822/149‑1 (MQ=255)
                                                                                                                                          ccGACCCTTCAGGGTCGG‑‑‑CTCGGGGCCAGCAACCAATCACTTGATCA‑‑GGCTGAGGAA‑‑‑CTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGACACCGCGcatcatca          >  3:85303/1‑149 (MQ=255)
                                                                                                                                          ccGACCCTGAAGGGTCGG‑‑‑CTCGGGGCCAGCAACCAATCACTTGATCA‑‑GGCTGAGGAA‑‑‑CTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGACACCGCGcatcatca          >  6:287872/1‑149 (MQ=255)
                                                                                                                                            gACCCTGAAGGGTCGG‑‑‑CTCGGGGCCAGCAACCAATCACTTGATCA‑‑GGCTGAGGAA‑‑‑CTCGCTGCGGGTGGCGGCGTTGTCGCGGACCTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGACACCGCGCATCATCATg        >  6:174876/1‑149 (MQ=255)
                                                                                                                                             aCCCTGAAGGGTCGG‑‑‑CTCGGGGCCAGCAACCAATCACTTGATCA‑‑GGCTGAGGAA‑‑‑CTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGACACCGCGCATCATCATGc       <  5:19960/149‑1 (MQ=255)
                                                                                                                                             aCCCTGAAGGGTCGG‑‑‑CTCGGGGCCAGCAACCAATCACTTGATCA‑‑GGCTGAGGAA‑‑‑CTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGACACCGAGCATCATCATg        >  7:344836/1‑148 (MQ=255)
                                                                                                                                                  gAAGGGTCGG‑‑‑CTCGGGGCCAGCAACCAATCACTTGATCA‑‑GGCTGAGGAA‑‑‑CTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGACACCGCGCATCATCATGCACAtg  >  2:312185/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                             
CTTGCTCGATGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGG‑‑‑CTCGGGGCCAGCAACCAATCACTTGATCA‑‑GGCTGAGGAA‑‑‑CTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGACACCGCGCATCATCATGCACATG  >  NC_002947/2049157‑2049451

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: