Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,074,707 A→T 11.5% L50Q (CTG→CAG)  gatA ← glutamyl‑tRNA(Gln) amidotransferase subunit A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,074,7070AT11.5% 164.3 / 9.8 60L50Q (CTG→CAG) gatAglutamyl‑tRNA(Gln) amidotransferase subunit A
Reads supporting (aligned to +/- strand):  ref base A (20/33);  new base T (3/4);  total (23/37)
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.37e-01

GCGCCTTGAAGTTGTCGAGCATCTTCGAGCCGCAGCTGGTGCGTACGCCGTTGGTGCAGAACAGGTCTTTGTGGGCGATCGGTGCACCCAGCAGCGCGCTGGTTTCGCCCGCGGCGCGACGGGCATCGGCGGCACGCGCCTGGCCCAGGGCCAGGTCTTCGGTGACGCTGATGAAGCTGTTGATCTGTGGGTCGAGCTGCTTGATGCGCGCCAGCAGGGCGCCGGTCAGCTCTTCGGAAGAAAACGATTTGTCGGCGAGTCCGCGGGCGATCTCGGCCAGGGTCAATTGATG  >  NC_002947/1074561‑1074852
                                                                                                                                                  |                                                                                                                                                 
gcgcCTTGAAGTTGTCGAGCATCTTCGAGCCGCAGCTGGTGCGTACGCCGTTGGTGCAGAACAGGTCTTTGTGGGCGATCGGTGCACCCAGCAGCGCGCTGGTTTCGCCCGCGGCGCGACGGGCATCGGCGGCACGCGCCTGGCCCAgg                                                                                                                                                 <  8:190836/149‑1 (MQ=255)
gcgcCTTGAAGTTGTCGAGCATCTTCGAGCCGCAGCTGGTGCGTACGCCGTTGGTGCAGAACAGGTCTTTGTGGGCGATCGGTGCACCCAGCAGCGCGCTGGTTTCGCCCGCGGCGCGACGGGCATCGGCGGCACGCGCCTGGCCCAgg                                                                                                                                                 >  1:265881/1‑149 (MQ=255)
gcgcCTTGAAGTTGTCGAGCATCTTCGAGCCGCAGCTGGTGCGTACGCCGTTGGTGCAGAACAGGTCTTTGTGGGCGATCGGTGCACCCAGCAGCGCGCTGGTTTCGCCCGCGGCGCGACGGGCATCGGCGGCACGCGCCTGGCCCAgg                                                                                                                                                 <  2:265881/149‑1 (MQ=255)
    cTTGAAGTTGTCGAGCATCTTCGAGCCGCAGCTGGTGCGTACGCCGTTGGTGCAGAACAGGTCTTTGTGGGCGATCGGTGCACCCAGCAGCGCGCTGGTTTCGCCCGCGGCGCGACGGGCATCGGCGGCACGCGCCTGGCCCAGGGCCa                                                                                                                                             <  3:61413/149‑1 (MQ=255)
         aGTTGTCGAGCATCTTCGAGCCGCAGCTGGTGCGTACGCCGTTGGTGCAGAACAGGTCTTTGTGGGCGATCGGTGCACCCAGCAGCGCGCTGGTTTCGCCCGCGGCGCGACGGGCATCGGCGGCACGCGCCTGGCCCTGGGCCAggcgc                                                                                                                                        <  7:156835/149‑4 (MQ=255)
          gTTGTCGAGCATCTTCGAGCCGCAGCTGGTGCGTACGCCGTTGGTGCAGAACAGGTCTTTGTGGGCGATCGGTGCACCCAGCAGCGCGCTGGTTTCGCCCGCGGCGCGACGGGCATCGGCGGCACGCGCCTGGCCCTGGGCCAggc                                                                                                                                          >  1:17064/1‑145 (MQ=255)
          gTTGTCGAGCATCTTCGAGCCGCAGCTGGTGCGTACGCCGTTGGTGCAGAACAGGTCTTTGTGGGCGATCGGTGCACCCAGCAGCGCGCTGGTTTCGCCCGCGGCGCGACGGGCATCGGCGGCACGCGCCTGGCCCTGGGCCAggc                                                                                                                                          <  2:17064/146‑2 (MQ=255)
          gTTGTCGAGCATCTTCGAGCCGCAGCTGGTGCGTACGCCGTTGGTGCAGAACAGGTCTTTGTGGGCGATCGGTGCACCCAGCAGCGCGCTGGTTTCGCCCGCGGCGCGACGGGCATCGGCGGCACGCGCCTGGCCCTGGGCCAggc                                                                                                                                          <  4:228632/146‑2 (MQ=255)
          gTTGTCGAGCATCTTCGAGCCGCAGCTGGTGCGTACGCCGTTGGTGCAGAACAGGTCTTTGTGGGCGATCGGTGCACCCAGCAGCGCGCTGGTTTCGCCCGCGGCGCGACGGGCATCGGCGGCACGCGCCTGGCCCTGGGCCAggc                                                                                                                                          >  3:228632/1‑145 (MQ=255)
             gTCGAGCATCTTCGAGCCGCAGCTGGTGCGTACGCCGTTGGTGCAGAACAGGTCTTTGTGGGCGATCGGTGCACCCAGCAGCGCGCTGGTTTCGCCCGCGGCGCGACGGGCATCGGCGGCACGCGCCTGGCCCAGGGCCAGGTCTTCgg                                                                                                                                    >  4:71877/1‑149 (MQ=255)
             gTCGAGCATCTTCGAGCCGCAGCTGGTGCGTACGCCGTTGGTGCAGAACAGGTCTTTGTGGGCGATCGGTGCACCCAGCAGCGCGCTGGTTTCGCCCGCGGCGCGACGGGCATCGGCGGCACGCGCCTGGCCCAGGGCCAGGTCTTCg                                                                                                                                     <  5:92141/148‑1 (MQ=255)
                gAGCATCTTCGAGCCGCAGCTGGTGCGTACGCCGTTGGTGCAGAACAGGTCTTTGTGGGCGATCGGTGCACCCAGCAGCGCGCTGGTTTCGCCCGCGGCGCGACGGGCATCGGCGGCACGCGCCTGGCCCAGGGCCAGGTCTTCGGTg                                                                                                                                  <  2:196841/148‑1 (MQ=255)
                gAGCATCTTCGAGCCGCAGCTGGTGCGTACGCCGTTGGTGCAGAACAGGTCTTTGTGGGCGATCGGTGCACCCAGCAGCGCGCTGGTTTCGCCCGCGGCGCGACGGGCATCGGCGGCACGCGCCTGGCCCAGGGCCAGGTCTTCGGTGa                                                                                                                                 <  5:333034/149‑1 (MQ=255)
                   cATCTTCGAGCCGCAGCTGGTGCGTACGCCGTTGGTGCAGAACAGGTCTTTGTGGGCGATCGGTGCACCCAGCAGCGCGCTGGTTTCGCCCGCGGCGCGACGGGCATCGGCGGCACGCGCCTGGCCCAGGGCCAGGTc                                                                                                                                         >  8:302772/1‑138 (MQ=255)
                   cATCTTCGAGCCGCAGCTGGTGCGTACGCCGTTGGTGCAGAACAGGTCTTTGTGGGCGATCGGTGCACCCAGCAGCGCGCTGGTTTCGCCCGCGGCGCGACGGGCATCGGCGGCACGCGCCTGGCCCAGGGCCAGGTc                                                                                                                                         <  7:302772/138‑1 (MQ=255)
                   cATCTTCGAGCCGCAGCTGGTGCGTACGCCGTTGGTGCAGAACAGGTCTTTGTGGGCGATCGGTGCACCCAGCAGCGCGCTGGTTTCGCCCGCGGCGCGACGGGCATCGGCGGCACGCGCCTGGCCCAGGGCCAGGTCTTCGGTGACGc                                                                                                                              >  4:285283/1‑149 (MQ=255)
                    aTCTTCGAGCCGCAGCTGGTGCGTACGCCGTTGGTGCAGAACAGGTCTTTGTGGGCGATCGGTGCACCCAGCAGCGCGCTGGTTTCGCCCGCGGCGCGACGGGCATCGGCGGCACGCGCCTGGCCCAGGGCCAGGTCTTCGGTGACGCt                                                                                                                             >  4:292312/1‑149 (MQ=255)
                      cTTCGAGCCGCAGCTGGTGCGTACGCCGTTGGTGCAGAACAGGTCTTTGTGGGCGATCGGTGCACCCAGCAGCGCGCTGGTTTCGCCCGCGGCGCGACGGGCATCGGCGGCACGCGCCTGGCCCAGGGCCAGGTCTTCGGTGACGCtga                                                                                                                           >  6:25189/1‑149 (MQ=255)
                         cGAGCCGCAGCTGGTGCGTACGCCGTTGGTGCAGAACAGGTCTTTGTGGGCGATCGGTGCACCCAGCAGCGCGCTGGTTTCGCCCGCGGCGCGACGGGCATCGGCGGCACGCGCCTGGCCCAGGGCCAGGTCTTCGGTGACGCtgat                                                                                                                          <  7:89357/147‑1 (MQ=255)
                         cGAGCCGCAGCTGGTGCGTACGCCGTTGGTGCAGAACAGGTCTTTGTGGGCGATCGGTGCACCCAGCAGCGCGCTGGTTTCGCCCGCGGCGCGACGGGCATCGGCGGCACGCGCCTGGCCCAGGGCCAGGTCTTCGGTGACGCtgat                                                                                                                          <  7:13237/147‑1 (MQ=255)
                         cGAGCCGCAGCTGGTGCGTACGCCGTTGGTGCAGAACAGGTCTTTGTGGGCGATCGGTGCACCCAGCAGCGCGCTGGTTTCGCCCGCGGCGCGACGGGCATCGGCGGCACGCGCCTGGCCCAGGGCCAGGTCTTCGGTGACGCtga                                                                                                                           <  2:222845/146‑1 (MQ=255)
                            gCCGCAGCTGGTGCGTACGCCGTTGGTGCAGAACAGGTCTTTGTGGGCGATCGGTGCACCCAGCAGCGCGCTGGTTTCGCCCGCGGCGCGACGGGCATCGGCGGCACGCGCCTGGCCCAGGGCCAGGTCTTCGGTGACGCTGATGAAGc                                                                                                                     >  3:126466/1‑149 (MQ=255)
                             ccGCAGCTGGTGCGTACGCCGTTGGTGCAGAACAGGTCTTTGTGGGCGATCGGTGCACCCAGCAGCGCGCTGGTTTCGCCCGCGGCGCGACGGGCATCGGCGGCACGCGCCTGGCCCAGGGCCAGGTCTTCGGTGACGCTGATGAAGCt                                                                                                                    <  3:292312/149‑1 (MQ=255)
                               gCAGCTGGTGCGTACGCCGTTGGTGCAGAACAGGTCTTTGTGGGCGATCGGTGCACCCAGCAGCGCGCTGGTTTCGCCCGCGGCGCGACGGGCATCGGCGGCACGCGCCTGGCCCAGGGCCAGGTCTTCGGTGACGCTGATGAAGCTGt                                                                                                                  <  7:36565/149‑1 (MQ=255)
                                cAGCTGGTGCTTACGCCGTTGGTGCAGAACAGGTCTTTGTGGGCGATCGGTGCACCCAGCAGCGCGCTGGTTTCGCCCGCGGCGCGACGGGCATCGGCGGCACGCGCCTGGCCCAGGGCCAGGTCTTCGGTGACGCTGATGAAGCTGtt                                                                                                                 <  1:160073/149‑1 (MQ=255)
                                cAGCTGGTGCGTACGCCGTTGGTGCAGAACAGGTCTTTGTGGGCGATCGGTGCACCCAGCAGCGCGCTGGTTTCGCCCGCGGCGCGACGGGCATCGGCGGCACGCGCCTGGCCCAGGGCCAGGTCTTCGGTGACGCTGATGAAGCTGtt                                                                                                                 <  7:241797/149‑1 (MQ=255)
                                           tACGCCGTTGGTGCAGAACAGGTCTTTGTGGGCGATCGGTGCACCCAGCAGCGCGCTGGTTTCGCCCGCGGCGCGACGGGCATCGGCGGCACGCGCCTGGCCCAGGGCCAGGTCTTCGGTGACGCTGATGAAGCTGTTGATCTGTGGGt                                                                                                      >  6:223613/1‑149 (MQ=255)
                                              gCCGTTGGTGCAGAACAGGTCTTTGTGGGCGATCGGTGCACCCAGCAGCGCGCTGGTTTCGCCCGCGGCGCGACGGGCATCGGCGGCACGCGCCTGGCCCAGGGCCAGGTCTTCGGTGACGCTGATGAAGCTGTTGATCTGTGGGTc                                                                                                     >  4:169305/1‑147 (MQ=255)
                                                cGTTGGTGCAGAACAGGTCTTTGTGGGCGATCGGTGCACCCAGCAGCGCGCTGGTTTCGCCCGCGGCGCGACGGGCATCGGCGGCACGCGCCTGGCCCAGGGCCAGGTCTTCGGTGACGCTGATGAAGCTGTTGATCTGTGGGTCGAgc                                                                                                 >  3:199409/1‑149 (MQ=255)
                                                        cAGAACAGGTCTTTGTGGGCGATCGGTGCACCCAGCAGCGCGCTGGTTTCGCCCGCGGCGCGACGGGCATCGGCGGCACGCGCCTGGCCCAGGGCCAGGTCTTCGGTGACGCTGATGAAGCTGTTGATCTGTGGGTCGAGCTGCTTGAt                                                                                         <  8:42617/149‑1 (MQ=255)
                                                            aCAGGTCTTTGTGGGCGATCGGTGCACCCAGCAGCGCGCTGGTTTCGCCCGCGGCGCGACGGGCATCGGCGGCACGCGCCTGGCCCAGGGCCAGGTCTTCGGTGACGCTGATGAAGCTGTTGATCTGTGGGTCGAGCTGCTTGATgcgc                                                                                     <  3:71877/149‑1 (MQ=255)
                                                            aCAGGTCTTTGTGGGCGATCGGTGCACCCAGCAGCGCGCTGGTTTCGCCCGCGGCGCGACGGGCATCGGCGGCACGCGCCTGGCCCAGGGCCAGGTCTTCGGTGACGCTGATGAAGCTGTTGATCTGTGGGTCGAGCTGCTTGATgcgc                                                                                     <  4:126466/149‑1 (MQ=255)
                                                              aGGTCTTTGTGGGCGATCGGTGCACCCAGCAGCGCGCTGGTTTCGCCCGCGGCGCGACGGGCATCGGCGGCACGCGCCTGGCCCAGGGCCAGGTCTTCGGTGACGCTGATGAAGCTGTTCATCTGTGGGTCGAGCTGCTTGATgcgcgc                                                                                   >  3:198351/1‑149 (MQ=255)
                                                                gTCTTTGTGGGCGATCGGTGCACCCAGCAGCGCGCTGGTTTCGCCCGCGGCGCGACGGGCATCGGCGGCACGCGCCTGGCCCAGGGCCAGGTCTTCGGTGACGCTGATGAAGCTGTTGATCTGTGGGTCGAGCTGCTTGATgcgcgc                                                                                   >  1:339945/1‑147 (MQ=255)
                                                                   tttGTGGGCGATCGGTGCACCCAGCAGCGCGCTGGTTTCGCCCGCGGCGCGACGGGCATCGGCGGCACGCGCCTGGCCCAGGGCCAGGTCTTCGGTGACGCTGATGAAGCTGTTGATCTGTGGGTCGAGCTGCTTGATGCGCGCcag                                                                                >  5:196995/1‑147 (MQ=255)
                                                                             aTCGGTGCACCCAGCAGCGCGCTGGTTTCGCCCGCGGCGCGACGGGCATCGGCGGCACGCGCCTGGCCCAGGGCCAGGTCTTCGGTGACGCTGATGAAGCTGTTGATCTGTGGGTCGAGCTGCTTGATGCGCGCCAGCAGGGCGCCGGt                                                                    >  6:237117/1‑149 (MQ=255)
                                                                                ggTGCACCCAGCAGCGCGCTGGTTTCGCCCGCGGCGCGACGGGCATCGGCGGCACGCGCCTGGCCCAGGGCCAGGTCTTCGGTGACGCTGATGAAGCTGTTGATCTGTGGGTCGAGCTGCTTGATGCGCGCCAGCAGGGCGCCGGCCAg                                                                 <  5:223613/149‑1 (MQ=255)
                                                                                 gTGCACCCAGCAGCGCGCTGGTTTCGCCCGCGGCGCGACGGGCATCGGCGGCACGCGCCTGGCCCAGGGCCAGGTCTTCGGTGACGCTGATGAAGCTGTTGATc                                                                                                             <  8:301496/104‑1 (MQ=255)
                                                                                 gTGCACCCAGCAGCGCGCTGGTTTCGCCCGCGGCGCGACGGGCATCGGCGGCACGCGCCTGGCCCAGGGCCAGGTCTTCGGTGACGCTGATGAAGCTGTTGATc                                                                                                             >  7:301496/1‑104 (MQ=255)
                                                                                 gTGCACCCAGCAGCGCGCTGGTTTCGCCCGCGGCGCGACGGGCATCGGCGGCACGCGCCTGGCCCAGGGCCAGGTCTTCGGTGACGCTGATGAAGCTGTTGATCTGTGGGTCGAGCTGCTTGATGCGCGCCAGCAGGGCGCCGGTCAGc                                                                <  6:161336/149‑1 (MQ=255)
                                                                                        cagcagCGCGCTGGTTTCGCCCGCGGCGCGACGGGCATCGGCGGCACGCGCCTGGCCCAGGGCCAGGTCTTCGGTGACGCTGATGAAGCTGTTGATCTCTGGGTCGAGCTGCTTGATGCGCGCCAGCAGGGCGCCGGTCAGCTCTTCg                                                          >  7:94493/1‑148 (MQ=255)
                                                                                         aggagcGCGCTGGTTTGGCCCGCGGCGCGACGGGCATCGGCGGCACGCGCCTGGCCCAGGGCCAGGTCTTCGGTGACGCTGATGAAGCTGTTGATCTGTGGGTCGAGCTGCTTGATGCGCGCCAGCAGGGCGCCGGTCAGCTCTTCGga                                                        <  1:17211/146‑1 (MQ=255)
                                                                                          gcagcGCGCTGGTTTCGCCCGCGGCGCGACGGGCATCGGCGGCACGCGCCTGGCCCAGGGCCAGGTCTTCGGTGACGCTGATGAAGCTGTTGATCTGTGGGTCGAGCTGCTTGATGCGCGCCAGCAGGGCGCCGGTCAGCTCTTCGgaa                                                       <  2:116051/149‑1 (MQ=255)
                                                                                          gcagcGCGCTGGTTTCGCCCGCGGCGCGACGGGCATCGGCGGCACGCGCCTGGCCCAGGGCCAGGTCTTCGGTGACGCTGATGAAGCTGTTGATCTGTGGGTCGAGCTGCTTGATGCGCGCCAGCAGGGCGCCGGTCAGCTCTTCGgaa                                                       <  3:285283/149‑1 (MQ=255)
                                                                                               gcgcTGGTTTCGCCCGCGGCGCGACGGGCATCGGCGGCACGCGCCTGGCCCAGGGCCAGGTCTTCGGTGACGCTGATGAAGCTGTTGATCTGTGGGTCGAGCTGCTTGATGCGCGCCAGCAGGGCGCCGGTCAGCTCTTCGGAAGaaaa                                                  <  2:119331/149‑1 (MQ=255)
                                                                                                  cTGGTTTCGCCCGCGGCGCGACGGGCATCGGCGGCACGCGCCTGGCCCAGGGCCAGGTCTTCGGTGACGCTGATGAAGCTGTTGATCTGTGGGTCGAGCTGCTTGATGCGCGCCAGCAGGGCGCCGGTCAGCTCTTCGGAAGAAAACGa                                               <  5:25189/149‑1 (MQ=255)
                                                                                                   tGGTTTCGCCCGCGGCGCGACGGGCATCGGCTGCACGCGCCTGGCCCAGGGCCAGGTCTTCGGTGACGCTGATGAAGCTGTTGATCTGTGGGTCGAGCTGCTTGATGCGCGCCAGCAGGGCGCCGGTCAGCTCTTCGGAAGAAAACGAt                                              <  6:196995/149‑1 (MQ=255)
                                                                                                        tCGCCCGCGGCGCGACGGGCATCGGCGGCACGCGCCTGGCCCAGGGCCAGGTCTTCGGTGACGCTGATGAAGCTGTTGATCTGTGGGTCGAGCTGCTTGATGCGCGCCAGCAGGGCGCCGGTCAGCTCTTCGGAAGAAAACGATTTGt                                          <  6:10222/148‑1 (MQ=255)
                                                                                                          gCCCGCGGCGCGACGGGCATCGGCGGCACGCGCCTGGCCCAGGGCCAGGTCTTCGGTGACGCTGATGAAGCTGTTGATCTGTGGGTCGAGCTGCTTGATGCGCGCCAGCAGGGCGCCGGTCAGCTCTTCGGAAGAAAACGATTTGTCgg                                       <  4:199409/149‑1 (MQ=255)
                                                                                                            ccGCGGCGCGACGGGCATCGGCGGCACGCGCCTGGCCCAGGGCCAGGTCTTCGGTGACGCTGATGAAGCTGTTGATCTGTGGGTCGAgctgct                                                                                             >  3:307862/1‑93 (MQ=255)
                                                                                                            ccGCGGCGCGACGGGCATCGGCGGCACGCGCCTGGCCCAGGGCCAGGTCTTCGGTGACGCTGATGAAGCTGTTGATCTGTGGGTCGAgctgct                                                                                             <  4:307862/93‑1 (MQ=255)
                                                                                                                         ggCATCGGCGGCACGCGCCTGGCCCAGGGCCAGGTCTTCGGTGACGCTGATGAAGCTGTTGATCTGTGGGTCGAGCTGCTTGATGCGCGCCAGCAGGGCGCCGGTCAGCTCTTCGGAAGAAAACGATTTGTCGGCGAGTCCGCGGGCGa                        >  7:279280/1‑149 (MQ=255)
                                                                                                                                   ccgaagaccTGGCCCTGGGCCAGGTCTTCGGTGACGCTGATGAAGCTGTTGATCTGTGGGTCGAGCTGCTTGATGCGCGCCAGCAGGGCGCCGGTCAGCTCTTCGGAAGAAAACGATTTGTCGGCGAGTCCGCGGGCGATCTCg                   >  2:250725/8‑144 (MQ=255)
                                                                                                                                   ccgaagaccTGGCCCTGGGCCAGGTCTTCGGTGACGCTGATGAAGCTGTTGATCTGTGGGTCGAGCTGCTTGATGCGCGCCAGCAGGGCGCCGGTCAGCTCTTCGGAAGAAAACGATTTGTCGGCGAGTCCGCGGGCGATCTCg                   <  1:250725/137‑1 (MQ=255)
                                                                                                                                     aCGCGCCTGGCCCAGGGCCAGGTCTTCGGTGACGCTGATGAAGCTGTTGATCTGTGGGTCGAGCTGCTTGATGCGCGCCAGCAGGGCGCCGGTCAGCTCTTCGGAAGAAAACGATTTGTCGGCGAGTCCGCGGGCGATCTCGGCCAggg            >  6:228410/1‑149 (MQ=255)
                                                                                                                                       gcgcCTGGCCCAGGGCCAGGTCTTCGGTGACGCTGATGAAGCTGTTGATCTGTGGGTCGAGCTGCTTGATGCGCGCCAGCAGGGCGCCGGTCAGCTCTTCGGAAGAAAACGATTTGTCGGCGAGTCCGCGGGCGATCTCGGCCAGGGTc          <  4:198351/149‑1 (MQ=255)
                                                                                                                                       gcgcCTGGCCCAGGGCCAGGTCTTCGGTGACGCTGATGAAGCTGTTGATCTGTGGGTCGAGCTGCTTGATGCGCGCCAGCAGGGCGCCGGTCAGCTCTTCGGAAGAAAACGATTTGTCGGCGAGTCCGCGGGCGATCTCGGCCAGGGTc          <  2:339945/149‑1 (MQ=255)
                                                                                                                                         gcCTGGCCCAGGGCCAGGTCTTCGGTGACGCTGATGAAGCTGTTGATCTGTGGGTCGAGCTGCTTGATGCGCGCCAGCAGGGCGCCGGTCAGCTCTTCGGAAGAAAACGATTTGTCGGCGAGTCCGCGGGCGATCTCGGCCAGGGTCaa        <  3:169305/149‑1 (MQ=255)
                                                                                                                                           cTGGCCCAGGGCCAGGTCTTCGGTGACGCTGATGAAGCTGTTGATCTGTGGGTCGAGCTGCTTGATGCGCGCCAGCAGGGCGCCGGTCAGCTCTTCGGAAGAAAACTATTTGTCGGCGAGTCCGCGGGCGATCTCGGCCAGGGGCAAtt      <  5:237117/149‑1 (MQ=255)
                                                                                                                                                ccAGGGCCAGGTCTTCGGTGACGCTGATGAAGCTGTTGATCTGTGGGTCGAGCTGCTTGATGCGCGCCAGCAGGGCGCCGGTCAGCTCTTCGGAAGAAAACGATTTGTCGGCGAGTCCGCGGGCGATCTCGGCCAGGGTCAATTGATg  >  4:187136/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
GCGCCTTGAAGTTGTCGAGCATCTTCGAGCCGCAGCTGGTGCGTACGCCGTTGGTGCAGAACAGGTCTTTGTGGGCGATCGGTGCACCCAGCAGCGCGCTGGTTTCGCCCGCGGCGCGACGGGCATCGGCGGCACGCGCCTGGCCCAGGGCCAGGTCTTCGGTGACGCTGATGAAGCTGTTGATCTGTGGGTCGAGCTGCTTGATGCGCGCCAGCAGGGCGCCGGTCAGCTCTTCGGAAGAAAACGATTTGTCGGCGAGTCCGCGGGCGATCTCGGCCAGGGTCAATTGATG  >  NC_002947/1074561‑1074852

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