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.T99.4% 654.5 / ‑3.8 177intergenic (‑119/‑742)ppiC/reppeptidyl‑prolyl cis‑trans isomerase C/rep helicase
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base T (77/99);  total (78/99)
Fisher's exact test for biased strand distribution p-value = 4.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.16e-01

Only 100 of 187 total aligned reads displayed.

TAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCTGCTCGCAGGTCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTG  >  NC_000913/3957470‑3957626
                                                                                      |                                                                                
tAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTcc                                                                     <  2:969172/93‑1 (MQ=255)
tAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTcc                                                                     >  1:3915776/1‑93 (MQ=255)
 aGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTccc                                                                    <  1:1155389/93‑2 (MQ=255)
 aGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTccc                                                                    <  1:3983779/93‑2 (MQ=255)
    aTGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTccctca                                                                 >  2:1048148/1‑89 (MQ=255)
    aTGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGGACATATATGGCAAAAGCCCGGAAAAACCGGGCTTCGCAGGTccctca                                                                 >  1:1354221/1‑89 (MQ=255)
      gATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAg                                                               >  1:2021806/1‑93 (MQ=255)
      gATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAg                                                               >  2:759493/1‑93 (MQ=255)
      gATCCGCAGAGGCTACCGTAATTGTCAGTCAGGGAGTGGCTGGTACATATAGGGCAAAAGCCCGGAAAAGCCGGGCTTCGCAGGTCCCTCAAg                                                               <  2:3701423/93‑1 (MQ=255)
       aTCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAgg                                                              >  2:1247691/1‑93 (MQ=255)
         ccGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAgggg                                                            >  2:865979/1‑93 (MQ=255)
         ccGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAgggg                                                            <  2:4565138/93‑1 (MQ=255)
         ccGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAgggg                                                            >  2:3143894/1‑93 (MQ=255)
          cGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGga                                                           <  1:1627949/93‑1 (MQ=255)
           gCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGgag                                                          >  1:1817352/1‑93 (MQ=255)
             aGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGaa                                                        <  1:2513073/93‑1 (MQ=255)
             aGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGGACATATATGGCAAAAGGCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAggggggta                                                        >  1:5871785/1‑89 (MQ=255)
              gATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGaaa                                                       <  2:1386428/93‑1 (MQ=255)
                  cTACCGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCGGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACtt                                                   <  2:110304/93‑1 (MQ=255)
                   tAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGGCCCTCAAGGGGAGAAAACTTa                                                  >  2:1320988/1‑93 (MQ=255)
                   tAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTa                                                  >  2:6646005/1‑93 (MQ=255)
                   tAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTa                                                  <  1:5144386/93‑1 (MQ=255)
                    aaCGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGTGAAAACTTAg                                                 <  1:4930401/93‑1 (MQ=255)
                    aaCGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAg                                                 >  2:808080/1‑93 (MQ=255)
                    aaCGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAg                                                 <  1:4410629/93‑1 (MQ=255)
                    aaCGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAg                                                 <  2:3972214/93‑1 (MQ=255)
                    aaCGTAATTGTCAGTCAGGGAGTGGCTGGGACATATATGGCAAAAACCCGGAAAATCCCGGCTTCGCAGGGCCCTCAAGGGGAGGAAACCTAg                                                 >  1:4331500/1‑93 (MQ=255)
                     aCGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAgg                                                >  2:6653436/1‑93 (MQ=255)
                     aCGTAATTGTAAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCGGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAgg                                                <  2:3423881/93‑1 (MQ=255)
                       gTAATTGTCAGTCAGGGAGTGGCTGGTACATATCTGGCAACAGCCCCGAAAATCCGGGCTTCGCAGGGCGCTCAAGGGGAGAAAACTTAAGGc                                              >  1:4687706/1‑93 (MQ=255)
                       gTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGc                                              >  1:3481259/1‑93 (MQ=255)
                       gTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCAGGAAAAACAGGGCTTCGCAGGTGCCTCAAGGGGAGTCAACTTAGGGc                                              >  2:6321610/1‑93 (MQ=255)
                        tAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGcc                                             <  2:1515009/93‑1 (MQ=255)
                          aTTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCctc                                           <  1:3218724/93‑1 (MQ=255)
                          aTTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCctc                                           <  2:362953/93‑1 (MQ=255)
                          aTTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCctc                                           <  2:1989253/93‑1 (MQ=255)
                            tGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTa                                         >  1:5683762/1‑93 (MQ=255)
                             gtcagtcaGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTAt                                        <  2:699866/93‑1 (MQ=255)
                              tcagtcagGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATg                                       >  2:4454262/1‑93 (MQ=255)
                              tcagtcagGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATg                                       <  2:6693520/93‑1 (MQ=255)
                               tagtcagGCAGTGGCTGGTACATATAAGGCAAAAACCCGTAAAATCCGGGCTTCGCAGGTCCCTCAAGAGGAGAAAACTTAGGGCCTCTATgg                                      <  1:5029746/92‑1 (MQ=255)
                               cagtcagGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCATCAAGGGGAGAAAACTTAGGGCCTCTATgg                                      <  2:2120300/93‑1 (MQ=255)
                                agtcagGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATggg                                     <  2:4294687/93‑1 (MQ=255)
                                 ggcagGGAGTGGCTGGTACATAAATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGc                                    <  1:4158017/91‑1 (MQ=255)
                                  tcagGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCt                                   >  1:1038896/1‑93 (MQ=255)
                                   cagGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTa                                  <  2:4107296/93‑1 (MQ=255)
                                    agGGAGTGGCTGGTACATATATGGCAAAAGCCCGGGAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTAc                                 >  1:6448225/1‑93 (MQ=255)
                                    agGGAGTGGCTGGGACATATATGGCAAAAGCCCGGAAAATCCTGGCTTCGCAGGTCCCTCAAGGGGGGCAAACTTAGGGCCTCTATGGGCTAc                                 >  2:3300764/1‑93 (MQ=255)
                                     gggAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACa                                >  1:3781667/1‑93 (MQ=255)
                                      ggAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACaa                               >  2:3617409/1‑93 (MQ=255)
                                      ggAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACaa                               <  2:418332/93‑1 (MQ=255)
                                      ggAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACaa                               >  2:4863574/1‑93 (MQ=255)
                                       gAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACaaa                              <  2:5109465/93‑1 (MQ=255)
                                       gAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACaaa                              <  2:6489889/93‑1 (MQ=255)
                                       gAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACaaa                              <  2:153785/93‑1 (MQ=255)
                                       gAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACaaa                              >  2:6298493/1‑93 (MQ=255)
                                        aGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAg                             >  1:3079269/1‑93 (MQ=255)
                                         gTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGt                            <  1:4039461/93‑1 (MQ=255)
                                         gTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGt                            >  1:5409305/1‑93 (MQ=255)
                                         gTGGCTGGGACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGt                            >  1:6611024/1‑93 (MQ=255)
                                          tGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTc                           >  1:6345480/1‑93 (MQ=255)
                                          tGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTc                           >  1:2157087/1‑93 (MQ=255)
                                             cTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCAtt                        <  1:6265014/93‑1 (MQ=255)
                                              tGGTACATATATGGCAAAAGCCGGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTa                       <  2:6158199/93‑1 (MQ=255)
                                              tGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTa                       >  1:5335640/1‑93 (MQ=255)
                                              tGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTa                       <  2:2427798/93‑1 (MQ=255)
                                              tGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTa                       <  1:4620272/93‑1 (MQ=255)
                                               ggTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTAt                      <  1:563064/93‑1 (MQ=255)
                                               ggTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTAt                      <  1:4504490/93‑1 (MQ=255)
                                               ggTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTAt                      <  1:2707573/93‑1 (MQ=255)
                                                 tACATATATGGCAAAAGCCGGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCg                    <  1:1314297/93‑1 (MQ=255)
                                                 tACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCg                    <  2:2910676/93‑1 (MQ=255)
                                                 tACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCg                    <  2:3684660/93‑1 (MQ=255)
                                                  aCATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCgg                   <  2:4838305/93‑1 (MQ=255)
                                                  aCATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCgg                   <  1:1076434/93‑1 (MQ=255)
                                                   cATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCggg                  <  1:4948776/93‑1 (MQ=255)
                                                   cATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCggg                  >  2:463073/1‑93 (MQ=255)
                                                     tatatGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCggggt                >  2:2214145/1‑92 (MQ=255)
                                                     tatatGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGa                >  1:4841110/1‑93 (MQ=255)
                                                      atatGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGAt               >  1:3395212/1‑93 (MQ=255)
                                                      atatGGCAAAAGCCCGGAAAATCCGGTCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGAt               <  1:842305/93‑1 (MQ=255)
                                                       tatGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATg              <  2:2575804/93‑1 (MQ=255)
                                                       tatGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATg              <  1:2072800/93‑1 (MQ=255)
                                                       tatGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATg              <  1:376323/93‑1 (MQ=255)
                                                        atGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATgg             >  1:5286981/1‑93 (MQ=255)
                                                         tGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGt            >  2:6009949/1‑93 (MQ=255)
                                                         tGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGt            <  2:5606876/93‑1 (MQ=255)
                                                         tGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGt            >  2:4964314/1‑93 (MQ=255)
                                                           gCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGtta          >  1:3624222/1‑93 (MQ=255)
                                                            cAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGttat         <  2:5518218/93‑1 (MQ=255)
                                                            cAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGttat         <  2:5078529/93‑1 (MQ=255)
                                                             aaaaGCCCGGAAAATCGGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGttatt        <  1:3583311/93‑1 (MQ=255)
                                                             aaaaGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGttatt        >  1:4783063/1‑93 (MQ=255)
                                                              aaaGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGttattt       >  2:1325332/1‑93 (MQ=255)
                                                              aaaGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGttattt       >  1:3299876/1‑93 (MQ=255)
                                                              aaaGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGttattt       >  1:4344103/1‑93 (MQ=255)
                                                                  cccGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGttattttatt   <  1:6431476/93‑1 (MQ=255)
                                                                   ccGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTg  >  2:2761726/1‑93 (MQ=255)
                                                                   ccGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTg  >  2:5415787/1‑93 (MQ=255)
                                                                   ccGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTg  >  1:3249338/1‑93 (MQ=255)
                                                                                      |                                                                                
TAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCTGCTCGCAGGTCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTG  >  NC_000913/3957470‑3957626

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