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.G94.7% 854.4 / 38.8 285intergenic (‑116/‑745)ppiC/reppeptidyl‑prolyl cis‑trans isomerase C/rep helicase
Reads supporting (aligned to +/- strand):  ref base . (9/6);  new base G (132/138);  total (141/144)
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.36e-02
*NC_0009133,957,5480TG100.0% 902.8 / NA 285intergenic (‑117/‑744)ppiC/reppeptidyl‑prolyl cis‑trans isomerase C/rep helicase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (141/144);  total (141/144)

Only 100 of 299 total aligned reads displayed.

AGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCTGCTCGCAGGTCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGT  >  NC_000913/3957471‑3957627
                                                                               |||                                                                                 
aGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTccc                                                                    <  2:3198887/93‑2 (MQ=255)
  cATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGATAATCCGGGCTTCGCAGGTccctc                                                                  >  2:5644126/1‑90 (MQ=255)
     gATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAg                                                                <  2:3852705/93‑1 (MQ=255)
      aTCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAgg                                                               <  1:3947108/93‑1 (MQ=255)
      aTCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGGACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAgg                                                               >  2:3378234/1‑93 (MQ=255)
       tCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAggg                                                              <  1:3178266/93‑1 (MQ=255)
       tCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGGACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAggg                                                              >  2:5975651/1‑93 (MQ=255)
         cGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGga                                                            >  1:6675462/1‑93 (MQ=255)
          gCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGgag                                                           >  2:364921/1‑93 (MQ=255)
           cAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGgaga                                                          <  2:1012947/93‑1 (MQ=255)
            aGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGaa                                                         <  2:7488492/93‑1 (MQ=255)
              aTGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGaaaa                                                       <  2:4461026/93‑1 (MQ=255)
               tGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAAc                                                      >  2:2373865/1‑93 (MQ=255)
               tGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAAc                                                      >  1:7340250/1‑93 (MQ=255)
                gCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACt                                                     >  1:4931482/1‑93 (MQ=255)
                gCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACt                                                     >  1:4734322/1‑93 (MQ=255)
                 cTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATTGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACtt                                                    >  2:5280072/1‑93 (MQ=255)
                 cTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACtt                                                    >  1:1016958/1‑93 (MQ=255)
                 cTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACtt                                                    <  2:6397098/93‑1 (MQ=255)
                  tAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTa                                                   >  2:6885328/1‑93 (MQ=255)
                  tAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTa                                                   <  2:1915729/93‑1 (MQ=255)
                  caaCGTAATTGTCAGACAAGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTCCGCAGGTCCCTCAAGGGGAGAAAACTTa                                                   <  1:6379891/92‑1 (MQ=255)
                   agcGTAATTGTCAGTCAGGGAGTGGCTGTTACATATATGGCAAAAGCCCGGAAAAACCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAg                                                  <  2:4695108/91‑1 (MQ=255)
                   aaCGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAg                                                  <  2:6156502/93‑1 (MQ=255)
                   aaCGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAg                                                  <  1:7862306/93‑1 (MQ=255)
                   aaCGTAATTGTCAGTCAGGGAGCGGCTGGGGCATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAg                                                  >  2:331943/1‑93 (MQ=255)
                    aCGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAgg                                                 <  2:2526323/93‑1 (MQ=255)
                    aCGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAgg                                                 >  1:2266842/1‑93 (MQ=255)
                    aCGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAgg                                                 >  2:414246/1‑93 (MQ=255)
                     cGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAAATTAggg                                                >  2:3024452/1‑93 (MQ=255)
                      gTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGc                                               >  2:6954633/1‑93 (MQ=255)
                      gTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGc                                               <  1:2834348/93‑1 (MQ=255)
                       tAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGcc                                              <  1:823182/93‑1 (MQ=255)
                        aaTTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCct                                             >  1:4623059/1‑93 (MQ=255)
                        aaTTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCct                                             >  1:5152088/1‑93 (MQ=255)
                         aTTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCctc                                            <  1:1187606/93‑1 (MQ=255)
                          ttGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCCGGCTTCGCAGGTCCCTCCAGGGGTGAAAACTTAGGGCctct                                           >  2:904281/1‑93 (MQ=255)
                          ttGTCAGTCAGGGAGTGGCTGGGACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCctct                                           >  1:1395017/1‑93 (MQ=255)
                           tGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTa                                          <  1:2973049/93‑1 (MQ=255)
                            gtcagtcaGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTAt                                         <  1:6884597/93‑1 (MQ=255)
                            gtcagtcaGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTAt                                         <  2:1025249/93‑1 (MQ=255)
                             tcagtcagGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCCAGGGGAGAAAACTTAGGGCCTCTTTg                                        >  1:6958774/1‑93 (MQ=255)
                             tcagtcagGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATg                                        <  2:4414582/93‑1 (MQ=255)
                             tcagtcagGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATg                                        <  2:705691/93‑1 (MQ=255)
                             tcagtcagGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATg                                        <  1:3865050/93‑1 (MQ=255)
                             tcagtcagGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATg                                        <  2:6714198/93‑1 (MQ=255)
                             tcagtcagGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATg                                        <  2:5383523/93‑1 (MQ=255)
                               agtcagGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATggg                                      >  2:2477588/1‑93 (MQ=255)
                                gtcagGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGc                                     <  1:5519262/93‑1 (MQ=255)
                                 tcagGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCt                                    >  1:1697314/1‑93 (MQ=255)
                                 tcagGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCt                                    <  2:6523225/93‑1 (MQ=255)
                                 tcagGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCt                                    >  1:598675/1‑93 (MQ=255)
                                  cagGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTa                                   <  2:7535113/93‑1 (MQ=255)
                                   agGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTAc                                  >  1:63430/1‑93 (MQ=255)
                                   agGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTAc                                  >  2:7179247/1‑93 (MQ=255)
                                    gggAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACa                                 <  2:226591/93‑1 (MQ=255)
                                    gggAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACa                                 >  2:2224101/1‑93 (MQ=255)
                                    gggAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACa                                 >  2:32648/1‑93 (MQ=255)
                                    gggAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACa                                 >  1:2042940/1‑93 (MQ=255)
                                     ggAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACaa                                <  1:7429987/93‑1 (MQ=255)
                                      gAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACaaa                               <  2:1215628/93‑1 (MQ=255)
                                      gAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACaaa                               <  2:2873017/93‑1 (MQ=255)
                                         tGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTc                            >  1:5772427/1‑93 (MQ=255)
                                          ggCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCa                           <  2:1539195/93‑1 (MQ=255)
                                           gCTGGTACATATATGGTAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCAt                          <  2:5179412/93‑1 (MQ=255)
                                           gCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCAt                          <  2:5701422/93‑1 (MQ=255)
                                           gCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCAt                          >  1:5268562/1‑93 (MQ=255)
                                            cTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCAtt                         <  1:1307521/93‑1 (MQ=255)
                                            cTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCAtt                         <  1:448682/93‑1 (MQ=255)
                                             tGGTACATATATGGCAAAAGCCCGGGAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTa                        <  2:5135937/93‑1 (MQ=255)
                                             tGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTa                        <  1:2390388/93‑1 (MQ=255)
                                              ggTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTAt                       <  1:7142401/93‑1 (MQ=255)
                                              ggTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTAt                       <  1:6165937/93‑1 (MQ=255)
                                                tACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCg                     <  2:6243400/93‑1 (MQ=255)
                                                tACATATACGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCg                     <  1:5086729/93‑1 (MQ=255)
                                                 aCATATATGGAAAAAGCCGGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCgg                    <  2:2836171/93‑1 (MQ=255)
                                                  cATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCggg                   >  1:7687587/1‑93 (MQ=255)
                                                   atatatGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCgggg                  >  2:2760747/1‑93 (MQ=255)
                                                    tatatGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGa                 >  2:1704657/1‑93 (MQ=255)
                                                    tatatGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGa                 >  1:3093287/1‑93 (MQ=255)
                                                      tatGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATg               >  1:1103232/1‑93 (MQ=255)
                                                      tatGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATg               >  1:3033570/1‑93 (MQ=255)
                                                        tGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGt             <  2:5548878/93‑1 (MQ=255)
                                                        tGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGt             >  2:7609268/1‑93 (MQ=255)
                                                          gCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGtta           >  2:4933352/1‑93 (MQ=255)
                                                           cAAAAGCCGGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGttat          <  1:5569400/93‑1 (MQ=255)
                                                           cAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGttat          <  1:4368706/93‑1 (MQ=255)
                                                            aaaaGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGttatt         <  2:2131257/93‑1 (MQ=255)
                                                            aaaaGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGttatt         >  1:7402266/1‑93 (MQ=255)
                                                            aaaaGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGttatt         >  1:1146561/1‑93 (MQ=255)
                                                            aaaaGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCACTATCGGGGATGGttatt         >  2:7013300/1‑93 (MQ=255)
                                                             aaaGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGttattt        >  1:4037987/1‑93 (MQ=255)
                                                             aaaGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGttattt        >  2:3104555/1‑93 (MQ=255)
                                                               aGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGttatttta      <  2:4901037/93‑1 (MQ=255)
                                                                gCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGttattttat     >  2:6476595/1‑93 (MQ=255)
                                                                 cccGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGttattttatt    <  2:1144544/93‑1 (MQ=255)
                                                                  ccGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTg   >  2:4046988/1‑93 (MQ=255)
                                                                  ccGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTg   >  1:2567450/1‑93 (MQ=255)
                                                                  ccGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTg   >  1:1568677/1‑93 (MQ=255)
                                                                  ccGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTg   <  2:7837778/93‑1 (MQ=255)
                                                                               |||                                                                                 
AGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCTGCTCGCAGGTCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGT  >  NC_000913/3957471‑3957627

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