Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,710,722 T→C 6.8% N274D (AAC→GAC)  gpsA ← glycerol‑3‑phosphate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,710,7220TC 6.8% 212.1 / 4.7 73N274D (AAC→GAC) gpsAglycerol‑3‑phosphate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (38/30);  new base C (2/3);  total (40/33)
Fisher's exact test for biased strand distribution p-value = 6.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.11e-01

GGTGCTCGATCACCTGGTTCAGCGTGCGCCCTTCGAACAGGATGGCATGCAGGCCAGCCACCAGCGGCATGTACACCTGCACTTGCTGGGCCTTGGTCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGCGTAGCCAACCTGATAGTTGCGGCTCTTGGGCGAGGAGCAGGTGACGATCAGGTCACCGACACCGGCCAGGCCGAGGAACGTCATGGGGTTGGCGCCCTGGCTGACGGCG  >  NC_002947/4710575‑4710867
                                                                                                                                                    |                                                                                                                                                 
ggTGCTCGATCACCTGGTTCAGCGTGCGCCCTTCGAACAGGATGGCATGCAGGCCAGCCACCAGCGGCATGTACACCTGCACTTGCTGGGCCTTGGTCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGtt                                                                                                                                                   >  1:200847/1‑148 (MQ=255)
ggTGCTCGATCACCTGGTTCAGCGTGCGCCCTTCGAACAGGATGGCATGCAGGCCAGCCACCAGCGGCATGTACACCTGCACTTGCTGGGCCTTGGTCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTg                                                                                                                                                  >  6:98342/1‑149 (MQ=255)
       gATCACCTGGTTCAGCGTGCGCCCATCGAACAGGATGGCATGCAGGCCAGCCACCAGCGGCATGTACACCTGCACTTGCTGGGCCTTGGTCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCct                                                                                                                                           >  7:56337/1‑149 (MQ=255)
        aTCACCTGGTTCAGCGTGCGCCCTTCGAACAGGATGGCATGCAGGCCAGCCACCAGCGGCATGTACACCTGCACTTGCTGGGCCTTGGTCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCctc                                                                                                                                          >  5:385152/1‑149 (MQ=255)
        aTCACCTGGTTCAGCGTGCGCCCTTCGAACAGGATGGCATGCAGGCCAGCCACCAGCGGCATGTACACCTGCACTTGCTGGGCCTTGGTCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCctc                                                                                                                                          <  2:200847/149‑1 (MQ=255)
         tCACCTGGTTCAGCGTGCGCCCTTCGAACAGGATGGCATGCAGGCCAGCCACCAGCGGCATGTACACCTGCACTTGCTGGGCCTTGGTCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTCAACCGCctgg                                                                                                                                         <  5:107801/149‑3 (MQ=255)
         tCACCTGGTTCAGCGTGCGCCCTTCGAACAGGATGGCATGCAGGCCAGCCACCAGCGGCATGTACACCTGCACTTGCTGGGCCTTGGTCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTCAACCGCctgg                                                                                                                                         <  5:7376/149‑3 (MQ=255)
             cTGGTTCAGCGTGCGCCCTTCGAACAGGATGGCATGCAGGCCAGCCACCAGCGGCATGTACACCTGCACTTGCTGGGCCTTGGTCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTc                                                                                                                                       >  8:203122/1‑147 (MQ=255)
             cTGGTTCAGCGTGCGCCCTTCGAACAGGATGGCATGCAGGCCAGCCACCAGCGGCATGTACACCTGCACTTGCTGGGCCTTGGTCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCa                                                                                                                                     <  6:367317/149‑1 (MQ=255)
                   cAGCGTGCGCCCTTCGAACAGGATGGCATGCAGGCCAGCCACCAGCGGCATGTACACCTGCACTTGCTGGGCCTTGGTCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTcc                                                                                                                                      >  2:92434/1‑142 (MQ=255)
                   cAGCGTGCGCCCTTCGAACAGGATGGCATGCAGGCCAGCCACCAGCGGCATGTACACCTGCACTTGCTGGGCCTTGGTCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTcc                                                                                                                                      <  1:92434/142‑1 (MQ=255)
                    aGCGTGCGCCCTTCGAACAGGATGGCATGCAGGCCAGCCACCAGCGGCATGTACACCTGCACTTGCTGGGCCTTGGTCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTg                                                                                                                              <  2:100981/149‑1 (MQ=255)
                      cGTGCGCCCTTCGAACAGGATGGCATGCAGGCCAGCCACCAGCGGCATGTACACCTGCACTTGCTGGGCCTTGGTCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGc                                                                                                                            >  1:368332/1‑149 (MQ=255)
                            cccTTCGAACAGGATGGCATGCAGGCCAGCCACCAGCGGCATGTACACCTGCACTTGCTGGGCCTTGGTCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGcc                                                                                                                      >  1:193734/1‑149 (MQ=255)
                                  gAACAGGATGGCATGCAGGCCAGCCACCAGCGGCATGTACACCTGCCACTTGCTGGGCCTTGGTCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGAcc                                                                                                                                               >  1:275887/1‑119 (MQ=255)
                                   aaCAGGATGGCATGCAGGCCAGCCACCAGCGGCATGTACACCTGCACTTGCTGGGCCTTGGTCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGc                                                                                                               <  6:385152/149‑1 (MQ=255)
                                        gATGGCATGCAGGCCAGCCACCAGCGGCATGTACACCTGCACTTGCTGGGCCTTGGTCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGAcc                                                                                                                                               >  8:63123/1‑112 (MQ=255)
                                        gATGGCATGCAGGCCAGCCACCAGCGGCATGTACACCTGCACTTGCTGGGCCTTGGTCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGCGTAGc                                                                                                          >  8:46477/1‑149 (MQ=255)
                                        gATGGCATGCAGGCCAGCCACCAGCGGCATGTACACCTGCACTTGCTGGGCCTTGGTCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGCGTAGc                                                                                                          >  1:283635/1‑149 (MQ=255)
                                               tGCAGGCCAGCCACCAGCGGCATGTACACCTGCACTTGCTGGGCCTTGGTCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGCGTAGCCAACCt                                                                                                    >  8:30808/1‑148 (MQ=255)
                                                 cAGGCCAGCCACCAGCGGCATGTACACCTGCACTTGCTGGGCCTTGGTCTTGAGCACCTTGCGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGCGTAGCCAACCTGa                                                                                                  <  7:203122/148‑1 (MQ=255)
                                                 cAGGCCAGCCACCAGCGGCATGTACACCTGCACTTGCTGGGCCTTGGTCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGCGTAGCCAACCTGa                                                                                                  <  7:30808/148‑1 (MQ=255)
                                                  aGGCCAGCCACCAGCGGCATGTACACCTGCACTTGCTGGGCCTTGGTCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGCGTAGCCAACCTGa                                                                                                  >  3:224474/1‑147 (MQ=255)
                                                  aGGCCAGCCACCAGCGGCATGTACACCTGCACTTGCTGGGCCTTGGTCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGCGTAGCCAACCTGATa                                                                                                <  4:313315/149‑1 (MQ=255)
                                                    gccagccaCCAGCGGCATGTACACCTGCACTTGCTGGGCCTTGGTCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGCGTAGCCAACCTGATAGt                                                                                              >  1:200396/1‑149 (MQ=255)
                                                      cagccaCCAGCGGCATGTACACCTGCACTTGCTGGGCCTTGGTCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGCGTAGCCAACCTGATAGTTg                                                                                            <  1:244143/149‑1 (MQ=255)
                                                        gccaCCAGCGGCATGTACACCTGCACTTGCTGGGCCTTGGTCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGCGTAGCCAACCTGATAGTTGCg                                                                                          <  2:368332/149‑1 (MQ=255)
                                                        gccaCCAGCGGCATGTACACCTGCACTTGCTGGGCCTTGGTCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGCGTAGCCAACCTGATAGTTGCg                                                                                          <  3:365504/149‑1 (MQ=255)
                                                         ccacccGCGGCAGGTACACCTGCACTTGCTGGGCCTTGGGCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGCGTAGCCAACCTGATAGTTGCgg                                                                                         <  4:312241/149‑1 (MQ=255)
                                                              aGCGTCATGTACACCTGCACTTGCTGGGCCTTGGTCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGCGTAGCCAACCTGAAAGTTGCGGCTCtt                                                                                    <  8:56337/149‑1 (MQ=255)
                                                               gCGGCATGTACACCTGCACTTGCTGGGCCTTGGTCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGCGTAGCCAACCTGATAGTTGCGGCTCTTg                                                                                   >  5:88343/1‑149 (MQ=255)
                                                                         cacCTGCACTTGCTGGGCCTTGGTCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGTCCCTGCCCCAGTGCGTAGCCAACCTGATAGTTGCGGCTCTTGGGCGAGGAGc                                                                         <  6:251561/149‑1 (MQ=255)
                                                                         cacCTGCACTTGCTGGGCCTTGGTCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGCGTAGCCAACCTGATAGTTGCGGCTCTTGGGCgaggag                                                                          >  5:123154/1‑148 (MQ=255)
                                                                          acCTGCACTTGCTGGGCCTTGGTCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTcc                                                                                                                                      >  2:259062/1‑87 (MQ=255)
                                                                          acCTGCACTTGCTGGGCCTTGGTCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTcc                                                                                                                                      <  1:259062/87‑1 (MQ=255)
                                                                          acCTGCACTTGCTGGGCCTTGGTCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGCGTAGCCAACCTGATAGTTGCGGCTCTTGGGCgaggag                                                                          >  7:16935/1‑147 (MQ=255)
                                                                           ccTGCACTTGCTGGGCCTTGGTCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGCGTAGCCAACCTGATAGTTGCGGCTCTTGGGCGAGGAGCAg                                                                       <  6:66231/149‑1 (MQ=255)
                                                                           ccTGCACTTGCTGGGCCTTGGTCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGCGTAGCCAACCTGATAGTTGCGGCTCTTGGGCGAGGAGCAg                                                                       >  3:109859/1‑149 (MQ=255)
                                                                            cTGCACTTGCTGGGCCTTGGTCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGCGTAGCCAACCTGATAGTTGCGGCTCTTGGGCGAGGAGCAgg                                                                      >  6:349691/1‑149 (MQ=255)
                                                                                cACTTGCTGGGCCTTGGTCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGCGTAGCCAACCTGATAGTTGCGGCTCTTGGGCGAGGAGCAGGt                                                                     <  6:88343/147‑1 (MQ=255)
                                                                                  cTTGCTGGGCCTTGGTCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCAGGCTTTGGCCCTGCCCCAGTGCGTAGCCAACCTGATAGTTGCGGCTCTTGGGCGAGGAGCAGGTGACGa                                                                <  5:98342/149‑1 (MQ=255)
                                                                                      cTGGGCCTTGGTCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGCGTAGCCAACCTGATAGTTGCGGCTCTTGGGCGAGGAGCAGGTGACGATc                                                              <  6:52880/148‑1 (MQ=255)
                                                                                      cTGGGCCTTGGTCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGCGTAGCCAACCTGATAGTTGCGGCTCTTGGGCGAGGAGCAGGTGACGATc                                                              <  7:102908/148‑1 (MQ=255)
                                                                                        gggCCTTGGTCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGCGTAGCCAACCTGATAGTTGCGGCTCTTGGGCGAGGAGCAGGTGACGATCAgg                                                           >  3:27840/1‑149 (MQ=255)
                                                                                         ggCCTTGGTCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGCCCCATTGCGTAGCCAACCTGATAGTTGCGGCTCTTGGGCGAGGAGCAGGTGACGAt                                                               >  6:282078/1‑144 (MQ=255)
                                                                                         ggCCTTGGTCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGCGTAGCCAACCTGATAGTTGCTGCTCTTGGGCGAGGAGCAGGTGACGAt                                                               <  5:282078/144‑1 (MQ=255)
                                                                                           ccTTGGTCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGCGTAGCCAACCTGATAGTTGCGGCTCTTGGGCGAGGAGCAGGTGACGATCAGGTCa                                                        >  3:211280/1‑149 (MQ=255)
                                                                                             ttGGTCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGCGTAGCCAACCTGATAGTTGCGGCTCTTGGGCGAGGAGCAGGTGACGATCAGGTCAcc                                                      >  5:141304/1‑149 (MQ=255)
                                                                                               ggTCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGCGTAGCCAACCTGATAGTTGCGGCTCTTGGGCGAGGAGCAGGTGACGATCAGGTCAcc                                                      <  7:46477/147‑1 (MQ=255)
                                                                                               ggTCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGCGTAGCCAACCTGATAGTTGCGGCTCTTGGGCGAGGAGCAGGTGACGATCAGGTCACCGa                                                    >  6:216030/1‑149 (MQ=255)
                                                                                                gTCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGCGTAGCCAACCTGATAGTTGCGGCTCTTGGGCGAGGAGCAGGTGACGATCAGGTCACCGa                                                    >  8:130446/1‑148 (MQ=255)
                                                                                                     gAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGCGTAGCCAACCTGATAGTTGCGGCTCTTGGGCGAGGAGCAGGTGACGATCAGGTCACCGacac                                                 <  6:123154/146‑1 (MQ=255)
                                                                                                       gCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGCGTAGCCAACCTGATAGTTGCGGCTCTTGGGCGAGGAGCAGGTGACGATCAGGTCACCGACACCGGc                                             >  6:49710/1‑148 (MQ=255)
                                                                                                       gCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGCGTAGCCAACCTGATAGTTGCGGCTCTAGGGCGAGGAGCAGGTGACGATCAGGTCACCGACACCGGcc                                            <  2:283635/149‑1 (MQ=255)
                                                                                                            ttGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGCGTAGCCAACCTGATAGTTGCGGCTCTTGGGCGAGGAGCAGGTGACGATCAGGTCACCGACACCGGCCAGGcc                                       >  6:381432/1‑149 (MQ=255)
                                                                                                                       gACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGCGTAGCCAACCTGATAGTTGCGGCTCTTGGGCGAGGAGCAGGTGACGATCAGGTCACCGACACCGGCCAGGCCGAGGAACGTc                             <  6:141304/148‑1 (MQ=255)
                                                                                                                       gACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGCGTAGCCAACCTGATAGTTGCGGCTCTTGGGCGAGGAGCAGGTGACGATCAGGTCACCGACACCGGCCAGGCCGAGGAACGTCa                            >  6:73064/1‑149 (MQ=255)
                                                                                                                          ccccTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGCGTAGCCAACCTGATAGTTGCGGCTCTTGGGCGAGGAGCAGGTGACGATCAGGTCACCGACACCGGCCAGGCCGAGGAACGTCATgg                         >  5:337164/1‑149 (MQ=255)
                                                                                                                          ccccTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGCGTAGCCAACCTGATAGTTGCGGCTCTTGGGCGAGGAGCAGGTGACGATCAGGTCACCGACACCGGCCAGGCCGAGGAACGTCATgg                         >  6:63794/1‑149 (MQ=255)
                                                                                                                            ccTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGCGTAGCCAACCTGATAGTTGCGGCTCTTGGGCGAGGAGCAGGTGACGATCAGGTCACCGACACCGGCCAGGCCGAGGAACGTCATgggg                       >  3:337669/1‑149 (MQ=255)
                                                                                                                            ccTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGCGTAGCCAACCTGATAGTTGCGGCTCTTGGGCGAGGAGCAGGTGACGATCAGGTCACCGACACCGGCCAGGCCGAGGAACGTCATgggg                       >  1:387599/1‑149 (MQ=255)
                                                                                                                            ccTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGCGTAGCCAACCTGATAGTTGCGGCTCTTGGGCGAGGAGCAGGTGACGATCAGGTCACCGACACCGGCCAGGCCGAGGAACGTCATgggg                       >  6:329039/1‑149 (MQ=255)
                                                                                                                                ggCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGCGTAGCCAACCTGATAGTTGCGGCTCTTGGGCGAGGAGCAGGTGACGATCAGGTCACCGACACCGGCCAGGCCGAGGAACGTCATGGGGTTg                    <  2:387599/148‑1 (MQ=255)
                                                                                                                                 gCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGCGTAGCCAACCTGATAGTTGCGGCTCTTGGGCGAGGAGCAGGTGACGATCAGGTCACCGACACCGGc                                             >  7:203699/1‑122 (MQ=255)
                                                                                                                                 gCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGCGTAGCCAACCTGATAGTTGCGGCTCTTGGGCGAGGAGCAGGTGACGATCAGGTCACCGACACCGGc                                             <  8:203699/122‑1 (MQ=255)
                                                                                                                                      ctggaagaGGCGGTCAACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGCGTAGCCAACCTGATAGTTGCGGCTCTTGGGCGAGGAGCAGGTGACGATCAGGTCACCGACACCGGCCAGGCCGAGGAACGTCATGGGGTTGGCGcc               <  3:344100/140‑1 (MQ=255)
                                                                                                                                       tggaagaGGCGGTCAACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGCGTAGCCAACCTGATAGTTGCGGCTCTTGGGCGAGGAGCAGGTGACGATCAGGTCACCGACACCGGCCAGGCCGAGGAACGTCATGGGGTTGGCGCCCTg            >  6:44906/7‑149 (MQ=255)
                                                                                                                                         gCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGCGTAGCCAACCTGATAGTTGCGGCTCTTGGGCGAGGAGCAGGTGACGATCAGGTCACCGACACCGGCCAGGCCGAGGAACGTCATGGGGTTGGCGCCCTgg           <  2:193734/148‑1 (MQ=255)
                                                                                                                                         gCCCAGGCGGGTGACCGCCTCTTCCAGGCTGTGGCCCTGCCCCAGTGCGTAGCCAACCTGATAGTTGCGGCGCTTGGGCGAGTAGCAGGTGACGATCAGGTCACCGACACCGGCCAGGCCGAGGAACGTCATGGGGTTGGCGCCCTgg           <  4:224474/148‑1 (MQ=255)
                                                                                                                                              ggCGGTCAACCGCCTCTTCCAGGCTTTGGCCCTGCGCCAGTGCGTAGCCAACCTGATAGTTGCGGCTCTTGGGCGAGGAGCAGGTGACGATCAGGTCACCGACACCGGCCAGGCCGAGGAACGTCATGGGGTTGGCGCCCTGGCTGACg     >  3:142470/1‑149 (MQ=255)
                                                                                                                                               gCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGCGTAGCCAACCTGATAGTTGCGGCTCTTGGGCGAGGAGCAGGTGACGATCAGGTCACCGACACCGGCCAGGCCGAGGAACGTCATGGGGTTGGCGCCCTGGCTGACgg    <  4:27840/149‑1 (MQ=255)
                                                                                                                                                 ggTTGACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGCGTAGCCAACCTGATAGTTGCGGCTCTTGGGCGAGGAGCAGGTGACGATCAGGTCACCGACACCGGCCAGGCCGAGGAACGTCATGGGGTTGGCGCCCTGGCTGACGGc   >  2:110124/1‑148 (MQ=255)
                                                                                                                                                  gTTGACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGCGTAGCCAACCTGATAGTTGCGGCTCTTGGGCGAGGAGCAGGTGACGATCAGGTCACCGACACCGGCCAGGCCGAGGAACGTCATGGGGTTGGCGCCCTGGCTGACGGCg  >  6:148944/1‑148 (MQ=255)
                                                                                                                                                  gTTGACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGCGTAGCCAACCTGATAGTTGCGGCTCTTGGGCGAGGAGCAGGTGACGATCAGGTCACCGACACCGGCCAGGCCGAGGAACGTCATGGGGTTGGCGCCCTGGCTGACGGCg  >  3:381112/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
GGTGCTCGATCACCTGGTTCAGCGTGCGCCCTTCGAACAGGATGGCATGCAGGCCAGCCACCAGCGGCATGTACACCTGCACTTGCTGGGCCTTGGTCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGCGTAGCCAACCTGATAGTTGCGGCTCTTGGGCGAGGAGCAGGTGACGATCAGGTCACCGACACCGGCCAGGCCGAGGAACGTCATGGGGTTGGCGCCCTGGCTGACGGCG  >  NC_002947/4710575‑4710867

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: