Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,958,647 A→G 10.6% L33P (CTG→CCG)  fliG ← flagellar motor switch protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,958,6470AG10.6% 177.9 / 9.9 66L33P (CTG→CCG) fliGflagellar motor switch protein
Reads supporting (aligned to +/- strand):  ref base A (28/31);  new base G (4/3);  total (32/34)
Fisher's exact test for biased strand distribution p-value = 7.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.94e-01

GCATCAGAACCGACGCCCAGGCTGGTCTGGTCGCCGACGATGTCGACGAACTCGCTCATCACCTGTTCGACCTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAACGGCTCGGTTATCACTCATCGGCGTTGATCCAGTCTTTCACGACCTGAGCCACACGGCCCGGG  >  NC_002947/4958500‑4958788
                                                                                                                                                   |                                                                                                                                             
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                                                                                                                                       >  2:119013/1‑140 (MQ=255)
                ccAGGCTGGTCTGGTCGCCGACGATGTCGACGAACTCGCTCATCACCTGTTCGACCTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCAGGTGCTGCGGCACATgg                                                                                                                                       <  1:119013/140‑1 (MQ=255)
                 cAGGCTGGTCTGGTCGCCGACGATGTCGACGAACTCGCTCATCACCTGTTCGACCTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGt                                                                                                                            <  2:229277/150‑1 (MQ=255)
                 cAGGCTGGTCTGGTCGCCGACGATGTCGACGAACTCGCTCATCACCTGTTCGACCTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCAGGTGCTGCGGCACATgg                                                                                                                                       <  1:850910/139‑1 (MQ=255)
                 cAGGCTGGTCTGGTCGCCGACGATGTCGACGAACTCGCTCATCACCTGTTCGACCTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCAGGTGCTGCGGCACATgg                                                                                                                                       >  2:850910/1‑139 (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                                                                      >  2:656291/1‑150 (MQ=255)
                                                                       cTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGc                                                                      >  1:971744/1‑150 (MQ=255)
                                                                       cTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGc                                                                      >  1:971896/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                                       >  1:472777/1‑140 (MQ=255)
                                                                                                                ccacccGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAACGGCTCGGTTATCACTCATCGGCGtt                                       <  2:472777/140‑1 (MQ=255)
                                                                                                                 cacccGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAACGGCTCgg                                                         >  1:1103659/1‑121 (MQ=255)
                                                                                                                 cacccGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAACGGCTCgg                                                         <  2:1103659/121‑1 (MQ=255)
                                                                                                                  acccGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAAc                                                                >  2:1348429/1‑113 (MQ=255)
                                                                                                                  acccGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAAc                                                                <  1:1348429/113‑1 (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     >  2:1287117/1‑150 (MQ=255)
                                                                                                                                        ccAGGTGCTGCGGCACATGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAACGGCTCGGTTATCACTCATCGGCGTTGATCCAGTCTTTCACGACCTGAGCCACACGGccc     >  1:870106/1‑150 (MQ=255)
                                                                                                                                          aggtgCTGCGGCACATGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAACGGCTCGGTTATCACTCATCGGCGTTGATCCAGTCTTTc                          <  1:1215051/125‑1 (MQ=255)
                                                                                                                                          aTGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAACGGCTCGGTTATCACTCATCGGCGTTGATCCAGTCTTTCACGACCTGAGCCACACGGCCCgg   <  1:1209600/150‑1 (MQ=255)
                                                                                                                                           tgtgCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGAc                                                                                     >  2:559013/1‑67 (MQ=255)
                                                                                                                                           tgtgCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGAc                                                                                     <  1:559013/67‑1 (MQ=255)
                                                                                                                                           tgtgCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAACGGCTCGGTTATCACTCATCGGCGTTGATCCAGTCTTTCACGACCTGAGCCACACGGCCCggg  <  2:223204/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                             
GCATCAGAACCGACGCCCAGGCTGGTCTGGTCGCCGACGATGTCGACGAACTCGCTCATCACCTGTTCGACCTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAACGGCTCGGTTATCACTCATCGGCGTTGATCCAGTCTTTCACGACCTGAGCCACACGGCCCGGG  >  NC_002947/4958500‑4958788

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