Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,990,293 T→A 8.3% L193Q (CTG→CAG) ‡ PP_3516 → AraC family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,990,2930TA 8.3% 168.9 / 5.3 60L193Q (CTG→CAG) ‡PP_3516AraC family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (23/32);  new base A (2/3);  total (25/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.75e-01

CGGCATTGCTTCGACGCCCGTCAGGGCCTGGGGCGCATCTTCGTCGAGCAACTCGACCTGTGTGCCCGCCATTTCGACACCAGCCCCCCGGCAGCCCGGCATTTACTGCTGGAACAGGCCTCCGCCACCCTGCTGCTGGCCTTGCAGCAGGATGAGCGCGTGCTGGGCAGCGAAGGCAGCAACCTCAGCACCCTGCACCTGACGCGGGTCGAGCAGTACGTGCAGCAAAACCTGGGTGACCCCGAGCTGTCGCCACAAACCATCGCCGGCGCCTGTGGG  >  NC_002947/3990157‑3990435
                                                                                                                                        |                                                                                                                                              
cGGCATTGCTTCGACGCCCGTCAGGTCCTGGGGCGCATCTTCGTCGAGCAACTCGACCTGTGTGCCCGCCATTTCGACACCAGCCCCCCGGCAGCCCGGCATTTACTGCTGGAACAGGCCTCCGCCACCCTGCTGCTGGCCTTgcagca                                                                                                                                    <  3:65454/149‑1 (MQ=255)
  gCATTGCTTCGACGCCCGTCAGGGCCTGGGGCGCATCTTCGTCGAGCAACTCGACCTGTGTGCCCGCCATTTCGACACCAGCCCCCCGGCAGCCCGGCATTTACTGCTGGAACAGGCCTCCGCCACCCTGCTGCTGGCCTTGCAGCAgg                                                                                                                                  <  1:339899/149‑1 (MQ=255)
  gCATTGCTTCGACGCCCGTCAGGGCCTGGGGCGCATCTTCGTCGAGCAACTCGACCTGTGTGCCCGCCATTTCGACACCAGCCCCCCGGCAGCCCGGCATTTACTGCTGGAACAGGCCTCCGCCACCCTGCTGCTGGCCTTGCAGCAgg                                                                                                                                  <  2:56046/149‑1 (MQ=255)
        cTTCGACGCCCGTCAGGGCCTGGGGCGCATCTTCGTCGAGCAACTCGACCTGTGTGCCCGCCATTTCGACACCAGCCCCCCGCCAGCCCGGCATTTACTGCTGGAACAGGCCTCCGCCACCCTGCTGCAAGGCCAGCAGCAGGGTggc                                                                                                                             <  7:180872/148‑3 (MQ=255)
         ttCGACGCCCGTCAGGGCCTGGGGCGCATCTTCGTCGAGCAACTCGACCTGTGTGCCCGCCATTTCGACACCAGCCCCCCGGCAGCCCGGCATTTACTGCTGGAACAGGCCTCCGCCACCCTGCTGCTGGCCTTGCAGCAGGATGAgcg                                                                                                                           >  8:287713/1‑149 (MQ=255)
                cccGTCAGGGCCTGGGGCGCATCTTCGTCGAGCAACTCGACCTGTGTGCCCGGCATTTCGACACCAGCCCCCCGGCAGCCCGGCATTTACTGCTGGAACAGGCCTCCGCCACCCTGCTGCTGGCCTTGCAGCAGGATGAGCGCGTGCTg                                                                                                                    <  5:121708/149‑1 (MQ=255)
                 ccGTCAGGGCCTGGGGCGCATCTTCGTCGAGCAACTCGACCTGTGTGCCCGCCATTTCGACACCAGCCCCCCGGCAGCCCGGCATTTACTGCTGGAACAGGCCTCCGCCACCCTGCTGCTGGCCTTGCAGCAGGATGAGCGCGTGCTgg                                                                                                                   >  1:296186/1‑149 (MQ=255)
                  cGTCAGGGCCTGGGGCGCATCTTCGTCGAGCAACTCGACCTGTGTGCCCGCCATTTCGACACCAGCCCCCCGGCAGCCCGGCATTTACTGCTGGAACAGGCCTCCGCCACCCTGCTGCTGGCCTTGCAGCAGGATGAGCGCGTGCTggg                                                                                                                  >  5:117155/1‑149 (MQ=255)
                  cGTCAGGGCCTGGGGCGCATCTTCGTCGAGCAACTCGACCTGTGTGCCCGCCATTTCGACACCAGCCCCCCGGCAGCCCGGCATTTACTGCTGGAACAGGCCTCCGCCACCCTGCTGCAAGGCCAGCAGCAgg                                                                                                                                  >  5:98628/1‑133 (MQ=255)
                  cGTCAGGGCCTGGGGCGCATCTTCGTCGAGCAACTCGACCTGTGTGCCCGCCATTTCGACACCAGCCCCCCGGCAGCCCGGCATTTACTGCTGGAACAGGCCTCCGCCACCCTGCTGCAAGGCCAGCAGCAgg                                                                                                                                  <  6:98628/133‑1 (MQ=255)
                        ggCCTGGGGCGCATCTTCGTCGAGCAACTCGACCTGTGTGCCCGCCATTTCGACACCAGCCCCCCGGCAGCCCGGCATTTACTGCTGGAACAGGCCTCCGCCACCCTGCTGCTGGCCTTGCAGCAGGATGAGCGCGTGCTGGGCAGCg                                                                                                             >  8:351679/1‑148 (MQ=255)
                         gCCTGGGGCGCATCTTCGTAGAGCAACTCGACCTGTGAGCCCTCCATTTCGACACCAGCCCCCCGGCAGCCCGGCATTTACTGCTGGAACAGGCCTCCGCCACCCTGCTTCTGGCCTTGCAGCATGATGAGCTCGTGCTGGGCAGCGaa                                                                                                           <  3:233344/149‑1 (MQ=255)
                           cTGGGGCGCATCTTCGTCGAGCAACTCGACCTGTGTGCCCGCCATTTCGACACCAGCCCCCCGTCAGCCCGGCATTTACTGCTGGAACAGGCCTCCGCCACCCTGCTGCTGGCCTTGCAGCAGGATGAGCGCGTGCTGGGCAGCGAAgg                                                                                                         <  8:63386/149‑1 (MQ=255)
                                  gcATCTTCGTCGAGCAACTCGACCTGTGTGCCCGCCATTTCGACACCAGCCCCCCGGCAGCCCGGCATTTACTGCTGGAACAGGCCTCCGCCACCCTGCTGCTGGCCTTGCAGCAGGATGAGCGCGTGCTGGGCAGCGAAGGCAGCAAc                                                                                                  <  4:356013/149‑1 (MQ=255)
                                             gaacAACTCGACCTGTGTGCCCGCCATTTCGACACCAGCCCCCCGGCAGCCCGGCATTTACTGCTGGAACAGGCCTCCGCCACCCTGCTGCTGGCCTTGCAGCAGGATGAGCGCGTGCTGGGCAGCGAAGGCAGCAACCTCAGCACCCt                                                                                       <  7:215491/146‑1 (MQ=255)
                                              aGCAACTCGACCTGTGTGCCCGCCATTTCGACACCAGCCCCCCGGCAGCCCGGCATTTACTGCTGGAACAGGCCTCCGCCACCCTGCTGCTGGCCTTGCAGCAGGATGAGCGCGTGCTGGGCAGCGAAGGCAGCAACCTCAGCACCCTg                                                                                      >  6:195449/1‑149 (MQ=255)
                                               gCAACTCGACCTGTGTGCCCGCCATTTCGACACCAGCCCCCCGGCAGCCCGGCATTTACTGCTGGAACAGGCCTCCGCCACCCTGCTGCTGGCCTTGCAGCAGGATGAGCGCGTGCTGGGCAGCGAAGGCAGCAACCTCAGCACCCTg                                                                                      >  4:336039/1‑148 (MQ=255)
                                               gCAACTCGACCTGTGTGCCCGCCATTTCGACACCAGCCCCCCGGCAGCCCGGCATTTACTGCTGGAACAGGCCTCCGCCACCCTGCTGCTGGCCTTGCAGCAGGATGAGCGCGTGCTGGGCAGCGAAGGCAGCAACCTCAGCACCCTGc                                                                                     >  1:74603/1‑149 (MQ=255)
                                                  aCTCGACCTGTGTGCCCGCCATTTCGACACCAGCCCCCCGGCAGCCCGGCATTTACTGCTGGAACAGGCCTCCGCCACCCTGCTGCTGGCCTTGCAGCAGGATGAGCGCGTGCTGGGCAGCGAAGGCAGCAACCTCAGCACCCTGCAcc                                                                                  >  8:177890/1‑149 (MQ=255)
                                                   cTCGACCTGTGTGCCCGCCATTTCGACACCAGCCCCCCGGCAGCCCGGCATTTACTGCTGGAACAGGCCTCCGCCACCCTGCTGCAAGGCCAGCAGCAgggt                                                                                                                                <  2:320525/102‑3 (MQ=255)
                                                   cTCGACCTGTGTGCCCGCCATTTCGACACCAGCCCCCCGGCAGCCCGGCATTTACTGCTGGAACAGGCCTCCGCCACCCTGCTGCAAGGCCAGCAGCAgggt                                                                                                                                >  1:320525/1‑100 (MQ=255)
                                                       aCCTGTGTGCCCGCCATTTCGACACCAGCCCCCCGGCAGCCCGGCATTTACTGCTGGAACAGGCCTCCGCCACCCTGCTGCTGGCCTTGCAGCAGGATGAGCGCGTGCTGGGCAGCGAAGGCAGCAACCTCAGCACCCTGCACCTGAcg                                                                             >  3:15924/1‑149 (MQ=255)
                                                        ccTGTGTGCCCGCCATTTCGACACCAGCCCCCCGGCAGCCCGGCATTTACTGCTGGAACAGGCCTCCGCCACCCTGCTGCTGGCCTTGCAGCAGGATGAGCGCGTGCTGGGCAGCGAAGGCAGCAACCTCAGCACCCTGCACCTGAcgc                                                                            <  2:74603/149‑1 (MQ=255)
                                                            tgtgCCCGCCATTTGGACACCAGCCCCCCGGCAGCCCGGCATTTACTGCTGGAACAGGCCTCCGCCACCCTGCTGCTGGCCTTGCAGCAGGATGAGCGCGTGCTGGGCAGCGAAGGCAGCAACCTCAGCACCCTGCACCTGACGCgg                                                                          <  5:362554/147‑1 (MQ=255)
                                                            tgtgCCCGCCATTTCGACACCAGCCCCCCGGCAGCCCGGCATTTACTGCTGGAACAGGCCTCCGCCACCCTGCTGCTGGCCTTGCAGCAGGATGAGCGCGTGCTGGGCAGCGAAGGCAGCAACCTCAGCACCCTGCACCTGACGCgg                                                                          <  6:117155/147‑1 (MQ=255)
                                                               gCCCGCCATTTCGACACCAGCCCCCCGGCAGCCCGGCATTTACTGCTGGAACAGGCCTCCGCCACCCTGCTGCTGGCCTTGCAGCAGGATGAGCGCGTGCTGGGCAGCGAAGGCAGCAACCTCAGCACCCTGCACCTGACGCGGGTCGa                                                                     <  5:41127/149‑1 (MQ=255)
                                                                cccGCCATTTCGACACCAGCCCCCCGGCAGCCCGGCATTTACTGCTGGAACAGGCCTCCGCCACCCTGCTGCTGGCCTTGCAGCAGGATGAGCGCGTGCTGGGCAGCGAAGGCAGCAACCTCAGCACCCTGCAGCTGACGCGGGTCGa                                                                     <  8:88249/148‑1 (MQ=255)
                                                                cccGCCATTTCGACACCAGCCCCCCGGCAGCCCGGCATTTACTGCTGGAACAGGCCTCCGCCACCCTGCTGCTGGCCTTGCAGCAGGATGAGCGCGTGCTGGGCAGCGAAGGCAGCAACCTCAGCACCCTGCACCTGACGCGGGTCGAg                                                                    <  7:287713/149‑1 (MQ=255)
                                                                cccGCCAGTTCGACACCAGCCCCCCGGCAGCCCGGCATATACTGCTGGAACAGGCCTCCGCCACCCTGCTGCTGGCCTTGCAGCAGGATGAGCGCGTGCTGGGCAGCGAAGGCAGCAACCTCAGCACCCTGCACCTGACGCGGGTCGAg                                                                    <  4:267526/149‑1 (MQ=255)
                                                                 ccGCCATTTCGACACCAGCCCCCCGGCAGCCCGGCATTTACTGCTGGAACAGGCCTCCGCCACCCTGCTGCTGGCCTTGCAGCAGGATGAGCGCGTGCTGGGCAGCGAAGGCAGCAACCTCAGCACCCTGCACCTGACGCGGGTCGAGc                                                                   <  4:170260/149‑1 (MQ=255)
                                                                  cGCCATTTCGACACCAGCCCCCCGGCAGCCCGGCATTTACTGCTGGAACAGGCCTCCGCCACCCTGCTGCTGGCCTTGCAGCAGGATGAGCGCGTGCTGGGCAGCGAAGGCAGCAACCTCAGCACCCTGCACCTGACGCGGGTCgcg                                                                    >  6:340407/1‑145 (MQ=255)
                                                                   gCCATTTCGACACCAGCCCCCCGGCAGCCCGGCATTTACTGCTGGAACAGGCCTCCGCCACCCTGCTGCTGGCCTTGCAGCAGGATGAGCGCGTGCTGGGCAGCGAAGGCAGCAACCTCAGCACCCTGCACCTGACGCGGGTCGAGTAg                                                                 >  5:84077/1‑149 (MQ=255)
                                                                    ccATTTCGAGACCAGCCCCCCGGCAGCCCGGCATTTACTGCTGGAACAGGCCTCCGCCACCCTGCTGGTGGCCTTGCAGCAGGATGAGCGCGTGCTGGGCAGCGAAGGCAGCAACCTCAGCACCCTGCACCTGACGCGGGTCGAGCAGt                                                                <  3:317236/149‑1 (MQ=255)
                                                                           gACACCAGCCCCCCGGCAGCCCGGCATTTACTGCTGGAACAGGCCTCCGCCACCCTGCTGCTGGCCTTGCAGCAGGATGAGCGCGTGCTGGGCAGCGAAGGCAGCAACCTCAGCACCCTGCACCTGACGCGGGTCGAGCAGTACGTgc                                                          <  5:195449/148‑1 (MQ=255)
                                                                             cacCAGCCCCCCGGCAGCCCGGCATTTACTGCTGGAACAGGCCTCCGCCACCCTGCTGCTGGCCTTGCAGCAGGATGAGCGCGTGCTGGGCAGCGAAGGCAGCAACCTCa                                                                                              >  8:293407/1‑110 (MQ=255)
                                                                             cacCAGCCCCCCGGCAGCCCGGCATTTACTGCTGGAACAGGCCTCCGCCACCCTGCTGCTGGCCTTGCAGCAGGATGAGCGCGTGCTGGGCAGCGAAGGCAGCAACCTCAGCACCCTGCACCTGACGCGGGTCGAGCAGTACGTgcag                                                        >  3:24280/1‑148 (MQ=255)
                                                                                cAGCCCCCCGGCAGCCCGGCATTTACTGCTGGAACAGGCCTCCGCCACCCTGCTGCTGGCCTTGCAGCAGGATGAGCGCGTGCTGGGCAGCGAAGGCAGCAACCTCAGCAcc                                                                                         >  8:304524/1‑112 (MQ=255)
                                                                                cAGCCCCCCGGCAGCCCGGCATTTACTGCTGGAACAGGCCTCCGCCACCCTGCTGCTGGCCTTGCAGCAGGATGAGCGCGTGCTGGGCAGCGAAGGCAGCAACCTCAGCAcc                                                                                         <  7:304524/112‑1 (MQ=255)
                                                                                  gCCCCCCGGCAGCCCGGCATTTACTGCTGGAACAGGCCTCCGCCACCCTGCTGCTGGCCTTGCAGCAGGATGAGCGCGTGCTGGGCAGCGAAGGCAGCAACCTCAGCACCCTGCACCTGACGCGGGTCGAGCAGTACGTGCAGCAAAAc                                                  >  4:139524/1‑149 (MQ=255)
                                                                                     ccccGGCAGCCCGGCATTTACTGCTGGAACAGGCCTCCGCCACCCTGCTGCTGGCCTTGCAGCAGGATGAGCGCGTGCTGGGCAGCGAAGGCAGCAACCTCAGCACCCTGCACCTGACGCGGGTCGAGCAGTACGTGCAGCAAAACCTg                                               <  3:139524/149‑1 (MQ=255)
                                                                                      cccGGCAGCCCGGCATTTACTGCTGGAACAGGCCTCCGCCACCCTGCTGCTGGCCTTGCAGCAGGATGAGCGCGTGCTGGGCAGCGAAGGCAGCAACCTCAGCACCCTGCACCTGACGCGGGTCGAGTAGTACGTGCAGCAAAACCTgg                                              <  6:84077/149‑1 (MQ=255)
                                                                                      cccGGCAGCCCGGCATTTACTGCTGGAACAGGCCTCCGCCACCCTGCTGCTGGCCTTGCAGCAGGATGAGCGCGTGCTGGGCAGCGAAGGCAGCAACCTCAGCACCCTGCACCTGACGCGGGTCGAGCAGTACGTGCAGCAAAACCTgg                                              <  8:208818/149‑1 (MQ=255)
                                                                                         ggCAGCCCGGCATTTACTGCTGGAACAGGCCTCCGCCACCCTGCTGCTGGCCTTGCAGCAGGATGAGCGCGTGCTGGGCAGCGAAGGCAGCAACCTCAGCACCCTGCACCTGACGCGGGTCGAGCAGTACGTGAAGCAAAACCTGGGt                                            >  5:308654/1‑148 (MQ=255)
                                                                                                cGGCATTTACTGCTGGAACAGGCCTCCGCCACCCTGCTGCTGGCCTTGCAGCAGGATGAGCGCGTGCTGGGCAGCGAAGGCAGCAACCTCAGCACCCTGCACCTGACGCGGGTCGAGCAGTACGTGCAGCAAAACCTGGGTGAcccc                                      >  5:22389/1‑147 (MQ=255)
                                                                                                cGGCATTTACTGCTGGAACAGGCCTCCGCCACCCTGCTGCTGGCCTTGCAGCAGGATGAGCGCGTGCTGGGCAGCGAAGGCAGCAACCTCAGCACCCTGCACCTGACGCGGGTCGAGCAGTACGTGCAGCAAAACCTGGGTGACCCCGa                                    <  7:351679/149‑1 (MQ=255)
                                                                                                       tACTGCTGGAACAGGCCTCCGCCACCCTGCTGCTGGCCTTGCAGCAGGATGAGCGCGTGCTGGGCAGCGAAGGCAGCAACCTCAGCACCCTGCACCTGACGCGGGTCGAGCAGTACGTGCAGCAAAACCTGGGTGACCCCGAGCTGTc                              <  4:15924/148‑1 (MQ=255)
                                                                                                        aCTGCTGGAACAGGCCTCCGCCACCCTGCTGCTGGCCTTGCAGCAGGATGAGCGCGTGCTGGGCAGCGAAGGCAGCAACCTCAGCACCCTGCACCTGACGCGGGTCGAGCAGTACGTGCAGCAAAACCTGGGTGACCCCGAGCTGTCg                             <  4:24280/148‑1 (MQ=255)
                                                                                                         ctgctgGAACAGGCCTCCGCCACCCTGCTGCTGGCCTTGCAGCAGGATGAGCGCGTGCTGGGCAGCGAAGGCAGCAACCTCAGCACCCTGCACCTGACGCGGGTCGAGCAGTACGTGCAGCAAAACCTGGGTGACCCCGAGCTGTCGc                            <  1:159199/148‑1 (MQ=255)
                                                                                                            ctgGAACAGGCCTCCGCCACCCTGCTGCTGGCCTTGCAGCAGGATGAGCGCGTGCTGGGCAGCGAAGGCAGCAACCTCAGCACCCTGCACCTGACGCGGGTCGAGCAGTAc                                                              >  2:75076/1‑111 (MQ=255)
                                                                                                            ctgGAACAGGCCTCCGCCACCCTGCTGCTGGCCTTGCAGCAGGATGAGCGCGTGCTGGGCAGCGAAGGCAGCAACCTCAGCACCCTGCACCTGACGCGGGTCGAGCAGTAc                                                              <  1:75076/111‑1 (MQ=255)
                                                                                                                    ggCCTCCGCCACCCTGCTGCTGGCCTTGCAGCAGGATGAGCGCGTGCTGGGCAGCGAAGGCAGCAACCTCAGCACCCTGCACCTGACGCGGGTCGAGCAGTACGTGCAGCAAAACCTGGGTGACCCCGAGCTGTCGCCACAAACCATCg                <  7:177890/149‑1 (MQ=255)
                                                                                                                     gCCTCCGCCACCCTGCTGCTGGCCTTGCAGCAGGATGAGCGCGTGCTGGGCAGCGAAGGCAGCAACCTCAGCACCCTGCACCTGACGCGGGTCGAGCAGTACGTGCAGCAAAACCTGGGTGa                                          >  8:206318/1‑122 (MQ=255)
                                                                                                                      ccTCCGCCACCCTGCTGCTGGCCTTGCAGCAGGATGAGCGCGTGCTGGGCAGCGAAGGCAGCAACCTCAGCACCCTGCACCTGACGCGGGTCGAGCAGTACGTGCAGCAAAACCTGGGTGACCCCGAg                                   >  2:81536/1‑128 (MQ=255)
                                                                                                                      ccTCCGCCACCCTGCTGCTGGCCTTGCAGCAGGATGAGCGCGTGCTGGGCAGCGAAGGCAGCAACCTCAGCACCCTGCACCTGACGCGGGTCGAGCAGTACGTGCAGCAAAACCTGGGTGACCCCGAg                                   <  1:81536/128‑1 (MQ=255)
                                                                                                                           gCCACCCTGCTGCTGGCCTTGCAGCAGGATGAGCGCGTGCTGGGCAGCGAAGGCAGCAACCTCAGCACCCTGCACCTGACGCGGGTCGAGCAGTACGTGCAGCAAAACCTGGGTGACCCCGAGCTGTCGCCACAAACCATCGCCgg            <  6:247013/146‑1 (MQ=255)
                                                                                                                           gCCACCCTGCTGCTGGCCTTGCAGCAGGATGAGCGCGTGCTGGGCAGCGAAGGCAGCAACCTCAGCACCCTGCACCTGACGCGGGTCGAGCAGTACGTGCAGCAAAACCTGGGTGACCCCGAGCTGTCGCCACAAACCATCGCCgg            >  5:247013/1‑146 (MQ=255)
                                                                                                                              aCCCTGCTGCTGGCCTTGCAGCAGGATGAGCGCGTGCTGGGCAGCGAAGGCAGCAACCTCAGCACCCTGCACCTGACGCGGGTCGAGCAGTACGTGCAGCAAAACCTGGGTGACCCCGAGCTGTCGCCACAAACCATCGCCGGCGCCtg      >  6:8183/1‑149 (MQ=255)
                                                                                                                               cccTGCTGCTGGCCTTGCAGCAGGATGAGCGCGTGCTGGGCAGCGAAGGCAGCAACCTCAGCACCCTGCACCTGACGCGGGTCGAGCAGTACGTGCAGCAAAACCTGGGTGACCCCGAGCTGTCGCCACAAACCATCGCCGGCGCCtgt     <  3:336039/149‑1 (MQ=255)
                                                                                                                                ccTGCTGCTGGCCTTGCAGCAGGATGAGCGCGTGCTTGGCAGCGAAGGCAGCAAACTCAGCACCCTGCACCTGACGCGGGTCGAGCAGTACGTGCAGCAAAACCTGGGTGACCCCGAGCTGTCGCCACAAACCATCGCCGGCGCCtgtg    <  5:8183/149‑1 (MQ=255)
                                                                                                                                ccTGCTGCTGGCCTTGCAGCAGGATGAGCGCGTGCTGGGCAGCGAAGGCAGCAACCTCAGCACCCTGAACCTGACGCGGGTCGAGCAGTACGTCCAGCAAAACCTGGGTGACCCCGAGCTGTCGCCACAAACCAACGCCGGCGCCtcag    >  3:337071/1‑146 (MQ=255)
                                                                                                                                  tgctgctgGCCTTGCAGCAGGATGAGCGCGTGCTGGGCAGCGAAGGCAGCAACCTCAGCACCCTGCACCTGACGCGGGTCGAGCAGTACGTGCAGCAAAACCTGGGTGACCCCGAGCTGTCGCCACAAACCATCGCCGGCGCCTGTggg  >  7:10335/1‑149 (MQ=255)
                                                                                                                                    ctgctgGGCTTGCAGCAGGATGAGCGCGGGCTGGGCAGCGTAGGCAGCAACCTCAGCACCCTGCAGCTGACGCGGGTCGAGCAGTACGTGCAGCAATACCTGGGTGACCCCGAGCTGTCGCCACAAACCATCGCCGGCGCCTGTggg  <  5:340407/147‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                              
CGGCATTGCTTCGACGCCCGTCAGGGCCTGGGGCGCATCTTCGTCGAGCAACTCGACCTGTGTGCCCGCCATTTCGACACCAGCCCCCCGGCAGCCCGGCATTTACTGCTGGAACAGGCCTCCGCCACCCTGCTGCTGGCCTTGCAGCAGGATGAGCGCGTGCTGGGCAGCGAAGGCAGCAACCTCAGCACCCTGCACCTGACGCGGGTCGAGCAGTACGTGCAGCAAAACCTGGGTGACCCCGAGCTGTCGCCACAAACCATCGCCGGCGCCTGTGGG  >  NC_002947/3990157‑3990435

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