Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,951,839 G→C 11.6% P568A (CCT→GCT) ‡ mnmC ← methyltransferase/FAD‑dependent oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,951,8390GC11.6% 168.4 / 9.6 59P568A (CCT→GCT) ‡mnmCmethyltransferase/FAD‑dependent oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (24/28);  new base C (4/3);  total (28/31)
Fisher's exact test for biased strand distribution p-value = 6.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTCAACGGTGCCGTGATCAACCCGCGTGAGCCATGCCCGCTGTTCACATACAACCCGCCCAGCCAAGGGCAGGGCACATCCGGTACTTGCCGCGCATCCCGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCAGGTAATCCGGGCTGGTGCAACGGAACGCAGCACGTCCCTGCAATTGCTCGGGGTCCAGTTCAGCCACAGCCAGCCGTTGCGCCAGATCGACGGAAATTTCGTCCAGCAGCGCCAGGTTGCCCTGGTGCTC  >  NC_002947/1951691‑1951983
                                                                                                                                                    |                                                                                                                                                
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                                                                                                     <  7:210110/149‑1 (MQ=255)
                                             acaTACAACCCGCCCAGCCAAGGGCAGGGCACATCCGGTACTTGCCGCGCATCCCGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCAGGTAATCCGGGCTGGTGCAACGGAACgcagc                                                                                                     <  8:251042/149‑1 (MQ=255)
                                               aTACAACCCGCCCAGCCAAGGGCAGGGCACATCCGGTACTTGCCGCGCATCCCGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCAg                                                                                                                                   <  2:266762/117‑1 (MQ=255)
                                               aTACAACCCGCCCAGCCAAGGGCAGGGCACATCCGGTACTTGCCGCGCATCCCGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCAg                                                                                                                                   >  1:266762/1‑117 (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                                                                               <  3:249897/149‑1 (MQ=255)
                                                                   ggcagggcaCATCCGGTACTTGCCGCGCATCCCGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCAGGTAATCCGGGCTGGTGCAACGGAACGCAGCACGTCCCTGCAATTGCTCgggg                                                                               >  4:249897/1‑149 (MQ=255)
                                                                          cacaTCCGGTACTTGCCGCGCATCCCGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCAGGTAATCCGGGCTGGTGCAACGGAACGCAGCACGTCCCTGCAATTGCTCGGGGTCCAGtt                                                                        <  5:206443/149‑1 (MQ=255)
                                                                             aTCCGGTACTTGCCGCGCATCCCGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCAGGTAATCCGGGCTGGTGCAACGGAACGCAGCACGTCCCTGCAATTGCTc                                                                                   >  4:231855/1‑135 (MQ=255)
                                                                             aTCCGGTACTTGCCGCGCATCCCGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCAGGTAATCCGGGCTGGTGCAACGGAACGCAGCACGTCCCTGCAATTGCTc                                                                                   <  3:231855/135‑1 (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)
                                                                                                                                                aTAGGGCCGACGATCGGCAGGTAATCCGGGCTGGTGCAACGGAACGCAGCACGTCCCTGCAATTGCTCGGGGTCCAGTTCAGCCACAGCCAGCCGTTGCGCCAGATCGACGGAAATTTCGTCCAGCAGCGCCAGGTTGCCCTGGTGCTc  <  4:45942/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
CTCAACGGTGCCGTGATCAACCCGCGTGAGCCATGCCCGCTGTTCACATACAACCCGCCCAGCCAAGGGCAGGGCACATCCGGTACTTGCCGCGCATCCCGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCAGGTAATCCGGGCTGGTGCAACGGAACGCAGCACGTCCCTGCAATTGCTCGGGGTCCAGTTCAGCCACAGCCAGCCGTTGCGCCAGATCGACGGAAATTTCGTCCAGCAGCGCCAGGTTGCCCTGGTGCTC  >  NC_002947/1951691‑1951983

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