Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,653,195 C→A 8.2% Q361K (CAG→AAG) ‡ PP_3218 → NtaA/SnaA/SoxA family monooxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,653,1950CA 8.2% 185.2 / 5.2 61Q361K (CAG→AAG) ‡PP_3218NtaA/SnaA/SoxA family monooxygenase
Reads supporting (aligned to +/- strand):  ref base C (32/24);  new base A (2/3);  total (34/27)
Fisher's exact test for biased strand distribution p-value = 6.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.65e-01

ACCTGGATGCCCCGTTCCCGCAAGAGAGCCTGGAGTTTGCACGGAGCAGCTTTTTCACGAATGCCAAGCGTATTGCCGACCAGGCGCGGGACAAGGCCTTCAG‑TCTGCGT‑‑CA‑GGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGATGGGTTCAATATCTACATTGGTCATCCTGGCCAGTTCCGGCGGTTTACCGAAGAAGTGGTACC  >  NC_002947/3653064‑3653340
                                                                                                                                       |                                                                                                                                                 
aCCTGGATGCCCCGTTCCCGCAAGAGAGCCTGGAGTTTGCACGGAGCAGCTTTTTCACGAATGCCAAGCGTATTGCCGACCAGGCGCGGGACAAGGCCTTCAG‑TCTGCGT‑‑CA‑GGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGcc                                                                                                                                  >  7:273311/1‑149 (MQ=255)
 ccTGGATGCCCCGTTCCCGCAAGAGAGCCTGGAGTTTGCACGGAGCAGCTTTTTCACGAATGCCAAGCGTATTGCCGACCAGGCGCGGGACAAGGCCTTCAG‑TCTGCGT‑‑CA‑GGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCt                                                                                                                                 <  8:149695/149‑1 (MQ=255)
 ccTGGATGCCCCGTTCCCGCAAGAGAGCCTGGAGTTTGCACGGAGCAGCTTTTTCACGAATGCCAAGCGTATTGCCGACCAGGCGCGGGACAAGGCCTTCAG‑TCTGCGT‑‑CA‑GGCTGTTGAGTTTGGCCGTAACTGACGGCCaaactcaa                                                                                                                                 <  8:287047/149‑8 (MQ=255)
    ggATGCCCCGTTCCCGCAAGAGAGCCTGGAGTTTGCACGGAGCAGCTTTTTCACGAATGCCAAGCGTATTGCCGACCAGGCGCGGGACAAGGCCTTCAG‑TCTGCGT‑‑CA‑GGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCg                                                                                                                              >  2:237382/1‑149 (MQ=255)
     gATGCCCCGTTCCCGCAAGAGAGCCTGGAGTTTGCACGGAGCAGCTTTTTCACGAATGCCAAGCGTATTGCCGACCAGGCGCGGGACAAGGCCTTCAG‑TCTGCGT‑‑CA‑GGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGt                                                                                                                             <  3:116447/149‑1 (MQ=255)
         ccccGTTCCCGCAAGAGAGCCTGGAGTTTGCACGGAGCAGCTTTTTCACGAATGCCAAGCGTATTGCCGACCAGGCGCGGGACAAGGCCTTCAG‑TCTGCGT‑‑CA‑GGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTgg                                                                                                                          <  4:339232/148‑1 (MQ=255)
          cccGTTCCCGCAAGAGAGCCTGGAGTTTGCACGGAGCAGCTTTTTCACGAATGCCAAGCGTATTGCCGACCAGGCGCGGGACAAGGCCTTCAG‑TCTGCGT‑‑CA‑GGCTGTTGAGTTTGGCCGTCAGTTAc                                                                                                                                             >  7:267376/1‑128 (MQ=255)
          cccGTTCCCGCAAGAGAGCCTGGAGTTTGCACGGAGCAGCTTTTTCACGAATGCCAAGCGTATTGCCGACCAGGCGCGGGACAAGGCCTTCAG‑TCTGCGT‑‑CA‑GGCTGTTGAGTTTGGCCGTCAGTTAc                                                                                                                                             <  8:267376/128‑1 (MQ=255)
          cccGTTCCCGCAAGAGAGCCTGGAGTTTGCACGGAGCAGCTTTTTCACGAATGCCAAGCGTATTGCCGACCAGGCGCGGGACAAGGCCTTCAG‑TCTGCGT‑‑CA‑GGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGtt                                                                                                                        >  5:123034/1‑149 (MQ=255)
              ttCCCCCAAGACCGCCTGGATGTTGCACGGATCAGCGTTTTCACCTATTCCAACCGTTATGGCGACCAGGCGCGGGACAAGGCCTTCAG‑TCTGCGT‑‑CA‑GGCTGTTGAGTTTGGCCGTAAGTTACGGCCAGGGGCCTTCGTTGGTTCGGc                                                                                                                    <  2:329584/149‑1 (MQ=255)
                cccGCAAGAGAGCCTGGAGTTTGCACGGAGCAGCTTTTTCACGAATGCCAAGCGTATTGCCGACCAGGCGCGGGACAAGGCCTTCAG‑TCTGCGT‑‑CA‑GGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGcc                                                                                                                   >  8:304998/1‑148 (MQ=255)
                    cAAGAGAGCCTGGAGTTTGCACGGAGCAGCTTTTTCACGAATGCCAAGCGTATTGCCGACCAGGCGCGGGACAAGGCCTTCAG‑TCTGCGT‑‑CA‑GGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAAc                                                                                                              >  4:214241/1‑149 (MQ=255)
                             cTGGAGTTTGCACGGAGCAGCTTTTTCACGAATGCCAAGCGTATTGCCGACCAGGCGCGGGACAAGGCCTTCAG‑TCTGCGT‑‑CA‑GGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTggcggc                                                                                                     >  6:275670/1‑149 (MQ=255)
                                gAGTTTGCACGGAGCAGCTTTTTCACGAATGCCAAGCGTATTGCCGACCAGGCGCGGGACAAGGCCTTCAG‑TCTGCGT‑‑CA‑GGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTg                                                                                                                           >  1:191807/1‑124 (MQ=255)
                                gAGTTTGCACGGAGCAGCTTTTTCACGAATGCCAAGCGTATTGCCGACCAGGCGCGGGACAAGGCCTTCAG‑TCTGCGT‑‑CA‑GGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTg                                                                                                                           <  2:191807/124‑1 (MQ=255)
                                gAGTTTGCACGGAGCAGCTTTTTCACGAATGCCAAGCGTATTGCCGACCAGGCGCGGGACAAGGCCTTCAG‑TCTGCGT‑‑CA‑GGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGt                                                                                                                         >  3:221728/1‑126 (MQ=255)
                                gAGTTTGCACGGAGCAGCTTTTTCACGAATGCCAAGCGTATTGCCGACCAGGCGCGGGACAAGGCCTTCAG‑TCTGCGT‑‑CA‑GGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGt                                                                                                                         <  4:221728/126‑1 (MQ=255)
                                       cACGGAGCAGCTTTTTCACGAATGCCAAGCGTATTGCCGACCAGGCGCGGGACAAGGCCTTCAG‑TCTGCGT‑‑CA‑GGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGc                                                                                            <  8:273311/148‑1 (MQ=255)
                                         cGGAGCAGCTTTTTCACGAATGCCAAGCGTATTGCCGACCAGGCGCGGGACAAGGCCTTCAG‑TCTGCGT‑‑CA‑GGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGa                                                                                         <  3:214241/149‑1 (MQ=255)
                                                gcTTTTTCACGAATGCCAAGCGTATTGCCGACCAGGCGCGGGACAAGGCCTTCAGCTCTGTCTACCATTGCTTTTGTGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTg                                                                                      <  1:115210/149‑1 (MQ=255)
                                                      tCACGAATGCCAAGCGTATTGCCGACCAGGCGCGGGACAAGGCCTTCAG‑TCTGCGT‑‑CA‑GGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAgcg                                                                            >  6:215837/1‑149 (MQ=255)
                                                      tCACGAATGCCAAGCGTATTGCCGACCAGGCGCGGGACAAGGCCTTCAG‑TCTGCGT‑‑CA‑GGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAgcg                                                                            >  5:62316/1‑149 (MQ=255)
                                                       cACGAATGCCAAGCGTATTGCCGACCAGGCGCGGGACAAGGCCTTCAG‑TCTGCGT‑‑CA‑GGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAgcg                                                                            <  2:72874/148‑1 (MQ=255)
                                                        aCGAATGCCAAGCGTATTGCCGACCAGGCGCGGGACAAGGCCTTCAG‑TCTGCGT‑‑CA‑GGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGcgcg                                                                          >  5:343707/1‑149 (MQ=255)
                                                        aCGAATGCCAAGCGTATTGCCGACCAGGCGCGGGACAAGGCCTTCAG‑TCTGCGT‑‑CA‑GGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGcgcg                                                                          >  2:24010/1‑149 (MQ=255)
                                                              gCCAAGCGTATTGCCGACCAGGCGCGGGACAAGGCCTTCAG‑TCTGCGT‑‑CA‑GGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGa                                                                                         >  3:112859/1‑128 (MQ=255)
                                                              gCCAAGCGTATTGCCGACCAGGCGCGGGACAAGGCCTTCAG‑TCTGCGT‑‑CA‑GGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGa                                                                                         <  4:112859/128‑1 (MQ=255)
                                                                 aaGCGAATTGCCGACCAGGCGCGGGACAAGGCCGTCAG‑TCTGCGT‑‑CA‑GGCAGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGAGGCGGCGAAGATGGCCGAGTGGATCGAAGCGCGGGCGCTgca                                                                 >  2:34475/1‑147 (MQ=255)
                                                                       aTTGCCGACCAGGCGCGGGACAAGGCCTTCAG‑TCTGCGT‑‑CA‑GGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTgg                                                                  >  6:137411/1‑142 (MQ=255)
                                                                       aTTGCCGACCAGGCGCGGGACAAGGCCTTCAG‑TCTGCGT‑‑CA‑GGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTgg                                                                  <  5:137411/142‑1 (MQ=255)
                                                                           ccGACCAGGCGCGGGACAAGGCCTTCAG‑TCTGCGT‑‑CA‑GGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGATGGGTTCaa                                                        <  6:123034/148‑1 (MQ=255)
                                                                            cGACCAGGCGCGGGACAAGGCCTTCAG‑TCTGCGT‑‑CA‑GGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGATGGGTTCAata                                                      >  4:29542/1‑149 (MQ=255)
                                                                              aCCAGGCGCGGGACAAGGCCTTCAG‑TCTGCGT‑‑CA‑GGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGATGGGTTCAATATc                                                    >  7:5679/1‑149 (MQ=255)
                                                                                cAGGCGCGGGACAAGGCCTTCAG‑TCTGCGT‑‑CA‑GGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGATGGGTTCAATATCTa                                                  >  1:287746/1‑149 (MQ=255)
                                                                                cAGGCGCGGGACAAGGCCTTCAG‑TCTGCGT‑‑CA‑GGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGATGGGTTCAATATCTa                                                  >  8:231940/1‑149 (MQ=255)
                                                                                  ggCGCGGGACAAGGCCTTCAG‑TCTGCGT‑‑CA‑GGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGAt                                                                <  8:182819/133‑1 (MQ=255)
                                                                                  ggCGCGGGACAAGGCCTTCAG‑TCTGCGT‑‑CA‑GGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGAt                                                                >  7:182819/1‑133 (MQ=255)
                                                                                   gcgcGGGACAAGGCCTTCAG‑TCTGCGT‑‑CA‑GGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGATGGGGTCAATATCTACAt                                               <  5:215837/149‑1 (MQ=255)
                                                                                    cgcgGGACAAGGCCTTCAG‑TCTGCGT‑‑CA‑GGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGATGGGCTCAATATCTACAtt                                              <  7:304998/149‑1 (MQ=255)
                                                                                           cAAGGCCTTCAG‑TCTGCGT‑‑CA‑GGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGATGGGTTCAATATCTACATTGGTCATc                                       <  5:275670/149‑1 (MQ=255)
                                                                                              ggCCTTCAG‑TCTGCGT‑‑CA‑GGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGt                                                                                                           <  1:123129/78‑1 (MQ=255)
                                                                                              ggCCTTCAG‑TCTGCGT‑‑CA‑GGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGt                                                                                                           >  2:123129/1‑78 (MQ=255)
                                                                                                  ttCAG‑TCTGCGT‑‑CA‑GGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGATGGGTTCAATATCTACATTGGTCATCCTGGc                                  >  4:168208/1‑147 (MQ=255)
                                                                                                    cAG‑TCTGCGT‑‑CA‑GGCTGTTGAGTTTGGCCGTCAGTTATGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGATGGGTTCAATATCTACATTGGTCATCCTGGCCAGt                              >  4:219389/1‑149 (MQ=255)
                                                                                                        tCTGCGT‑‑CA‑GGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAAGCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGATGGGtt                                                           <  6:44675/117‑1 (MQ=255)
                                                                                                        tCTGCGT‑‑CA‑GGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGATGGGtt                                                           >  5:44675/1‑117 (MQ=255)
                                                                                                          tGCGT‑‑CA‑GGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGATGGGTTCAATATCTACATTGGTCATCCTGGCCAGTTCcgg                         <  8:5679/149‑1 (MQ=255)
                                                                                                            cGT‑‑CA‑GGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGATGGGTTCAATATCTACATTGGTCATCCTGGCCAGTTCcgg                         <  2:287746/147‑1 (MQ=255)
                                                                                                             gT‑‑CA‑GGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGATGGGTTCAATATCTACATTGGTCATCCTGGCCAGTTCcggcgg                      >  4:28308/1‑149 (MQ=255)
                                                                                                                       tGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGATGGGTTCAATATCTACATTGGTCATCCTGGCCAGTTCCGGCGGTTTAcc                >  1:52700/1‑148 (MQ=255)
                                                                                                                         ttGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGATGGGTTCAATATCTACATTGGTCATCCTGGCCAGTTCCGGCGGTTTAc                 >  5:13812/1‑145 (MQ=255)
                                                                                                                         ttGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGATGGGTTCAATATCTACATTGGTCATCCTGGCCAGTTCCGGCGGTTTAc                 <  6:13812/145‑1 (MQ=255)
                                                                                                                         ttGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGATGGGTTCAATATCTACATTGGTCATCCTGGCCAGTTCCGGCGGTTTACCg               >  2:340169/1‑147 (MQ=255)
                                                                                                                             ccctGGCCGTAACTGACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGATGGGTTCAATATCTTCATTGGTCATCCTGGCCAGTTCcgg                         >  5:308866/4‑133 (MQ=255)
                                                                                                                             ccctGGCCGTAACTGACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGATGGGTTCAATATCTACATTGGTCATCCTGGCCAGTTCcgg                         <  6:308866/130‑1 (MQ=255)
                                                                                                                             ccctGGCCGTAACTGACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGATGGGTTCAATATCTACATTGGTCATCCTGGCCAGTTCCGGCGGTTTACCGaagaag         >  5:209824/4‑149 (MQ=255)
                                                                                                                              tttGGCTGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGATGGGTTCAATATCTACATTGGTCATCCTGGCCAGTTCCGGCGGTTTACCGAAGAAGt        >  8:305150/1‑149 (MQ=255)
                                                                                                                                 ggCCGTCAGTTACGGCCA‑GGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGATGGGTTCAATATCTACATTGGTCATCCTGGCCAGTTCCGGCGGTTTACCGAAGAAGTGGTa    >  7:33796/1‑149 (MQ=255)
                                                                                                                                    cGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGATGGGTTCAATATCTACATTGGTCATCCTGGCCAGTTCCGGCGGTTTACCGAAGAAGTGGTAcc  <  6:343707/149‑1 (MQ=255)
                                                                                                                                    cGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGATGGGTTCAATATCCACATAGGTCATCCTGGCCAGTTCCGGCGGTTTACCGAAGAAGt        >  8:117516/1‑143 (MQ=255)
                                                                                                                                    cGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGATGGGTTCAATATCCACATAGGTCATCCTGGCCAGTTCCGGCGGTTTACCGAAGAAGt        <  7:117516/143‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                                 
ACCTGGATGCCCCGTTCCCGCAAGAGAGCCTGGAGTTTGCACGGAGCAGCTTTTTCACGAATGCCAAGCGTATTGCCGACCAGGCGCGGGACAAGGCCTTCAG‑TCTGCGT‑‑CA‑GGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGATGGGTTCAATATCTACATTGGTCATCCTGGCCAGTTCCGGCGGTTTACCGAAGAAGTGGTACC  >  NC_002947/3653064‑3653340

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