Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,957,548 1 bp→GG intergenic (‑117/‑744) ppiC ← / → rep peptidyl‑prolyl cis‑trans isomerase C/rep helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,957,5471.G95.1% 375.9 / 11.9 122intergenic (‑116/‑745)ppiC/reppeptidyl‑prolyl cis‑trans isomerase C/rep helicase
Reads supporting (aligned to +/- strand):  ref base . (1/5);  new base G (54/62);  total (55/67)
Fisher's exact test for biased strand distribution p-value = 2.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.05e-01
*NC_0009133,957,5480TG99.2% 389.5 / NA 123intergenic (‑117/‑744)ppiC/reppeptidyl‑prolyl cis‑trans isomerase C/rep helicase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  major base G (54/66);  minor base A (0/1);  total (55/68)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 128 total aligned reads displayed.

CAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCTGCTCGCAGGTCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTG  >  NC_000913/3957457‑3957626
                                                                                            |||                                                                                  
cAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATTGCCAAAGCCCGGAAAATCCTg                                                                                    >  2:4557639/1‑93 (MQ=255)
              aGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTccc                                                                   <  2:596648/93‑2 (MQ=255)
              aGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTccc                                                                   <  1:3528767/93‑2 (MQ=255)
              aGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTccc                                                                   >  1:2983036/1‑92 (MQ=255)
                 aTGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGGACATATATGGCAAAAGCCCGGAAAATCCGCGCTTCGCAGGTccctca                                                                >  2:2461867/1‑89 (MQ=255)
                   gATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAg                                                               <  2:4308653/93‑1 (MQ=255)
                   gATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAg                                                               >  2:1798487/1‑93 (MQ=255)
                     tCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAggg                                                             >  1:719859/1‑93 (MQ=255)
                     tCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAggg                                                             <  2:3044564/93‑1 (MQ=255)
                     tCCGAAGATGCTAACGTATTTGTCAGTCAGGGAGTGGCTGGAACATATATGGAAAAACCCGGGAAATTCCGGGCTTCGCAGGTCCCTCAAggg                                                             <  2:4234829/93‑1 (MQ=255)
                      ccGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAgggg                                                            <  1:4260472/93‑1 (MQ=255)
                       cGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGga                                                           <  1:3895607/93‑1 (MQ=255)
                       cGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGga                                                           <  2:3104340/93‑1 (MQ=255)
                        gCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGgag                                                          <  2:4484044/93‑1 (MQ=255)
                          aGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGaa                                                        <  1:1396783/93‑1 (MQ=255)
                          aGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGaa                                                        >  1:2237251/1‑93 (MQ=255)
                           gATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATTGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGaaa                                                       >  2:2592414/1‑93 (MQ=255)
                           gATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCCGGTCCCTCAAGGGGAGaaa                                                       >  1:4350201/1‑93 (MQ=255)
                           gATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGaaa                                                       >  1:2899088/1‑93 (MQ=255)
                             tGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGGAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCCCCAGGGGAGAAAAc                                                     >  2:1450562/1‑93 (MQ=255)
                             tGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAAc                                                     >  2:2159224/1‑93 (MQ=255)
                                tAACGTAATTGTCAGACAGGAAATGGCTAGTACATATATGGCAAAAGCCAGGAAAATCCGAGCTTCGCAGGTCACTCAAGGGAAGAAAACTTa                                                  <  1:4196628/93‑1 (MQ=255)
                                 agcGTAATGGTAAGTCAGGGAGTGGCTGGTACATATAGGGCAAAAGCCGGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAg                                                 <  2:1742699/91‑1 (MQ=255)
                                 aaCGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAg                                                 >  1:3451878/1‑93 (MQ=255)
                                 aaCGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAg                                                 <  2:4744306/93‑1 (MQ=255)
                                  aCGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAgg                                                <  2:1288017/93‑1 (MQ=255)
                                   cGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAggg                                               <  2:3673431/93‑1 (MQ=255)
                                    gTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGc                                              >  2:4378360/1‑93 (MQ=255)
                                      aaTTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCct                                            >  2:2989376/1‑93 (MQ=255)
                                       aTTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCctc                                           >  1:4702288/1‑93 (MQ=255)
                                       aTTGGCAGTCAGGGAGTGGCTGGGACATATATTTCAAAAGCCCGGGAAATCCGGGCTTCGCAGGGCACTCACGGGGAGAAAAATTAGGGCctc                                           >  2:1737049/1‑93 (MQ=255)
                                         tGTCAGTCAGGGAGTGGCTGGTACATGTATGGCAAAAGCCGGGAAAATCCGGGCTTCGCAGGGCCCTCAAGGGGAGAAAACTTAGGGCCTCTa                                         <  2:282603/93‑1 (MQ=255)
                                         tGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTa                                         <  2:2779441/93‑1 (MQ=255)
                                          gtcagtcaGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTAt                                        <  2:440679/93‑1 (MQ=255)
                                          gtcagtcaGGGAGTGCCGGGTACATATAGGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTAt                                        <  1:3541469/93‑1 (MQ=255)
                                           tcagtcagGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATg                                       <  2:1909374/93‑1 (MQ=255)
                                            cagtcagGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTAtcg                                      >  2:232461/1‑91 (MQ=255)
                                              gtcagGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCGAGGGGAGAAAACTTAGGGCCTCTATGGGc                                    >  2:4689933/1‑93 (MQ=255)
                                               tcagGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCt                                   >  2:754053/1‑93 (MQ=255)
                                                cagGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTa                                  >  1:1629085/1‑93 (MQ=255)
                                                 agGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTAc                                 >  2:4616872/1‑93 (MQ=255)
                                                  gggAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACa                                >  2:2705324/1‑93 (MQ=255)
                                                  gggAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACa                                <  1:992955/93‑1 (MQ=255)
                                                   ggAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACaa                               <  2:3403299/93‑1 (MQ=255)
                                                   ggAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACaa                               <  2:3403298/93‑1 (MQ=255)
                                                   ggAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACaa                               >  1:2598603/1‑93 (MQ=255)
                                                    gAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACaaa                              >  1:4074900/1‑93 (MQ=255)
                                                    gAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACaaa                              >  2:2403508/1‑93 (MQ=255)
                                                    gAGTGACTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACaaa                              <  1:1126539/93‑1 (MQ=255)
                                                      gTGGCTGGTCCATATATGGCAAAAGCCCGGAAAAACCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAAGGCCTCTATGGGATACAAAGt                            >  2:2012928/1‑93 (MQ=255)
                                                      gTGGCTGGTACATATATGGCAAAAGCCGGGAAAATCCGGTCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGt                            <  1:4255012/93‑1 (MQ=255)
                                                      gTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGt                            <  1:1397864/93‑1 (MQ=255)
                                                       tGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTc                           <  1:1024811/93‑1 (MQ=255)
                                                       tGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTc                           <  2:2284363/93‑1 (MQ=255)
                                                       tGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTc                           >  1:777776/1‑93 (MQ=255)
                                                          cTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCAtt                        >  1:1481620/1‑93 (MQ=255)
                                                          cTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCAtt                        <  1:2086089/93‑1 (MQ=255)
                                                          cTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCAtt                        >  2:1520120/1‑93 (MQ=255)
                                                           tGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTa                       >  2:3300711/1‑93 (MQ=255)
                                                           tGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTa                       <  1:953571/93‑1 (MQ=255)
                                                           tGGAACATATATGGCAAACGCCAGAAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTa                       <  2:63370/93‑1 (MQ=255)
                                                            ggTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTAt                      >  1:3301713/1‑93 (MQ=255)
                                                            ggTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTAt                      >  1:3766076/1‑93 (MQ=255)
                                                            ggTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTAt                      <  2:4749242/93‑1 (MQ=255)
                                                            ggTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTAt                      <  1:3480991/93‑1 (MQ=255)
                                                            ggTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTAt                      <  2:3307059/93‑1 (MQ=255)
                                                            ggTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTAt                      >  1:508741/1‑93 (MQ=255)
                                                             gTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATc                     <  2:4680445/93‑1 (MQ=255)
                                                             gTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATc                     <  2:2348704/93‑1 (MQ=255)
                                                             gTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATc                     <  1:2352706/93‑1 (MQ=255)
                                                              tACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGTGCTACAAAGTCATTATCg                    >  2:2970955/1‑93 (MQ=255)
                                                              tACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGTGCTACAAAGTCATTATCg                    >  2:2970956/1‑93 (MQ=255)
                                                              tACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCg                    <  1:2119259/93‑1 (MQ=255)
                                                              tACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCg                    >  1:105707/1‑93 (MQ=255)
                                                                cATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCggg                  <  1:837688/93‑1 (MQ=255)
                                                                cATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCggg                  <  1:2074688/93‑1 (MQ=255)
                                                                  tatatGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGa                >  1:875851/1‑93 (MQ=255)
                                                                  tatatGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGa                >  1:1444537/1‑93 (MQ=255)
                                                                  tatatGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGa                >  2:4734310/1‑93 (MQ=255)
                                                                  tatatGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGa                >  2:1542904/1‑93 (MQ=255)
                                                                  tatatGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGa                <  2:2460505/93‑1 (MQ=255)
                                                                   atatGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGAt               <  1:1144018/93‑1 (MQ=255)
                                                                   atatGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGAt               >  1:1822942/1‑93 (MQ=255)
                                                                    tatGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATg              <  1:744283/93‑1 (MQ=255)
                                                                     atGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATgg             <  2:2005597/93‑1 (MQ=255)
                                                                     atGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATgg             <  1:2281113/93‑1 (MQ=255)
                                                                      tGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGCAAACTTAGGGCCTCTATTGGCTACAAAGTCATTATCGGGGATGGt            >  1:4656259/1‑93 (MQ=255)
                                                                      tGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGt            <  1:164173/93‑1 (MQ=255)
                                                                       ggCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGtt           <  1:1461593/93‑1 (MQ=255)
                                                                       ggCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGtt           <  2:326793/93‑1 (MQ=255)
                                                                       ggCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGtt           <  2:1647124/93‑1 (MQ=255)
                                                                          aaaaGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGttatt        >  1:1413870/1‑93 (MQ=255)
                                                                           aaaGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGttattt       >  1:2033833/1‑93 (MQ=255)
                                                                            aaGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGttatttt      >  1:2911528/1‑93 (MQ=255)
                                                                             aGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGttatttta     <  2:537388/93‑1 (MQ=255)
                                                                              gCCTGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGttattttat    >  1:1179011/1‑93 (MQ=255)
                                                                              gCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGttattttat    <  1:2275320/93‑1 (MQ=255)
                                                                              gCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGttattttat    >  1:3174822/1‑93 (MQ=255)
                                                                               cccGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGttattttatt   <  2:3251590/93‑1 (MQ=255)
                                                                                ccGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTg  >  1:327859/1‑93 (MQ=255)
                                                                                            |||                                                                                  
CAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCTGCTCGCAGGTCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTG  >  NC_000913/3957457‑3957626

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