Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,951,834 G→C 10.9% I569M (ATC→ATG) ‡ mnmC ← methyltransferase/FAD‑dependent oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,951,8340GC10.9% 179.2 / 13.0 64I569M (ATC→ATG) ‡mnmCmethyltransferase/FAD‑dependent oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (29/28);  new base C (4/3);  total (33/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.83e-01

CACTCAACGGTGCCGTGATCAACCCGCGTGAGCCATGCCCGCTGTTCACATACAACCCGCCCAGCCAAGGGCAGGGCACATCCGGTACTTGCCGCGCATCCCGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCAGGTAATCCGGGCTGGTGCAACGGAACGCAGCACGTCCCTGCAATTGCTCGGGGTCCAGTTCAGCCACAGCCAGCCGTTGCGCCAGATCGACGGAAATTTCGTCCAGCAGCGCCAGGTTGCCCTGGT  >  NC_002947/1951689‑1951979
                                                                                                                                                 |                                                                                                                                                 
cACTCAACGGTGCCGTGATCAACCCGCGTGAGCCATGCCCGCTGTTCACATACAACCCGCCCAGCCAAGGGCAGGGCACATCCGGTACTTGCCGCGCATCCCGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGa                                                                                                                                                  >  8:252799/1‑147 (MQ=255)
 aCTCAACGGTGCCGTGATCAACCCGCGTGAGCCATGCCCGCTGTTCACATACAACCCGCCCAGCCAAGGGCAGGGCACATCCGGTACTTGCCGCGCATCCCGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAg                                                                                                                                               >  8:239044/1‑149 (MQ=255)
  cTCAACGGTGCCGTGATCAACCCGCGTGAGCCATGCCCGCTGTTCACATACAACCCGCCCAGCCAAGGGCAGGGCACATCCGGTACTTGCCGCGCATCCCGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAgg                                                                                                                                              >  7:251042/1‑149 (MQ=255)
   tCAACGGTGCCGTGATCAACCCGCGTGAGCCATGCCCGCTGTTCACATACAACCCGCCCAGCCAAGGGCAGGGCACATCCGGTACTTGCCGCGCATCCCGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAggg                                                                                                                                             >  1:6338/1‑149 (MQ=255)
   tCAACGGTGCCGTGATCAACCCGCGTGAGCCATGCCCGCTGTTCACATACAACCCGCCCAGCCAAGGGCAGGGCACATCCGGTACTTGCCGCGCATCCCGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAggg                                                                                                                                             <  4:1685/149‑1 (MQ=255)
   tCAACGGTGCCGTGATCAACCCGCGTGAGCCATGCCCGCTGTTCACATACAACCCGCCCAGCCAAGGGCAGGGCACATCCGGTACTTGCCGCGCATCCCGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAggg                                                                                                                                             >  2:311829/1‑149 (MQ=255)
     aaCGGTGCCGTGATCAACCCGCGTGAGCCATGCCCGCTGTTCACATACAACCCGCCCAGCCAAGGGCAGGGCACATCCGGTACTTGCCGCGCATCCCGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGcc                                                                                                                                           <  1:263908/149‑1 (MQ=255)
     aaCGGTGCCGTGATCAACCCGCGTGAGCCATGCCCGCTGTTCACATACAACCCGCCCAGCCAAGGGCAGGGCACATCCGGTACTTGCCGCGCATCCCGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGc                                                                                                                                            <  1:311829/148‑1 (MQ=255)
       ctgTGCCGTGATCAACCCGCGTGAGCCATGCCCGCTGTTCACATACAACCCGCCCAGCCAAGGGCAGGGCACATCCGGTACTTGCCGCGCATCCCGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCcga                                                                                                                                         <  1:229090/147‑1 (MQ=255)
        ggTGCCGTGATCAACCCGCGTGAGCCATGCCCGCTGTTCACATACAACCCGCCCAGCCAAGGGCAGGGCACATCCGGTACTTGCCGCGCATCCCGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCcgac                                                                                                                                        <  2:229193/149‑1 (MQ=255)
        ggTGCCGTGATCAACCCGCGTGAGCCATGCCCGCTGTTCACATACAACCCGCCCAGCCAAGGGCAGGGCACATCCGGTACTTGCCGCGCATCCCGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCCCTATCGCcgac                                                                                                                                        <  5:112777/149‑1 (MQ=255)
            ccGTGATCAACCCGCGTGAGCCATGCCCGCTGTTCACATACAACCCGCCCAGCCAAGGGCAGGGCACATCCGGTACTTGCCGCGCATCCCGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATc                                                                                                                                    <  7:239044/149‑1 (MQ=255)
                   cAACCCGCGTGAGCCATGCCCGCTGTTCACATACAACCCGCCCAGCCAAGGGCAGGGCACATCCGGTACTTGCCGCGCATCCCGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCAgg                                                                                                                              >  4:119019/1‑148 (MQ=255)
                    aaCCCGCGTGAGCCATGCCCGCTGTTCACATACAACCCGCCCAGCCAAGGGCAGGGCACATCCGGTACTTGCCGCGCATCCCGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCcga                                                                                                                                         >  3:309340/1‑136 (MQ=255)
                    aaCCCGCGTGAGCCATGCCCGCTGTTCACATACAACCCGCCCAGCCAAGGGCAGGGCACATCCGGTACTTGCCGCGCATCCCGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCcga                                                                                                                                         <  4:309340/136‑1 (MQ=255)
                    aaCCCGCGTGAGCCATGCCCGCTGTTCACATACAACCCGCCCAGCCAAGGGCAGGGCACATCCGGTACTTGCCGCGCATCCCGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCCCTATCGCcgacgc                                                                                                                                      >  3:19408/1‑138 (MQ=255)
                    aaCCCGCGTGAGCCATGCCCGCTGTTCACATACAACCCGCCCAGCCAAGGGCAGGGCACATCCGGTACTTGCCGCGCATCCCGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCCCTATCGCcgacgc                                                                                                                                      <  4:19408/139‑2 (MQ=255)
                        cgcgTGAGCCATGCCCGCTGTTCACATACAACCCGCCCAGCCAAGGGCAGGGCACATCCGGTACTTGCCGCGCATCCCGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCAg                                                                                                                               >  2:322732/1‑142 (MQ=255)
                        cgcgTGAGCCATGCCCGCTGTTCACATACAACCCGCCCAGCCAAGGGCAGGGCACATCCGGTACTTGCCGCGCATCCCGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCAg                                                                                                                               <  1:322732/142‑1 (MQ=255)
                           gTGAGCCATGCCCGCTGTTCACATACAACCCGCCCAGCCAAGGGCAGGGCACATCCGGTACTTGCCGCGCATCCCGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCAGGTAATCCggg                                                                                                                     <  7:252799/149‑1 (MQ=255)
                              aGCCATGCCCGCTGTTCACATACAACCCGCCCAGCCAAGGGCAGGGCACATCCGGTACTTGCCGCGCATCCCGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCAGGTAATCCGGGCTg                                                                                                                  <  3:119019/149‑1 (MQ=255)
                               gCCATGCCCGCTGTTCACATACAACCCGCCCAGCCAAGGGCAGGGCACATCCGGTACTTGCCGCGCATCCCGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCAGGTAATCCGGGCTg                                                                                                                  <  2:217573/148‑1 (MQ=255)
                               gCCATGCCCGCTGTTCACATACAACCCGCCCAGCCAAGGGCAGGGCACATCCGGTACTTGCCGCGCATCCCGGCCCAGTACTGCATAGGCTTCGGCAAAAGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCAGGTAATCCGGGCTg                                                                                                                  >  1:217573/1‑148 (MQ=255)
                                ccATGCCCGCTGTTCACATACAACCCGCCCAGCCAAGGGCAGGGCACATCCGGTACTTGCCGCGCATCCCGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCCCTATCGCcgac                                                                                                                                        >  4:105245/1‑125 (MQ=255)
                                ccATGCCCGCTGTTCACATACAACCCGCCCAGCCAAGGGCAGGGCACATCCGGTACTTGCCGCGCATCCCGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCCCTATCGCcgac                                                                                                                                        <  3:105245/125‑1 (MQ=255)
                                    gCCCGCTGTTCACATACAACCCGCCCAGCCAAGGGCAGGGCACATCCGGTACTTGCCGCGCATCCCGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCAGGTAATCCGGGCTGGTGCaa                                                                                                            >  8:210110/1‑149 (MQ=255)
                                        gCTGTTCACATACAACCCGCCCAGCCAAGGGCAGGGCACATCCGGTACTTGCCGCGCATCCCGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCAGGTAATCCGGGCTGGTGCAACgg                                                                                                         >  3:155525/1‑148 (MQ=255)
                                          tGTTCACATACAACCCGCCCAGCCAAGGGCAGGGCACATCCGGTACTTGCCGCGCATCCCGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCAGGTAATCCGGGCTGGTGCAACGGAAc                                                                                                      >  6:1681/1‑149 (MQ=255)
                                            ttCACATACAACCCGCCCAGCCAAGGGCAGGGCACATCCGGTACTTGCCGCGCATCCCGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCAGGTAATCCGGGCTGGTGCAACGGAACgc                                                                                                    <  1:311126/149‑1 (MQ=255)
                                              cacaTACAACCCGCCCAGCCAAGGGCAGGGCACATCCGGTACTTGCCGCGCATCCCGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCAGGTAATCCGGGCTGGTGCAACGGAACgca                                                                                                   <  5:1681/148‑1 (MQ=255)
                                               acaTACAACCCGCCCAGCCAAGGGCAGGGCACATCCGGTACTTGCCGCGCATCCCGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCAGGTAATCCGGGCTGGTGCAACGGAACgcagc                                                                                                 <  8:251042/149‑1 (MQ=255)
                                               acaTACAACCCGCCCAGCCAAGGGCAGGGCACATCCGGTACTTGCCGCGCATCCCGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCAGGTAATCCGGGCTGGTGCAACGGAACgcagc                                                                                                 <  7:210110/149‑1 (MQ=255)
                                                 aTACAACCCGCCCAGCCAAGGGCAGGGCACATCCGGTACTTGCCGCGCATCCCGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCAg                                                                                                                               >  1:266762/1‑117 (MQ=255)
                                                 aTACAACCCGCCCAGCCAAGGGCAGGGCACATCCGGTACTTGCCGCGCATCCCGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCAg                                                                                                                               <  2:266762/117‑1 (MQ=255)
                                                  tACAACCCGCCCAGCCAAGGGCAGGGCACATCCGGTACTTGCCGCGCATCCCGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCAGGTAATCCGGGCTGGTGCAACGGAACGCAGCACg                                                                                              >  1:134535/1‑149 (MQ=255)
                                                    cAACCCGCCCAGCCAAGGGCAGGGCACATCCGGTACTTGCCGCGCATCCCGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCAGGTAATCCGGGCTGGTGCAACGGAACGCAGCACGTc                                                                                            >  7:107048/1‑149 (MQ=255)
                                                     aaCCCGCCCAGCCAAGGGCAGGGCACATCCGGTACTTGCCGCGCATCCCGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCAGGTAATCCGGGCTGGTGCAACGGAACGCAGCACGTcc                                                                                           >  7:46132/1‑149 (MQ=255)
                                                        ccGCCCAGCCAAGGGCAGGGCACATCCGGTACTTGCCGCGCATCCCGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCAGGTAATCCGGGCTGGTGCAACGGAACGCAGCACGTCCCTg                                                                                        <  2:134535/149‑1 (MQ=255)
                                                            ccAGCCAAGGGCAGGGCACATCCGGTACTTGCCGCGCATCCCGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCAGGTAATCCGGGCTGGTGCAACGGAACGCAGCACGTCCCTGCAAt                                                                                    >  1:58698/1‑149 (MQ=255)
                                                                 cAAGGGCAGGGCACATCCGGTACTTGCCGCGCATCCCGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCAGGTAATCCGGGCTGGTGCAACGGAACGCAGCACGTCCCTGCAATTGCTc                                                                               >  6:206443/1‑149 (MQ=255)
                                                                     ggcagggcaCATCCGGTACTTGCCGCGCATCCCGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCAGGTAATCCGGGCTGGTGCAACGGAACGCAGCACGTCCCTGCAATTGCTCgggg                                                                           >  4:249897/1‑149 (MQ=255)
                                                                     ggcagggcaCATCCGGTACTTGCCGCGCATCCCGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCAGGTAATCCGGGCTGGTGCAACGGAACGCAGCACGTCCCTGCAATTGCTCgggg                                                                           <  3:249897/149‑1 (MQ=255)
                                                                            cacaTCCGGTACTTGCCGCGCATCCCGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCAGGTAATCCGGGCTGGTGCAACGGAACGCAGCACGTCCCTGCAATTGCTCGGGGTCCAGtt                                                                    <  5:206443/149‑1 (MQ=255)
                                                                               aTCCGGTACTTGCCGCGCATCCCGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCAGGTAATCCGGGCTGGTGCAACGGAACGCAGCACGTCCCTGCAATTGCTc                                                                               <  3:231855/135‑1 (MQ=255)
                                                                               aTCCGGTACTTGCCGCGCATCCCGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCAGGTAATCCGGGCTGGTGCAACGGAACGCAGCACGTCCCTGCAATTGCTc                                                                               >  4:231855/1‑135 (MQ=255)
                                                                                tCCGGTACTTGCCGCGCATCCCGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCAGGTAATCCGGGCTGGTGCAACGGAACGCAGCACGTCCCTGCAATTGCTCGGGGTCCAGTTcagc                                                                <  6:25253/149‑1 (MQ=255)
                                                                                             gcgcATCCCGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCAGGTAATCCGGGCTGGTGCAACGGAACGCAGCACGTCCCTGCAATTGCTCggg                                                                            >  5:67494/1‑124 (MQ=255)
                                                                                             gcgcATCCCGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCAGGTAATCCGGGCTGGTGCAACGGAACGCAGCACGTCCCTGCAATTGCTCggg                                                                            <  6:67494/124‑1 (MQ=255)
                                                                                                   cccGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCAGGTAATCCGGGCTGGTGCAACGGAACGCAGCACGTCCCTGCAATTGCTCGGGGTCCAGTTCAGCCACAGCCAGCCGTTgcgc                                              >  6:261635/1‑148 (MQ=255)
                                                                                                           aGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCAGGTAATCCGGGCTGGTGCAACGGAACGCAGCACGTCCCTGCAATTGCTCGGGGTCCAGTTCAGCCACAGCCAGCCGTTGCGCCAGATCGAc                                     <  8:107048/149‑1 (MQ=255)
                                                                                                                  cATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCAGGTAATCCGGGCTGGTGCAACGGAACGCAGCACGTCCCTGCAATTGCTCGGGGTCCAGTTCAGCCACAGCCAGCCGTTGCGCCAGATCGACGGAAAtt                              >  7:73528/1‑149 (MQ=255)
                                                                                                                  cATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCAGGTAATCCGGGCTGGTGCAACGGAACGCAGCACGTCCCTGCAATTGCTCGGGGTCCAGTTCAGCCACAGCCAGCCGTTGCGCCAGATCGACGGAAAtt                              >  7:213166/1‑149 (MQ=255)
                                                                                                                      ggCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCAGGTAATCCGGGCTGGTGCAACGGAACGCAGCACGTCCCTGCAATTGCTCGGGGTCCAGTTCAGCCACAGCCAGCCGTTGCGCCAGATCGACGGAAATTTCg                           <  2:6338/148‑1 (MQ=255)
                                                                                                                       gCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCAGGTAATCCGGGCTGGTGCAACGGAACGCAGCACGTCCCTGCAATTGCTCGGGGTCCAGTTCAGCCACAGCCAGCCGTTGCGCCAGATCGACGGAAATTTCGTc                         >  3:127099/1‑149 (MQ=255)
                                                                                                                             gCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCAGGTAATCCGGGCTGGTGCAACGGAACGCAGCACGTCCCTGCAATTGCTCg                                                                              >  5:76960/1‑90 (MQ=255)
                                                                                                                             gCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCAGGTAATCCGGGCTGGTGCAACGGAACGCAGCACGTCCCTGCAATTGCTCg                                                                              <  6:76960/90‑1 (MQ=255)
                                                                                                                              cAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCAGGTAATCCGGGCTGGTGCAACGGAACGCAGCACGTCCCTGCAATTGCTCGGGGTCCAGTTCAGCCACAGCCAGCCGTTGCGCCAGATCGACGGAAATTTCGTCCagcagc                  <  8:46132/149‑1 (MQ=255)
                                                                                                                                   gCCTGCGCGTCGGCGATAGGGCCGACGATCGGCAGGTAATCCGGGCTGGTGCAACGGAACGCAGCACGTCCCTGCAATTGCTCGGGGTCCAGTTCAGCCACAGCCAGCCGTTGCGCCAGATCGACGGAAATTTCGTCCAGCAGCGCCAg             <  2:58698/149‑1 (MQ=255)
                                                                                                                                     cTGCGCGTCGGCGATAGGGCCGACGATCGGCAGGTAATCCGGGCTGGTGCAACGGAACGCAGCACGTCCCTGCAATTGCTCGGGGTCCAGTTCAGCCACAGCCAGCCGTTGCGCCAGATCGACGGAAATTTCGTCCAGCAGCGCCAGGt           >  5:339982/1‑149 (MQ=255)
                                                                                                                                     cTGCGCGTCGGCGATAGGGCCGACGATCGGCAGGTAATCCGGGCTGGTGCAACGGAACGCAGCACGTCCCTGCAATTGCTCGGGGTCCAGTTCAGCCACAGCCAGCCGTTGCGCCAGATCGACGGAAATTTCGTCCAGCAGCGCCAGGt           >  3:45942/1‑149 (MQ=255)
                                                                                                                                       gatcgTCGGCCCTATCGCCGACGATCGGCAGGTAATCCGGGCTGGTGCAACGGAACGCAGCACGTCCCTGCAATTGCTCGGGGTCCAGTTCAGCCACAGCCAGCCGTTGCGCCAGATCGACGGAAATTTCGTCCAGCAGCGCCAGGTTg         >  4:292532/4‑149 (MQ=255)
                                                                                                                                        atcgTCGGCCCTATCGCCGACGATCGGCAGGTAATCCGGGCTGGTGCAACGGAACGCAGCACGTCCCTGCAATTGCTCGGGGTCCAGTTCAGCCACAGCCAGCCGTTGCGCCAGATCGACGGAAATTTCGTCCAGCAGCGCCAGGTTGc        >  5:177977/3‑149 (MQ=255)
                                                                                                                                            tCGGCGATAGGGCCGACGATCGGCAGGTAATCCGGGCTGGTGCAACGGAACGCAGCACGTCCCTGCAATTGCTCGGGGTCCAGTTCAGCCACAGCCAGCCGTTGCGCCAGATCGACGGAAATTTCGTCCAGCAGCGCCAGGTTGCCCTg    >  3:15152/1‑149 (MQ=255)
                                                                                                                                              ggCGATAGGGCCGACGATCGGCAGGTAATCCGGGCTGGTGCAACGGAACGCAGCACGTCCCTGCAATTGCTCGGGGTCCAGTTCAGCCACAGCCAGCCGTTGCGCCAGATCGACGGAAATTTCGTCCAGCAGCGGCAGGTTGCCCTGGt  <  4:15152/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
CACTCAACGGTGCCGTGATCAACCCGCGTGAGCCATGCCCGCTGTTCACATACAACCCGCCCAGCCAAGGGCAGGGCACATCCGGTACTTGCCGCGCATCCCGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCAGGTAATCCGGGCTGGTGCAACGGAACGCAGCACGTCCCTGCAATTGCTCGGGGTCCAGTTCAGCCACAGCCAGCCGTTGCGCCAGATCGACGGAAATTTCGTCCAGCAGCGCCAGGTTGCCCTGGT  >  NC_002947/1951689‑1951979

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: ‑