Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,957,551 +T intergenic (‑120/‑741) 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,5501.T98.5% 498.5 / ‑0.3 138intergenic (‑119/‑742)ppiC/reppeptidyl‑prolyl cis‑trans isomerase C/rep helicase
Reads supporting (aligned to +/- strand):  ref base . (2/0);  new base T (60/75);  total (63/75)
Fisher's exact test for biased strand distribution p-value = 2.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

Only 100 of 139 total aligned reads displayed.

TAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCTG‑‑CTCGCAGGTCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGT  >  NC_000913/3957470‑3957627
                                                                                     |                                                                               
tAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTcc                                                                     <  1:4646966/93‑1 (MQ=255)
   cATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTccctc                                                                  >  2:4558835/1‑90 (MQ=255)
    aTGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTccctca                                                                 <  1:3594481/93‑5 (MQ=255)
       aTCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAgg                                                               <  1:826261/93‑1 (MQ=255)
       aTCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAgg                                                               <  1:151896/93‑1 (MQ=255)
        tCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAggg                                                              >  2:4325677/1‑93 (MQ=255)
        tCCGCAGATGCTAACATAATTGTCAGTCAGGGAGTGGCTGGGACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAggg                                                              >  2:2895915/1‑93 (MQ=255)
         ccGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAgggg                                                             <  2:5923120/93‑1 (MQ=255)
          cGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGga                                                            <  2:3388717/93‑1 (MQ=255)
               aTGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGaaaa                                                       >  1:308986/1‑93 (MQ=255)
                tGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAAc                                                      >  2:4544786/1‑93 (MQ=255)
                tGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGGACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAAc                                                      >  1:2605310/1‑93 (MQ=255)
                 gCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACt                                                     >  1:2938992/1‑93 (MQ=255)
                 gCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CGTCGCAGGTCCCTCAAGGGGAGAAAACt                                                     >  2:1293094/1‑93 (MQ=255)
                    aaCGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAg                                                  <  1:1018440/93‑1 (MQ=255)
                    aaCGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAg                                                  >  1:1900277/1‑93 (MQ=255)
                    aaCGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAg                                                  <  2:1925912/93‑1 (MQ=255)
                     aCGTAATTGTCAGTCAGGGAGTCGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAgg                                                 <  1:3201179/93‑1 (MQ=255)
                       gTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGc                                               >  1:3249213/1‑93 (MQ=255)
                       gTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGc                                               <  2:1798259/93‑1 (MQ=255)
                       gTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCTCTCAAGGGGAGAAAACTTAGGGc                                               >  1:3443540/1‑93 (MQ=255)
                        tAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGcc                                              <  2:5958304/93‑1 (MQ=255)
                        tAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGcc                                              >  2:1670348/1‑93 (MQ=255)
                        tAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGcc                                              <  1:1056185/93‑1 (MQ=255)
                         aaTTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCct                                             <  2:6071234/93‑1 (MQ=255)
                          aTTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCctc                                            <  2:4959993/93‑1 (MQ=255)
                          aTTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCctc                                            <  2:1815027/93‑1 (MQ=255)
                            tGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTa                                          <  2:3506523/93‑1 (MQ=255)
                              tcagtcagGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATg                                        <  1:1129571/93‑1 (MQ=255)
                                agtcagGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATggg                                      >  2:828683/1‑93 (MQ=255)
                                agtcagGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATggg                                      >  2:1298989/1‑93 (MQ=255)
                                agtcagGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATggg                                      >  2:1197914/1‑93 (MQ=255)
                                agtcagGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATggg                                      <  1:521336/93‑1 (MQ=255)
                                agtcagGGAGTGGCTGACACATATATGGCAAAACCCCTGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATggg                                      <  2:2739448/93‑1 (MQ=255)
                                  tcagGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCt                                    >  2:1566262/1‑93 (MQ=255)
                                  tcagGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCt                                    >  1:4757245/1‑93 (MQ=255)
                                    agGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGGGAAAACTTAGGGCCCCCATggggtcc                                  >  2:1508563/1‑89 (MQ=255)
                                    agGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTAc                                  >  1:4052741/1‑93 (MQ=255)
                                    agGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTAc                                  <  2:2358848/93‑1 (MQ=255)
                                    agGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTAc                                  <  1:619340/93‑1 (MQ=255)
                                    agGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTAc                                  >  2:1508564/1‑93 (MQ=255)
                                     gggAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACa                                 >  2:1727812/1‑93 (MQ=255)
                                     gggAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACa                                 >  2:4555456/1‑93 (MQ=255)
                                      ggAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACaa                                <  2:3183266/93‑1 (MQ=255)
                                      ggAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCCTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACaa                                >  2:5165397/1‑93 (MQ=255)
                                       gAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACaaa                               <  2:585246/93‑1 (MQ=255)
                                       gAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACaaa                               <  1:5480111/93‑1 (MQ=255)
                                        aGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAg                              <  1:4886800/93‑1 (MQ=255)
                                        aGTGGCTGGTACATATAAGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAg                              <  1:1996586/93‑1 (MQ=255)
                                         gTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGt                             <  1:3299215/93‑1 (MQ=255)
                                          tGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTc                            >  1:4814745/1‑93 (MQ=255)
                                            gCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCAt                          <  1:2671362/93‑1 (MQ=255)
                                             cTGGTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCAtt                         <  1:907817/93‑1 (MQ=255)
                                              tGGTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTa                        >  1:1533937/1‑93 (MQ=255)
                                              tGGTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTa                        >  1:4211893/1‑93 (MQ=255)
                                              tGGTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTa                        <  1:5584220/93‑1 (MQ=255)
                                               ggTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTTTATGGGCTACAAAGTCATTAt                       <  1:4971884/93‑1 (MQ=255)
                                               ggTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTAt                       <  1:5958066/93‑1 (MQ=255)
                                               ggTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTAt                       <  2:2162068/93‑1 (MQ=255)
                                               ggTACATATATGGCAAAAGCCCGAAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTAt                       <  1:2561476/93‑1 (MQ=255)
                                                gTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATc                      >  1:255486/1‑93 (MQ=255)
                                                gTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATc                      >  2:1131429/1‑93 (MQ=255)
                                                gTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATc                      <  2:3238432/93‑1 (MQ=255)
                                                gTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATc                      >  1:3562583/1‑93 (MQ=255)
                                                gTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGAGCCTCTATGGGCTACAAAGTCATTATc                      >  2:235248/1‑93 (MQ=255)
                                                 tACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCg                     <  1:4456508/93‑1 (MQ=255)
                                                 tACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCg                     >  2:4773745/1‑93 (MQ=255)
                                                 tACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCg                     <  2:4320497/93‑1 (MQ=255)
                                                   cATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCggg                   >  2:2286091/1‑93 (MQ=255)
                                                   cATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCggg                   >  1:5961297/1‑93 (MQ=255)
                                                    atatatGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCgggg                  >  1:5810927/1‑93 (MQ=255)
                                                    atatatGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCgggg                  >  1:3787/1‑93 (MQ=255)
                                                    atatatGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCgggg                  >  2:802404/1‑93 (MQ=255)
                                                     tatatGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGa                 <  1:686916/93‑1 (MQ=255)
                                                     tatatGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGa                 <  1:5811190/93‑1 (MQ=255)
                                                     tatatGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGa                 <  1:3881312/93‑1 (MQ=255)
                                                     tatatGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGa                 >  2:260567/1‑93 (MQ=255)
                                                      atatGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGAt                <  2:3942384/93‑1 (MQ=255)
                                                      atagGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGAt                <  2:158199/93‑1 (MQ=255)
                                                       tatGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAACGTCATTATCGGGGATg               >  2:307469/1‑93 (MQ=255)
                                                       tatGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATg               <  2:5639335/93‑1 (MQ=255)
                                                         tGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGt             <  1:5699089/93‑1 (MQ=255)
                                                         tGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGt             <  2:651666/93‑1 (MQ=255)
                                                         tGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGt             <  1:6090494/93‑1 (MQ=255)
                                                           gCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGtta           <  1:870907/93‑1 (MQ=255)
                                                           gCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGtta           <  1:3366842/93‑1 (MQ=255)
                                                             aaaaGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGttatt         >  1:2762356/1‑93 (MQ=255)
                                                              aaaGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGttattt        <  2:417131/93‑1 (MQ=255)
                                                               aaGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGttatttt       <  2:3318719/93‑1 (MQ=255)
                                                               aaGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGttatttt       >  1:3605180/1‑93 (MQ=255)
                                                                aGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGttatttta      >  2:5262722/1‑93 (MQ=255)
                                                                 gCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGttattttat     >  2:3086145/1‑93 (MQ=255)
                                                                 gCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGttattttat     >  2:2094205/1‑93 (MQ=255)
                                                                 gCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGttattttat     <  2:1168370/93‑1 (MQ=255)
                                                                  cccGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGttattttatt    >  1:469397/1‑93 (MQ=255)
                                                                  cccGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGttattttatt    <  1:1553052/93‑1 (MQ=255)
                                                                   ccGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTg   >  1:560076/1‑93 (MQ=255)
                                                                   ccGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTg   <  2:1736995/93‑1 (MQ=255)
                                                                   ccGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTg   <  2:3791210/93‑1 (MQ=255)
                                                                    cGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGt  <  1:2959575/93‑1 (MQ=255)
                                                                                     |                                                                               
TAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCTG‑‑CTCGCAGGTCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGT  >  NC_000913/3957470‑3957627

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