Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,300,903 A→G G262G (GGT→GGC yjcR ← putative membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,300,9030AG99.7% 1223.5 / NA 386G262G (GGT→GGCyjcRputative membrane protein
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (168/217);  total (169/217)
Fisher's exact test for biased strand distribution p-value = 4.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 398 total aligned reads displayed.

GTCATCCGGTAGCACGCCGTAGCCAATCGAATCCACTTTACCCTCGAAGGTTTTGCCGCTGTCACTCATCAGGCGAATCGTTGCGGGTGTACCTGAGCGAATATTTTTCAGATCGGTTTCGCGGAAGTTGGCGATCACATACCAGTGACGAGTGTCGATTAGGGTAAAAATAGGGCGCAT  >  NC_000913/4300813‑4300992
                                                                                          |                                                                                         
gTCATCCGGTAGCACGCCGTAGCCAATCGAATCCACTTTACCCTCGAAGGTTTTGCCGCTGTCACTCATCAGGCGAATCGTTGCGGGTGTGcc                                                                                         <  1:1097798/93‑1 (MQ=255)
gTCATCCGGTAGCACGCCGTAGCCAATCGAATCCACTTTACCCTCGAAGGTTTTGCCGCTGTCACTCATCAGGCGAATCGTTGCGGGTGTGcc                                                                                         >  1:4584743/1‑93 (MQ=255)
 tCATCCGGTAGCACGCCGTAGCCAATCGAATCCACTTTACCCTCGAAGGTTTTGCCGCTGTCACTCATCAGGCGAATCGTTGCGGGTGTGCCt                                                                                        <  1:3383432/93‑1 (MQ=255)
    tCCGGTAGCACGCCGTAGCCAATCGAATCCACTTTACCCTCGAAGGTTTTGCCGCTGTCACTCATCAGGCGAATCGTTGCGGGTGTGCCTGAg                                                                                     >  1:4884166/1‑93 (MQ=255)
     ccGGTAGCACGCCGTAGCCAATCGAATCCACTTTACCCTCGAAGGTTTTGCCGCTGTCACTCATCAGGCGAATCGTTGCGGGTGTGCCTGAGc                                                                                    <  1:5698763/93‑1 (MQ=255)
      cGGTAGCACGCCGTAGCCAATCGAATCCACTTTACCCTCGAAGGTTTTGCCGCTGTCACTCATCAGGCGAATCGTTGCGGGTGTGCCTGAGCg                                                                                   >  2:1835528/1‑93 (MQ=255)
        gTAGCACGCCGTAGCCAATCGAATCCACTTTACCCTCGAAGGTTTTGCCGCTGTCACTCATCAGGCGAATCGTTGCGGGTGTGCCTGAGCGaa                                                                                 <  2:769691/93‑1 (MQ=255)
        gTAGCACGCCGTAGCCAATCGAATCCACTTTACCCTCGAAGGTTTTGCCGCTGTCACTCATCAGGCGAATCGTTGCGGGTGTGCCTGAGCGaa                                                                                 >  1:5458273/1‑93 (MQ=255)
          aGCACGCCGTAGCCAATCGAATCCACTTTACCCTCGAAGGTTTTGCCGCTGTCACTCATCAGGCGAATCGTTGCGGGTGTGCCTGAGCGAata                                                                               <  2:20324/93‑1 (MQ=255)
           gCACGCCGTAGCCAATCGAATCCACTTTACCCTCGAAGGTTTTGCCGCTGTCACTCATCAGGCGAATCGTTGCGGGTGTGCCTGAGCGAatat                                                                              <  1:1702584/93‑1 (MQ=255)
             aCGCCGTAGCCAATCGAATCCACTTTACCCTCGAAGGTTTTGCCGCTGTCACTCATCAGGCGAATCGTTGCGGGTGTGCCTGAGCGAATAttt                                                                            >  1:5146038/1‑93 (MQ=255)
                ccGTAGCCAATCGAATCCACTTTACCCTCGAAGGTTTTGCCGCTGTCACTCATCAGGCGAATCGTTGCGGGTGTGCCTGAGCGAATATTTTTc                                                                         <  2:713653/93‑1 (MQ=255)
                 cGTAGCCAATCGAATCCACTTTACCCTCGAAGGTTTTGCCGCTGTCACTCATCAGGCGAATCGTTGCGGGTGTGCCTGAGCGAATATTTTTCa                                                                        <  2:3316873/93‑1 (MQ=255)
                 cGTAGCCAATCGAATCCACTTTACCCTCGAAGGTTTTGCCGCTGTCACTCATCAGGCGAATCGTTGCGGGTGTGCCTGAGCGAATATTTTTCa                                                                        <  2:3972647/93‑1 (MQ=255)
                 cGTAGCCAATCGAATCCACTTTACCCTCGAAGGTTTTGCCGCTGTCACTCATCAGGCGAATCGTTGCGGGTGTGCCTGAGCGAATATTTTTCa                                                                        <  1:7419895/93‑1 (MQ=255)
                  gTAGCCAATCGAATCCACTTTACCCTCGAAGGTTTTGCCGCTGTCACTCATCAGGCGAATCGTTGCGGGTGTGCCTGAGCGAATATTTTTCAg                                                                       >  1:7163959/1‑93 (MQ=255)
                   tAGCCAATCGAATCCACTTTACCCTCGAAGGTTTTGCCGCTGTCACTCATCAGGCGAATCGTTGCGGGTGTGCCTGAGCGAATATTTTTCAGa                                                                      >  2:7909985/1‑93 (MQ=255)
                   tAGCCAATCGAATCCACTTTACCCTCGAAGGTTTTGCCGCTGTCACTCATCAGGCGAATCGTTGCGGGTGTGCCTGAGCGAATATTTTTCAGa                                                                      <  1:2083774/93‑1 (MQ=255)
                    aGCCAATCGAATCCACTTTACCCTCGAAGGTTTTGCCGCTGTCACTCATCAGGCGAATCGTTGCGGGTGTGCCTGAGCGAATATTTTTCAGAt                                                                     >  1:2786107/1‑93 (MQ=255)
                    aGCCAATCGAATCCACTTTACCCTCGAAGGTTTTGCCGCTGTCACTCATCAGGCGAATCGTTGCGGGTGTGCCTGAGCGAATATTTTTCAGAt                                                                     <  1:7378850/93‑1 (MQ=255)
                      ccAATCGAATCCACTTTACCCTCGAAGGTTTTGCCGCTGTCACTCATCAGGCGAATCGTTGCGGGTGTGCCTGAGCGAATATTTTTCAGATCg                                                                   <  1:4156451/93‑1 (MQ=255)
                        aaTCGAATCCACTTTACCCTCGAAGGTTTTGCCGCTGTCACTCATCAGGCGAATCGTTGCGGGTGTGCCTGAGCGAATATTTTTCAGATCGGt                                                                 >  1:7688689/1‑93 (MQ=255)
                         aTCGAATCCACTTTACCCTCGAAGGTTTTGCCGCTGTCACTCATCAGGCGAATCGTTGCGGGTGTGCCTGAGCGAATATTTTTCAGATCGGtt                                                                <  2:33675/93‑1 (MQ=255)
                         aTCGAATCCACTTTACCCTCGAAGGTTTTGCCGCTGTCACTCATCAGGCGAATCGTTGCGGGTGTGCCTGAGCGAATATTTTTCAGATCGGtt                                                                >  2:658051/1‑93 (MQ=255)
                          tCGAATCCACTTTACCCTCGAAGGTTTTGCCGCTGTCACTCATCAGGCGAATCGTTGCGGGTGTGCCTGAGCGAATATTTTTCAGATCGGttt                                                               <  2:1375651/93‑1 (MQ=255)
                          tCGAATCCACTTTACCCTCGAAGGTTTTGCCGCTGTCACTCATCAGGCGAATCGTTGCGGGTGTGCCTGAGCGAATATTTTTCAGATCGGttt                                                               <  2:3756256/93‑1 (MQ=255)
                          tCGAATCCACTTTACCCTCGAAGGTTTTGCCGCTGTCACTCATCAGGCGAATCGTTGCGGGTGTGCCTGAGCGAATATTTTTCAGATCGGttt                                                               <  1:7823352/93‑1 (MQ=255)
                          tCGAATCCACTTTACCCTCGAAGGTTTTGCCGCTGTCACTCATCAGGCGAATCGTTGCGGGTGTGCCTGAGCGAATATTTTTCAGATCGGttt                                                               <  1:6507540/93‑1 (MQ=255)
                             aaTCCACTTTACCCTCGAAGGTTTTGCCGCTGTCACTCATCAGGCGAATCGTTGCGGGTGTGCCTGAGCGAATATTTTTCAGATCGGTTTcgc                                                            <  1:462975/93‑1 (MQ=255)
                             aaTCCACTTTACCCTCGAAGGTTTTGCCGCTGTCACTCATCAGGCGAATCGTTGCGGGTGTGCCTGAGCGAATATTTTTCAGATCGGTTTcgc                                                            >  2:4447436/1‑93 (MQ=255)
                              aTCCACTTTACCCTCGAAGGTTTTGCCGCTGTCACTCATCAGGCGAATCGTTGCGGGTGTGCCTGAGCGAATATTTTTCAGATCGGTTTcgcg                                                           <  1:5627336/93‑1 (MQ=255)
                               tCCACTTTACCCTCGAAGGTTTTGCCGCTGTCACTCATCAGGCGAATCGTTGCGGGTGTGCCTGAGCGAATATTTTTCAGATCGGTTTCGCgg                                                          >  1:3248979/1‑93 (MQ=255)
                                ccACTTTACCCTCGAAGGTTTTGCCGCTGTCACTCATCAGGCGAATCGTTGCGGGTGTGCCTGAGCGAATATTTTTCAGATCGGTTTCGCGGa                                                         >  1:1275317/1‑93 (MQ=255)
                                  aCTTTACCCTCGAAGGTTTTGCCGCTGTCACTCATCAGGCGAATCGTTGCGGGTGTGCCTGAGCGAATATTTTTCAGATCGGTTTCGCGGAAg                                                       >  2:5070620/1‑93 (MQ=255)
                                  aCTTTACCCTCGAAGGTTTTGCCGCTGTCACTCATCAGGCGAATCGTTGCGGGTGTGCCTGAGCGAATATTTTTCAGATCGGTTTCGCGGAAg                                                       <  2:221665/93‑1 (MQ=255)
                                  aCTTTACCCTCGAAGGTTTTGCCGCTGTCACTCATCAGGCGAATCGTTGCGGGTGTGCCTGAGCGAATATTTTTCAGATCGGTTTCGCGGAAg                                                       >  2:5831105/1‑93 (MQ=255)
                                   cTTTACCCTCGAAGGTTTTGCCGCTGTCACTCATCAGGCGAATCGTTGCGGGTGTGCCTGAGCGAATATTTTTCAGATCGGTTTCGCGGAAGt                                                      <  1:6656183/93‑1 (MQ=255)
                                   cTTTACCCTCGAAGGTTTTGCCGCTGTCACTCATCAGGCGAATCGTTGCGGGTGTGCCTGAGCGAATATTTTTCAGATCGGTTTCGCGGAAGt                                                      >  2:2907951/1‑93 (MQ=255)
                                    tttACCCTCGAAGGTTTTGCCGCTGTCACTCATCAGGCGAATCGTTGCGGGTGTGCCTGAGCGAATATTTTTCAGATCGGTTTCGCGGAAGtt                                                     >  1:2518603/1‑93 (MQ=255)
                                    tttACCCTCGAAGGTTTTGCCGCTGTCACTCATCAGGCGAATCGTTGCGGGTGTGCCTGAGCGAATATTTTTCAGATCGGTTTCGCGGAAGtt                                                     >  1:2154042/1‑93 (MQ=255)
                                     ttACCCTCGAAGGTTTTGCCGCTGTCACTCATCAGGCGAATCGTTGCGGGTGTGCCTGAGCGAATATTTTTCAGATCGGTTTCGCGGAAGTTg                                                    <  2:7227819/93‑1 (MQ=255)
                                      tACCCTCGAAGGTTTTGCCGCTGTCACTCATCAGGCGAATCGTTGCGGGTGTGCCTGAGCGAATATTTTTCAGATCGGTTTCGCGGAAGTTgg                                                   <  1:1968886/93‑1 (MQ=255)
                                      tACCCTCGAAGGTTTTGCCGCTGTCACTCATCAGGCGAATCGTTGCGGGTGTGCCTGAGCGAATATTTTTCAGATCGGTTTCGCGGAAGTTgg                                                   <  2:144882/93‑1 (MQ=255)
                                       aCCCTCGAAGGTTTTGCCGCTGTCACTCATCAGGCGAATCGTTGCGGGTGTGCCTGAGCGAATATTTTTCAGATCGGTTTCGCGGAAGTTGGc                                                  <  2:1099968/93‑1 (MQ=255)
                                        cccTCGAAGGTTTTGCCGCTGTCACTCATCAGGCGAATCGTTGCGGGTGTGCCTGAGCGAATATTTTTCAGATCGGTTTCGCGGAAGTTGGCg                                                 >  1:6363827/1‑93 (MQ=255)
                                         ccTCGAAGGTTTTGCCGCTGTCACTCATCAGGCGAATCGTTGCGGGTGTGCCTGAGCGAATATTTTTCAGATCGGTTTCGCGGAAGTTGGCGa                                                >  2:3444104/1‑93 (MQ=255)
                                         ccTCGAAGGTTTTGCCGCTGTCACTCATCAGGCGAATCGTTGCGGGTGTGCCTGAGCGAATATTTTTCAGATCGGTTTCGCGGAAGTTGGCGa                                                <  1:4769664/93‑1 (MQ=255)
                                         ccTCGAAGGTTTTGCCGCTGTCACTCATCAGGCGAATCGTTGCGGGTGTGCCTGAGCGAATATTTTTCAGATCGGTTTCGCGGAAGTTGGCGa                                                <  2:5919284/93‑1 (MQ=255)
                                          cTCGAAGGTTTTGCCGCTGTCACTCATCAGGCGAATCGTTGCGGGTGTGCCTGAGCGAATATTTTTCAGATCGGTTTCGCGGAAGGTGGCGAt                                               >  1:1879585/1‑93 (MQ=255)
                                           tCGAAGGTTTTGCCGCTGTCACTCATCAGGCGAATCGTTGCGGGTGTGCCTGAGCGAATATTTTTCAGATCGGTTTCGCGGAAGTTGGCGATc                                              <  2:3083744/93‑1 (MQ=255)
                                            cGAAGGTTTTGCCGCTGTCACTCATCAGGCGAATCGTTGCGGGTGTGCCTGAGCGAATATTTTTCAGATCGGTTTCGCGGAAGTTGGCGATca                                             <  1:5847628/93‑1 (MQ=255)
                                            cGAAGGTTTTGCCGCTGTCACTCATCAGGCGAATCGTTGCGGGTGTGCCTGAGCGAATATTTTTCAGATCGGTTTCGCGGAAGTTGGCGATca                                             <  1:1214569/93‑1 (MQ=255)
                                             gAAGGTTTTGCCGCTGTCACTCATCAGGCGAATCGTTGCGGGTGTGCCTGAGCGAATATTTTTCAGATCGGTTTCGCGGAAGTTGGCGATcac                                            >  2:4578879/1‑93 (MQ=255)
                                             gAAGGTTTTGCCGCTGTCACTCATCAGGCGAATCGTTGCGGGTGTGCCTGAGCGAATATTTTTCAGATCGGTTTCGCGGAAGTTGGCGATcac                                            >  2:6348886/1‑93 (MQ=255)
                                             gAAGGTTTTGCCGCTGTCACTCATCAGGCGAATCGTTGCGGGTGTGCCTGAGCGAATATTTTTCAGATCGGTTTCGCGCAAGTTGGCGATcac                                            >  2:5698424/1‑93 (MQ=255)
                                             gAAGGTTTTGCCGCTGTCACTCATCAGGCGAATCGTTGCGGGTGTGCCTGAGCGAATATTTTTCAGATCGGGTTCGCGGAAGTTGGCGATcac                                            >  2:3220155/1‑93 (MQ=255)
                                              aaGGTTTTGCCGCTGTCACTCATCAGGCGAATCGTTGCGGGTGTGCCTGAGCGAATATTTTTCAGATCGGTTTCGCGGAAGTTGGCGATcaca                                           <  1:3622912/93‑1 (MQ=255)
                                              aaGGTTTTGCCGCTGTCACTCATCAGGCGAATCGTTGCGGGTGTGCCTGAGCGAATATTTTTCAGATCGGTTTCGCGGAAGTTGGCGATcaca                                           <  1:7472741/93‑1 (MQ=255)
                                                agTTTTGCCGCTGTCACTCATCAGGCGAATCGTTGCGGGTGTGCCTGAGCGAATATTTTTCAGATCGGTTTCGCGGAAGTTGGCGATCACATa                                         <  1:2974895/92‑1 (MQ=255)
                                                 gTTTTGCCGCTGTCACTCATCAGGCGAATCGTTGCGGGTGTGCCTGAGCGAATATTTTTCAGATCGGTTTCGCGGAAGTTGGCGATCACATAc                                        >  2:325758/1‑93 (MQ=255)
                                                    ttGCCGCTGTCACTCATCAGGCGAATCGTTGCGGGTGTGCCTGAGCGAATATTTTTCAGATCGGTTTCGCGGAAGTTGGCGATCACATACCAg                                     <  1:368126/93‑1 (MQ=255)
                                                    ttGCCGCTGTCACTCATCAGGCGAATCGTTGCGGGTGTGCCTGAGCGAATATTTTTCAGATCGGTTTCGCGGAAGTTGGCGATCACATACCAg                                     <  1:2706774/93‑1 (MQ=255)
                                                    ttGCCGCTGTCACTCATCAGGCGAATCGTTGCGGGTGTGCCTGAGCGAATATTTTTCAGATCGGTTTCGCGGAAGTTGGCGATCACATACCAg                                     <  2:5369636/93‑1 (MQ=255)
                                                     tGCCGCTGTCACTCATCAGGCGAATCGTTGCGGGTGTGCCTGAGCGAATATTTTTCAGATCGGTTTCGCGGAAGTTGGCGATCACATACCAGt                                    <  2:5202834/93‑1 (MQ=255)
                                                     tGCCGCTGTCACTCATCAGGCGAATCGTTGCGGGTGTGCCTGAGCGAATATTTTTCAGATCGGTTTCGCGGAAGTTGGCGATCACATACCAGt                                    <  1:7677722/93‑1 (MQ=255)
                                                     tGCCGCTGTCACTCATCAGGCGAATCGTTGCGGGTGTGCCTGAGCGAATATTTTTCAGATCGGTTTCGCGGAAGTTGGCGATCACATACCAGt                                    <  1:2304859/93‑1 (MQ=255)
                                                     tGCCGCTGTCACTCATCAGGCGAATCGTTGCGGGTGTGCCTGAGCGAATATTTTTCAGATCGGTTTCGCGGAAGTTGGCGATCACATACCAGt                                    <  1:4293249/93‑1 (MQ=255)
                                                      gCCGCTGTCACTCATCAGGCGAATCGTTGCGGGTGTGCCTGAGCGCATATTTTTCAGATCGGTTTTGCGGAAGTTGGCGATCACATACCAGTg                                   >  1:7121974/1‑93 (MQ=255)
                                                      gCCGCTGTCACTCATCAGGCGAATCGTTGCGGGTGTGCCTGAGCGAATATTTTTCAGATCGGTTTCGCGGAAGTTGGCGATCACATACCAGTg                                   <  2:5777853/93‑1 (MQ=255)
                                                       gcGCTGTCACTCATGAGGCGAATCGTTGCGGGGGTGCCTGAGCGAATATTTTTCAGATCGTTTTCGCGGAAGTTGGCGATCACATACCAGTGa                                  <  1:2445059/92‑1 (MQ=255)
                                                       ccGCTGTCACTCATCAGGCGAATCGTTGCGGGTGTGCCTGAGCGAATATTTTTCAGATCGGTTTCGCGGAAGTTGGCGATCACATACCAGTGa                                  <  2:3643448/93‑1 (MQ=255)
                                                       ccGCTGTCACTCATCAGGCGAATCGTTGCGGGTGTGCCTGAGCGAATATTTTTCAGATCGGTTTCGCGGAAGTTGGCGATCACATACCAGTGa                                  <  2:555683/93‑1 (MQ=255)
                                                       ccGCTGTCACTCATCAGGAGAATCGTTGCGGGTGTGCCTGAGCGAATATTTTTCAGATCGGTTTCGCGGAAGTTGGCGATCACATACCAGTGa                                  <  1:6907971/93‑1 (MQ=255)
                                                          cTGTCACTCATCAGGCGAATCGTTGCGGGTGTGCCTGAGCGAATATTTTTCAGATCGGTTTCGCGGAAGTTGGCGATCACATACCAGTGACGa                               <  2:6099827/93‑1 (MQ=255)
                                                           tGTCACTCATCAGGCGAATCGTTGCGGGTGTGCCTGAGCGAATATTTTTCAGATCGGTTTCGCGGAAGTTGGCGATCACATACCAGTGACGAg                              <  1:390448/93‑1 (MQ=255)
                                                           tGTCACTCATCAGGCGAATCGTTGCGGGTGTGCCTGAGCGAATATTTTTCAGATCGGTTTCGCGGAAGTTGGCGATCACATACCAGTGACGAg                              >  1:4748816/1‑93 (MQ=255)
                                                            gTCACTCATCAGGCGAATCGTTGCGGGTGTGCCTGAGCGAATATTTTTCAGATCGGTTTCGCGGAAGTTGGCGATCACATACCAGTGACGAgt                             <  1:4089861/93‑1 (MQ=255)
                                                             tCACTCATCAGGCGAATCGTTGCGGGTGTGCCTGAGCGAATATTTTTCAGATCGGTTTCGCGGAAGTTGGCGATCACATACCAGTGACGAgtg                            <  1:1562990/93‑1 (MQ=255)
                                                              cACTCATCAGGCGAATCGTTGCGGGTGTGCCTGAGCGAATATTTTTCAGATCGGTTTCGCGGAAGTTGGCGATCACATACCAGTGACGAgtgt                           >  1:3734603/1‑93 (MQ=255)
                                                               aCTCATCAGGCGAATCGTTGCGGGTGTGCCTGAGCGAATATTTTTCAGATCGGTTTCGCGGAAGTTGGCGATCACATACCAGTGACGAGTGTc                          >  1:87807/1‑93 (MQ=255)
                                                               aCTCATCAGGCGAATCGTTGCGGGTGTGCCTGAGCGAATATTTTTCAGATCGGTTTCGCGGAAGTTGGCGATCACATACCAGTGACGAGTGTc                          <  2:75868/93‑1 (MQ=255)
                                                               aCTCATCAGGCGAATCGTTGCGGGTGTGCCTGAGCGAATATTTTTCAGATCGGTTTCGCGGAAGTTGGCGATCACATACCAGTGACGAGTGTc                          <  2:436922/93‑1 (MQ=255)
                                                                cTCATCAGGCGAATCGTTGCGGGTGTGCCTGAGCGAATATTTTTCAGATCGGTTTCGCGGAAGTTGGCGATCACATCCCAGTTACGAGTGTcc                         >  2:6439434/1‑92 (MQ=255)
                                                                 ccatcaGGCGAATCGTTGCGGGTGTGCCTGAGCGAATATTTTTCAGATCGGTTTCGCGGAAGTTGGCGATCACATACCAGTGACGAGTGTCGa                        <  1:1392715/92‑1 (MQ=255)
                                                                  catcaGGCGAATCGTTGCGGGTGTGCCTGAGCGAATATTTTTCAGATCGGTTTCGCGGAAGTTGGCGATCACATACCAGTGACGAGTGTCGAt                       <  2:4015506/93‑1 (MQ=255)
                                                                   atcaGGCGAATCGTTGCGGGTGGCCCTGAGCGAATATTTTTCAGATCGGTTTCGCGGAAGTTGGCGATCACATACCAGTGACGAGGGTCGAtt                      <  2:6818003/93‑1 (MQ=255)
                                                                      aGGCGAATCGTTGCGGGTGTGCCTGAGCGAATATTTTTCAGATCGGTTTCGCGGAAGTTGGCGATCACATACCAGTGACGAGTGTCGATTAgg                   <  2:4712942/93‑1 (MQ=255)
                                                                      aGGCGAATCGTTGCGGGTGTGCCTGAGCGAATATTTTTCAGATCGGTTTCGCGGAAGTTGGCGATCACATACCAGTGACGAGTGTCGATTAgg                   >  2:6670551/1‑93 (MQ=255)
                                                                       ggCGAATCGTTGCGGGTGTGCCTGAGCGAATATTTTTCAGATCGGTTTCGCGGAAGTTGGCGATCACATACCAGTGACGAGTGTCGATTAggg                  <  2:7444359/93‑1 (MQ=255)
                                                                           aaTCGTTGCGGGTGTGCCTGAGCGAATATTTTTCAGATCGGTTTCGCGGAAGTTGGCGATCACATACCAGTGACGAGTGTCGATTAGGGTaaa              <  1:1510108/93‑1 (MQ=255)
                                                                            aTCGTTGCGGGTGTGCCTGAGCGAATATTTTTCAGATCGGTTTCGCGGAAGTTGGCGATCACATACCAGTGACGAGTGTCGATTAGGGTaaaa             >  1:4372216/1‑93 (MQ=255)
                                                                             tCGTTGCGGGTGTGCCTGAGCGAATATTTTTCAGATCGGTTTCGCGGAAGTTGGCGATCACATACCAGTGACGAGTGTCGATTAGGGTaaaaa            <  1:454587/93‑1 (MQ=255)
                                                                             tCGTTGCGGGTGTGCCTGAGCGAATATTTTTCAGATCGGTTTCGCGGAAGTTGGCGATCACATACCAGTGACGAGTGTCGATTAGGGTaaaaa            >  1:6128619/1‑93 (MQ=255)
                                                                               gTTGCGGGTGTGCCTGAGCGAATATTTTTCAGATCGGTTTCGCGGAAGTTGGCGATCACATACCAGTGACGAGTGTCGATTAGGGTAAAAATa          <  1:1007391/93‑1 (MQ=255)
                                                                                ttGCGGGTGTGCCTGAGCGAATATTTTTCAGATCGGTTTCGCGGAAGTTGGCGATCACATACCAGTGACGAGTGTCGATTAGGGTAAAAATAg         >  2:905258/1‑93 (MQ=255)
                                                                                     ggTGTGCCTGAGCGAATATTTTTCAGATCGGTTTCGCGGAAGTTGGCGATCACATACCAGTGACGAGTGTCGATTAGGGTAAAAATAGGgcgc    >  1:5558115/1‑93 (MQ=255)
                                                                                      gtgtGCCTGAGCGAATATTTTTCAGATCGGTTTCGCGGAAGTTGGCGATCACATACCAGTGACGAGTGTCGATTAGGGTAAAAATAGGGCGCa   <  2:6937240/93‑1 (MQ=255)
                                                                                      gtgtGCCTGAGCGAATATTTTTCAGATCGGTTTCGCGGAAGTTGGCGATCACATACCAGTGACGAGTGTCGATTAGGGTAAAAATAGGGCGCa   <  2:4896828/93‑1 (MQ=255)
                                                                                       tgtgCCTGAGCGAATATTTTTCAGATCGGTTTCGCGGAAGTTGGCGATCACATACCAGTGACGAGTGTCGATTAGGGTAAAAATAGGGCGCAt  >  1:80719/1‑93 (MQ=255)
                                                                                       tgtgCCTGAGCGAATATTTTTCAGATCGGTTTCGCGGAAGTTGGCGATCACATACCAGTGACGAGTGTCGATTAGGGTAAAAATAGGGCGCAt  <  2:1263487/93‑1 (MQ=255)
                                                                                          |                                                                                         
GTCATCCGGTAGCACGCCGTAGCCAATCGAATCCACTTTACCCTCGAAGGTTTTGCCGCTGTCACTCATCAGGCGAATCGTTGCGGGTGTACCTGAGCGAATATTTTTCAGATCGGTTTCGCGGAAGTTGGCGATCACATACCAGTGACGAGTGTCGATTAGGGTAAAAATAGGGCGCAT  >  NC_000913/4300813‑4300992

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