Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,224,501 C→A 8.4% intergenic (‑613/+208) PP_3703 ← / ← PP_3704 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,224,5010CA 8.4% 166.1 / 6.3 60intergenic (‑613/+208)PP_3703/PP_3704hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (29/26);  new base A (3/2);  total (32/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 = 4.17e-01

GCCCTATCAATCCAAACGTAAATGTCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATACAGAACCTTAGGATTGTCCATGTAAACCTGGAAAAGCGCCAAGGAAGGTTCTGTGGAGAAGGGAT  >  NC_002947/4224356‑4224649
                                                                                                                                                   |                                                                                                                                                    
gCCCTATCAATCCAAACGTAAATGTCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCtgt                                                                                                                                                   >  6:274134/1‑149 (MQ=255)
  ccTATCAATCCAAACGTAAATGTCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCAGTAACAAAGACGGCCAGATAAGTGATAAGCGGAGCGATGCGGGAAAACCCCGAGGAAAATAAACATAGGTGGC‑‑CCAGCAGTGCACCtgtg                                                                                                                                                  >  5:376390/1‑148 (MQ=255)
   cTATCAATCCAAACGTAAATGTCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCa                                                                                                                                                <  7:99428/149‑1 (MQ=255)
           ccAAACGTAAATGTCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGt                                                                                                                                        >  5:380035/1‑149 (MQ=255)
            cacaCGTAAATGTCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGtt                                                                                                                                       <  4:149606/146‑1 (MQ=255)
             aaaCGTAAATGTCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTTGGTGGCAACCAGCAATGCACAGGTGCACTGCTggg                                                                                                                                        <  1:56768/149‑2 (MQ=255)
             aaaCGTAAATGTCCCTAAGCCTGGAGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTg                                                                                                                                      >  5:341750/1‑149 (MQ=255)
              aaCGTAAATGTCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGc                                                                                                                                     <  6:341750/149‑1 (MQ=255)
                    aaTGTCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACgtg                                                                                                                               <  6:380035/149‑1 (MQ=255)
                        tCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTaaa                                                                                                                           <  2:343809/149‑1 (MQ=255)
                        tCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTaaa                                                                                                                           <  8:148637/149‑1 (MQ=255)
                        tCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTa                                                                                                                             >  8:219496/1‑147 (MQ=255)
                          ccTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAgg                                                                                                                         >  2:314578/1‑149 (MQ=255)
                          ccTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAgg                                                                                                                         >  1:12876/1‑149 (MQ=255)
                           cTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGc                                                                                                                        >  7:75532/1‑149 (MQ=255)
                                ccTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTAt                                                                                                                    >  6:23403/1‑148 (MQ=255)
                                ccTGGCGGCGGGGAACAACCACGCTCACGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTCGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTAt                                                                                                                    >  7:207682/1‑148 (MQ=255)
                                         ggggAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGa                                                                                                          >  8:143083/1‑149 (MQ=255)
                                         ggggAACAACCACGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTCGTTGCCACGTGTAAAGGCTTATGAGCTAAAGa                                                                                                          >  4:269513/1‑149 (MQ=255)
                                          gggAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCGCCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGAc                                                                                                         <  7:104219/149‑1 (MQ=255)
                                            gAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAg                                                                                                       <  6:215066/149‑1 (MQ=255)
                                                aacCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACt                                                                                                   >  8:160251/1‑149 (MQ=255)
                                                 cccTCGCTCCCGTAACAAAGACGGCCAGATTCGTGATAAGCGGAGCGATGCGGGAAAACCGCGCGGAAAATAAACCTAGGGGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGTCTg                                                                                                  <  8:207682/148‑1 (MQ=255)
                                                  ccTCTCTCCATTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTgg                                                                                                 <  3:180631/149‑1 (MQ=255)
                                                  ccTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTgg                                                                                                 >  8:41941/1‑149 (MQ=255)
                                                    tCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTgg                                                                                                 >  2:220209/1‑147 (MQ=255)
                                                    tCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTgg                                                                                                 <  1:220209/147‑1 (MQ=255)
                                                     cGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTa                                                                                              >  6:294043/1‑149 (MQ=255)
                                                           cGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCgcgc                                                                                        <  1:314578/149‑1 (MQ=255)
                                                            gTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGcc                                                                                       >  4:146132/1‑149 (MQ=255)
                                                               aCAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGa                                                                                                                  >  6:165708/1‑119 (MQ=255)
                                                               aCAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCa                                                                                     <  7:143083/148‑1 (MQ=255)
                                                                cAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATc                                                                                   >  6:297122/1‑149 (MQ=255)
                                                                cAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAATGCACAGGTGCActgctg                                                                                                                                          <  4:104799/94‑1 (MQ=39)
                                                                cAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAATGCACAGGTGCActgctg                                                                                                                                          >  3:104799/1‑94 (MQ=39)
                                                                 aaaGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCt                                                                                  <  5:23403/149‑1 (MQ=255)
                                                                  aaGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTACAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTg                                                                                 >  4:241675/1‑149 (MQ=255)
                                                                       ggCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCGCATCTGTAGCg                                                                            >  4:203990/1‑149 (MQ=255)
                                                                          cAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCg                                                                            >  5:81418/1‑146 (MQ=255)
                                                                            gATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAgat                                                                       <  5:294043/149‑1 (MQ=255)
                                                                                cGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAgatgat                                                                    <  3:241675/148‑1 (MQ=255)
                                                                                 gTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATAc                                                                  <  8:75532/149‑1 (MQ=255)
                                                                                   gATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATa                                                                   >  7:112477/1‑146 (MQ=255)
                                                                                          ggAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATACAGAACCTTa                                                         <  7:160251/149‑1 (MQ=255)
                                                                                           gAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATACAGAACCTTAg                                                        <  5:274134/149‑1 (MQ=255)
                                                                                                     ggAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATACAGAACCTTAGGATTGTCCAt                                              <  7:41941/149‑1 (MQ=255)
                                                                                                            cccGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAATACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATACAGAACCTTAGGATTGTCCATGTAATcc                                       >  8:224194/1‑149 (MQ=255)
                                                                                                                 ggAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATACAGAACCTTAGGATTGTCCATGTAAACCTGGaa                                  <  8:112477/149‑1 (MQ=255)
                                                                                                                  gAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATACAGAACCTTAGGATTGTCCATGTAAACCTGGaa                                  >  2:6283/1‑148 (MQ=255)
                                                                                                                     aaTAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATACAGAACCTTAGGATTGTCCATGTAAACCTGGAAAAgc                              <  1:357449/149‑1 (MQ=255)
                                                                                                                     aaTAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATACAGAACCTTAGGATTGTCCATGTAAACCTGGAAAAgc                              >  2:106714/1‑149 (MQ=255)
                                                                                                                            cTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATACAGAACCTTAGGATTGTCCATGTAAACCTGGAAAAGCGCCaag                        >  7:111535/1‑148 (MQ=255)
                                                                                                                              aGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCAACTGTAGCGGAGATGATACAGAACCTTAGGATTGTCCATGTAAACCTGGAAAAGCGCCaagg                       >  3:5171/1‑147 (MQ=255)
                                                                                                                                 tggcaaccAGCAATGCACAGGTGCACTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATACAGAACCTTAGGATTGTCCATGTAAACCTGGAAAAGCGCCaaggaag                    >  1:56796/7‑149 (MQ=255)
                                                                                                                                   gC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACATGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATACAGAACCTTAGGATTGTCCATGTAAACCTGGAAAAGCGCCAAGGAAGGtt                 <  8:217880/148‑1 (MQ=255)
                                                                                                                                        cagcagTGCACCTGTGCATTGGTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATACAGAACCTTAGGATTGTCCATGTAAACCTGGAAAAGCGCCAAGGAAGGTTCtgt             <  3:269513/149‑1 (MQ=255)
                                                                                                                                        cagcaaTGCACAGGTGCACTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATACAGAACCTTAGGATTGTCCATGTAAACCTGGAAAAGCGCCAAGGAAGGTTCtgt             >  5:98112/1‑149 (MQ=255)
                                                                                                                                         ccgatctcACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATACAGAACCTTAGGATTGTCCATGTAAACCTGGAAAAGCGCCAAGGAAGGTTCtgt             <  3:106509/141‑1 (MQ=255)
                                                                                                                                                cACCTGTGCTTTGCTGGTTGCCACGTGTAAAGGATTATTAGCTAAAGACAGGACTGTGTATCGCGCCCATCTGTAGCGGAGATGATACAGAACCTTAGGATTGTCCATGTAAACCTTGAAAATCGCCAAGGAAGGTTCTGTGGAGAAg      <  6:376390/148‑1 (MQ=255)
                                                                                                                                                 aCCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATACAGAACCTTAGGATTGTCCATGTAAACCTGGAAAAGCGCCAAGGAAGGTTCtgt             >  4:106509/1‑140 (MQ=255)
                                                                                                                                                   cTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATACAGAACCTTAGGATTGTCCATGTAAACCTGGAAAAGCGCCAAGGAAGGTTCTGTGGAGAAGGGat  <  6:81418/149‑1 (MQ=255)
                                                                                                                                                   cTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATACAGAACCTTAGGATTGTCCATGTAAACCTGGAAAAGCGCCAAGGAAGGTTCTGTGGAGAAGGGa   >  5:217061/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GCCCTATCAATCCAAACGTAAATGTCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATACAGAACCTTAGGATTGTCCATGTAAACCTGGAAAAGCGCCAAGGAAGGTTCTGTGGAGAAGGGAT  >  NC_002947/4224356‑4224649

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑