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.C97.6% 1372.8 / 15.5 411intergenic (‑127/‑734)ppiC/reppeptidyl‑prolyl cis‑trans isomerase C/rep helicase
Reads supporting (aligned to +/- strand):  ref base . (9/1);  new base C (185/216);  total (194/217)
Fisher's exact test for biased strand distribution p-value = 7.69e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 426 total aligned reads displayed.

TAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCTGCTCGCAGGTCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATG  >  NC_000913/3957470‑3957645
                                                                                              |                                                                                       
tAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTcc                                                                                       >  1:3519080/1‑93 (MQ=255)
 aGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTccc                                                                                      <  2:1064051/93‑2 (MQ=255)
  gCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTccct                                                                                     <  2:4730922/93‑3 (MQ=255)
  gCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGGACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTccct                                                                                     >  2:181928/1‑91 (MQ=255)
       aTCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAgg                                                                                 <  1:5444915/93‑1 (MQ=255)
        tCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCGAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAggg                                                                                >  2:1697448/1‑93 (MQ=255)
        tCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAggg                                                                                <  2:1164613/93‑1 (MQ=255)
         ccGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAgggg                                                                               >  2:353311/1‑93 (MQ=255)
          cGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGga                                                                              <  2:4881286/93‑1 (MQ=255)
           gCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGgag                                                                             >  2:4566382/1‑93 (MQ=255)
           gCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGgag                                                                             >  2:6941401/1‑93 (MQ=255)
            cAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGgaga                                                                            >  2:7186010/1‑93 (MQ=255)
             aGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGaa                                                                           <  2:6658278/93‑1 (MQ=255)
             aGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGaa                                                                           >  2:1413406/1‑93 (MQ=255)
                tGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAAc                                                                        >  2:6072275/1‑93 (MQ=255)
                tGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGGACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAAc                                                                        >  1:6240871/1‑93 (MQ=255)
                 gCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACt                                                                       <  1:1858189/93‑1 (MQ=255)
                   tAACGTAATTGTCAGTCAGGGAGTGGCTGGGACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTa                                                                     >  1:769790/1‑93 (MQ=255)
                    aaCGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAg                                                                    <  2:3930457/93‑1 (MQ=255)
                    aaCGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAg                                                                    <  1:7242985/93‑1 (MQ=255)
                    aaCGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAg                                                                    >  1:4539520/1‑93 (MQ=255)
                    aaCGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAg                                                                    <  1:2142060/93‑1 (MQ=255)
                     aCGTAATTGTCAGTCAGGGAGTGGCTGGGACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGGCCCTCAAGGGGAGAAAACTTAgg                                                                   >  1:16063/1‑93 (MQ=255)
                     aCGTAATTGTCAGTCAGGGAGTGGCTGGGACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAgg                                                                   >  1:1213898/1‑93 (MQ=255)
                        tAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGcc                                                                >  2:4290977/1‑93 (MQ=255)
                        tAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGcc                                                                <  2:5047279/93‑1 (MQ=255)
                        tAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGcc                                                                >  2:5215680/1‑93 (MQ=255)
                        tAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGcc                                                                <  2:1145423/93‑1 (MQ=255)
                         aaTTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCct                                                               <  1:4167935/93‑1 (MQ=255)
                          aTTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCctc                                                              >  2:1955670/1‑93 (MQ=255)
                          aTTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCctc                                                              <  2:4005295/93‑1 (MQ=255)
                          aTTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCctc                                                              <  1:4588165/93‑1 (MQ=255)
                          aTTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCctc                                                              <  1:2804720/93‑1 (MQ=255)
                          aTTGTCAGTCAGGGAGTGGCTGGGACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCctc                                                              >  1:4908951/1‑93 (MQ=255)
                          aTTGTCAGTCAGGGAGTGGCTGGCACATATATGGCAAAAGCCCGGAAAAACCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAATGCctc                                                              >  2:6733661/1‑93 (MQ=255)
                           ttGTCAGTCAGGGAGTGGCTGGGACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCctct                                                             >  1:5725900/1‑93 (MQ=255)
                            tGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTa                                                            <  2:1246648/93‑1 (MQ=255)
                            tGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTa                                                            <  2:237985/93‑1 (MQ=255)
                            tGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTa                                                            <  2:2562609/93‑1 (MQ=255)
                            tGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTa                                                            <  2:368421/93‑1 (MQ=255)
                             gtcagtcaGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTAt                                                           <  1:5248299/93‑1 (MQ=255)
                             gtcagtcaGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTAt                                                           <  2:2664266/93‑1 (MQ=255)
                             gtcagtcaGGGAGTGGCTGGGACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTAt                                                           >  1:5038275/1‑93 (MQ=255)
                              tcagtcagGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATg                                                          <  1:1784803/93‑1 (MQ=255)
                              tcagtcagGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATg                                                          <  1:1297772/93‑1 (MQ=255)
                               cagtcagGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATgg                                                         <  1:2012237/93‑1 (MQ=255)
                                 gtcagGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGc                                                       <  1:4038328/93‑1 (MQ=255)
                                  tcagGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCt                                                      >  2:2827864/1‑93 (MQ=255)
                                  tcagGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCt                                                      >  1:3801659/1‑93 (MQ=255)
                                  tcagGGAGTGGCTGGGACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGGGAAAACTTAGGGGCTCTATGGGCt                                                      >  1:6527962/1‑93 (MQ=255)
                                    agGGAGTGGCTGGTCCATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTAc                                                    >  2:6441802/1‑93 (MQ=255)
                                    agGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTAc                                                    >  1:3701852/1‑93 (MQ=255)
                                    agGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCCAc                                                    >  1:4709057/1‑93 (MQ=255)
                                     gggAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACa                                                   >  2:3648854/1‑93 (MQ=255)
                                     gggAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACa                                                   >  2:3157907/1‑93 (MQ=255)
                                     gggAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACa                                                   >  1:5612325/1‑93 (MQ=255)
                                       gAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCCGGCTTCGCAGGTCCCTCAAGGGGGGAAAACTTAGGGCCTCTATGGGCTACaaa                                                 >  1:6648032/1‑93 (MQ=255)
                                       gAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACaaa                                                 <  1:3921276/93‑1 (MQ=255)
                                       gAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACaaa                                                 <  1:946153/93‑1 (MQ=255)
                                       gAGTGGCTGGGACATATATGGCAAAAGCCCGGAAAAACCGGGCTTCGCAGGGCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACaaa                                                 >  1:5829414/1‑93 (MQ=255)
                                        aGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAg                                                >  2:3458358/1‑93 (MQ=255)
                                             cTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCAtt                                           <  2:6353838/93‑1 (MQ=255)
                                              tGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTa                                          <  1:2435657/93‑1 (MQ=255)
                                              tGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTa                                          <  2:5844472/93‑1 (MQ=255)
                                               ggTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTAt                                         <  2:741027/93‑1 (MQ=255)
                                                  aCATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCgg                                      >  2:569953/1‑93 (MQ=255)
                                                  aCATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCgg                                      <  2:6223907/93‑1 (MQ=255)
                                                   cATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCggg                                     <  2:6993193/93‑1 (MQ=255)
                                                     tatatGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGa                                   <  1:2244573/93‑1 (MQ=255)
                                                     tatatGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGa                                   <  1:3214330/93‑1 (MQ=255)
                                                       tatGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATg                                 >  1:3147509/1‑93 (MQ=255)
                                                          ggCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGtt                              <  1:6866911/93‑1 (MQ=255)
                                                          ggCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGtt                              <  1:7080039/93‑1 (MQ=255)
                                                           gCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGtta                             <  2:5476560/93‑1 (MQ=255)
                                                            cAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGttat                            <  2:5388990/93‑1 (MQ=255)
                                                             aaaaGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGttatt                           >  1:5535530/1‑93 (MQ=255)
                                                             aaaaGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGttatt                           >  2:4167343/1‑93 (MQ=255)
                                                             aaaaGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGttatt                           >  2:2454424/1‑93 (MQ=255)
                                                               aaGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGttatttt                         <  2:804487/93‑1 (MQ=255)
                                                                aGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGttatttta                        <  2:4809044/93‑1 (MQ=255)
                                                                 gCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGttattttat                       >  1:6119491/1‑93 (MQ=255)
                                                                  cccGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGttattttatt                      <  1:5375366/93‑1 (MQ=255)
                                                                   ccGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTg                     <  2:2133588/93‑1 (MQ=255)
                                                                   ccGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTg                     <  2:6880842/93‑1 (MQ=255)
                                                                   ccGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTg                     >  2:6103386/1‑93 (MQ=255)
                                                                    cGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGt                    <  2:6050192/93‑1 (MQ=255)
                                                                    cGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGt                    >  1:6459262/1‑93 (MQ=255)
                                                                    cGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGt                    <  2:3300779/93‑1 (MQ=255)
                                                                                 ggcttCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAg         >  1:2677529/5‑93 (MQ=255)
                                                                                 ggcttCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAg         >  1:5150964/5‑93 (MQ=255)
                                                                                 ggcttCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAg         <  1:1154996/89‑1 (MQ=255)
                                                                                  gcttCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGc        <  2:3048222/90‑1 (MQ=255)
                                                                                  gcttCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGc        >  2:7306014/4‑93 (MQ=255)
                                                                                   cttCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCg       <  1:1482368/91‑1 (MQ=255)
                                                                                     tCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGat     <  1:835016/93‑1 (MQ=255)
                                                                                       gCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGatat   >  1:191796/1‑93 (MQ=255)
                                                                                       gCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGatat   <  2:4415005/93‑1 (MQ=255)
                                                                                       gCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGatat   <  1:2548976/93‑1 (MQ=255)
                                                                                       gCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGatat   >  1:2301319/1‑93 (MQ=255)
                                                                                        cAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATg  >  1:3059417/1‑93 (MQ=255)
                                                                                              |                                                                                       
TAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCTGCTCGCAGGTCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATG  >  NC_000913/3957470‑3957645

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