Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,653,199 T→G 6.4% L362* (TTA→TGA)  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,1990TG 6.4% 194.5 / 4.0 62L362* (TTA→TGA) PP_3218NtaA/SnaA/SoxA family monooxygenase
Reads supporting (aligned to +/- strand):  ref base T (32/26);  new base G (2/2);  total (34/28)
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.98e-01

ACCTGGATGCCCCGTTCCCGCAAGAGAGCCTGGAGTTTGCACGGAGCAGCTTTTTCACGAATGCCAAGCGTATTGCCGACCAGGCGCGGGACAAGGCCTTCAG‑TCTGCGT‑‑CA‑GGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGATGGGTTCAATATCTACATTGGTCATCCTGGCCAGTTCCGGCGGTTTACCGAAGAAGTGGTACCGTTG  >  NC_002947/3653064‑3653344
                                                                                                                                           |                                                                                                                                                 
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                                                                                                                               <  2:191807/124‑1 (MQ=255)
                                gAGTTTGCACGGAGCAGCTTTTTCACGAATGCCAAGCGTATTGCCGACCAGGCGCGGGACAAGGCCTTCAG‑TCTGCGT‑‑CA‑GGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTg                                                                                                                               >  1:191807/1‑124 (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                                                                                >  5:62316/1‑149 (MQ=255)
                                                      tCACGAATGCCAAGCGTATTGCCGACCAGGCGCGGGACAAGGCCTTCAG‑TCTGCGT‑‑CA‑GGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAgcg                                                                                >  6:215837/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                                                                                             <  4:112859/128‑1 (MQ=255)
                                                              gCCAAGCGTATTGCCGACCAGGCGCGGGACAAGGCCTTCAG‑TCTGCGT‑‑CA‑GGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGa                                                                                             >  3:112859/1‑128 (MQ=255)
                                                                 aaGCGAATTGCCGACCAGGCGCGGGACAAGGCCGTCAG‑TCTGCGT‑‑CA‑GGCAGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGAGGCGGCGAAGATGGCCGAGTGGATCGAAGCGCGGGCGCTgca                                                                     >  2:34475/1‑147 (MQ=255)
                                                                       aTTGCCGACCAGGCGCGGGACAAGGCCTTCAG‑TCTGCGT‑‑CA‑GGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTgg                                                                      <  5:137411/142‑1 (MQ=255)
                                                                       aTTGCCGACCAGGCGCGGGACAAGGCCTTCAG‑TCTGCGT‑‑CA‑GGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTgg                                                                      >  6:137411/1‑142 (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)
                                                                                                                                        aGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGATGGGTTCAATATCTACATTGGTCATCCTGGCCAGTTCCGGCGGTTTACCGAAGAAGTGGTACCGTtg  <  1:237382/149‑1 (MQ=255)
                                                                                                                                          ttACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGATGGGTTCAATATCTACATTGGTCATCCTGGCCAGTTCCGGCGGTTTa                      >  7:253354/1‑127 (MQ=255)
                                                                                                                                          ttACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGATGGGTTCAATATCTACATTGGTCATCCTGGCCAGTTCCGGCGGTTTa                      <  8:253354/127‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                                 
ACCTGGATGCCCCGTTCCCGCAAGAGAGCCTGGAGTTTGCACGGAGCAGCTTTTTCACGAATGCCAAGCGTATTGCCGACCAGGCGCGGGACAAGGCCTTCAG‑TCTGCGT‑‑CA‑GGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGATGGGTTCAATATCTACATTGGTCATCCTGGCCAGTTCCGGCGGTTTACCGAAGAAGTGGTACCGTTG  >  NC_002947/3653064‑3653344

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