Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,252,715 A→G S821P (TCC→CCC)  b1202 ← putative adhesion and penetration protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,252,7150AG99.5% 682.2 / NA 201S821P (TCC→CCC) b1202putative adhesion and penetration protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (78/122);  minor base T (0/1);  total (78/123)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 204 total aligned reads displayed.

GCCATTACTCGCCTTCATAGTCATCTCATTCTGGTGGCTGTATGTAAGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGACAGGTTGAACCTCTGCCCGGCTTCCAGGGAGATGCTCATTCCATTCGCAGTGCCGTTGGCGTTAACGCCGTTGTTCTGACTGTCCAGTA  >  NC_000913/1252626‑1252804
                                                                                         |                                                                                         
gCCATTACTCGCCTTCATAGTCATCTCATTCTGGTGGCTGTATGTAAGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAg                                                                                        >  2:291948/1‑93 (MQ=255)
gCCATTACTCGCCTTCATAGTCATCTCATTCTGGTGGCTGTATGTAAGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAg                                                                                        >  1:3908615/1‑93 (MQ=255)
gCCATTACTCGCCTTCATAGTCATCTCATTCTGGTGGCTGTATGTAAGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAg                                                                                        <  1:5461509/93‑1 (MQ=255)
   aTTACTCGCCTTCATAGTCATCTCATTCTGGTGGCTGTATGTAAGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAGGtt                                                                                     <  1:1967737/93‑1 (MQ=255)
    ttACTCGCCTTCATAGTCATCTCATTCTGGTGGCTGTATGTAAGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAGGTTg                                                                                    >  1:5944001/1‑93 (MQ=255)
      aCTCGCCTTCATAGTCATCTCATTCTGGTGGCTGTATGTAAGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAGGTTGaa                                                                                  >  2:6148789/1‑93 (MQ=255)
       cTCGCCTTCATAGTCATCTCATTCTGGTGGCTGTATGTAAGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAGGTTGAAc                                                                                 <  2:4306222/93‑1 (MQ=255)
       cTCGCCTTCATAGTCATCTCATTCTGGTGGCTGTATGTAAGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAGGTTGAAc                                                                                 <  1:4833088/93‑1 (MQ=255)
         cGCCTTCATAGTCATCTCATTCTGGTGGCTGTATGTAAGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAGGTTGAACct                                                                               >  1:125907/1‑93 (MQ=255)
          gCCTTCATAGTCATCTCATTCTGGTGGCTGTATGTAAGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAGGTTGAACctc                                                                              >  2:5967339/1‑93 (MQ=255)
          gCCTTCATAGTCATCTCATTCTGGTGGCTGTATGTAAGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAGGTTGAACctc                                                                              <  2:5630269/93‑1 (MQ=255)
           ccTTCATAGTCATCTCATTCTGGTGTCTGTATGTAAGCTGGGTGTGCGGCTCTATATAGAACCCATAACCAGTAGGGGTCAGGTTGAACctct                                                                             <  1:6082655/93‑1 (MQ=255)
           ccTTCATAGTCATCTCATTCTGGTGGCTGTATGTAAGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAGGTTGAACctct                                                                             <  1:4681111/93‑1 (MQ=255)
           ccTTCATAGTCATCTCATTCTGGTGGCTGTATGTAAGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAGGTTGAACctct                                                                             <  1:2733868/93‑1 (MQ=255)
            cTTCATAGTCATCTCATTCTGGTGGCTGTATGTAAGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAGGTTGAACCTCTg                                                                            >  2:4182638/1‑93 (MQ=255)
            cTTCATAGTCATCTCATTCTGGTGGCTGTATGTAAGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAGGTTGAACCTCTg                                                                            >  2:1807205/1‑93 (MQ=255)
            cTTCATAGTCATCTCATTCTGGTGGCTGTATGTAAGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAGGTTGAACCTCTg                                                                            >  1:5150921/1‑93 (MQ=255)
             ttCATAGTCATCTCATTCTGGTGGCTGTATGTAAGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAGGTTGAACCTCTGc                                                                           >  1:1090948/1‑93 (MQ=255)
             ttCATAGTCATCTCATTCTGGTGGCTGTATGTAAGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAGGTTGAACCTCTGc                                                                           <  2:478164/93‑1 (MQ=255)
              tCATAGTCATCTCATTCTGGTGGCTGTATGTAAGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAGGTTGAACCTCTGcc                                                                          <  2:3113057/93‑1 (MQ=255)
              tCATAGTCATCTCATTCTGGTGGCTGTATGTAAGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAGGTTGAACCTCTGcc                                                                          <  1:3223893/93‑1 (MQ=255)
              tCATAGTCATCTCATTCTGGTGGCTGTATGTAAGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAGGTTGAACCTCTGcc                                                                          <  2:3382530/93‑1 (MQ=255)
               taTAGTCATCTCATTCTGGTGGCTGTATGTAAGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAGGTTGAACCTCTGccc                                                                         <  2:1393474/92‑1 (MQ=255)
                aTAGTCATCTCATTCTGGTGGCTGTATGTAAGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAGGTTGAACCTCTGCCCg                                                                        >  1:1019250/1‑93 (MQ=255)
                aTAGTCATCTCATTCTGGTGGCTGTATGTAAGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAGGTTGAACCTCTGCCCg                                                                        >  2:6282827/1‑93 (MQ=255)
                  aGTCATCTCATTCTGGTGGCTGTATGTAAGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAGGTTGAACCTCTGCCCGGc                                                                      <  2:2355087/93‑1 (MQ=255)
                   gTCATCTCATTCTGGTGGCTGTATGTAAGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAGGTTGAACCTCTGCCCGGCt                                                                     <  2:3375625/93‑1 (MQ=255)
                   gTCATCTCATTCTGGTGGCTGTATGTAAGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAGGTTGAACCTCTGCCCGGCt                                                                     <  1:3446064/93‑1 (MQ=255)
                   gTCATCTCATTCTGGTGGCTGTATGTAAGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAGGTTGAACCTCTGCCCGGCt                                                                     <  2:784812/93‑1 (MQ=255)
                     cacctcATTCTGGTGGTTGTATGTAACCTGGGTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAGGTTGAACCTCTGCCCGGCTTcc                                                                  <  2:2428385/90‑1 (MQ=255)
                     cATCTCATTCTGGTGGCTGTATGTAAGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAGGTTGAACCTCTGCCCGGCTTc                                                                   <  2:3617954/93‑1 (MQ=255)
                     cATCTCATTCTGGTGGCTGTATGTAAGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAGGTTGAACCTCTGCCCGGCTTc                                                                   >  2:2745266/1‑93 (MQ=255)
                     cATCTCATTCTGGTGGCTGTATGTAAGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAGGTTGAACCTCTGCCCGGCTTc                                                                   >  2:646231/1‑93 (MQ=255)
                       tctcATTCTGGTGGCTGTATGTAAGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAGGTTGAACCTCTGCCCGGCTTCCa                                                                 <  1:5196641/93‑1 (MQ=255)
                        ctcATTCTGGTGGCTGTATGTAAGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAGGTTGAACCTCTGCCCTGCTTCCAg                                                                <  1:6364476/93‑1 (MQ=255)
                        ctcATTCTGGTGGCTGTATGTAAGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAGGTTGAACCTCTGCCCGGCTTCCAg                                                                <  1:3835200/93‑1 (MQ=255)
                        ctcATTCTGGTGGCTGTATGTAAGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAGGTTGAACCTCTGCCCGGCTTCCAg                                                                <  1:828273/93‑1 (MQ=255)
                         tcATTCTGGTGGCTGTATGTAAGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAGGTTGAACCTCTGCCCGGCTTCCAgg                                                               <  1:5907087/93‑1 (MQ=255)
                         tcATACTGGTGGCTGTATGTAAGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAGGTTGAACCTCTGCCCGGCTTCCAgg                                                               <  1:5795782/93‑1 (MQ=255)
                          cATTCTGGTGGCTGTATGTAAGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAGGTTGAACCTCTGCCCGGCTTCCCggg                                                              <  1:3194031/93‑1 (MQ=255)
                          cATTCTGGTGGCTGTATGTAAGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAGGTTGAACCTCTGCCCGGCTTCCAggg                                                              >  2:4317738/1‑93 (MQ=255)
                           aTTCTGGTGGCTGTATGTAAGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAGGTTGAACCTCTGCCCGGCTTCCAGGga                                                             >  1:2186658/1‑93 (MQ=255)
                                ggtggCTGTATGTAAGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAGGTTGAACCTCTGCCCGGCTTCCAGGGAGATGc                                                        >  1:2814315/1‑93 (MQ=255)
                                 gtggCTGTATGTAAGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAGGTTGAACCTCTGCCCGGCTTCCAGGGAGATGCt                                                       <  2:87552/93‑1 (MQ=255)
                                   ggCTGTATGTAAGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAGGTTGAACCTCTGCCCGGCTTCCAGGGAGATGCTca                                                     >  2:1949650/1‑93 (MQ=255)
                                        tatgtaAGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAGGTTGAACCTCTGCCCGGCTTCCAGGGAGATGCTcattcca                                                >  1:3544308/1‑93 (MQ=255)
                                          tgtaAGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAGGTTGAACCTCTGCCCGGCTTCCAGGGAGATGCTcattccatt                                              <  1:5317716/93‑1 (MQ=255)
                                           gtaAGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAGGTTGAACCTCTGCCCGGCTTCCAGGGAGATGCTcattccattc                                             <  1:1837511/93‑1 (MQ=255)
                                            taAGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAGGTTGAACCTCTGCCCGGCTTCCAGGGAGATGCTCATTCCATTCg                                            <  1:1809310/93‑1 (MQ=255)
                                               gCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAGGTTGAACCTCTGCCCGGCTTCCAGGGAGATGCTCATTCCATTCGCAg                                         >  1:3493133/1‑93 (MQ=255)
                                                  gggTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAGGTTGAACCTCTGCCCGGCTTCCAGGGAGATGCTCATTCCATTCGCAGTGc                                      >  1:1533979/1‑93 (MQ=255)
                                                  gggTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAGGTTGAACCTCTGCCCGGCTTCCAGGGAGATGCTCATTCCATTCGCAGTGc                                      <  1:296603/93‑1 (MQ=255)
                                                    gTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAGGTTGAACCTCTGCCCGGCTTCCAGGGAGATGCTCATTCCATTCGCAGTGCCg                                    >  1:3048529/1‑93 (MQ=255)
                                                     tttGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAGGTTGAACCTCTGCCCGGCTTCCAGGGAGATGCTCATTCCATTCGCAGTGCCGt                                   >  1:5715565/1‑93 (MQ=255)
                                                     tttGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAGGTTGAACCTCTGCCCGGCTTCCAGGGAGATGCTCATTCCATTCGCAGTGCCGt                                   <  1:4068538/93‑1 (MQ=255)
                                                      ttGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAGGTTGAACCTCTGCCCGGCTTCCAGGGAGATGCTCATTCCATTCGCAGTGCCGtt                                  <  1:5247263/93‑1 (MQ=255)
                                                       tGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAGGTTGAACCTCTGCCCGGCTTCCAGGGAGATGCTCATTCCATTCGCAGTGCCGTTg                                 <  1:5230345/93‑1 (MQ=255)
                                                       tGCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAGGTTGAACCTCTGCCCGGCTTCCAGGGAGATGCTCATTCCATTCGCAGTGCCGTTg                                 <  2:2211725/93‑1 (MQ=255)
                                                        gCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAGGTTGAACCTCTGCCCGGCTTCCAGGGAGATGCTCATTCCATTCGCAGTGCCGTTgg                                >  2:1129899/1‑93 (MQ=255)
                                                        gCGGCTCTATATAGAACCCATAACCAGTAGGGGGCAGGTTGAACCTCTGCCCGGCTTCCAGGGAGATGCTCATTCCATTCGCAGTGCCGTTgg                                <  2:2965506/93‑1 (MQ=255)
                                                           gCTCTATATAGAACCCATAACCAGTAGGGGGCTGGTTGACCCTCTGCCCGGCTTCCAGGGAGATGCTCATTCCATTCGCAGTGCCGTTGGCGt                             <  1:4287428/93‑1 (MQ=255)
                                                           gCTCTATATAGAACCCATAACCAGTAGGGGGCAGGTTGAACCTCTGCCCGGCTTCCAGGGAGATGCTCATTCCATTCGCAGTGCCGTTGGCGt                             <  1:4911500/93‑1 (MQ=255)
                                                                atataGAACCCATAACCAGTAGGGGGCAGGTTGAACCTCTGCCCGGCTTCCAGGGAGATGCTCATTCCATTCGCAGTGCCGTTGGCGTTAACg                        <  1:1727950/93‑1 (MQ=255)
                                                                 tataGAACCCATAACCAGTAGTGGGCAGGTTGAACCTCTGCCCGGCTTCCAGGGAGATGCTCATTCCATTCGCAGTGCCGTTGGCGTTAACGc                       <  1:226694/93‑1 (MQ=255)
                                                                 tataGAACCCATAACCAGTAGGGGGCAGGTTGAACCTCTGCCCGGCTTCCAGGGAGATGCTCATTCCATTCGCAGTGCCGTTGGCGTTAACGc                       <  2:4940966/93‑1 (MQ=255)
                                                                 tataGAACCCATAACCAGTAGGGGGCAGGTTGAACCTCTGCCCGGCTTCCAGGGAGATGCTCATTCCATTCGCAGTGCCGTTGGCGTTAACGc                       <  2:5508278/93‑1 (MQ=255)
                                                                  ataGAACCCATAACCAGTAGGGGGCAGGTTGAACCTCTGCCCGGCTTCCAGGGAGATGCTCATTCCATTCGCAGTGCCGTTGGCGTTAACGcc                      <  1:4975535/93‑1 (MQ=255)
                                                                  ataGAACCCATAACCAGTAGGGGGCAGGTTGAACCTCTGCCCGGCTTCCAGGGAGATGCTCATTCCATTCGCAGTGCCGTTGGCGTTAACGcc                      >  2:3418931/1‑93 (MQ=255)
                                                                    aGAACCCATAACCAGTAGGGGGCAGGTTGAACCTCTGCCCGGCTTCCAGGGAGATGCTCATTCCATTCGCAGTGCCGTTGGCGTTAACGCCgt                    <  2:3203761/93‑1 (MQ=255)
                                                                    aGAACCCATAACCAGTAGGGGGCAGGTTGAACCTCTGCCCGGCTTCCAGGGAGATGCTCATTCCATTCGCAGTGCCGTTGGCGTTAACGCCgt                    >  2:3308462/1‑93 (MQ=255)
                                                                    aGAACCCATAACCAGTAGGGGGCAGGTTGAACCTCTGCCCGGCTTCCAGGGAGATGCTCATTCCATTCGCAGTGCCGTTGGCGTTAACGCCgt                    <  2:3073866/93‑1 (MQ=255)
                                                                    aGAACCCATAACCAGTAGGGGGCAGGTTGAACCTCTGCCCGGCTTCCAGGGAGATGCTCATTCCATTCGCAGTGCCGTTGGCGTTAACGCCgt                    <  2:5319755/93‑1 (MQ=255)
                                                                     aaaCCCATAACCAGTAGGGGGCAGGTTGAACCTCTGCCCGGCTTCCAGGGAGATGCTCATTCCATTCGCAGTGCCGTTGGCGTTAACGCCgtt                   <  1:4144846/92‑1 (MQ=255)
                                                                       aCCCATAACCAGTAGGGGGCAGGTTGAACCTCTGCCCGGCTTCCAGGGAGATGCTCATTCCATTCGCAGTGCCGTTGGCGTTAACGCCgttgt                 >  2:1529317/1‑93 (MQ=255)
                                                                         ccATAGTCATTAGGGGGCAGGTTGAACCTCTGCCCGGCTTCCAGGGAGATGCTCATTCCATTCGCAGTGCCGTTGGCGTTAACGCCGTTGTTc               <  1:1396998/93‑1 (MQ=255)
                                                                         ccATAACCAGTAGGGGGCAGGTTGAACCTCTGCCCGGCTTCCAGGGAGATGCTCATTCCATTCGCAGTGCCGTTGGCGTTAACGCCgttgttt               >  2:1456051/1‑92 (MQ=255)
                                                                         ccATAACCAGTAGGGGGCAGGTTGAACCTCTGCCCGGCTTCCAGGGAGATGCTCATTCCATTCGCAGTGCCGTTGGCGTTAACGCCGTTGTTc               <  2:3761594/93‑1 (MQ=255)
                                                                         ccATAACCAGTAGGGGGCAGGTTGAACCTCTGCCCGGCTTCCAGGGAGATGCTCATTCCATTCGCAGTGCCGTTGCCGTTAACGCCGTTGTTc               >  1:458156/1‑93 (MQ=255)
                                                                          cATAACCAGTAGGGGGCAGGTTGAACCTCTGCCCGGCTTCCAGGGAGATGCTCATTCCATTCGCAGTGCCGTTGGCGTTAACGCCGTTGTTCt              >  2:5418050/1‑93 (MQ=255)
                                                                          cATAACCAGTAGGGGGCAGGTTGAACCTCTGCCCGGCTTCCAGGGAGATGCTCATTCCATTCGCAGTGCCGTTGGCGTTAACGCCGTTGTTCt              <  2:5068906/93‑1 (MQ=255)
                                                                          cATAACCAGTAGGGGGCAGGTTGAACCTCTGCCCGGCTTCCAGGGAGATGCTCATTCCATTCGCAGTGCCGTTGGCGTTAACGCCGTTGTTCt              <  1:3605137/93‑1 (MQ=255)
                                                                           aTAACCAGTAGGGGGCAGGTTGAACCTCTGCCCGGCTTCCAGGGAGATGCTCATTCCATTCGCAGTGCCGTTGGCGTTAACGCCGTTGTTCTg             <  1:4661964/93‑1 (MQ=255)
                                                                           aTAACCAGTAGGGGGCAGGTTGAACCTCTGCCCGGCTTCCAGGGAGATGCTCATTCCATTCGCAGTGCCGTTGGCGTTAACGCCGTTGTTCTg             <  1:4053243/93‑1 (MQ=255)
                                                                            tAACCAGTAGGGGGCAGGTTGAACCTCTGCCCGGCTTCCAGGGAGATGCTCATTCCATTCGCAGTGCCGTTGGCGTTAACGCCGTTGTTCTGa            <  2:4615801/93‑1 (MQ=255)
                                                                             aaCCAGTAGGGGGCAGGTTGAACCTCTGCCCGGCTTCCAGGGAGATGCTCATTCCATTCGCAGTGCCGTTGGCGTTAACGCCGTTGTTCTGAc           >  1:61988/1‑93 (MQ=255)
                                                                             aaCCAGTAGGGGGCAGGTTGAACCTCTGCCCGGCTTCCAGGGAGATGCTCATTCCATTCGCAGTGCCGTTGGCGTTAACGCCGTTGTTCTGAc           <  2:6401419/93‑1 (MQ=255)
                                                                              aCCAGTAGGGGGCAGGTTGAACCTCTGCCCGGCTTCCAGGGAGATGCTCATTCCATTCGCAGTGCCGTTGGCGTTAACGCCGTTGTTCTGACt          >  2:4749623/1‑93 (MQ=255)
                                                                               ccAGTAGGGGGCAGGTTGAACCTCTGCCCGGCTTCCAGGGAGATGCTCATTCCATTCGCAGTGCCGTTGGCGTTAACGCCGTTGTTCTGACTg         <  1:1248842/93‑1 (MQ=255)
                                                                               ccAGTAGGGGGCAGGTTGAACCTCTGCCCGGCTTCCAGGGAGATGCTCATTCCATTCGCAGTGCCGTTGGCGTTAACGCCGTTGTTCTGACTg         <  2:1788975/93‑1 (MQ=255)
                                                                               ccAGTAGGGGGCAGGTTGAACCTCTGCCCGGCTTCCAGGGAGATGCTCATTCCATTCGCAGTGCCGTTGGCGTTAACGCCGTTGTTCTGACTg         <  2:5578341/93‑1 (MQ=255)
                                                                               ccAGTAGGGGGCAGGTTGAACCTCTGCCCGGCTTCCAGGGAGATGCTCATTCCATTCGCAGTGCCGTTGGCGTTAACGCCGTTGTTCTGACTg         <  1:5826517/93‑1 (MQ=255)
                                                                                cAGTAGGGGGCAGGTTGAACCTCTGCCCGGCTTGCAGGGAGATGCTCATTCCATTCGCAGTGCCGTTGGCGTTAACGCCGTTGTTCTGACTGt        >  2:4218427/1‑93 (MQ=255)
                                                                                cAGTAGGGGGCAGGTTGAACCTCTGCCCGGCTTCCAGGGAGATGCTCATTCCATTCGCAGTGCCGTTGGCGTTAACGCCGTTGTTCTGACTGt        <  2:3017865/93‑1 (MQ=255)
                                                                                 aGTAGGGGGCAGGTTGAACCTCTGCCCGGCTTCCAGGGAGATGCTCATTCCATTCGCAGTGCCGTTGGCGTTAACGCCGTTGTTCTGACTGTc       <  2:224072/93‑1 (MQ=255)
                                                                                  gTAGGGGGCAGGTTGAACCTCTGCCCGGCTTCCAGGGAGATGCTCATTCCATTCGCAGTGCCGTTGGCGTTAACGCCGTTGTTCTGACTGTcc      <  1:3637678/93‑1 (MQ=255)
                                                                                   tAGGGGGCAGGTTGAACCTCTGCCCGGCTTCCAGGGAGATGCTCATTCCATTCGCAGTGCCGTTGGCGTTAACGCCGTTGTTCTGTCTGTCCa     >  2:3781878/1‑93 (MQ=255)
                                                                                   tAGGGGGCAGGTTGAACCTCTGCCCGGCTTCCAGGGAGATGCTCATTCCATTCGCAGTGCCGTTGGCGTTAACGCCGTTGTTCTGACTGTCCa     <  1:6149687/93‑1 (MQ=255)
                                                                                    aGGGGGCAGGTTGAACCTCTGCCCGGCTTCCAGGGAGATGCTCATTCCATTCGCAGTGCCGTTGGCGTTAACGCCGTTGTTCTGACTGTCCAg    >  2:4328375/1‑93 (MQ=255)
                                                                                     ggggGCAGGTTGAACCTCTGCCCGGCTTCCAGGGAGATGCTCATTCCATTCGCAGTGCCGTTGGCGTTAACGCCGTTGTTCTGACTGTCCAGt   <  1:555176/93‑1 (MQ=255)
                                                                                     ggggGCAGGTTGAACCTCTGCCCGGCTTCCAGGGAGATGCTCATTCCATTCGCAGTGCCGTTGGCGTTAACGCCGTTGTTCTGACTGTCCAGt   >  2:404227/1‑93 (MQ=255)
                                                                                         |                                                                                         
GCCATTACTCGCCTTCATAGTCATCTCATTCTGGTGGCTGTATGTAAGCTGGGTTTGCGGCTCTATATAGAACCCATAACCAGTAGGGGACAGGTTGAACCTCTGCCCGGCTTCCAGGGAGATGCTCATTCCATTCGCAGTGCCGTTGGCGTTAACGCCGTTGTTCTGACTGTCCAGTA  >  NC_000913/1252626‑1252804

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