Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,957,560 +C intergenic (‑129/‑732) 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,5581.C99.7% 1269.5 / ‑5.6 386intergenic (‑127/‑734)ppiC/reppeptidyl‑prolyl cis‑trans isomerase C/rep helicase
Reads supporting (aligned to +/- strand):  ref base . (1/0);  major base C (188/196);  minor base G (1/0);  total (190/196)
Fisher's exact test for biased strand distribution p-value = 4.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 423 total aligned reads displayed.

TAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCTGCTCGCAGGTCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATG  >  NC_000913/3957470‑3957645
                                                                                                |                                                                                         
 aGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTccc                                                                                       <  1:3522547/93‑2 (MQ=255)
  gCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTccct                                                                                      <  2:4651424/93‑3 (MQ=255)
  gCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTccct                                                                                      <  1:7895890/93‑3 (MQ=255)
  gCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGGACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTccct                                                                                      >  1:8714947/1‑91 (MQ=255)
   cATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTccctc                                                                                     <  1:4617034/93‑4 (MQ=255)
    aTGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTccctca                                                                                    >  2:3388512/1‑89 (MQ=255)
      gATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAg                                                                                  >  1:9662376/1‑93 (MQ=255)
       aTCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAgg                                                                                 <  2:10023356/93‑1 (MQ=255)
       aTCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAgg                                                                                 >  1:8183801/1‑93 (MQ=255)
        tCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGTCCGGGAAATCCGGGCTTCGCAGGTCCCTCAAggg                                                                                >  1:2050685/1‑93 (MQ=255)
        tCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAggg                                                                                >  2:7401021/1‑93 (MQ=255)
         ccGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCCGGCTTCGCCGGTCCCTCAAgggg                                                                               >  2:758161/1‑93 (MQ=255)
         ccGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAgggg                                                                               <  2:5127500/93‑1 (MQ=255)
          cGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGga                                                                              <  1:686728/93‑1 (MQ=255)
          cGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGGACATATATGGCAAAAGACCGGAAAATTCGGGCTTCGCAGGTCCCTCAAGGGga                                                                              >  2:3920551/1‑93 (MQ=255)
              gATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGGAAAAGCCCCGAAAATCCGGGCTTCGCAGGGCCCTCAAGGGGGGaaa                                                                          >  2:1940500/1‑93 (MQ=255)
              gATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGaaa                                                                          >  2:4050144/1‑93 (MQ=255)
              gATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGaaa                                                                          >  1:7651234/1‑93 (MQ=255)
                tGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCTCAGGTCCCTCAAGGGGAGAAAAc                                                                        >  2:7259791/1‑93 (MQ=255)
                tGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAAc                                                                        <  2:2569141/93‑1 (MQ=255)
                tGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGGACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGCAAAc                                                                        >  1:3730247/1‑93 (MQ=255)
                 gCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCCAGGGGAGAAAACt                                                                       >  1:7822958/1‑93 (MQ=255)
                 gCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACt                                                                       >  1:5967259/1‑93 (MQ=255)
                  cTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGGAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACtt                                                                      >  1:3096040/1‑93 (MQ=255)
                  cTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACtt                                                                      <  1:693321/93‑1 (MQ=255)
                   tAACGTAATTGTCAGTCGGGGAGTGGCGGGTACATATATAGCAAAAGCCAGGAAAACCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTa                                                                     <  1:1918711/93‑1 (MQ=255)
                   tAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTa                                                                     <  1:6626596/93‑1 (MQ=255)
                   tAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTa                                                                     <  1:6056807/93‑1 (MQ=255)
                   tAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTa                                                                     <  1:1108621/93‑1 (MQ=255)
                   tAACGTAATTGTCAGTCAGGGAGTGGCTGGGACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTa                                                                     >  1:10004361/1‑93 (MQ=255)
                   tAACGTAATTGTCAGTCAGGGAGTGGCTGGGACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTa                                                                     >  2:9078832/1‑93 (MQ=255)
                      cGTAATTGTCAGTCAGGGAGTGGCTGGGACCTTTATGGCAAAAGCCCGGGAAATCCGGGCTTCGCCGGGCTCTCAAGGGGAGCAAACCTCggg                                                                  >  1:7285211/1‑93 (MQ=255)
                       gTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGc                                                                 >  1:5285254/1‑93 (MQ=255)
                         aaTTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCct                                                               <  2:2927154/93‑1 (MQ=255)
                          aTTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCctc                                                              <  1:4040385/93‑1 (MQ=255)
                          aTTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCctc                                                              >  2:6701435/1‑93 (MQ=255)
                           ttGTCAGTCAGGGAGTGGCTGGGACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCctct                                                             >  1:8818324/1‑93 (MQ=255)
                            tGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTa                                                            <  1:9429408/93‑1 (MQ=255)
                            tGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTa                                                            <  1:5501688/93‑1 (MQ=255)
                             gtcagtcaGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTAt                                                           >  2:5358060/1‑93 (MQ=255)
                             gtcagtcaGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTAt                                                           >  2:7870572/1‑93 (MQ=255)
                             gtcagtcaGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTAt                                                           <  2:1473900/93‑1 (MQ=255)
                             gtcagtcaGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTAt                                                           <  1:3964854/93‑1 (MQ=255)
                               cagtcagGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATgg                                                         <  2:6958507/93‑1 (MQ=255)
                                 gtcagGGAGTGGCTGGTACATATATGGCAAAAGCCCGGGAAAACCGGGCTTCGCAGGTCCCTCAAGGGGAGGAAACTTAGGGCCTCTATGGGc                                                       >  2:1262187/1‑93 (MQ=255)
                                  tcagGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCt                                                      >  1:2795195/1‑93 (MQ=255)
                                  tcagGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCt                                                      >  1:7426141/1‑93 (MQ=255)
                                    agGGAGTGGCTGGTACATATATTGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCCAAGGGAGAAAACTTAGGGCCTCTATGGGCTAc                                                    >  2:9589144/1‑93 (MQ=255)
                                    agGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTAc                                                    <  1:9195212/93‑1 (MQ=255)
                                     gggAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACa                                                   >  1:7201277/1‑93 (MQ=255)
                                        aGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAg                                                >  2:6343641/1‑93 (MQ=255)
                                         gTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGt                                               >  2:9851498/1‑93 (MQ=255)
                                          tGGCTGGTACATATATGGCAAAAGCGCGGAAAATCCGGGCTTCGCAGGTCCCTCAGGGGGAGAAACCTTAGGGCCTCTATGGGCTACAAAGTc                                              <  1:8475989/93‑1 (MQ=255)
                                          tGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTc                                              >  1:8584216/1‑93 (MQ=255)
                                           ggCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCa                                             <  1:6353632/93‑1 (MQ=255)
                                           ggCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCa                                             <  1:5045612/93‑1 (MQ=255)
                                           ggCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCa                                             >  2:4872432/1‑93 (MQ=255)
                                             cTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCAtt                                           >  1:1408765/1‑93 (MQ=255)
                                               ggTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTAt                                         >  2:6039658/1‑93 (MQ=255)
                                               ggTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTAt                                         <  2:7739029/93‑1 (MQ=255)
                                               ggTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTAt                                         >  1:2965123/1‑93 (MQ=255)
                                                gTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATc                                        <  1:7494616/93‑1 (MQ=255)
                                                 tACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCg                                       >  1:9564880/1‑93 (MQ=255)
                                                 tACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCg                                       >  2:1718760/1‑93 (MQ=255)
                                                  aCATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCgg                                      <  1:9814052/93‑1 (MQ=255)
                                                  aCATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCgg                                      <  1:2141926/93‑1 (MQ=255)
                                                   cATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCggg                                     >  1:4792785/1‑93 (MQ=255)
                                                     tatatGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCggggt                                   >  2:8522725/1‑92 (MQ=255)
                                                      atatGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGAt                                  <  2:6280978/93‑1 (MQ=255)
                                                        atGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATgg                                >  2:810853/1‑93 (MQ=255)
                                                        atGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATgg                                >  2:1207331/1‑93 (MQ=255)
                                                         tGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGt                               <  2:8432336/93‑1 (MQ=255)
                                                           gCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGtta                             <  2:4357849/93‑1 (MQ=255)
                                                            cAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGttat                            >  2:5567846/1‑93 (MQ=255)
                                                             aaaaGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGttatt                           >  2:367124/1‑93 (MQ=255)
                                                              aaaGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGttattt                          >  2:8001115/1‑93 (MQ=255)
                                                              aaaGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGttattt                          >  2:5855670/1‑93 (MQ=255)
                                                              aaaGCCCGGAAAATCCGGGCTTCGCAGGTCCATCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGttattt                          >  2:4995059/1‑93 (MQ=255)
                                                               aaGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGttatttt                         >  2:2952948/1‑93 (MQ=255)
                                                                 gCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGttattttat                       <  2:6148915/93‑1 (MQ=255)
                                                                  cccGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGttattttatt                      <  2:4268401/93‑1 (MQ=255)
                                                                   ccGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTg                     <  2:5775253/93‑1 (MQ=255)
                                                                   ccGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTg                     <  1:3258836/93‑1 (MQ=255)
                                                                   ccGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTg                     >  1:7140836/1‑93 (MQ=255)
                                                                                   ggcttCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAg         <  1:4814537/89‑1 (MQ=255)
                                                                                   ggcttCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAg         >  2:4389236/5‑93 (MQ=255)
                                                                                     cttCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCg       >  1:170937/3‑93 (MQ=255)
                                                                                      ttCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGa      <  1:4928346/92‑1 (MQ=255)
                                                                                      ttCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGa      <  1:2556450/92‑1 (MQ=255)
                                                                                       tCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGat     <  2:9319051/93‑1 (MQ=255)
                                                                                       tCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGat     >  1:1012311/1‑93 (MQ=255)
                                                                                        cGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGTGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGata    <  1:1618346/93‑1 (MQ=255)
                                                                                        cGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGata    <  1:5567847/93‑1 (MQ=255)
                                                                                        cGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGata    <  1:265884/93‑1 (MQ=255)
                                                                                        cGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGata    <  2:9892381/93‑1 (MQ=255)
                                                                                         gCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGatat   >  2:6465832/1‑93 (MQ=255)
                                                                                         gCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGatat   >  2:2478643/1‑93 (MQ=255)
                                                                                         gCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGatat   <  2:8899620/93‑1 (MQ=255)
                                                                                         gCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGatat   <  1:3348081/93‑1 (MQ=255)
                                                                                          cAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATg  <  1:4473386/93‑1 (MQ=255)
                                                                                                |                                                                                         
TAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCTGCTCGCAGGTCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATG  >  NC_000913/3957470‑3957645

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