Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,958,639 T→G 9.4% M36L (ATG→CTG)  fliG ← flagellar motor switch protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,958,6390TG 9.4% 165.9 / 9.1 64M36L (ATG→CTG) fliGflagellar motor switch protein
Reads supporting (aligned to +/- strand):  ref base T (28/30);  new base G (4/2);  total (32/32)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.28e-01

GGATGTACGCATCAGAACCGACGCCCAGGCTGGTCTGGTCGCCGACGATGTCGACGAACTCGCTCATCACCTGTTCGACCTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAACGGCTCGGTTATCACTCATCGGCGTTGATCCAGTCTTTCACGACCTGAGCCACACGGCCCGGG  >  NC_002947/4958492‑4958788
                                                                                                                                                   |                                                                                                                                                     
ggATGTACGCATCAGAACCGACGCCCAGGCTGGTCTGGTCGCCGACGATGTCGACGAACTCGCTCATCACCTGTTCGACCTGGTCACGGTGCACATTGCCCATCTGCGCCCTGGCCACACCCTCCCGCTGCACTTCCTTCGGGCCCAtgt                                                                                                                                                     >  1:163532/1‑150 (MQ=255)
ggATGTACGCATCAGAACCGACGCCCAGGCTGGTCTGGTCGCCGACGATGTCGACGAACTCGCTCATCACCTGTTCGACCTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCAtgt                                                                                                                                                     <  2:880386/150‑1 (MQ=255)
     tACGCATCAGAACCGACGCCCAGGCTGGTCTGGTCGCCGACGATGTCGACGAACTCGCTCATCACCTGTTCGACGTGTGCACGGTGCACATTGCTCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCgc                                                                                                                                                <  1:884392/150‑1 (MQ=255)
        gCATCAGAACCGACGCCCAGGCTGGTCTGGTCGCCGACGATGTCGACGAACTCGCTCATCACCTGTTCGACCTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCgcagc                                                                                                                                             >  1:288687/1‑150 (MQ=255)
              gAACCGACGCCCAGGCTGGTCTGGTCGCCGACGATGTCGACGAACTCGCTCATCACCTGTTCGACCTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTgg                                                                                                                                       <  2:175055/150‑1 (MQ=255)
               aaCCGACGCCCAGGCTGGTCTGGTCGCCGACGATGTCGACGAACTCGCTCATCACCTGTTCGACCTGGTCACGATGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTggg                                                                                                                                      <  2:1206312/150‑1 (MQ=255)
                 ccGACGCCCAGGCTGGTCTGGTCGCCGACGATGTCGACGAACTCGCTCATCACCTGTTCGACCTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGcc                                                                                                                                    <  2:221461/150‑1 (MQ=255)
                   gACGCCCAGGCTGGTCTGGTCGCCGACGATGTCGACGAACTCGCTCATCACCTGTTCGACCTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGc                                                                                                                                  >  2:1343712/1‑150 (MQ=255)
                        ccAGGCTGGTCTGGTCGCCGACGATGTCGACGAACTCGCTCATCACCTGTTCGACCTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCAGGTGCTGCGGCACATgg                                                                                                                                       <  1:119013/140‑1 (MQ=255)
                        ccAGGCTGGTCTGGTCGCCGACGATGTCGACGAACTCGCTCATCACCTGTTCGACCTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCAGGTGCTGCGGCACATgg                                                                                                                                       >  2:119013/1‑140 (MQ=255)
                         cAGGCTGGTCTGGTCGCCGACGATGTCGACGAACTCGCTCATCACCTGTTCGACCTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGt                                                                                                                            <  2:229277/150‑1 (MQ=255)
                         cAGGCTGGTCTGGTCGCCGACGATGTCGACGAACTCGCTCATCACCTGTTCGACCTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCAGGTGCTGCGGCACATgg                                                                                                                                       >  2:850910/1‑139 (MQ=255)
                         cAGGCTGGTCTGGTCGCCGACGATGTCGACGAACTCGCTCATCACCTGTTCGACCTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCAGGTGCTGCGGCACATgg                                                                                                                                       <  1:850910/139‑1 (MQ=255)
                                      tcGCCGACGATTTCGACGAACTCACTCATTACCTGTTCGACCTGGGCACGGTGCACAGTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCgag                                                                                                               <  2:163532/150‑1 (MQ=255)
                                      tcGCCGACGATGTCGACGAACTCGCTCATCACCTGTTCGACCTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCgag                                                                                                               >  1:314715/1‑150 (MQ=255)
                                       cGCCGACGATGTCGACGAACTCGCTCATCACCTGTTCGACCTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCgaga                                                                                                              <  2:1316044/150‑1 (MQ=255)
                                        gCCGACGATGTCGACGAACTCGCTCATCACCTGTTCGACCTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGagag                                                                                                             <  1:8083/150‑1 (MQ=255)
                                            acgaTGTCGACGAACTCGCTCATCACCTGTTCGACCTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCagg                                                                                                         >  1:708635/1‑150 (MQ=255)
                                                tGTCGACGAACTCGCTCATCACCTGTTCGACCTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAggagga                                                                                                     <  1:1343712/150‑1 (MQ=255)
                                                           tCGCTCATCACCTGTTCGACCTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTg                                                                                          >  2:1253307/1‑150 (MQ=255)
                                                            cGCTCATCACCTGTTCGACCTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTACTTCGGGCCCATGTTCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGt                                                                                         <  2:983879/150‑1 (MQ=255)
                                                               tcatcaCCTGTTCGACCTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGa                                                                                      >  2:17516/1‑150 (MQ=255)
                                                                catcaCCTGTTCGACCTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTg                                                                                                                                        <  1:271312/99‑1 (MQ=255)
                                                                catcaCCTGTTCGACCTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTg                                                                                                                                        >  2:271312/1‑99 (MQ=255)
                                                                 atcaCCTGTTCGACCTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCCACCTCCCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGAcg                                                                                    >  1:1084984/1‑150 (MQ=255)
                                                                     ccTGTTCGACCTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATcgccgc                                                                                              >  2:17840/1‑136 (MQ=255)
                                                                     ccTGTTCGACCTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATcgccgc                                                                                              <  1:17840/136‑1 (MQ=255)
                                                                     ccTGTTCGACCTGGTCACGGTGCACATTGCCCACCAGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGc                                                                                >  1:337186/1‑150 (MQ=255)
                                                                         ttCCACCTGGTCACGGTGCACATTGCCCATCCGCGCCATGGCCACCCCCCCCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGAATCGCCGCCTTGTCGACGCGGCGCAg                                                                            <  2:1028557/150‑1 (MQ=255)
                                                                            gACCTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGGGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCtt                                                                         >  1:809989/1‑150 (MQ=255)
                                                                               cTGGTCACGGTGCACATTGCCCATCTGCGCCATGGTCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGc                                                                      <  2:337186/150‑1 (MQ=255)
                                                                               cTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGc                                                                      >  1:971896/1‑150 (MQ=255)
                                                                               cTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGc                                                                      >  1:971744/1‑150 (MQ=255)
                                                                               cTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGc                                                                      >  2:656291/1‑150 (MQ=255)
                                                                                tGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCg                                                                     <  2:924370/150‑1 (MQ=255)
                                                                                tGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGAAGGATCGCCGCCTTGTGGACGCGGCTCAGCTTGGCg                                                                     >  2:182400/1‑150 (MQ=255)
                                                                                       ggTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAACgg                                                              >  1:223204/1‑150 (MQ=255)
                                                                                          gCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCagga                                                                                                        <  2:1106062/105‑1 (MQ=255)
                                                                                          gCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCagga                                                                                                        >  1:1106062/1‑105 (MQ=255)
                                                                                             caTTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAACGGCTCGGt                                                        <  2:1204471/150‑1 (MQ=255)
                                                                                                 gCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTACGGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAACGGCTCGGTTATc                                                    <  1:728137/150‑1 (MQ=255)
                                                                                                   ccATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAACGGCTCGGTTATCAc                                                  >  2:1209600/1‑150 (MQ=255)
                                                                                                   ccATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAACGGCTCGGTTATCAc                                                  >  1:747470/1‑150 (MQ=255)
                                                                                                      tCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAACGGCTCGGTTATCACTCa                                               >  2:902070/1‑150 (MQ=255)
                                                                                                           gcCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAACGGCTCGGTTATCACTCATCGGc                                          <  1:967455/150‑1 (MQ=255)
                                                                                                            ccATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAACGGCTCGGTTATCACTCATCGGCg                                         >  2:1166968/1‑150 (MQ=255)
                                                                                                              aTGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGCGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTACCGGCTCGGTTATCACTCCTCGGCGtt                                       >  2:586321/1‑150 (MQ=255)
                                                                                                              aTGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAACGGCTCGGTTATCACTCATCGGCGtt                                       <  2:747470/150‑1 (MQ=255)
                                                                                                                  ccACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAACGGCTCGGTTATCACTCATCGGCGTTGATc                                   <  1:17516/150‑1 (MQ=255)
                                                                                                                      acccacccGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAACGGCTCGGTTATCACTCATCGGCGTTGATCCAGt                               <  2:1270197/150‑1 (MQ=255)
                                                                                                                       cccacccGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAACGGCTCGGTTATCACTCATCGGCGTTGTTCCAGTc                              <  1:902070/150‑1 (MQ=255)
                                                                                                                        ccacccGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAACGGCTCGGTTATCACTCATCGGCGtt                                       <  2:472777/140‑1 (MQ=255)
                                                                                                                        ccacccGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAACGGCTCGGTTATCACTCATCGGCGtt                                       >  1:472777/1‑140 (MQ=255)
                                                                                                                         cacccGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAACGGCTCgg                                                         >  1:1103659/1‑121 (MQ=255)
                                                                                                                         cacccGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAACGGCTCgg                                                         <  2:1103659/121‑1 (MQ=255)
                                                                                                                          acccGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAAc                                                                <  1:1348429/113‑1 (MQ=255)
                                                                                                                          acccGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAAc                                                                >  2:1348429/1‑113 (MQ=255)
                                                                                                                           cccGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAACGGCTCGGTTATCACTCATCGGCGTTGATCCAGTCTTTc                          <  2:314715/150‑1 (MQ=255)
                                                                                                                            ccGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGACCAGGAGGATCACCACCTTGTCGACGCGGCTCAGCTTGGCGGTAACGGCTCGGTTATGACTCCGCGGTGTTGATCCAGTCTTTCa                         >  2:601175/1‑150 (MQ=255)
                                                                                                                                tGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAACGGCTCGGTTATCACTCATCGGCGTTGATCCAGTCTTTCACGAc                     <  1:1253307/150‑1 (MQ=255)
                                                                                                                                        cttcGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAACGGCTCGGTTATCACTCATCGGCGTTGATCCa                                 <  1:1259404/130‑1 (MQ=255)
                                                                                                                                        cttcGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAACGGCTCGGTTATCACTCATCGGCGTTGATCCa                                 >  2:1259404/1‑130 (MQ=255)
                                                                                                                                               cccATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAACGGCGCGGTTATCACTCATCGGCGTTGATCCAGTCTTTCAAGACCTGAGCCACACGGcc      >  1:759864/1‑150 (MQ=255)
                                                                                                                                                ccAGGTGCTGCGGCACATGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAACGGCTCGGTTATCACTCATCGGCGTTGATCCAGTCTTTCACGACCTGAGCCACACGGccc     >  1:870106/1‑150 (MQ=255)
                                                                                                                                                ccAGGTGCTGCGGCACATGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAACGGCTCGGTTATCACTCATCGGCGTTGATCCAGTCTTTCACGACCTGAGCCACACGGccc     >  2:1287117/1‑150 (MQ=255)
                                                                                                                                                  aTGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAACGGCTCGGTTATCACTCATCGGCGTTGATCCAGTCTTTCACGACCTGAGCCACACGGCCCgg   <  1:1209600/150‑1 (MQ=255)
                                                                                                                                                   tgtgCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGAc                                                                                     <  1:559013/67‑1 (MQ=255)
                                                                                                                                                   tgtgCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGAc                                                                                     >  2:559013/1‑67 (MQ=255)
                                                                                                                                                   tgtgCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAACGGCTCGGTTATCACTCATCGGCGTTGATCCAGTCTTTCACGACCTGAGCCACACGGCCCggg  <  2:223204/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GGATGTACGCATCAGAACCGACGCCCAGGCTGGTCTGGTCGCCGACGATGTCGACGAACTCGCTCATCACCTGTTCGACCTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAACGGCTCGGTTATCACTCATCGGCGTTGATCCAGTCTTTCACGACCTGAGCCACACGGCCCGGG  >  NC_002947/4958492‑4958788

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: