Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,924,590 A→T 6.2% L75Q (CTG→CAG)  peaB ← quinohemoprotein amine dehydrogenase modification protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,924,5900AT 6.2% 196.5 / 3.0 64L75Q (CTG→CAG) peaBquinohemoprotein amine dehydrogenase modification protein
Reads supporting (aligned to +/- strand):  ref base A (28/32);  new base T (2/2);  total (30/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.41e-01

AGGTCTTCCTTGTAGCAATAGGTGCAGCTGAGGTTGCAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTGGATATCCTCGCCGGCAAACTGCTGGCCAAGGCGCTGCATCAGCTCCTCGGCCGAGCAGCCTTGCTGGCGCAAGGTGTCGATGATGCCGCCGGTCACCGCGTCAGCGGTGAACAGCGAACTGCTGGG  >  NC_002947/3924443‑3924735
                                                                                                                                                   |                                                                                                                                                  
aGGTCTTCCTTGTAGCAATAGGTGCAGCTGAGGTTGCAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAg                                                                                                                                                   >  6:336386/1‑149 (MQ=255)
 ggTCTTCCTTGTAGCAATAGGTGCAGCTGAGGTTGCAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGt                                                                                                                                                  >  7:260933/1‑149 (MQ=255)
    cTTCCTTGTAGCAATAGGTGCAGCTGAGGTTGCAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTcg                                                                                                                                               >  8:20536/1‑149 (MQ=255)
        cTTGTAGCAATAGGTGCAGCTGAGGTTGCAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTcgcg                                                                                                                                             >  2:259127/1‑147 (MQ=255)
         ttGTAGCAATAGGTGCAGCTGAGGTTGCAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAg                                                                                                                                          >  4:83590/1‑149 (MQ=255)
           gTAGCAATAGGTGCAGCTGAGGTTGCAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAgg                                                                                                                                         <  5:60799/148‑1 (MQ=255)
           gTAGCAATAGGTGCAGCTGAGGTTGCAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAg                                                                                                                                          <  1:85638/147‑1 (MQ=255)
               cAATAGGTGCAGCTGAGGTTGCAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTc                                                                                                                                    >  6:36221/1‑149 (MQ=255)
                 aTAGGTGCAGCTGAGGTTGCAACCGGTGTTGACGGTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCt                                                                                                                                  <  3:79176/149‑1 (MQ=255)
                  tAGGTGCAGCTGAGGTTGCAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTg                                                                                                                                 <  5:36221/149‑1 (MQ=255)
                     gTGCAGCTGAGGTTGCAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGGGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCTCGCAGGGTGTCCTggtt                                                                                                                              <  2:14575/149‑3 (MQ=255)
                     gTGCAGCTGAGGTTGCAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTGGat                                                                                                                              <  6:219937/149‑1 (MQ=255)
                     gTGCAGCTGAGGTTGCAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTGGa                                                                                                                               >  7:176027/1‑148 (MQ=255)
                      tGCAGCTGAGGTTGCAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTGGata                                                                                                                             <  1:259127/149‑1 (MQ=255)
                        cAGCTGAGGTTGCAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTGGATATc                                                                                                                           <  8:260933/149‑1 (MQ=255)
                         aGCTGAGGTTGCAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTGGATATcc                                                                                                                          >  2:277759/1‑149 (MQ=255)
                         aGCTGAGGTTGCAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTGGATATcc                                                                                                                          <  4:1535/149‑1 (MQ=255)
                           cTGAGGTTGCAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTGGATATCCTc                                                                                                                        >  3:153678/1‑149 (MQ=255)
                                 ttGCAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTGGATATCCTCGCCGGc                                                                                                                  <  4:270978/149‑1 (MQ=255)
                                   gCAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTGGATATCCTCGCCGGCa                                                                                                                 <  7:78656/148‑1 (MQ=255)
                                     aaCCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTGGATATc                                                                                                                           <  4:20048/136‑1 (MQ=255)
                                     aaCCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTGGATATCCTCGCCGGCAAAc                                                                                                              >  5:15342/1‑149 (MQ=255)
                                     aaCCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTGGATATCCTCGCCGGCAAAc                                                                                                              >  4:24274/1‑149 (MQ=255)
                                     aaCCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCCGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTGGATATc                                                                                                                           >  3:20048/1‑136 (MQ=255)
                                        cGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAgg                                                                                                                                         >  5:99717/1‑119 (MQ=255)
                                        cGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAgg                                                                                                                                         <  6:99717/119‑1 (MQ=255)
                                         ggTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTGGATATCCTCGCCGGCAAActgct                                                                                                          <  7:20536/149‑1 (MQ=255)
                                          gtgtTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTGGATATCCTCGCCGGCAAActgctg                                                                                                        <  5:101794/149‑1 (MQ=255)
                                               gACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTGGATATCCTCGCCGGCAAACTGCTGGcc                                                                                                     >  6:335779/1‑147 (MQ=255)
                                                 cGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTGGATATCCTCGCCGGCAAACTGCTGGCCAAgg                                                                                                 <  5:336386/149‑1 (MQ=255)
                                                  gTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTGGATATCCTCGCCGGCAAACTGCTGGCCAAGgc                                                                                                <  8:226596/149‑1 (MQ=255)
                                                      aGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTGGATATCCTCGCCGGCAAACTGCTGGCCAAGGCGCTg                                                                                            <  5:335779/149‑1 (MQ=255)
                                                      aGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTGGATATCCTCGCCGGCAAACTGCTGGCCAAGGCGCTg                                                                                            >  7:87649/1‑149 (MQ=255)
                                                       gCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTGGATATCCTCGCCGGCAAACTGCTGGCCAAGGCGCTGc                                                                                           <  4:153678/149‑1 (MQ=255)
                                                           cacGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTGGATATCCTCGCCGGCAAACTGCTGGCCAAGGCGCTGCAt                                                                                         >  1:12042/1‑147 (MQ=255)
                                                                     aGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTGGATATCCTCGCCGGCAAACTGCTGGCCAAGGCGCTGCATCa                                                                                       >  5:199868/1‑139 (MQ=255)
                                                                     aGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTGGATATCCTCGCCGGCAAACTGCTGGCCAAGGCGCTGCATCa                                                                                       <  6:199868/139‑1 (MQ=255)
                                                                                   gACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTGGATATCCTCGCCGGCAAACTGCTCGGCCGAGGAGCTGATGCAGCGCCTTGGCcagcagtttgccg                                                               <  6:15342/149‑8 (MQ=25)
                                                                                          cTGAGGGCGATTTCCGGTGTCAGCGGCGAGCCGGCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTGGATATCCTCGCCGGCAAACTGCTGGCCAAGGCGCTGCATCAGCTCCTCGGCCGAGCAGCCTTGCTGGCGCaa                                                        <  8:157237/149‑1 (MQ=255)
                                                                                          cTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTGGATATCCTCGCCGGCAAACTGCTGGCCAAGGCGCTGCATCAGCTCCTCGGCCGAGCAGCCTTGCTGGCGCaa                                                        <  3:83590/149‑1 (MQ=255)
                                                                                              gggCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGAGATCAGTTCGCGCAGGGTGTCCTGGATATCCTCGCCGGCAAACTGCTGGCCAAGGCGCTGCATCAGCTCCTCGGCCGAGCAGCCTTGCTGGCGCAAGgtg                                                    <  8:176027/149‑1 (MQ=255)
                                                                                                gCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTGGATATCCTCGCCGGCAAACTGCTGGCCAAGGCGCTGCATCAGCTCCTCGGCCGAGCAGCCTTGCTGGCGCaa                                                        >  7:301276/1‑143 (MQ=255)
                                                                                                gCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTGGATATCCTCGCCGGCAAACTGCTGCCCAAGGCGCTGCATCAGCTCCTCGGCCGAGCAGCCTTGCTGGCGCaa                                                        <  8:301276/143‑1 (MQ=255)
                                                                                                    tttCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTGGATATCCTCGCCGGCAAACTGCTGGCCAAGGCGCTGCATCAGCTCCTCGGCCGAGCAGCCTTGCTGGCGCAAGGTGTCgat                                               >  8:61922/1‑148 (MQ=255)
                                                                                                    tttCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTGGATATCCTCGCCGGCAAACTGCTGGCCAAGGCGCTGCATCAGCTCCTCGGCCGAGCAGCCTTGCTGGCGCAAGGTGTCgat                                               >  6:106407/1‑148 (MQ=255)
                                                                                                       ccGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTGGATATCCTCGCCGGCAAACTGCTGGCCAAGGCGCTGCATCAGCTCCTCGGCCGAGCAGCCTTGCTGGCGCAAGGTGTCgatgat                                            <  3:24274/148‑1 (MQ=255)
                                                                                                           ggTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTGGATATGCTCGCCGGCAAACTGCTGGCCAAGGCGCTGCATCAGCTCCTCGGCCGAGGAGCCTTGCTGGCGCAAGGTGTCGATGATgccgc                                       <  2:275323/149‑1 (MQ=255)
                                                                                                             tCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTGGATATCCTCGCCGGCAAACTGCTGGCCAAGGCGCTGCATCAGCTCCTCGGCCGag                                                                         <  2:321778/113‑1 (MQ=255)
                                                                                                               aGCGGCGAGCCGTCGCTGACCAGCTCCAGAGCGATCAGTTCGCGCAGGGTGTCCTGGATATCCTCGCCGGCAAACTGCTGGCCAAGGCGCTGCATCAGCTCCTCGGCCGAGCAGCCTTGCTGGCGCAAGGTGTCGATGATGCCGCCGGt                                   >  3:173597/1‑149 (MQ=255)
                                                                                                                         cGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTGGATATCCTCGCCGGCAAACTGCTGGCCAAGGCGCTGCATCAGCTCCTCGGCCGAGCAGCCTTGCTGGCGCAAGGTGTCGATGATGCCGCCGGTCACcgc                             >  1:57531/1‑145 (MQ=255)
                                                                                                                         cGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTGGATATCCTCGCCGGCAAACTGCTGGCCAAGGCGCTGCATCAGCTCCTCGGCCGAGCAGCCTTGCTGGCGCAAGGTGTCGATGATGCCGCCGGTCACcgc                             <  2:57531/145‑1 (MQ=255)
                                                                                                                         cGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTGGATATCCTCGCCGGCAAACTGCTGGCCAAGGCGCTGCATCAGCTCCTCGGCCGAGCAGCCTTGCTGGCGCAAGGTGTCGATGATGCCGCCGGTCACCGCGTCa                         <  2:12042/149‑1 (MQ=255)
                                                                                                                                 aCCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTGGATATCCTCGCCGGCAAACTGCTGGCCAAGGCGCTGCATCAGCTCCTCGGCCGAGCAGCCTTGCTGGCGCAAGGTGTCGATGATGCCGCCGGTCACCGCGTCAGCGGTGaa                 >  3:58006/1‑149 (MQ=255)
                                                                                                                                  ccAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTGGATATCCTCGCCGGCAAACTGCTGGCCAAGGCGCTGCATCAGCTCCTCGGCCGAGCAGCCTTGCTGGCGCAAGGTGTCGATGATGCCGCCGGTCACCGCGTCAGCGGTGaa                 >  7:171553/1‑148 (MQ=255)
                                                                                                                                        cctgcgcGAACTGATCGCGCAGGGTGTCCTGGATATCCTCGCCGGCAAACTGCTGGCCAAGGCGCTGCATCAGCTCCTCGGCCGAGCAGCCTTGCTGGCGCAAGGTGTCGATGATGCCGCCGGTCACCGCGTCAGCGGTGAACAg              >  3:136940/4‑145 (MQ=255)
                                                                                                                                        cctgcgcGAACTGATCGCGCAGGGTGTCCTGGATATCCTCGCCGGCAAACTGCTGGCCAAGGCGCTGCATCAGCTCCTCGGCCGAGCAGCCTTGCTGGCGCAAGGTGTCGATGATGCCGCCGGTCACCGCGTCAGCGGTGAACAg              <  4:136940/142‑1 (MQ=255)
                                                                                                                                         cAGCGCGATCAGTTCGCGCAGGGTGTCCTGGATATCCTCGCCGGCAAACTGCTGGCCAAGGCGCTGCATCAGCTCCTCGGCCGAGCAGCCTTGCTGGCGCAAGGTGTCga                                                <  7:202681/110‑1 (MQ=255)
                                                                                                                                         cAGCGCGATCAGTTCGCGCAGGGTGTCCTGGATATCCTCGCCGGCAAACTGCTGGCCAAGGCGCTGCATCAGCTCCTCGGCCGAGCAGCCTTGCTGGCGCAAGGTGTCga                                                >  8:202681/1‑110 (MQ=255)
                                                                                                                                          tgcgcGAACTGATCGCGCAGGGTGTCCTGGATATCCTCGCCGGCAAACTGCTGGCCAAGGCGCTGCATCAGCTCCTCGGCCGAGCAGCCTTGCTGGCGCAAGGTGTCGGTGATGCCGCCGGTCACCGCGTCAGCGGTGAACa               >  5:330640/2‑142 (MQ=255)
                                                                                                                                          tgcgcGAACTGATCGCGCAGGGTGTCCTGGATATCCTCGCCGGCAAACTGCTGGCCAAGGCGCTGCATCAGCTCCTCGGCCGAGCAGCCTTGCTGGCGCAAGGTGTCGGTGATGCCGCCGGTCACCGCGTCAGCGGTGAACa               <  6:330640/141‑1 (MQ=255)
                                                                                                                                              cgATCAGTTCGCGCAGGGTGTCCTGGATATCCTCGCCGGCAAACTGCTGGCCAAGGCGCTGCATCAGCTCCTCGGCCGAGCAGCCTTGCTGGCGCAAGgtgt                                                   >  6:78455/1‑102 (MQ=255)
                                                                                                                                              cgATCAGTTCGCGCAGGGTGTCCTGGATATCCTCGCCGGCAAACTGCTGGCCAAGGCGCTGCATCAGCTCCTCGGCCGAGCAGCCTTGCTGGCGCAAGgtgt                                                   <  5:78455/102‑1 (MQ=255)
                                                                                                                                                aTCAGTTCGCGCAGGGTGTCCTGGATATCCTCGCCGGCAAACTGCTGGCCAAGGCGCTGCATCAGCTCCTCGGCCGAGCAGCCTTGCTGGCGCAAGGTGTCGATGATGCCGCCGGTCACCGCGTCAGCGGTGAACAGCGAACTGCTggg  >  5:89019/1‑149 (MQ=255)
                                                                                                                                                aTCAGTTCGCGCAGGGTGTCCTGGATATCCTCGCCGGCAAACTGCTGGCCAAGGCGCTGCATCAGCTCCTCGGCCGAGCAGCCTTGCTGGCGCAAGGTGTCGATGATGCCGCCGGTCACCGCGTCAGCGGTGAACAGCGAACTGCTggg  >  3:203198/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
AGGTCTTCCTTGTAGCAATAGGTGCAGCTGAGGTTGCAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTGGATATCCTCGCCGGCAAACTGCTGGCCAAGGCGCTGCATCAGCTCCTCGGCCGAGCAGCCTTGCTGGCGCAAGGTGTCGATGATGCCGCCGGTCACCGCGTCAGCGGTGAACAGCGAACTGCTGGG  >  NC_002947/3924443‑3924735

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: