Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 317,012 C→T 6.4% A307T (GCA→ACA)  dctD‑I ← C4‑dicarboxylate transport transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947317,0120CT 6.4% 184.3 / 2.1 62A307T (GCA→ACA) dctD‑IC4‑dicarboxylate transport transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (28/30);  new base T (2/2);  total (30/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.01e-01

GTTGCAGCGCATGCGGCGACAGGCCATGACGTTGGCTGGCAGCATCGGCGAAATGCTGGAACAACACCAATATGTCGTCGCCGCGTTCACGCAGTGGCGGAATACGCAGCGG‑CGCCACGTTCAGGCGGTAATACAGGTCTGCGCGGAAACGCCCTTGGTCGGCAGACTGACGCAGGTCTTCCTTGGTCGCGGCTATGATGCGGATGTCCAGCGGGATCAACTGGTTGCCGCCCAGGCGTTCGACCACTCGCTCCTGCAGCATGCGCAGCAACTTCACCTGTACATC  >  NC_002947/316871‑317156
                                                                                                                                              |                                                                                                                                                
gTTGCAGCGCATGCGGCGACAGGCCATGACGTTGGCTGGCAGCATCGGCGAAATGCTGGAACAACACCAATATGTCGTCGCCGCGTTCACGCAGTGGCGGAATACGCAGCGG‑CGCCACGTTCAGGCGGTAATACAGGTCTGTGCGGaaa                                                                                                                                           >  1:197067/1‑149 (MQ=255)
     aGCGCATGCGGCGACAGGCCATGACGTTGGCTGGCAGCATCGGCGAAATGCTGGAACAACACCAATATGTCGTCGCCGCGTTCACGCAGTGGCGGAATACGCAGCGG‑CGCCACGTTCAGGCGGTAATACAGGTCTGCGCGGAAACGccc                                                                                                                                      >  3:326068/1‑149 (MQ=255)
      gcgcATGCGGCGACAGGCCATGACGTTGGCTGGCAGCATCGGCGAAATGCTGGAACAACACCAATATGTCGTCGCCGCGTTCACGCAGTGGCGGAATACGCAGCGG‑CGCCACGTTCAGGCGGTAATACAGGTCTGTGCGGAAACGCCCt                                                                                                                                     <  2:197067/149‑1 (MQ=255)
       agcATGCGGCGACAGGCCATGACGTTGGCTGGCAGCATCGGCGAAATGCTGGAACAACACCAATATGTCGTCGCCTCGTTCACGCAGTTGCGGAATTCGCATCGG‑CGCCACGTTCAGGCGGTAATACAGGTCTGCGCGGAAACGCCCtt                                                                                                                                    <  7:138905/148‑1 (MQ=255)
        gcATGCGGCGACAGGCCATGACGTTGGCTGGCAGCATCGGCGAAATGCTGGAACAACACCAATATGTCGTCGCCGCGTTCACGCAGTGGCGGAATACGCAGCGG‑CGCCACGTTCAGGCGGTAATACAGGTCTGCGCGGAAACGCCCTTg                                                                                                                                   <  1:387722/149‑1 (MQ=255)
             cggcgACAGGCCATGACGTTGGCTGGCAGCATCGGCGAAATGCTGGAACAACACCAATATGTCGTCGCCGCGTTCACGCAGTGGCGGAATACGCAGCGG‑CGCCACGTTCAGGTGGTAATACAGGTCTGCGCGGAAACGCCCTTGGTCgg                                                                                                                              >  6:264603/1‑149 (MQ=255)
              ggcgACAGGCCATGACGTTGGCTGGCAGCATCGGCGAAATGCTGGAACAACACCAATATGTCGTCGCCGCGTTCACGCAGTGGCGGAATACGCAGCGG‑CGCCACGTTCAGGCGGTAATACAGGTCTGCGCGGAAACGCCCTTGGTCGGc                                                                                                                             <  4:326068/149‑1 (MQ=255)
               gcgACAGGCCATGACGTTGGCTGGCAGCATCGGCGAAATGCTGGAACAACACCAATATGTCGTCGCCGCGTTCACGCAGTGGCGGAATACGCAGCGG‑CGCCACGTTCAGGCGGTAATACAGGTCTGCGCGGTATCGCCCTTGGTCGGCa                                                                                                                            <  4:255268/149‑1 (MQ=255)
                 gACAGGCCATGACGTTGGCTGGCAGCATCGGCGAAATGCTGGAACAACACCAATATGTCGTCGCCGCGTTCACGCAGTGGCGGAATACGCAGCGG‑CGCCACGTTCAGGCGGTAATACAGGTCTGCGCGGAAACGCCCTTGGTCGGCAGa                                                                                                                          <  6:133188/149‑1 (MQ=255)
                  aCAGGCCATGACGTTGGCTGGCAGCATCGGCGAAATGCTGGAACAACACCAATATGTCGTCGCCGCGTTCACGCAGTGGCGGAATACGCAGCGG‑CGCCACGTTCAGGCGGTAATACAGGTCTGCGCGGAAACGCCCTTGGTCGGCAGAc                                                                                                                         <  4:25320/149‑1 (MQ=255)
                    aGGCCATGACGTTGGCTGGCAGCATCGGCGAAATGCTGGAACAACACCAATATGTCGTCGCCGCGTTCACGCAGTGGCGGAATACGCAGCGG‑CGCCACGTTCAGGCGGTAATACAGGTCTGTGCGGAAACGCCCTTGGTCGGCAGACt                                                                                                                        >  3:292402/1‑148 (MQ=255)
                      gCCATGACGTTGGCTGGCAGCATCGGCGAAATGCTGGAACAACACCAATATGTCGTCGCCGCGTTCACGCAGTGGCGGAATACGCAGCGG‑CGCCACGTTCAGGCGGTAATACAGGTCTGCGCGGAAACGCCCTTGGTCGGCAGACTGAc                                                                                                                     >  4:276164/1‑149 (MQ=255)
                       ccATGACGTTGGCTGGCAGCATCGGCGAAATGCTGGAACAACACCAATATGTCGTCGCCGCGTTCACGCAGTGGCGGAATACGCAGCGG‑CGCCACGTTCAGGCGGTAATACAGGTCTGCGCGGAAACGCCCTTGGTCGGCAGACTGACg                                                                                                                    <  6:110509/149‑1 (MQ=255)
                         aTGACGTTGGCTGGCAGCATCGGCGAAATGCTGGAACAACACCAATATGTCGTCGCCGCGTTCACGCAGTGGCGGAATACGCAGCGG‑CGCCACGTTCAGGCGGTAATACAGGTCTGCGCGGAAACGCCCTTGGTCGGCAGACTGACg                                                                                                                    >  8:47009/1‑147 (MQ=255)
                          tGACGTTGGCTGGCAGCATCGGCGAAATGCTGGAACAACACCAATATGTCGTCGCCGCGTTCACGCAGTGGCGGAATACGCAGCGG‑CGCCACGTTCAGGCGGTAATACAGGTCTGCGCGGAAACGCCCTTGGTCGGCAGACTGACGCAg                                                                                                                 <  4:9011/149‑1 (MQ=255)
                            aCGTTGGCTGGCAGCATCGGCGAAATGCTGGAACAACACCAATATGTCGTCGCCGCGTTCACGCAGTGGCGGAATACGCAGCGG‑CGCCACGTTCAGGCGGTAATACAGGTCTGCGCGGAAACGCCCTTGGTCGGCAGACTGACGCAGGt                                                                                                               <  3:276164/149‑1 (MQ=255)
                               ttGGCTGGCAGCATCGGCGAAATGCTGGAACAACACCAATATGTCGTCGCCGCGTTCACGCAGTGGCGGAATACGCAGCGG‑CGCCACGTTCAGGCGGTAATACAGGTCTGCGCGGAAACGCCCTTGGTCGGCAGACTGACGCAGGTCtt                                                                                                            >  3:42307/1‑149 (MQ=255)
                                 ggctggcAGCATCGGCGAAATGCTGGAACAACACCAATATGTCGTCGCCGCGTTCACGCAGTGGCGGAATACGCAGCGG‑CGCCACGTTCAGGCGGTAATACAGGTCTGCGCGGAAACGCCCTTGGTCGGCAGACTGACGCAGGTCTTc                                                                                                           <  2:335322/148‑1 (MQ=255)
                                            tCGGCGAAATGCTGGAACAACACCAATATGTCGTCGCCGCGTTCACGCAGTGGCGGAATACGCAGCGG‑CGCCACGTTCAGGCGGTAATACAGGTCTGCGCGGAAACGCCCTTGGTCGGCAGACTGACGCAGGTCTTCCTTGGTCGCGGc                                                                                               >  6:141178/1‑149 (MQ=255)
                                                cGAAATGCTGGAACAACACCAATATGTCGTCGCCGCGTTCACGCAGTGGCGGAATACGCAGCGG‑CGCCACGTTCAGGTGGTAATACAGGTCTGCGCGGAAACGCCCTTGGTCGGCAGACTGACGCAGGTCTTCCTTGGTCGCGGCTatg                                                                                           <  5:264603/149‑1 (MQ=255)
                                                cGAAATGCTGGAACAACACCAATATGTCGTCGCCGCGTTCACGCAGTGGCGGAATACGCAGCGG‑CGCCACGTTCAGGCGGTAATACAGGTCTGCGCGGAAACGCCCTTGGTCGGCAGACTGACGCAGGTCTTCCTTGGTCGCGGCTatg                                                                                           <  7:47009/149‑1 (MQ=255)
                                                        tGGAACAACACCAATATGTCGTCGCCGCGTTCACGCAGTGGCGGAATACGCAGCGG‑CGCCACGTTCAGGCGGTAATACAGGTCTGCGCGGAAACGCCCTTGGTCGGCAGACTGACGCAGGTCTTCCTTGGTCGCGGCTATGATGCgg                                                                                     >  8:109564/1‑147 (MQ=255)
                                                        tGGAACAACACCAATATGTCGTCGCCGCGTTCACGCAGTGGCGGAATACGCAGCGG‑CGCCACGTTCAGGCGGTAATACAGGTCTGCGCGGAAACGCCCTTGGTCGGCAGACTGACGCAGGTCTTCCTTGGTCGCGGCTATGATGCgg                                                                                     <  7:109564/147‑1 (MQ=255)
                                                        tGGAACAACACCAATATGTCGTCGCCGCGTTCACGCAGTGGCGGAATACGCAGCGG‑CGCCACGTTCAGGCGGTAATACAGGTCTGCGCGGAAACGCCCTTGGTCGGCAGACTGACGCAGGTCTTCCTTGGTAGCGGCTATGATGCGGAt                                                                                   >  7:201248/1‑149 (MQ=255)
                                                          gAACAACACCAATATGTCGTCGCCGCGTTCACGCAGTGGCGGAATACGCAGCGG‑CGCCACGTTCAGGCGGTAATACAGGTCTGCGCGGAAACGCCCTTGGTCGGCAGACTGACGCAGGTCTTCCTTGGTCGCGGCTATGATGCGGAt                                                                                   <  7:342573/147‑1 (MQ=255)
                                                            acaacaCCAATATGTCGTCGCCGCGTTCACGCAGTGGCGGAATACGCAGCGG‑CGCCACGTTCAGGCGGTAATACAGGTCTGCGCGGAAACGCCCTTGGTCGGCAGACTGACGCAGGTCTTCCTTGGTCGCGGCTATGATGCGGATGTc                                                                                >  2:361692/1‑148 (MQ=255)
                                                                cacCAATATGTCGTCGCCGCGTTCACGCAGTGGCGGAATACGCAGCGG‑CGCCACGTTCAGGCGGTAATACAGGTCTGCGCGGAAACGCCCTTGGTCGGCAGACTGACGCAGGTCTTCCTTGGTCGCGGCTATGATGCGGa                                                                                    >  6:204366/1‑140 (MQ=255)
                                                                cacCAATATGTCGTCGCCGCGTTCACGCAGTGGCGGAATACGCAGCGG‑CGCCACGTTCAGGCGGTAATACAGGTCTGCGCGGAAACGCCCTTGGTCGGCAGACTGACGCAGGTCTTCCTTGGTCGCGGCTATGATGCGGa                                                                                    <  5:204366/140‑1 (MQ=255)
                                                                cacCAATATGTCGTCGCCGCGTTCACGCAGTGGCGGAATACGCAGCGG‑CGCCACGTTCAGGCGGTAATACAGGTCTGCGCGGAAACGCCCTTGGTCGGCAGACTGACGCAGGTCTTCCTTGGTCGCGGCTATGATGCGGATGTCCAGCg                                                                           <  7:48461/149‑1 (MQ=255)
                                                                cacCAATATGTCGTCGCCGCGTTCACGCAGTGGCGGAATACGCAGCGG‑CGCCACGTTCAGGCGGTAATACAGGTCTGCGCGGAAACGCCCTTGGTCGGCAGACTGACGCAGGTCTTCCTTGGTCGCGGCTATGATGCGGATGTCCAGCg                                                                           >  6:260376/1‑149 (MQ=255)
                                                                   cAATATGTCGTCGCCGCGTTCACGCAGTGGCGGAATACGCAGCGG‑CGCCACGTTCAGGCGGTAATACAGGTCTGCGCGGAAACGCCCTTGGTCGGCAGACTGACGCAGGTCTTCCTTGGTCGCGGCTATGATGCGGATGTCCAGctctg                                                                        >  2:232415/1‑145 (MQ=255)
                                                                   cAATATGTCGTCGCCGCGTTCACGCAGTGGCGGAATACGCAGCGG‑CGCCACGTTCAGGCGGTAATACAGGTCTGCGCGGAAACGCCCTTGGTCGGCAGACTGACGCAGGTCTTCCTTGGTCGCGGCTATGATGCGGATGTCCAGCGGGa                                                                        >  7:363220/1‑149 (MQ=255)
                                                                   cAATATGTCGTCGCCGCGTTCACGCAGTGGCGGAATACGCAGCGG‑CGCCACGTTCAGGCGGTAATACAGGTCTGCGCGGAAACGCCCTTGGTCGGCAGACTGACGCAGGTCTTCCTTGGTCGCGGCTATGATGCGGATGTCCAGCGGGa                                                                        <  3:175438/149‑1 (MQ=255)
                                                                    aaTATGTCGTCGCCGCGTTCACGCAGTGGCGGAATACGCAGCGG‑CGCCACGTTCAGGCGGTAATACAGGTCTGCGCGGAAACGCCCTTGGTCGGCAGACTGACGCAGGTCTTCCTTGGTCGCGGCTATGATGCGGATGTCCAGCGGGAt                                                                       <  4:146864/149‑1 (MQ=255)
                                                                         gtcgtcGCCGCGTTCACGCAGTGGCGGAATACGCAGCGG‑CGCCACGTTCAGGCGGTAATACAGGTCTGCGCGGAAACGCCCTTGGTCGGCAGACTGACGCAGGTCTTCCTTGGTCGCGGCTATGATGCGGATGTCCAGCGGGATCAACt                                                                  >  2:24337/1‑149 (MQ=255)
                                                                              cgccgcGTTCACGCAGTGGCGGAATACGCAGCGG‑CGCCACGTTCAGGCGGTAATACAGGTCTGCGCGGAAACGCCCTTGGTCGGCAGACTGACGCAg                                                                                                                 >  4:232081/1‑97 (MQ=255)
                                                                              cgccgcGTTCACGCAGTGGCGGAATACGCAGCGG‑CGCCACGTTCAGGCGGTAATACAGGTCTGCGCGGAAACGCCCTTGGTCGGCAGACTGACGCAg                                                                                                                 <  3:232081/97‑1 (MQ=255)
                                                                               gccgcGTTCACGCAGTGGCGGAATACGCAGCGGACGCCACGTTCAGGCGGTAATACAGGTCTGCGCGGAAACGCCGTTGGTCGGCAGCCTGACGCAGGTCTTCCTTGGTCGCGGCTATGATGCGGGTGTCCAGCGGGATCAACTGGTTg                                                             >  1:258882/1‑149 (MQ=39)
                                                                                ccgcGTTCACGCAGTGGCGGAATACGCAGCGG‑CGCCACGTTCAGGCGGTAATACAGGTCTGCGCGGAAACGCCCTTGGTCGGCAGACTGACGCAGGTCTTCCTTGGTCGCGGCTATGATGCGGATGTCCAGCGGGATCAACTGGTTgcc                                                           >  7:17759/1‑149 (MQ=255)
                                                                                 cgcgTTCACGCAGTGGCGGAATACGCAGCGG‑CGCCACGTTCAGGCGGTAATACAGGTCTGCGCGGAAACGCCCTTGGTCGGCAGACTGACGCAGGTCTTCCTTGGTCGCGGCTATCATGCGGATGTCCAGCGGGATCAACTGGTTgacg                                                          >  6:219438/1‑149 (MQ=255)
                                                                                    gTTCACGCAGTGGCGGAATACGCAGCGG‑CGCCACGTTCAGGCGGTAATACAGGTCTGCGCGGAAACGCCCTTGGTCGGCAGACTGACGCAGGTCTTCCTTGGTCGCGGCTATGATGCGGATGTCCAGCGGGATCAACTGGTTgccgcc                                                        >  7:7474/1‑148 (MQ=255)
                                                                                        aCGCAGTGGCGGAATACGCAGCGG‑CGCCACGTTCAGGCGGTAATACAGGTCTGCGCGGAAACGCCCTTGGTCGGCAGACTGACGCAGGTCTTCCTTGGTCGCGGCTATGATGCGGATGTCCAGCGGGATCAACTGGTTGCCGCCCAGGc                                                   >  2:151420/1‑149 (MQ=255)
                                                                                          gCAGTGGCGGAATACGCAGCGG‑CGCCACGTTCAGGCGGTAATACAGGTCTGCGCGGAAACGCCCTTGGTCGGCAGACTGACGCAGGTCTTCCTTGGTCGCGGCTATGATGCGGATGTCCAGCGGGATCAACTGGTTGCCGCCCAGGCGt                                                 <  5:260376/149‑1 (MQ=255)
                                                                                            aGTGGCGGAATACGCAGCGG‑CGCCACGTTCAGGCGGTAATACAGGTCTGTGCGGAAACGCCCTTGGTCGGCAGACTGACGCAGGTCTTCCTTGGTCGCGGCTATGATGCGGATGTCCAGCGGGATCAACTGGTTGCCGCCCAGGCGTTc                                               <  4:292402/149‑1 (MQ=255)
                                                                                              tGGCGGAATACGCAGCGG‑CGCCACGTTCAGGCGGTAATACAGGTCTGCGCGGAAACGCCCTTGGTCGGCAGACTGACGCAGGTCTTCCTTGGTCGCGGCTATGATGCGGATGTCCAGCGGGATCAACTGGTTGCCGCCCAGGCGTTCGa                                             <  8:363220/149‑1 (MQ=255)
                                                                                               ggCGGAATACGCAGCGG‑CGCCACGTTCAGGCGGTAATACAGGTCTGCGCGGAAACGCCCTTGGTCGGCAGACTGACGCAGGTCTTCCTTGGTCGCGGCTATGATGCGGATGTCCAGCGGGATCAACTGGTTGCCGCCCAGGCGTTc                                               >  6:154151/1‑146 (MQ=255)
                                                                                               ggCGGAATACGCAGCGG‑CGCCACGTTCAGGCGGTAATACAGGTCTGCGCGGAAACGCCCTTGGTCGGCAGACTGACGCAGGTCTTCCTTGGTCGCGGCTATGATGCGGATGTCCAGCGGGATCAACTGGTTGCCGCCCAGGCGTTc                                               <  5:154151/146‑1 (MQ=255)
                                                                                                  ggAATACGCAGCGG‑CGCCACGTTCAGGCGGTAATACAGGTCTGCGCGGAAACGCCCTTGGTCGGCAGACTGACGCAGGTCTTCCTTGGTCGCGGCTATGATGCGGATGTCCAGCGGGATCAACTGGTTGCCGCCCAGGCGTTCGACCAc                                         <  5:141178/149‑1 (MQ=255)
                                                                                                           aGCGG‑CGCCACGTTCAGGCGGTAATACAGGTCTGCGCGGAAACGCCCTTGGTCGGCAGACTGACGCAGGTCTTCCTTGGTCGCGGCTATGATGCGGATGTCCAGCGGGATCAACTGGTTGCCGCCCAGGCGTTCGACCACTCGCTCCt                                 >  2:278062/1‑148 (MQ=255)
                                                                                                               g‑cgcCACGTTCAGGCGGTAATACAGGTCTGCGCGGAAACGCCCTTGGTCGGCAGACTGACGCAGGTCTTCCTTGGTCGCGGCTATGATGCGGATGTCCAGCGGGATCAACTGGTTGCCGCCCAGGCGTTCGACCACTCGCTCCTgcagc                            <  1:278062/149‑1 (MQ=255)
                                                                                                                   ccACGTTCAGGCGGTAATACAGGTCTGCGCGGAAACGCCCTTGGTCGGCAGACTGACGCAGGTCTTCCTTGGTCGCGGCTATGATGCGGATGTCCAGCGGGATCAACTGGTTGCCGCCCAGGCGTTCGACCACTCGCTCCTGCAGCATg                         >  4:129892/1‑149 (MQ=255)
                                                                                                                    cACGTTCAGGCGGTAATACAGGTCTGCGCGGAAACGCCCTTGGTCGGCAGACTGACGCAGGTCTTCCTTGGTCGCGGCTATGATGCGGATGTCCAGCGGGATCAACTGGTTGCCGCCCAGGCGTTCGAc                                            >  7:358337/1‑129 (MQ=255)
                                                                                                                    cACGTTCAGGCGGTAATACAGGTCTGCGCGGAAACGCCCTTGGTCGGCAGACTGACGCAGGTCTTCCTTGGTCGCGGCTATGATGCGGATGTCCAGCGGGATCAACTGGTTGCCGCCCAGGCGTTCGAc                                            <  8:358337/129‑1 (MQ=255)
                                                                                                                    cACGTTCAGGCGGTAATACAGGTCTGCGCGGAAACGCCCTTGGTCGGCAGACTGACACAGGTCTTCCTTGGTCGCGGCTATGATGCGGATGTCCAGCGGGATCAACTg                                                                 <  7:333652/108‑1 (MQ=255)
                                                                                                                    cACGTTCAGGCGGTAATACAGGTCTGCGCGGAAACGCCCTTGGTCGGCAGACTGACACAGGTCTTCCTTGGTCGCGGCTATGATGCGGATGTCCAGCGGGATCAACTg                                                                 >  8:333652/1‑108 (MQ=255)
                                                                                                                           aGGCGGTAATACAGGTCTGCGCGGAAACGCCCTTGGTCGGCAGACTGACGCAGGTCTTCCTTGGTCGCGGCTATGATGCGGATGTCCAGCGGGATCAACTGGTTGCCGCCCAGGCGTTCGACCACTCGCTCCTGCAGCATGCGCAGCaa                 <  3:129892/149‑1 (MQ=255)
                                                                                                                                      cAGGTCTGCGCGGAAACGCCCTTGGTCGGCAGACTGACGCAGGTCTTCCTTGGTCGCGGCTATGATGCGGATGTCCAGCGGGATCAACTGGTTGCCGCCCAGGCGTTCGACCACTCGCTCCTGCAGCATGCGCAGCAACTTCACCTGTa      <  1:151420/149‑1 (MQ=255)
                                                                                                                                      cAGGTCTGCGCGGAAACGCCCTTGGTCGGCAGACTGACGCAGGTCTTCCTTGGTCGCGGCTATGATGCGGATGTCCAGCGGGATCAACTGGTTGCCGCCCAGGCGTTCGACCACTCGCTCCTGCAGCATGCGCAGCAACTTCACCTGTa      >  8:26678/1‑149 (MQ=255)
                                                                                                                                          tCTGCGCGGAAACGCCCTTGGTCGGCAGACTGACGCAGGTCTTCCTTGGTCGCGGCTATGATGCGGATGTCCAGCGGGATCAACTGGTTGCCGCCCAGGCGTTCGACCACTCGCTCCTGCAGCATGCGCAGCAACTTCACCTGTACATc  >  8:223221/1‑149 (MQ=255)
                                                                                                                                          tCTGCGCGGAAACGCCCTTGGTCGGCAGACTGACGCAGGTCTTCCTTGGTCGCGGCTATGATGCGGATGTCCAGCGGGATCAACTGGTTGCCGCCCAGGCGTTCGACCACTCGCTCCTGCAGCATGCGCAGCAACTTCACCTGTACATc  <  4:380534/149‑1 (MQ=255)
                                                                                                                                          tCTGCGCGGAAACGCCCTTGGGCGGCAGACTGACGCAGGTCTTCCTTGGTCGCGGCTATGATGCGGATGTCCAGCGGGATCAACTGGTTGCCGCCCAGGCGTTCGACCACTCGCTCCTGCAGCATGCGCAGCAACTTCACCTGTACATc  <  4:255219/149‑1 (MQ=255)
                                                                                                                                            tGCGCGGAAACGCCCTTGGTCGGCAGACTGACGCAGGTCTTCCTTGGTCGCGGCTATGATGCGGATGTCCAGCGGGATCAACTGGTTGCCGCCCAGGCGTTCGACCACTCGCTCCTGCAGCATGCGCAGCAACTTCACCTGTACATc  >  5:346057/1‑147 (MQ=255)
                                                                                                                                              |                                                                                                                                                
GTTGCAGCGCATGCGGCGACAGGCCATGACGTTGGCTGGCAGCATCGGCGAAATGCTGGAACAACACCAATATGTCGTCGCCGCGTTCACGCAGTGGCGGAATACGCAGCGG‑CGCCACGTTCAGGCGGTAATACAGGTCTGCGCGGAAACGCCCTTGGTCGGCAGACTGACGCAGGTCTTCCTTGGTCGCGGCTATGATGCGGATGTCCAGCGGGATCAACTGGTTGCCGCCCAGGCGTTCGACCACTCGCTCCTGCAGCATGCGCAGCAACTTCACCTGTACATC  >  NC_002947/316871‑317156

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