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.C98.9% 582.7 / ‑5.1 184intergenic (‑127/‑734)ppiC/reppeptidyl‑prolyl cis‑trans isomerase C/rep helicase
Reads supporting (aligned to +/- strand):  ref base . (2/0);  major base C (84/96);  minor base A (1/1);  total (87/97)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 196 total aligned reads displayed.

TTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCC‑‑TGCTCGCAGGTCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATG  >  NC_000913/3957469‑3957645
                                                                                                 |                                                                                        
   gCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTccct                                                                                      >  2:4653394/1‑91 (MQ=255)
     aTGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTccctca                                                                                    >  2:3178309/1‑89 (MQ=255)
     aTGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGGACATATATGGCAAAAGCCCGGGAAATCCGGGCTTCGCAGGTccctcc                                                                                    >  2:2400370/1‑89 (MQ=255)
        aTCCGCAGATGCTAACGTAATTGTCAGTCGGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAgg                                                                                 <  1:1846983/93‑1 (MQ=255)
        aTCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAgg                                                                                 <  2:4503788/93‑1 (MQ=255)
        aTCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAgg                                                                                 <  1:4443476/93‑1 (MQ=255)
         tCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAggg                                                                                <  1:3730967/93‑1 (MQ=255)
         tCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAggg                                                                                >  2:4453096/1‑93 (MQ=255)
          ccGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAgggg                                                                               <  2:1534196/93‑1 (MQ=255)
            gCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAgggggg                                                                             >  1:1666371/1‑91 (MQ=255)
            gCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGGACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGgag                                                                             >  1:2855760/1‑93 (MQ=255)
             cAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTTCCTATATCGCAAAAGCCCGGGAAAAACG‑‑GGCGTCGCAGGTCCCTCAAAGGgaga                                                                            >  1:4587355/1‑93 (MQ=255)
              aGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGGAAAAGCCCGGAAAATCCC‑‑GGCTTCGCAAGTACCTCAAGGGgagca                                                                           >  1:4140935/1‑91 (MQ=255)
              aGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGgagga                                                                           >  1:3216590/1‑91 (MQ=255)
               gATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGaaa                                                                          >  2:367432/1‑93 (MQ=255)
               gATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGaaa                                                                          >  2:2703296/1‑93 (MQ=255)
                 tGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGGACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAAc                                                                        >  2:3318892/1‑93 (MQ=255)
                   cTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGACCCTCAAGGGGAGAAAACtt                                                                      >  1:1284556/1‑93 (MQ=255)
                   cTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACtt                                                                      >  1:4197072/1‑93 (MQ=255)
                    tAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCC‑‑GGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTa                                                                     <  2:2275360/93‑1 (MQ=255)
                     aaCGTAATTGTCAGTCAGGGAGTGGCTGTTACATATATGGCAAAATCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAg                                                                    <  2:4067072/93‑1 (MQ=255)
                     aaCGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAg                                                                    <  2:1949860/93‑1 (MQ=255)
                       cGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAggg                                                                  <  1:784265/93‑1 (MQ=255)
                        gTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGc                                                                 <  2:3910044/93‑1 (MQ=255)
                         tAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGGTTCGCAGGTCCCCCAAGGGGAGAAAACTTAggggc                                                                >  2:2508064/1‑91 (MQ=255)
                         tAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGcc                                                                >  1:3918564/1‑93 (MQ=255)
                         tAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGcc                                                                <  1:2078765/93‑1 (MQ=255)
                          aaTTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCct                                                               >  1:3329894/1‑93 (MQ=255)
                           atggtcagtcaGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCC‑‑GGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCctc                                                              <  2:677083/90‑1 (MQ=255)
                           aTTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCctc                                                              <  1:1020261/93‑1 (MQ=255)
                             tGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTa                                                            <  2:1273965/93‑1 (MQ=255)
                             tGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTa                                                            <  2:4150698/93‑1 (MQ=255)
                              gtcagtcaGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTAt                                                           <  2:130383/93‑1 (MQ=255)
                              gtcagtcaGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTAt                                                           <  2:3303971/93‑1 (MQ=255)
                               tcagtcagGGAGTGGCTGGTACATATATGGCAAAAGGCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCCATg                                                          >  2:1630245/1‑93 (MQ=255)
                               tcagtcagGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATg                                                          <  1:538367/93‑1 (MQ=255)
                                  gtcagGGAGTGGCTGGTACGTATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGc                                                       <  2:480700/93‑1 (MQ=255)
                                  gtcagGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGc                                                       <  2:2999988/93‑1 (MQ=255)
                                    cagGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTa                                                     <  1:2390802/93‑1 (MQ=255)
                                     agGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTAc                                                    <  1:3110353/93‑1 (MQ=255)
                                     agGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTAc                                                    >  1:3007812/1‑93 (MQ=255)
                                      gggAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACa                                                   >  2:1766907/1‑93 (MQ=255)
                                      gggAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACa                                                   >  2:241424/1‑93 (MQ=255)
                                      gggAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCC‑‑GGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACa                                                   >  1:1499344/1‑93 (MQ=255)
                                        gAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACaaa                                                 <  1:1416845/93‑1 (MQ=255)
                                          gTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGt                                               <  2:4116896/93‑1 (MQ=255)
                                          gTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGt                                               >  2:956580/1‑93 (MQ=255)
                                           tGGCTGGTACATATATGGCAAAAGCCCGGAAAATCC‑‑GGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTc                                              >  1:98210/1‑93 (MQ=255)
                                            ggCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCa                                             <  2:3399515/93‑1 (MQ=255)
                                            ggCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCa                                             >  2:4242645/1‑93 (MQ=255)
                                            ggCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCa                                             <  1:1675130/93‑1 (MQ=255)
                                            ggCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCa                                             <  1:2609307/93‑1 (MQ=255)
                                            ggCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCa                                             <  1:395835/93‑1 (MQ=255)
                                               tGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTa                                          <  2:1197357/93‑1 (MQ=255)
                                                ggTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTAt                                         <  2:4287924/93‑1 (MQ=255)
                                                 gTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATc                                        >  1:3559288/1‑93 (MQ=255)
                                                 gTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATc                                        >  1:3286258/1‑93 (MQ=255)
                                                   aCATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCgg                                      <  1:1031198/93‑1 (MQ=255)
                                                    cATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCggg                                     >  1:4676919/1‑93 (MQ=255)
                                                    cATATATGGCAAAAGCCCGGAAAATCC‑‑GGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCggg                                     >  2:3766126/1‑93 (MQ=255)
                                                     atatatGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCgggg                                    >  1:4603493/1‑93 (MQ=255)
                                                      tatatGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGa                                   >  1:1325728/1‑93 (MQ=255)
                                                       atatGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGAt                                  <  2:265810/93‑1 (MQ=255)
                                                        tatGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATg                                 >  2:510672/1‑93 (MQ=255)
                                                          tGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGt                               <  2:1462522/93‑1 (MQ=255)
                                                          tGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGt                               >  1:3142560/1‑93 (MQ=255)
                                                             cAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGttat                            >  1:4016532/1‑93 (MQ=255)
                                                              aaaaGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGCTGGttatt                           >  1:2239000/1‑93 (MQ=255)
                                                              aaaaGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGttatt                           >  1:2945734/1‑93 (MQ=255)
                                                              aaaaGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGttatt                           <  1:1949039/93‑1 (MQ=255)
                                                              aaaaGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGttatt                           >  1:3342342/1‑93 (MQ=255)
                                                              aaaaGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGttatt                           >  2:3647666/1‑93 (MQ=255)
                                                              aaaaGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGttatt                           <  2:3476585/93‑1 (MQ=255)
                                                                 aGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGttatttta                        <  2:1014414/93‑1 (MQ=255)
                                                                   cccGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGttattttatt                      <  2:1234236/93‑1 (MQ=255)
                                                                   cccGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGttattttatt                      <  1:3800201/93‑1 (MQ=255)
                                                                   cccGGAAAATCCGGGCGTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGttattttatt                      <  1:3258482/93‑1 (MQ=255)
                                                                    ccGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTg                     >  1:3605119/1‑93 (MQ=255)
                                                                    ccGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTg                     >  1:3630673/1‑93 (MQ=255)
                                                                    ccGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTg                     <  2:3996342/93‑1 (MQ=255)
                                                                                    ggcttCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAg         >  1:3440693/5‑93 (MQ=255)
                                                                                     gcttCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGc        >  1:1127069/4‑93 (MQ=255)
                                                                                     gcttCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGc        <  2:3855718/90‑1 (MQ=255)
                                                                                     gctgcGCAGGTCCCTCAAGGGGAGAAACCTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGc        <  1:972014/89‑1 (MQ=255)
                                                                                      tttCGCAGGTCCGTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTTCAAAGTCATTCTCGGGGATGCTTATTTTATTGTCAATTTATGAGCg       <  2:672912/91‑1 (MQ=255)
                                                                                      cttCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCg       >  1:3452303/3‑93 (MQ=255)
                                                                                       ttCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGa      >  1:377684/2‑93 (MQ=255)
                                                                                       ttCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGa      >  2:106642/2‑93 (MQ=255)
                                                                                       ttCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGa      <  2:3537196/92‑1 (MQ=255)
                                                                                       ttCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGa      >  1:3856433/2‑93 (MQ=255)
                                                                                        tCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGat     <  2:2853594/93‑1 (MQ=255)
                                                                                        tCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGat     <  1:3093857/93‑1 (MQ=255)
                                                                                         cGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGata    <  2:3088012/93‑1 (MQ=255)
                                                                                          gCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGatat   >  1:4312881/1‑93 (MQ=255)
                                                                                          gCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGatat   <  1:718913/93‑1 (MQ=255)
                                                                                          gCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGatat   >  1:3535182/1‑93 (MQ=255)
                                                                                          gCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGatat   <  2:82057/93‑1 (MQ=255)
                                                                                           cAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATg  <  1:2896696/93‑1 (MQ=255)
                                                                                           cAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATg  <  2:1820970/93‑1 (MQ=255)
                                                                                           cAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATg  >  2:3969323/1‑93 (MQ=255)
                                                                                                 |                                                                                        
TTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCC‑‑TGCTCGCAGGTCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATG  >  NC_000913/3957469‑3957645

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: