Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,778,991 G→A 7.2% I92I (ATC→ATT) ‡ pvdD ← non‑ribosomal peptide synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,778,9910GA 7.2% 205.0 / 4.4 69I92I (ATC→ATT) ‡pvdDnon‑ribosomal peptide synthetase
Reads supporting (aligned to +/- strand):  ref base G (32/32);  new base A (3/2);  total (35/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 = 9.91e-01

CTCGACGTGCTCGCGCAAGGCGTCGCTGGTCGCGCTCATGCCGATGAGCGGCTCGATGAGCAACTCGAACGGCATCTGCGGGTGGATGATCTGCACCACGCGGTCGCTGTCTTGTCGCAACGTGGTGCGCAGCACTTCGTGGCGCTGGATCAGCGCCTGGCAGGCACCGCGCAAGCCCTCGATATCCAGCCGCCCCTGCAGCTGCAGCGCCGCACCGATGTTGTAGGCTCCGCTGGCGGGGTCCAGTTGCCAGAGAAACCACTGCCGCTCCTGGGCGTAGGACACCGGGGTCGGGT  >  NC_002947/4778844‑4779139
                                                                                                                                                   |                                                                                                                                                    
cTCGACGTGCTCGCGCAAGGCGTCGCTGGTCGCGCTCATGCCGATGAGCGGCTCGATGAGCAACTCGAACGGCATCTGCGGGTGGATGATCTGCACCACGCGGTCGCTGTCTTGTCGCAACGTGGTGCGCAGCACTTCGTGGCGCTgg                                                                                                                                                      <  1:56629/148‑1 (MQ=255)
cTCGACGTGCTCGCGCAAGGCGTCGCTGGTCGCGCTCATGCCGATGAGCGGCTCGATGAGCAACTCGAACGGCATCTGCGGGTGGATGATCTGCACCACGCGGTCGCTGTCTTGTCGCAACGTGGTGCGCAGCACTTCGTGGCGCTgg                                                                                                                                                      >  2:56629/1‑148 (MQ=255)
    aCGTGCTCGCGCAAGGCGTCGCTGGTCGCGCTCATGCCGATGAGCGGCTCGATGAGCAACTCGAACGGCATCTGCGGGTGGATGATCTGCACCACGCGGTCGCTGTCTTGTCGCAACGTGGTGCGCAGCACTTCGTGGCGCTGGATCAg                                                                                                                                                 <  3:234595/149‑1 (MQ=255)
         cTCGCGCAAGGCGTCGCTGGTCGCGCTGATGCCGATGAGCGGCTCGATGAGCAACTCGAACGGCATCTGCGGGTGGATGATCTGCACCACGCGGTCGCTGTCTTGTCGCAACGTGGTGCGCAGCACTTCGTGGCGATGGATCAGCGcc                                                                                                                                             <  6:33621/148‑1 (MQ=255)
          tCGCGCAAGGCGTCGCTGGTCGCGCTCATGCCGATGAGCGGCTCGATGAGCAACTCGAACGGCATCTGCGGGTGGATGATCTGCACCACGCGGTCGCTGTCTTGTCGCAACGTGGTGCGCAGCACTTCGTGGCGCTGATCCAGCGcca                                                                                                                                            >  5:87623/1‑147 (MQ=255)
           cgcgCAAGGCGTCGCTGGTCGCGCTCATGCCGATGAGCGGCTCGATGAGCAACTCGAACGGCATCTGCGGGTGGATGATCTGCACCACGCGGTCGCTGTCTTGTCGCAACGTGGTGCGCAGCACTTCGTGGCGCTGGATCAGCGCCTgg                                                                                                                                          >  5:33845/1‑149 (MQ=255)
            gcgcAAGGCGTCGCTGGTCGCGCTCATGCCGATGAGCGGCTCGATGAGCAACTCGAACGGCATCTGCGGGTGGATGATCTGCACCACGCGGTCGCTGTCTTGTCGCAACGTGGTGCGCAGCACTTCGTGGCGCTGGATCAGCGCCTgg                                                                                                                                          <  6:33845/148‑1 (MQ=255)
            gcgcAAGGCGTCGCTGGTCGCGCTCATGCCGATGAGCGGCTCGATGAGCAACTCGAACGGCATCTGCGGGTGGATGATCTGCACCACGCGGTCGCTGTCTTGTCGCAACGTGGTGCGCAGCACTTCGTGGCGCTGATCCAGCGccacg                                                                                                                                          <  2:104546/148‑4 (MQ=255)
            gcgcAAGGCGTCGCTGGTCGCGCTCATGCCGATGAGCGGCTCGATGAGCAACTCGAACGGCATCTGCGGGTGGATGATCTGCACCACGCGGTCGCTGTCTTGTCGCAACGTGGTGCGCAGCACTTCGTGGCGCTGATCCAGCGcca                                                                                                                                            <  6:87623/146‑2 (MQ=255)
             cgcAAGGCGTCGCTGGTCGCGCTCATGCCGATGAGCGGCTCGATGAGCAACTCGAACGGCATCTGCGGGTGGATGATCTGCCCCACGCGGTCGCTGTCTTGTCGCAACGTGGTGCGCAGCACTTCGTGGCGCTGGATCAGCGCCTgg                                                                                                                                          >  8:279734/1‑147 (MQ=255)
                  ggCGTCGCTGGTCGCGCTCATGCCGATGAGCGGCTCGATGAGCAACTCGAACGGCATCTGCGGGTGGATGATCTGCACCACGCGGTCGCTGTCTTGTCGCAACGTGGTGCGCAGCACTTCGTGGCGCTGGATCAGCGCCTGGCAGGCAc                                                                                                                                   >  6:42029/1‑149 (MQ=255)
                     gTCGCTGGTCGCGCTCATGCCGATGAGCGGCTCGATGAGCAACTCGAACGGCATCTGCGGGTGGATGATCTGCACCACGCGGTCGCTGTCTTGTCGCAACGTGGTGCGCAGCACTTCGTGGCGCTGGATCAGCGCCTGGCAGGCACcgc                                                                                                                                <  5:42029/149‑1 (MQ=255)
                         cTGGTCGCGCTCATGCCGATGAGCGGCTCGATGAGCAACTCGAACGGCATCTGCGGGTGGATGATCTGCACCACGCGGTCGCTGTCTTGTCGCAACGTGGTGCGCAGCACTTCGTGGCGCTGGATCAGCGCCTGGCAGGCACCGCGCaa                                                                                                                            <  7:279734/149‑1 (MQ=255)
                            gTCGCGCTCATGCCGATGAGCGGCTCGATGAGCAACTCGAACGGCATCTGCGGGTGGATGATCTGCACCACGCGGTCGCTGTCTTGTCGCAACGTGGTGCGCAGCACTTCGTGGCGCTGGATCAGCGCCTGGCAGGCACCGCGCAAg                                                                                                                           <  8:8054/147‑1 (MQ=255)
                                 gcTCATGCCGATGAGCGGCTCGATGAGCAACTCGAACGGCATCTGCGGGTGGATGATCTGCACCACGCGGTCGCTGTCTTGTCGCAACGTGGTGCGCAGCACTTCGTGGCGCTGGATCAGCGCCTGGCAGGCACCGCGCAAGCCCTCGa                                                                                                                    <  1:297789/149‑1 (MQ=255)
                                      tGCCGATGAGCGGCTCGATGAGCAACTCGAACGGCATCTGCGGGTGGATGATCTGCACCACGCGGTCGCTGTCTTGTCGCAACGTGGTGCGCAGCACTTCGTGGCGCTGGATCAGCGCCTGGCAGGCACCGCGCAAGCCCTCGATATc                                                                                                                >  5:121004/1‑148 (MQ=255)
                                               gCGGCTCGATGAGCAACTCGAACGGCATCTGCGGGTGGATGATCTGCACCACGCGGTCGCTGTCTTGTCGCAACGTGGTGCGCAGCACTTCGTGGCGCTGGATCAGCGCCTGGCAGGCACCGCGCAAGCCCTCGATATCCAGCCGcccc                                                                                                      <  1:340838/149‑1 (MQ=255)
                                                  gCTCGATGAGCAACTCGAACGGCATCTGCGGGTGGATGATCTGCACCACGCGGTCGCTGTCTTGTCGCAACGTGGTGCGCAGCACTTCGTGGCGCTGGATCAGCGCCTGGCAGGCACCGCGCAAGCCCTCGATATCCAGCCGCCCctg                                                                                                    >  5:192215/1‑148 (MQ=255)
                                                      gATGAGCAACTCGAACGGCATCTGCGGGTGGATGATCTGCACCACGCGGTCGCTGTCTTGTCGCAACGTGGTGCGCAGCACTTCGTGGCGCTGGATCAGCGCCTGGCAGGCACCGCGCAAGCCCTCGATATCCAGCCGCCCctgcagct                                                                                               <  8:77967/149‑1 (MQ=255)
                                                       aTGAGCAACTCGAACGGCATCTGCGGGTGGATGATCTGCACCACGCGGTCGCTGTCTTGTCGCAACGTGGTGCGCAGCACTTCGTGGCGCTGGATCAGCGCCTGGCAGGCACCGCGCAAGCCCTCGATATCCAGCCGCCCctgcagctg                                                                                              <  4:338295/149‑1 (MQ=255)
                                                         gAGCAACTCGAACGGCATCTGCGGGTGGATGATCTGCACCACGCGGTCGCAGTCTTGTCGCAACGTGGTGCGCAGCACTTCGTGGCGCTGGATCAGCGCCTGGCAGGCACCGCGCAAGCCCTCGATATCCAGCCGCCCctgcagctg                                                                                              >  6:70643/1‑147 (MQ=255)
                                                            cAACTCGAACGGCATCTGCGGGTGGATGATCTGCACCACGCGGTCGCTGTCTTGTCGCAACGTGGTGCGCAGCACTTCGTGGCGCTGGATCAGCGCCTGGCAGGCACCGCGCAAGCCCTCGATATCCAGCCGCCCctg                                                                                                    <  6:79167/138‑1 (MQ=255)
                                                            cAACTCGAACGGCATCTGCGGGTGGATGATCTGCACCACGCGGTCGCTGTCTTGTCGCAACGTGGTGCGCAGCACTTCGTGGCGCTGGATCAGCGCCTGGCAGGCACCGCGCAAGCCCTCGATATCCAGCCGCCCctg                                                                                                    >  5:79167/1‑138 (MQ=255)
                                                             aaCTCGAACGGCATCTGCGGGTGGATGATCTGCACCACGCGGTCGCTGTCTTGTCGCAACGTGGTGCGCAGCACTTCGTGGCGCTGGATCAGCGCCTGGCAGGCACCGCGCAAGCCCTCGATATCCAGCCGCCCCTGCAGCTGCAgcgc                                                                                        >  6:239486/1‑149 (MQ=255)
                                                              aCTCGAACGGCATCTGCGGGTGGATGATCTGCACCACGCGGTCGCTGTCTTGTCGCAACGTGGTGCGCAGCACTTCGTGGCGCTGGATCAGCGCCTGGCAGGCACCGCGCAAGCCCTCGATATCCAGCCGCCCCTGCAGCTGCAGcgcc                                                                                       >  6:261060/1‑149 (MQ=255)
                                                               cTCGAACGGCATCTGCGGGTGGATGATCTGCACCACGCGGTCGCTGTCTTGTCGCAACGTGGTGCGCAGCACTTCGTGGCGCTGGATCAGCGCCTGGCAGGCACCGCGCAAGCCCTCGATATCCAGCCGCCCCTGCAGCTGCAGcgccg                                                                                      >  3:79726/1‑149 (MQ=255)
                                                               cTCGAACGGCATCTGCGGGTGGATGATCTGCACCACGCGGTCGCTGTCTTGTCGCAACGTGGTGCGCAGCACTTCGTGGCGCTGGATCAGCGCCTGGCAGGCACCGCGCAAGCCCTCGATATCCAGCCGCCCCTGCAGCTGCAGcgccg                                                                                      >  7:195957/1‑149 (MQ=255)
                                                               cTCGAACGGCATCTGCGGGTGGATGATCTGCACCACGCGGTAGCTGTCTTGTCGCAACGTGGTGCGCAGCACTTCGTGGCGCTGGATCAGCGCCTGGCAGGCACCGCGCAAGCCCTCGATATCCAGCCGCCCCTGCAGCTGCAGcgccg                                                                                      >  5:23329/1‑149 (MQ=255)
                                                                tCGAACGGCATCTGCGGGTGGATGATCTGCACCACGCGGTCGCTGTCTTGTCGCAACGTGGTGCGCAGCACTTCGTGGCGCTGGATCAGCGCCTGGCAGGCACCGCGCAAGCCCTCGATATCCAGCCGCCCCTGCAGCTGCAGcgccg                                                                                      >  6:296421/1‑148 (MQ=255)
                                                                 cGAACGGCATCTGCGGGTGGATGATCTGCACCACGCGGTCGCTGTCTTGTCGCAACGTGGTGCGCAGCACTTCGTGGCGCTGGATCAGCGCCTGGCAGGCAc                                                                                                                                   <  1:35737/102‑1 (MQ=255)
                                                                 cGAACGGCATCTGCGGGTGGATGATCTGCACCACGCGGTCGCTGTCTTGTCGCAACGTGGTGCGCAGCACTTCGTGGCGCTGGATCAGCGCCTGGCAGGCAc                                                                                                                                   >  2:35737/1‑102 (MQ=255)
                                                                             gCGGGTGGATGATCTGCACCACGCGGTCGCTGTCTTGTCGCAACGTGGTGCGCAGCACTTCGTGGCGCTGGATCAGCGCCTGGCAGGCACCGCGCAAGCCCTCGATATCCAGCCGCCCCTGCAGCTGCAGCGCCGCACCGATGTTGTAg                                                                        >  5:106380/1‑149 (MQ=255)
                                                                              cGGGTGGATGATCTGCACCACGCGGTCGCTGTCTTGTCGCAACGTGGTGCGCAGCACTTCGTGGCGCTGGATCAGCGCCTGGCAGGCACCGCGCAAGCCCTCGATATCCAGCCGCCCCTGCAGCTGCAGCGCCGCACCGATGTTGTa                                                                         <  8:166708/147‑1 (MQ=255)
                                                                              agggTGGATGATCTGCACCACGCGGTCGCTGGCTTGTCTCAACGTGGTGCGCAGCACTTCGTGGCTCTGGATCAGCGCCTGGCAGGCTCCGCGCAAGCCCTCGATATCCAGCCGCCCCTGCAGCTGCATCGCCGCACCGAtgttgt                                                                          <  2:227815/145‑1 (MQ=255)
                                                                                ggTGGATGATCTGCACCACGCGGTCGCTGTCTTGTCGCAACGTGGTGCGCAGCACTTCGTGGCGCTGGATCAGCGCCTGGCAGGCACCGCGCAAGCCCTCGATATCCAGCCGCCCCTGCAGCTGCAGCGCCGCACCGATGTTGTAGGCt                                                                     <  6:192215/149‑1 (MQ=255)
                                                                                 gTGGATGATCTGCACCACGCGGTCGCTGTCTTGTCGCAACGTGGTGCGCAGCACTTCGTGGCGCTGGATCAGCGCCTGGCAGGCACCGCGCAAGCCCTCGATATCCAGCCGCCCCTGCAGCTGCAGCGCCGCACCGAtgttg                                                                           <  6:36608/142‑1 (MQ=255)
                                                                                 gTGGATGATCTGCACCACGCGGTCGCTGTCTTGTCGCAACGTGGTGCGCAGCACTTCGTGGCGCTGGATCAGCGCCTGGCAGGCACCGCGCAAGCCCTCGATATCCAGCCGCCCCTGCAGCTGCAGCGCCGCACCGAtgttg                                                                           >  5:36608/1‑142 (MQ=255)
                                                                                   ggATGATCTGCACCACGCGGTCGCTGTCTTGTCGCAACGTGGTGCGCAGCACTTCGTGGCGCTGGATCAGCGCCTGGCAGGCACCGCGCAAGCCCTCGATATCCAGCCGCCCCTGCCGCTGCAGAGCCGCACAGATGTTGTAGGCTCCg                                                                  >  6:39929/1‑149 (MQ=255)
                                                                                    gatgatCTGCACCACGCGGTCGCTGTCTTGTCGCAACGTGGTGCGCAGCACTTCGTGGCGCTGGATCAGCGCCTGGCAGGCACCGCGCAAGCCCTCGATATCCAGCCGCCCCTGCAGCTGCAGCGCCGCACCGATGTTGTAGGCTCCGc                                                                 >  1:198336/1‑149 (MQ=255)
                                                                                    gatgatCTGCACCACGCGGTCGCTGTCTTGTCGCAACGTGGTGCGCAGCACTTCGTGGCGCTGGATCAGCGCCTGGCAGGCACCGCGCAAGCCCTCGATATCCAGCCGCCCCTGCAGCTGCAGCGCCGCACCGATGTTGTAGGCTCCGc                                                                 >  2:294631/1‑149 (MQ=255)
                                                                                      tgatCTGCACTACGCGGTCGCTGTCTTTTCGCAACGTGGTGCGCAGCACTTCGTGGCGCTGGATCAGCGCCTGGCAGGCACCGCGCAAGCCCTCTATATCCAGCCGCCCCTGCAGCTGCAGCGCCGCACCGATGTTGTAGGCTCCGCTg                                                               <  6:23329/149‑1 (MQ=255)
                                                                                             caccacGCGGTCGCTGTCTTGTCGCAACGTGGTGCGCAGCACTTCGTGGCGCTGGATCAGCGCCTGGCAGGCACCGCGCAAGCCCTCGATATCCAGCCGCCCCTGCAGCTGCAGCGCCGCACCGATGTTGTAGGCTCCGCTGGCgggg                                                         <  5:39929/148‑1 (MQ=255)
                                                                                                cacGCGGTCGCTGTCTTGTCGCAACGTGGTGCGCAGCACTTCGTGGCGCTGGATCAGCGCCTGGCAGGCACCGCGCAAGCCCTCGATATCCAGCCGCCCCTGCAGCTGCAGCGCCGCACCGATGTTGTAGGCTCCGCTGGCGGGGTCCa                                                     >  2:56209/1‑149 (MQ=255)
                                                                                                 acGCGGTCGCTGTCTTGTCGCAACGTGGTGCGCAGCACTTCGTGGCGCTGGATCAGCGCCTGGCAGGCACCGCGCAAGCCCTCGATATCCAGCCGCCCCTGCAGCTGCAGCGCCGCACCGATGTTGTAGGCTCCGCTGGCGGGGTCCAg                                                    >  7:126191/1‑149 (MQ=255)
                                                                                                   gcgGTCGCTGTCTTGTCGCAACGTGGTGCGCAGCACTTCGTGGCGCTGGATCAGCGCCTGGCAGGCACCGCGCAAGCCCTCGATATCCAGCCGCCCCTGCAGCTGCAGCGCCGCACCGATGTTGTAGGCTCCGCTGGCGGGGTCCAGtt                                                  <  5:261060/149‑1 (MQ=255)
                                                                                                   gcgGTCGCTGTCTTGTCGCAACGTGGTGCGCAGCACTTCGTGGCGCTGGATCAGCGCCTGGCAGGCACCGCGCAAGCCCTCGATATCCAGCCGCCCCTGCAGCTGCAGCGCCGCACCGATGTTGTAGGCTCCGCTGGCGGGGTCCAGtt                                                  <  4:315578/149‑1 (MQ=255)
                                                                                                     ggTCGCTGTCTTGTCGCAACGTGGTGCGCAGCACTTCGTGGCGCTGGATCAGCGCCTGGCAGGCACCGCGCAAGCCCTCGATATCCAGCCGCCCCTGCAGCTGCAGCGCCGCACCGATGTTGTAGGCTCCGCTGGCGGGGTCCAGTTGc                                                <  3:47218/149‑1 (MQ=255)
                                                                                                         gCTGTCTTGTCGCAACGTGGTGCGCAGCACTTCGTGGCGCTGGATCAGCGCCTGGCAGGCACCGCGCAAGCCCTCGATATCCAGCCGCCCCTGCAGCTGCAGCGCCGCACCGATGTTGTAGGCTCCGCTGGCGGGGTCCAGTTGCCaga                                            <  4:79726/149‑1 (MQ=255)
                                                                                                         gCTGTCTTGTCGCAACGTGGTGCGCAGCACTTCGTGGCGCTGGATCAGCGCCTGGCAGGCACCGCGCAAGCCCTCGATATCCAGCCGCCCCTGCAGCTGCAGCGCCGCACCGATGTTGTAGGCTCCGCTGGCGGGGTCCAGTTGCCa                                              >  5:198675/1‑147 (MQ=255)
                                                                                                          cTGTCTTGTCGCAACGTGGTGCGCAGCACTTCGTGGCGCTGGATCAGCGCCTGGCAGGCACCGCGCAAGCCCTCGATATCCAGCCGCCCCTGCAGCTGCAGCGCCGCACCGATGTTGTAGGCTCCGCTGGCGGGGTCCAGTTGCCagag                                           <  8:126191/149‑1 (MQ=255)
                                                                                                           tGTCTTGTCGCAACGTGGTGCGCAGCACTTCGTGGCGCTGGATCAGCGCCTGGCAGGCACCGCGCAAGCCCTCGATATCCAGCCGCCCCTGCAGCTGCAGCGCCGCACCGATGTTGTAGGCTCCGCTGGCGGGGTCCAGTTGCCagag                                           >  6:287576/1‑148 (MQ=255)
                                                                                                               ttGTCGCAACGTGGTGCGCAGCACTTCGTGGCGCTGGATCAGCGCCTGGCAGGCACCGCGCAAGCCCTCGATATCCAGCCGCCCCTGCAGCTGCAGCGCCGCACCGATGTTGTAGGCTcc                                                                   >  2:317920/1‑120 (MQ=255)
                                                                                                               ttGTCGCAACGTGGTGCGCAGCACTTCGTGGCGCTGGATCAGCGCCTGGCAGGCACCGCGCAAGCCCTCGATATCCAGCCGCCCCTGCAGCTGCAGCGCCGCACCGATGTTGTAGGCTcc                                                                   <  1:317920/120‑1 (MQ=255)
                                                                                                                      aaCGTGGTGCGCAGCACTTCGTGGCGCTGGATCAGCGCCTGGCAGGCACCGCGCAAGCCCTCGATATCCAGCCGCCCCTGCAGCTGCAGCGCCGCACCGATGTTGTAGGCTCCGCTGGCGGGGTCCAGTTGCCAGAGAAACCACTGc                                 >  7:218217/1‑147 (MQ=255)
                                                                                                                      aaCGTGGTGCGCAGCACTTCGTGGCGCTGGATCAGCGCCTGGCAGGCACCGCGCAAGCCCTCGATATCCAGCCGCCCCTGCAGCTGCAGCGCCGCACCGATGTTGTAGGCTCCGCTGGCGGGGTCCAGTTGCCAGAGAAACCACTGc                                 <  8:218217/147‑1 (MQ=255)
                                                                                                                           ggtgCGCAGCACTTCGTGGCGCTGGATCAGCGCCTGGCAGGCACCGCGCAAGCCCTCGATATCCAGCCGCCCCTGCAGCTGCAGCGCCGCACCGATGTTGTAGGCTCCGCTGGCGGGGTCCAGTTGCCAGAGAAACCACTGCCGCTCCt                          <  5:287576/149‑1 (MQ=255)
                                                                                                                            gtgCGCAGCACTTCGTGGCGCTGGATCAGCGCCTGGCAGGCACCGCGCAAGCCCTCGATATCCAGCCGCCCCTGCAGCTGCAGCGCCGCACCGATGTTGTAGGCTCCGCTGGCGGGGTCCAGTTGCCAGAGAAACCACTGCCGCTCCTg                         >  5:250642/1‑149 (MQ=255)
                                                                                                                              gcgcAGCACTTCGTGGCGCTGGATCAGCGCCTGGCAGGCACCGCGCAAGCCCTCGATATCCAGCCGCCCCTGCAGCTGCAGCGCCGCACCGATGTTGTAGGCTCCGCTGGCGGGGTCCAGTTGCCAGAGAAACCACTGCCGCTCCTggg                       >  8:305375/1‑149 (MQ=255)
                                                                                                                               cgcAGCACTTCGTGGCGCTGGATCAGCGCCTGGCAGGCACCGCGCAAGCCCTCGATATCCAGCCGCCCCTGCAGCTGCAGCGCCGCACCGATGTTGTAGGCTCCGCTGGCGGGGTCCAGTTGCCAGAGAAACCACTGCCGCTCCTggg                       <  7:179854/148‑1 (MQ=255)
                                                                                                                                gcagcaCTTCGTGGCGCTGGATCAGCGCCTGGCAGGCACCGCGCAAGCCCTCGATATCCAGCCGCCCCTGCAGCTGCAGCGCCGCa                                                                                    >  2:63267/1‑86 (MQ=255)
                                                                                                                                 cagcaCTTCGTGGCGCTGGATCAGCGCCTGGCAGGCACCGCGCAAGCCCTCGATATCCAGCCGCCCCTGCAGCTGCAGCGCCGCACCGATGTTGTAGGCTCCGCTGGCGGGGTCCAGTTGCCAGAGAAACCACTGCCGCTCCTGGGCg                     >  8:269621/1‑148 (MQ=255)
                                                                                                                                 cagcaCTTCGTGGCGCTGGATCAGCGCCTGGCAGGCACCGCGCAAGCCCTCGATATCCAGCCGCCCCTGCAGCTACAGCGCCGCACCGATGTTGTAGGCTCCGCTGGCGGGGTCCAGTTGcc                                               <  2:109803/122‑1 (MQ=255)
                                                                                                                                 cagcaCTTCGTGGCGCTGGATCAGCGCCTGGCAGGCACCGCGCAAGCCCTCGATATCCAGCCGCCCCTGCAGCTACAGCGCCGCACCGATGTTGTAGGCTCCGCTGGCGGGGTCCAGTTGcc                                               >  1:109803/1‑122 (MQ=255)
                                                                                                                                  agcaCTTCGTGGCGCTGGATCAGCGCCTGGCAGGCACCGCGCAAGCCCTCGATATCCAGCCGCCCCTGCAGCTGCAGCGCCGCACCGATGTTGTAGGCTCCGCTGGCGGGGTCCAg                                                    >  6:315568/1‑116 (MQ=255)
                                                                                                                                  agcaCTTCGTGGCGCTGGATCAGCGCCTGGCAGGCACCGCGCAAGCCCTCGATATCCAGCCGCCCCTGCAGCTGCAGCGCCGCACCGATGTTGTAGGCTCCGCTGGCGGGGTCCAg                                                    <  5:315568/116‑1 (MQ=255)
                                                                                                                                   gcaCTTCGTGGCGCTGGATCAGCGCCTGGCAGGCACCGCGCAAGCCCTCGATATCCAGCCGCCCCTGCAGCTGCAGCGCCGCACCGATGTTGTAGGCTCCGCTGGCGGGGTCCAGTTGCCAGAGAAACCACTGCCGCTCCTGGGCGTAg                  <  6:121004/149‑1 (MQ=255)
                                                                                                                                     aCTTCGTGGCGCTGGATCAGCGCCTGGCAGGCACCGCGCAAGCCCTCGATATCCAGCCGCCCCTGCAGCTGCAGCGCCGCACCGATGTTGTAGGCTCCGCTGGCGGGGTCCAGTTGCCAGAGAAACCACTGCCGCTCCTGGGCGTAGGa                <  6:106380/149‑1 (MQ=255)
                                                                                                                                       ttCGTGGCGCTGGATCAGCGCCTGGCAGGCACCGCGCAAGCCCTCGATATCCAGCCGCCCCTGCAGCTGCAGCGCCGCACCGATGTTGTAGGCTCCGCTGGCGGGGTCCAGTTGCCAGAGAAACCACTGCCGCTCCTGGGCGTAGGaca              >  8:192638/1‑149 (MQ=255)
                                                                                                                                         ccaggCGCTGATCCAGCGCCTGGCAGGCACCGCGCAAGCCCTCGATATCCAGCCGCCCCTGCAGCTGCAGCGCCGCACCGATGTTGTAGGCTCCGCTGGCGGGGTCCAGTTGCCAGAGAAACCACTGCCGCTCCTGGGCGTAGGACAcc            >  2:192491/4‑149 (MQ=255)
                                                                                                                                          gTGGCGCTGGATCAGCGCCTGGCAGGCACCGCGCAAGCCCTCGATATCCAGCCGCCCCTGCAGCTGCAGCGCCGCACCGATGTTGTAGGCTCCGCTGGCGGGGTCCAGTTGCCAGAGAAACCACTGCCGCTCCTGGGCGTAGGACACCg           <  5:239486/149‑1 (MQ=255)
                                                                                                                                          caggCGCTGATCCAGCGCCTGGCAGGCACCGCGCAAGCCCTCGATATCCAGCCGCCCCTGCAGCTGCAGCGCCGCACCGATGTTGTAGGCTCCGCTGGCGGGGTCCAGTTGCCAGAGAAACCACTGCCGCTCCTGGGCGTAGGACACCg           >  4:119370/3‑149 (MQ=255)
                                                                                                                                                   gATCAGCGCCTGGCAGGCACCGCGCAAGCCCTCGATATCCAGCCGCCCCTGCAGCTGCAGCGCCGCACCGATGTTGTAGGCTCCGCTGGCGGGGTCCAGTTGCCAGAGAAACCACTGCCGCTCCTGGGCGTAGGACACCGgggtcgggt  >  8:222480/1‑149 (MQ=255)
                                                                                                                                                   gATCAGCGCCTGGCAGGCACCGCGCAAGCCCTCGATATCCAGCCGCCCCTGCAGCTGCAGCGCCGCACCGATGTTGTAGGCTCCGCTGGCGGGGTCCAGTTGCCAGAGAAACCACTGCCGCTCCTGGGCGTAGGACACCGgggtcgg    <  2:198336/147‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CTCGACGTGCTCGCGCAAGGCGTCGCTGGTCGCGCTCATGCCGATGAGCGGCTCGATGAGCAACTCGAACGGCATCTGCGGGTGGATGATCTGCACCACGCGGTCGCTGTCTTGTCGCAACGTGGTGCGCAGCACTTCGTGGCGCTGGATCAGCGCCTGGCAGGCACCGCGCAAGCCCTCGATATCCAGCCGCCCCTGCAGCTGCAGCGCCGCACCGATGTTGTAGGCTCCGCTGGCGGGGTCCAGTTGCCAGAGAAACCACTGCCGCTCCTGGGCGTAGGACACCGGGGTCGGGT  >  NC_002947/4778844‑4779139

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: