Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,563,133 C→A 5.4% R243L (CGG→CTG)  vgrG‑III ← type VI secretion system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,563,1330CA 5.4% 239.8 / 3.4 74R243L (CGG→CTG) vgrG‑IIItype VI secretion system protein
Reads supporting (aligned to +/- strand):  ref base C (36/34);  new base A (2/2);  total (38/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.01e-01

AGGGCTTGCCCACCGCATCGCGCTTGTAGCGCCCCGGGTAATCAAAGCGCTCGTAGTCCAGTGACTGGTGCGTGACGCCACTGGTCTGGGCCTTGTGCTCCTGCCGGTAGGCCGGGTTGGTGAAGGTGTAGTCGCGCTGCACTTGCCGCGCAGTACGTACCTGCTCGCAATAGCGCAAGCGCCGCAGGCATGCGCGGGCCTGGTCGCCACCACTGCTGGCGTGGTACAGCACCTTGTCGGGGCCGACATCATCCTCATCCAGCCAAAGCGCCTCGTCGTCCTCTGCCTT  >  NC_002947/4562987‑4563275
                                                                                                                                                  |                                                                                                                                              
aGGGCTTGCCCACCGCATCGCGCTTGTAGCGCCCCGGGTAATCAAAGCGCTCGTAGTCCAGTGACTGGTGCGTGACGCCACTGGTCTGGGCCTTGTGCTCCTGCCGGTAGGCCGGGTTGGTGAAGGTGTAGTCGCGCTGCACTTGCcg                                                                                                                                               <  4:236653/148‑1 (MQ=255)
   gCTTGCCCACCGCATCGCGCTTGTAGCGCCCCGGGTAATCAAAGCGCTCGTAGTCCAGTGACTGGTGCGTGACGCCACTGGTCTGGGCCTTGTGCTCCTGCCGGTAGGCCGGGTTGGTGAAGGTGTAGTCGCGCTGCACTTGCCGCGCa                                                                                                                                           <  4:186768/149‑1 (MQ=255)
       gCCCACCGCATCGCGCTTGTAGCGCCCCGGGTAATCAAAGCGCTCGTAGTCCAGTGACTGGTGCGTGACGCCACTGGTCTGGGCCTTGTGCTCCTGCCGGTAGGCCGGGTTGGTGAAGGTGTAGTCGCGCTGCACTTGCCGCGCAgt                                                                                                                                         <  5:171150/147‑1 (MQ=255)
        cccACCGCATCGCGCTTGTAGCGCCCCGGGTAATCAAAGCGCTCGTAGTCCAGTGACTGGTGCGTGACGCCACTGGTCTGGGCCTTGTGCTCCTGCCGGTAGGCCGGGTTGGTGAAGGTGTAGTCGCGCGGCAAGTGCAgcgcgacta                                                                                                                                       <  4:257870/148‑6 (MQ=39)
        cccACCGCATCGCGCTTGTAGCGCCCCGGGTAATCAAAGCGCTCGTAGTCCAGTGACTGGTGCGTGACGCCACGGGTCTGGGCCTTGTGCTCCTGCCGGTAGGCCGGGTCGGTGAAGGTGTAGTCGCGCTGCACTTGCCGCGCAgtacg                                                                                                                                      >  1:78298/1‑149 (MQ=255)
         ccACCGCATCGCGCTTGTAGCGCCCCGGGTAATCAAAGCGCTCGTAGTCCAGTGACTGGTGCGTGACGCCACTGGTCTGGGCCTTGTGCTCCTGCCGGTAGGCCGGGTTGGTGAAGGTGTAGTCGCGCTGCACTTGCCGCGCa                                                                                                                                           <  2:357443/143‑1 (MQ=255)
         ccACCGCATCGCGCTTGTAGCGCCCCGGGTAATCAAAGCGCTCGTAGTCCAGTGACTGGTGCGTGACGCCACTGGTCTGGGCCTTGTGCTCCTGCCGGTAGGCCGGGTTGGTGAAGGTGTAGTCGCGCTGCACTTGCCGCGCa                                                                                                                                           >  1:357443/1‑143 (MQ=255)
          cACCGCATCGCGCTTGTAGCGCCCCGGGTAATCAAAGCGCTCGTAGTCCAGTGACTGGTGCGTGACGCCACTGGTCTGGGCCTTGTGCTCCTGCCGGTAGGCCGGGTTGGTGAAGGTGTAGTCGCGCTGCACTTGCCGCGCAgtacgta                                                                                                                                    <  2:31954/149‑1 (MQ=255)
           aCCGCATCGCGCTTGTAGCGCCCCGGGTAATCAAAGCGCTCGTAGTCCAGTGACTGGTGCGTGACGCCACTGGTCTGGGCCTTGTGCTCCTGCCGGTAGGCCGGGTTGGTGAAGGTGTAGTCGCGCTGCACTTGCCGCGCAgtacgtac                                                                                                                                   >  8:186352/1‑149 (MQ=255)
               cATCGCGCTTGTAGCGCCCCGGGTAATCAAAGCGCTCGTAGTCCAGTGACTGGTGCGTGACGCCACTGGTCTGGGCCTTGTGCTCCTGCCGGTAGGCCGGGTTGGTGAAGGTGTAGTCGCGCTGCACTTGCCGCGCAGTACGTACCTGc                                                                                                                               >  1:212247/1‑149 (MQ=255)
               cATCGCGCTTGTAGCGCCCCGGGTAATCAAAGCGCTCGTAGTCCAGTGACTGGTGCGTGACGCCACTGGTCTGGGCCTTGTGCTCCTGCCGGTAGGCCGGGTTGGTGAAGGTGTAGTCGCGCTGCACTTGCCGCGCAGTACGTACCTGc                                                                                                                               >  7:135511/1‑149 (MQ=255)
                     gcTTGTAGCGCCCCGGGTAATCAAAGCGCTCGTAGTCCAGTGACTGGTGCGTGACGCCACTGGTCTGGGCCTTGTGCTCCTGCCGGTAGGCCGGGTTGGTGAAGGTGTAGTCGCGCTGCACTTGCCGCGCAGTACGTACCTGCTCGCaa                                                                                                                         >  5:348780/1‑149 (MQ=255)
                            gcgcCCCGGGTAATCAAAGCGCTCGTAGTCCAGTGACTGGTGCGTGACGCCACTGGTCTGGGCCTTGTGCTCCTGCCGGTAGGCCGGGTTGGGGAAGGTGTAGTCGCGCTGCACTTGCCGCGCAGTACGTACCTGCTCGCAATAgcgc                                                                                                                   <  6:321738/148‑1 (MQ=255)
                              gcCCCGGGTAATCAAAGCGCTCGTAGTCCAGTGACTGGTGCGTGACGCCACTGGTCTGGGCCTTGTGCTCCTGCCGGTAGGCCGGGTTGGTGAAGGTGTAGTCGCGCTGCACTTGCCGCGCAGTACGTACCTGCTCGCAATAGCGCAAg                                                                                                                <  6:342902/149‑1 (MQ=255)
                                 ccGGGTAATCAAAGCGCTCGTAGTCCAGTGACTGGTGCGTGACGCCACTGGTCTGGGCCTTGTGCTCCTGCCGGTAGGCCGGGTTGGTGAAGGTGTAGTCGCGCTGCACTTGCCGCGCAGTACGTACCTGCTCGCAATAGCGCAAgcgc                                                                                                             >  8:110821/1‑149 (MQ=255)
                                     gTAATCAAAGCGCTCGTAGTCCAGTGACTGGTGCGTGACGCCACTGGTCTGGGCCTTGTGCTCCTGCCGGTAGGCCGGGTTGGTGAAGGTGTAGTCGCGCTGCACTTGCCGCGCAGTACGTACCTGCTCGCAATAGCGCAAGcgccg                                                                                                           >  5:181292/1‑147 (MQ=255)
                                      tAATCAAAGCGCTCGTAGTCCAGTGACTGGTGCGTGACGCCACTGGTCTGGGCCTTGTGCTCCTGCCGGTAGGCCGGGTTGGTGCAGGTGTAGTCGCGCTGCACTTGCCGCGCAGTACGTACCTGCTCGCAATAGCGCAAGCGCCGCAg                                                                                                        <  6:146300/149‑1 (MQ=255)
                                      tAATCAAAGCGCTCGTAGTCCAGTGACTGGTGCGTGACGCCACTGGTCTGGGCCTTGTGCTCCTGCCGGTAGGCCGGGTTGGTGAAGGTGTAGTCGCGCTGCACTTGCCGCGCAGTACGTACCTGCTCGCAATAGCGCAAGCGCCGCAg                                                                                                        <  5:150572/149‑1 (MQ=255)
                                        aTCAAAGCGCTCGTAGTCCAGTGACTGGTGCGTGACGCCACTGGTCTGGGCCTTGTGCTCCTGCCGGTAGGCCGGGTTGGTGAAGGTGTAGTCGCGCTGCACTTGCCGCGCAGTACGTACCTGCTCGCAATAGCGCAAGCGCCGCAGGc                                                                                                      <  6:348780/149‑1 (MQ=255)
                                          cAAAGCGCTCGTAGTCCAGTGACTGGTGCGTGACGCCACTGGTCTGGGCCTTGTGCTCCTGCCGGTAGGCCGGGTTGGTGAAGGTGTAGTCGCGCTGCACTTGCCGCGCAGTACATACCTGCTCGCAATAGCGCAAGCGCCGCAGGCAt                                                                                                    >  3:332007/1‑149 (MQ=255)
                                            aaGCGCTCGTAGTCCAGTGACTGGTGCGTGACGCCACTGGTCTGGGCCTTGTGCTCCTGCCGGTAGGCCGGGTTGGTGAAGGTGTAGTCGCGCTGCACTTGCCGCGCAGTACGTACCTGCTCGCAATAGCGCAAGCGCCGCAGGCATgc                                                                                                  >  7:217196/1‑149 (MQ=255)
                                                gcTCGTAGTCCAGTGACTGGTGCGTGACGCCACTGGTCTGGGCCTTGTGCTCCTGCCGGTAGGCCGGGTTGGTGAAGGTGTAGTCGCGCTGCACTTGCCGCGCAGTACGTACCTGCTCGCAATAGCGCAAGCGCCGCAGGCATGCGCgg                                                                                              <  1:80276/149‑1 (MQ=255)
                                                gcTCGTAGTCCAGTGACTGGTGCGTGACGCCACTGGTCTGGGCCTTGTGCTCCTGCCGGTAGGCCGGGTTGGTGAAGGTGTAGTCGCGCTGCACTTGCCGCGCAGTACGTACCTGCTCGCAATAGCGCAAGCGCCGCAGGCATGCGCgg                                                                                              <  5:366499/149‑1 (MQ=255)
                                                gcTCGTAGTCCAGTGACTGGTGCGTGACGCCACTGGTCTGGGCCTTGTGCTCCTGCCGGTAGGCCGGGCTGGTGAAGGTGTAGTCGCGCTGCACTTGCCGCGCAGTACGTACCTGCGCGCAATAGCGCAAGCGCCGCAGGCATGCGCgg                                                                                              <  6:317863/149‑1 (MQ=255)
                                                           aGTGACTGGTGCGTGACGCCACTGGTCTGGGCCTTGTGCTCCTGCCGGTAGGCCGGGTTGGTGAAGGTGTAGTCGCGCTGCACTTGCCGCGCAGTACGTACCTGCTCGCAATAGCGCAAGCGCCGCAGGCATGCGCGGGCCTGGTc                                                                                      <  8:224137/146‑1 (MQ=255)
                                                           aGTGACTGGTGCGTGACGCCACTGGTCTGGGCCTTGTGCTCCTGCCGGTAGGCCGGGTTGGTGAAGGTGTAGTCGCGCTGCACTTGCCGCGCAGTACGTACCTGCTCGCAATAGCGCAAGCGCCGCAGGCATGCGCGGGCCTGGTc                                                                                      >  7:224137/1‑146 (MQ=255)
                                                               aCTGGTGCGTGACGCCACTGGTCTGGGCCTTGTGCTCCTGCCGGTAGGCCGGGTTGGTGAAGGTGTAGTCGCGCTGCACTTGCCGCGCAGTACGTACCTGCTCGCAATAGCGCAAGCGCCGCAGGCATGCGCGGGCCTGGTCGccacca                                                                               <  7:186352/149‑1 (MQ=255)
                                                                cTGGTGCGTGACGCCACTGGTCTGGGCCTTGTGCTCCTGCCGGTAGGCCGGGTTGGTGAAGGTGTAGTCGCGCTGCACTTGCCGCGCAGTACGTACCTGCTCGCAATAGCGCAAGCGCCGCACGCATGCGCGGGCCTGGTCGCcaccac                                                                              <  2:212247/149‑1 (MQ=255)
                                                                 tGGTGCGTGACGCCACTGGTCTGGGCCTTGTGCTCCTGCCGGTAGGCCGGGTTGGTGAAGGTGTAGTCGCGCTGCACTTGCCGCGCAGTACGTACCTGCTCGCAATAGCGCAAGCGCCGCAGGCATGCGCGGGCCTGGTc                                                                                      >  1:149895/1‑140 (MQ=255)
                                                                 tGGTGCGTGACGCCACTGGTCTGGGCCTTGTGCTCCTGCCGGTAGGCCGGGTTGGTGAAGGTGTAGTCGCGCTGCACTTGCCGCGCAGTACGTACCTGCTCGCAATAGCGCAAGCGCCGCAGGCATGCGCGGGCCTGGTc                                                                                      <  2:149895/140‑1 (MQ=255)
                                                                    tGCGTGACGCCACTGGTCTGGGCCTTGTGCTCCTGCCGGTAGGCCGGGTTGGTGAAGGTGTAGTCGCGCTGCACTTGCCGCGCAGTACGTACCTGCTCGCAATAGCGCAAGCGCCGCAGGCATGCGCGGGCCTGGTCGCCACCActgct                                                                          >  4:34919/1‑149 (MQ=255)
                                                                      cGTGACGCCACTGGTCTGGGCCTTGTGCTCCTGACGGTAGGCCGGGTTGGTGAAGGTGTAGTCGCGCTGCACTTGCCGCGCAGTACGTACCTGCTCGCAATAGCGCAAGCGCCGCAGGCATGCGCGGGCCTGGTCGCCACCACTGCTgg                                                                        >  4:72080/1‑149 (MQ=255)
                                                                        tGACGCCACTGGTCTGGGCCTTGTGCTCCTGCCGGTAGGCCGGGTTGGTGAAGGTGTAGTCGCGCTGCACTTGCCGCGCAGTACGTACCTGCTCGCAATAGCGCAAGCGCCGCAGGCATGCGCGGGCCTGGTCGCCACCACTGCTGGCg                                                                      <  8:135511/149‑1 (MQ=255)
                                                                          aCGCCACTGGTCTGGGCCTTGTGCTCCTGCCGGTAGGCCGGGTTGGTGAAGGTGTAGTCGCGCTGCACTTGCCGCGCAGTACGTACCTGCTCGCAATAGCGCAAGCGCCGCAGGCATGCGCGGGCCTGGTCGCCACCACTGCTGGCGt                                                                     <  6:181292/148‑1 (MQ=255)
                                                                          aCGCCACTGGTCTGGGCCTTGTGCTCCTGCCGGTAGGCCGGGTTGGTGAAGGTGTAGTCGCGCTGCACTTGCCGCGCAGTACATACCTGCTCGCAATAGCGCAAGCGCCGCAGGCATGCGCGGGCCTGGTCGCCACCACTGCTGGCGTg                                                                    <  4:332007/149‑1 (MQ=255)
                                                                             ccACTGGTCTGGGCCTTGTGCTCCTGCCGGTAGGCCGGGTTGGTGAAGGTGTAGTCGCGCTGCACTTGCCGCGCAGTACGTACCTGCTCGCAATAGCGCAAGCGCCGCAGGCATGCGCGGGCCTGGTCGCCACCACTGCTGGCGTGGTa                                                                 <  7:110821/149‑1 (MQ=255)
                                                                                   gTCTGGGCCTTGTGCTCCTGCCGGTAGGCCGGGTTGGTGAAGGTGTAGTCGCGCTGCACTTGCCGCGCAGTACGTACCTGCTCGCAATAGCGCAAGCGCCGCAGGCATGCGCGGGCCTGGTCGCCACCACTGCTGGCGTGGTACAGCAc                                                           >  6:221505/1‑149 (MQ=255)
                                                                                       gggCCTTGTGCTCCTGCCGGTAGGCCGGGTTGGTGAAGGTGTAGTCGCGCTGCACTTGCCGCGCAGTACGTACCTGCTCGCAATAGCGCAAGCGCCGCAGGCATGCGCGGGCCTGGTCGCCACCACTGCTGGCGTGGTACAGCACCTTg                                                       <  3:272855/149‑1 (MQ=255)
                                                                                        ggCCTTGTGCTCCTGCCGGTAGGCCGGGTTGGTGAAGGTGTAGTCGCGCTGCACTTGCCGCGCAGTACGTACCTGCTCGCAATAGCGCAAGCGCCGCAGGCATGCGCGGGCCTGGTCGCCACCACTGCTGGCGTGGTACAGCACCTTGt                                                      >  3:1249/1‑149 (MQ=255)
                                                                                        ggCCTTGTGCTCCTGCCGGTAGGCCGGGTTGGTGAAGGTGTAGTCGCGCTGCACTTGCCGCGCAGTACGTACCTGCTCGCAATAGCGCAAGCGCCGCAGGCATGCGCGGGCCTGGTCGCCACCACTGCTGGCGTGGTACAGCACCTTGt                                                      >  5:200150/1‑149 (MQ=255)
                                                                                         gCCTTGTGCTCCTGCCGGTAGGCCGGGTTGGTGAAGGTGTAGTCGCGCTGCACTTGCCGCGCAGTACGTACCTGCTCGCAATAGCGCAAGCGCCGCAGGCATGCGCGGGCCTGGTCGCCACCACTGCTGGCGTGGTACAGCACCTTGTc                                                     <  3:34919/149‑1 (MQ=255)
                                                                                              gtgCTCCTGCCGGTAGGCCGGGTTGGTGAAGGTGTAGTCGCGCTGCACTTGCCGCGCAGTACGTACCTGCTCGCAATAGCGCAAGCGCCGCAGGCTTGCGCGGGCCTGGTCGCCACCACTGCTGGCGTGGTACAGc                                                             >  1:291638/1‑136 (MQ=255)
                                                                                              gtgCTCCTGCCGGTAGGCCGGGTTGGTGAAGGTGTAGTCGCGCTGCACTTGCCGCGCAGTACGTACCTGCTCGCAATAGCGCAAGCGCCGCAGGCATGCGCGGGCCTGGTCGCCACCACTGCTGGCGTGGTACAGc                                                             <  2:291638/136‑1 (MQ=255)
                                                                                                 cTCCTGCCGGTAGGCCGGGTTGGTGAAGGTGTAGTCGCGCTGCACTTGCCGCGCAGTACGTACCTGCTCGCAATAGCGCAAGCGCCGCAGGCATGCGCGGGCCTGGTCGCCACCACTGCTGGCGTGGTACAGCACCTTGTCGGGGCCGa                                             >  1:180283/1‑149 (MQ=255)
                                                                                                  tCCTGCCGGTAGGCCGGGTTGGTGAAGGTGTAGTCGCGCTGCACTTGCCGCGCAGTACGTACCTGCTCGCAATAGCGCAAGCGCCGCAGGCATGCGCGGGCCTGGTCGCCACCACTGCTGGCGTGGTACAGCACCTTGTCGGGGCCGa                                             >  4:336938/1‑148 (MQ=255)
                                                                                                    cTGCCGGTAGGCCGGGTTGGTGAAGGTGTAGTCGCGCTGCACTTGCCGCGCAGTACGTACCTGCTCGCAATAGCGCAAGCGCCGCAGGCATGCGCGGGCCTGGTCGCCACCACTGCTGGCGTGGTACAGCACCTTGTCGGGGCCGAcat                                          >  8:215997/1‑149 (MQ=255)
                                                                                                    cTGCCGGTAGGCCGGGTCGGTGAAGGTGTAGTCGCGCTGCACTTGCCGCGCAGTACGTACCTGCTCGCAATAGCGCAAGCGCCGCAGGCATGCGCGGGCCTGGTCGCCACCACTGCTGGCGTGGTACAGCACCTTGTCGGGGCCGAc                                            <  2:78298/147‑1 (MQ=255)
                                                                                                       ccGGTAGGCCGGGTTGGTGAAGGTGTAGTCGCGCTGCACTTGCCGCGCAGTACGTACCTGCTCGCAATAGCGCAAGCGCCGCAGGCATGCGCGGGCCTGGTCGCCACCACTGCTGGCGTGGTACAGCAc                                                           >  2:80679/1‑129 (MQ=255)
                                                                                                              gCCGGGTTGGTGAAGGTGTAGTCGCGCTGCACTTGCCGCGCAGTACGTACCTGCTCGCAATAGCGCAAGCGCCGCAGGCATGCGCGGGCCTGGTCGCCACCACTGCTGGCGTGGTACAGCACCTTGTCGGGGCCGAcatc                                         >  2:68045/1‑140 (MQ=255)
                                                                                                              gCCGGGTTGGTGAAGGTGTAGTCGCGCTGCACTTGCCGCGCAGTACGTACCTGCTCGCAATAGCGCAAGCGCCGCAGGCATGCGCGGGCCTGGTCGCCACCACTGCTGGCGTGGTACAGCACCTTGTCGGGGCCGAcatc                                         <  1:68045/140‑1 (MQ=255)
                                                                                                               ccGGGTTGGTGAAGGTGTAGTCGCGCTGCACTTGCCGCGCAGTACGTACCTGCTCGCAATAGCGCAAGCGCCGCAGGCATGCGCGGGCCTGGTCGCCACCACTGCTGGCGTGGTACAGCACCTTGTCGGGGCCGACAtcatcctcatc                                >  1:199930/1‑148 (MQ=255)
                                                                                                               ccGGGTTGGTGAAGGTGTAGTCGCGCTGCACTTGCCGCGCAGTACGTACCTGCTCGCAATAGCGCAAGCGCCGCAGGCATGCGCGGGCCTGGTCGCCACCACTGCTGGCGTGGTACAGCACCTTGTCGGGGCCGACAtcatcctcatc                                >  6:254587/1‑148 (MQ=255)
                                                                                                                cGGGTTGGTGAAGGTGTAGTCGCGCTGCACTTGCCGCGCAGTACGTACCTGCTCGCAATAGCGCAAGCGCCGCAGGCATGCGCGGGCCTGGTCGCCACCACTGCTGGCGTGGTACAGCACCTTGTCGGGGCCGACATCATCCTCATCCa                              <  6:48090/149‑1 (MQ=255)
                                                                                                                       gTGAAGGTGTAGTCGCGCTGCACTTGCCGCGCAGTACGTACCTGCTCGCAATAGCGCAAGCGCCGCAGGCATGCGCGGGCCTGGTCGCCACCACTGCTGGCGTGGTACAGCACCTTGTCggg                                                  >  1:157546/1‑122 (MQ=255)
                                                                                                                       gTGAAGGTGTAGTCGCGCTGCACTTGCCGCGCAGTACGTACCTGCTCGCAATAGCGCAAGCGCCGCAGGCATGCGCGGGCCTGGTCGCCACCACTGCTGGCGTGGTACAGCACCTTGTCggg                                                  <  2:157546/122‑1 (MQ=255)
                                                                                                                       gTGAAGGTGTAGTCGCGCTGCACTTGCCGCGCAGTACGTACCTGCTCGCAATAGCGCAAGCGCCGCAGGCATGCGCGGGCCTGGTCGCCACCACTGCTGGCGTGGTACAGCACCTTGTCGGGGCCGACATCATCCTCATCCAGCCaaa                        <  3:336938/148‑1 (MQ=255)
                                                                                                                       gTGAAGGTGTAGTCGCGCTGCACTTGCCGCGCAGTACGTACCTGCTCGCAATAGCGCAAGCGCCACAGGCATGCGCGGGCCTGGTCGCCACCACTGCTGGCGTGGTACAGCACCTTGTCGGGGCCGACATCATCCTCATCCAGCCAAAg                       >  6:146028/1‑149 (MQ=255)
                                                                                                                        tGAAGGTGTAGTCGCGCTGCACTTGCCGCGCAGTACGTACCTGCTCGCAATAGCGCAAGCGCCGCAGGCATGCGCGGGCCTGGTCGCCACCACTGCTGGCGTGGTACAGCACCTTGTCGGGGCCGACATCATCCTCATCCAGCCAAAgc                      <  2:4832/149‑1 (MQ=255)
                                                                                                                          aaGGTGTAGTCGCGCTGCACTTGCCGCGCAGTACGTACCTGCTCGCAATAGCGCAAGCGCCGCAGGCATGCGCGGGCCTGGTCGCCACCACTGCTGGCGTGGTACAGCACCTTGTCGGGGCCGACATCATCCTCATCCAGCCAAAgcgc                    >  7:99353/1‑149 (MQ=255)
                                                                                                                             gtgtAGTCGCGCTGCACTTGCCGCGCAGTACGTACCTGCTCGCAATAGCGCAAGCGCCGCAGGCATGCGCGGGCCTGGTCGCCACCACTGCTGGCGTTGTACAGCACCTTGTCGGGGCCGa                                             <  6:185941/121‑1 (MQ=255)
                                                                                                                             gtgtAGTCGCGCTGCACTTGCCGCGCAGTACGTACCTGCTCGCAATAGCGCAAGCGCCGCAGGCATGCGCGGGCCTGGTCGCCACCACTGCTGGCGTGGTACAGCACCTTGTCGGGGCCGa                                             >  5:185941/1‑121 (MQ=255)
                                                                                                                                 aGTCGCGCTGCACTTGCCGCGCAGTACGTACCTGCTCGCAATAGCGCAAGCGCCGCAGGCATGCGCGGGCCTGGTCGCCACCACTGCTGGCGTGGTACAGCACCTTGTCGGGGCCGACATCATCCTCATCCAGCCAAAGCGCCtcgtcg             >  5:258185/1‑149 (MQ=255)
                                                                                                                                  gTCGCGCTGCACTTGCCGCGCAGTACGTACCTGCTCGCAATAGCGCAAGCGCCGCAGGCATGCGCGGGCCTGGTCGCCACCACTGCTGGCGTGGTACAGCACCTTGTCGGGGCCGACatcatc                                      <  1:352027/123‑1 (MQ=255)
                                                                                                                                  gTCGCGCTGCACTTGCCGCGCAGTACGTACCTGCTCGCAATAGCGCAAGCGCCGCAGGCATGCGCGGGCCTGGTCGCCACCACTGCTGGCGTGGTACAGCACCTTGTCGGGGCCGACatcatc                                      >  2:352027/1‑123 (MQ=255)
                                                                                                                                  gTCGCGCTGCACTTGCCGCGCAGTACGTACCTGCTCGCAATAGCGCAAGCGCCGCAGGCATGCGCGGGCCTGGTCGCCACCACTGCTGGCGTGGTACAGCACCTTGTCGGGGCCGACATCATCCTCATCCAGCCAAAGCGCCtcgtcg             <  8:217196/148‑1 (MQ=255)
                                                                                                                                   ctgcgcGGCAAGTGCAGCGCAGTACGTACCTGCTCGCAATAGCGCAAGCGCCGCAGGCATGCGCGGGCCTGGTCGCCACCACTGCTGGCGTGGTACAGCACCTTGTCGGGGCCGACATCATCCTCATCCAGCCAAAGCGCCTCgtcgtc           >  8:310934/3‑149 (MQ=255)
                                                                                                                                   ctgcgcGGCAAGTGCAGCGCAGTACGTACCTGCTCGCAATAGCGCAAGCGCCGCAGGCATGCGCGGGCCTGGTCGCCACCACTGCTGGCGTGGTACAGCACCTTGTCGGGGCCGACATCATCCTCATCCAGCCAAAGCGCCTCgtcgtc           >  3:223261/3‑149 (MQ=255)
                                                                                                                                     gcgcTGCACTTGCCGCGCAGTACGTACCTGCTCGCAATAGCGCAAGCGCCGCAGGCATGCGCGGGCCTGGTCGCCACCACTGCTGGCGTGGTACAGCACCTTGTCGGGGCCGACATCATCCTCATCCAGCCAAAgc                      >  6:358090/1‑136 (MQ=255)
                                                                                                                                     gcgcTGCACTTGCCGCGCAGTACGTACCTGCTCGCAATAGCGCAAGCGCCGCAGGCATGCGCGGGCCTGGTCGCCACCACTGCTGGCGTGGTACAGCACCTTGTCGGGGCCGACATCATCCTCATCCAGCCAAAgc                      <  5:358090/136‑1 (MQ=255)
                                                                                                                                     gcgcTGCACTTGCCGCGCAGTACGTACCTGCTCGCAATAGCGCAAGCGCCGCAGGCATGCGCGGGCCTGGTCGCCACCACTGCTGGCGTGGTACAGCACCTTGTCGGGGCCGACATCATCCTCATCCAGCCAAAGCGCCTCGTCGTCct         >  1:140035/1‑149 (MQ=255)
                                                                                                                                     gcgcTGCACTTGCCGCGCAGTACGTACCTGCTCGCAATAGCGCAAGCGCCGCAGGCATGCGCGGGCCTGGTCGCCACCACTGCTGGCGTGGTACAGCACCTTGTCGGGGCCGACATCATCCTCATCCAGCCAAAGCGCCTCGTCGTCct         >  4:4548/1‑149 (MQ=255)
                                                                                                                                         tGCACTTGCCGCGCAGTACGTACCTGCTCGCAATAGCGCAAGCGACGCAGGCATGCGCGGGCCTGGTCGCCACCACTGCTGGCGTGGTACAGCACCTTGTCGGGGCCGACATCATCCTCATCCAGCCAAAGCGCCTCGTCGTCCTCTGc     >  6:152144/1‑149 (MQ=255)
                                                                                                                                          gcaagtGCAGCGCAGTACGTACCTGCTCGCAATAGCGCAAGCGCCGCAGGCATGCGCGGGCCTGGTCGCCACCACTGCTGGCGTGGTACAGCACCTTGTCGGGGCCGACATCATCCTCATCCAGCCAAAGCGCCTCGTCGTCCTCTGcc    <  4:223261/144‑1 (MQ=255)
                                                                                                                                            aCTTGCCGCGCAGTACGTACCTGCTCGCAATAGCGCAAGCGCCGCAGGCATGCGCGGGCCTGGTCGCCACCACTGCTGGCGTGGTACAGCACCTTGTCGGGGCCGACATCATCCTCATCCAGCCAAAGCGCCTCGTCGTCCTCTGCCtt  <  2:180283/149‑1 (MQ=255)
                                                                                                                                            aCTTGCCGCGCAGTACGTACCTGCTCGCAATAGCGCAAGCGCCGCAGGCATGCGCGGGCCTGGTCGCCACCACTGCTGGCGTGGTACAGCACCTTGTCGGGGCCGACATCATCCTCATCCAGCCAAAGCGCCTCGTCGTCCTCTGCCt   >  3:120795/1‑148 (MQ=255)
                                                                                                                                                  cgcgCAGTACGTACCTGCTCGCAATAGCGCAAGCGCCGCAGGCATGCGCGGGCCTGGTCGCCACCACTGCTGGCGTGGTACAGCACCTTGTCGGGGCCGa                                             <  3:200314/100‑1 (MQ=255)
                                                                                                                                                  cgcgCAGTACGTACCTGCTCGCAATAGCGCAAGCGCCGCAGGCATGCGCGGGCCTGGTCGCCACCACTGCTGGCGTGGTACAGCACCTTGTCGGGGCCGa                                             >  4:200314/1‑100 (MQ=255)
                                                                                                                                                  |                                                                                                                                              
AGGGCTTGCCCACCGCATCGCGCTTGTAGCGCCCCGGGTAATCAAAGCGCTCGTAGTCCAGTGACTGGTGCGTGACGCCACTGGTCTGGGCCTTGTGCTCCTGCCGGTAGGCCGGGTTGGTGAAGGTGTAGTCGCGCTGCACTTGCCGCGCAGTACGTACCTGCTCGCAATAGCGCAAGCGCCGCAGGCATGCGCGGGCCTGGTCGCCACCACTGCTGGCGTGGTACAGCACCTTGTCGGGGCCGACATCATCCTCATCCAGCCAAAGCGCCTCGTCGTCCTCTGCCTT  >  NC_002947/4562987‑4563275

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