Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,785,283 Δ1 bp 100% coding (1830/2379 nt) ppsA ← phosphoenolpyruvate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,785,2810C.67.0% 128.9 / 152.6 97coding (1832/2379 nt)ppsAphosphoenolpyruvate synthase
Reads supporting (aligned to +/- strand):  ref base C (15/17);  new base . (29/36);  total (44/53)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.52e-01

ATCTCAACGTTGGTCAGTCCCATGTCGTTGCGCACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGGCCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGCATATTCGTTCGATTTAAAATCAGAGAGAC  >  NC_000913/1785185‑1785380
                                                                                                |                                                                                                   
aTCTCAACGTTGGTCAGTCCCATGTCGTTGCGCACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGGccgc                                                                                                  >  2:3672058/1‑98 (MQ=255)
aTCTCAACGTTGGTCAGTCCCATGTCGTTGCGCACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGGCCcg                                                                                                  >  1:3525948/1‑100 (MQ=255)
  ctcAACGTTGGTCAGTCCCATGTCGTTGCGCACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGGccgcgc                                                                                                >  2:3070649/1‑96 (MQ=255)
   tcAACGTTGGTCAGTCCCATGTCGTTGCGCACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGG‑CCGCGcc                                                                                              >  2:1232396/1‑100 (MQ=255)
      aCGTTGGTCAGTCCCATGTCGTTGCGCACACGTTTGACTGCTTCACACTCCAGCGCGAAACAGGCGCGGAAGCTGTCGGAAACATAGCGGCCCGCGCCAc                                                                                            <  2:2593046/100‑1 (MQ=255)
      aCGTTGGTCAGTCCCATGTCGTTGCGCACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGACGCGGAAGCTGTCGGAAACATAGCGG‑CCGCGCCACg                                                                                           >  1:3937309/1‑100 (MQ=255)
          tGGTCAGTCCCATGTCGTTGCGCACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGG‑CCGCGCCACGGAAg                                                                                       <  1:2712109/100‑1 (MQ=255)
           ggTCAGTCCCATGTCGTTGCGCACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGG‑CCGCGCCACGGAAGc                                                                                      >  1:4054354/1‑100 (MQ=255)
           ggTCAGTCCCATGTCGTTGCGCACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGG‑CCGCGCCACGGAAGc                                                                                      >  1:2845731/1‑100 (MQ=255)
            gTCAGTCCCATGTCGTTGCGCACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGG‑CCGCGCCACGGAAGcc                                                                                     >  1:3107187/1‑100 (MQ=255)
            gTCAGTCCCATGTCGTTGCGCACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGG‑CCGCGCCACGGAAGcc                                                                                     >  1:2197493/1‑100 (MQ=255)
            gTCAGTCCCATGTCGTTGCGCACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGG‑CCGCGCCACGGAAGcc                                                                                     >  1:1658561/1‑100 (MQ=255)
             tCAGTCCCATGTCGTTGCGCACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGGCCCGCGCCACGGAAGcc                                                                                     >  1:2343663/1‑100 (MQ=255)
                  cccATGTCGTTGCGCACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGACCCGCGCCACGGAAGCCGAGCa                                                                                >  1:2226834/1‑100 (MQ=255)
                   ccATGTCGTTGCGCACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGGCCCGCGCCACGGAAGCCGAGCAt                                                                               <  1:3794187/100‑1 (MQ=255)
                    cATGTCGTTGCGCACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGGCCCGCGCCACGGAAGCCGAGCATc                                                                              <  1:3164239/100‑1 (MQ=255)
                     aTGTCGTTGCGCACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGGCCCGCGCCACGGAAGCCGAGCATCg                                                                             <  2:2764963/100‑1 (MQ=255)
                     aTGTCGTTGCGCACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGG‑CCGCGCCACGGAAGCCGAGCATCgg                                                                            <  1:3120934/100‑1 (MQ=255)
                      tGTCGTTGCGCACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGG‑CCGCGCCACGGAAGCCGAGCATCggg                                                                           <  2:2724978/100‑1 (MQ=255)
                      tGTCGTTGCGCACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGG‑CCGCGCCACGGAAGCCGAGCATCggg                                                                           <  1:1981846/100‑1 (MQ=255)
                       gTCGTTGCGCACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGGCCCGCGCCACGGAAGCCGAGCATCggg                                                                           <  2:616126/100‑1 (MQ=255)
                       gTCGTTGCGCACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGG‑CCGCGCCACGGAAGCCGAGCATCGGGt                                                                          <  2:2149638/100‑1 (MQ=255)
                       gTCGTTGCGCACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGG‑ACGCGCCACGGAAGCCGAGCATCggg                                                                           >  2:4051769/1‑99 (MQ=255)
                        tCGTTGCGCACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGG‑CCGCGCCACGGAAGCCGAGCATCGGGtt                                                                         <  1:2181645/100‑1 (MQ=255)
                        tCGGTGCGCACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGGCCCGCGCCACGGAAGCCGAGCATCGGGt                                                                          <  2:51494/100‑1 (MQ=255)
                             gcgcACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGG‑CCGCGCCACGGAAGCCGAGCATCGGGTTCTCtt                                                                    >  2:2320103/1‑100 (MQ=255)
                                cacaCGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGGCCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCa                                                                  >  2:889052/1‑100 (MQ=255)
                                cacaCGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGGCCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCa                                                                  >  1:3100940/1‑100 (MQ=255)
                                cacaCGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGG‑CCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCAt                                                                 >  1:2508545/1‑100 (MQ=255)
                                  cacGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGGCCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTc                                                                   >  1:1548531/1‑97 (MQ=255)
                                  cacGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGGCCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATc                                                                >  1:3554862/1‑100 (MQ=255)
                                    cGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGG‑CCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTgg                                                             >  1:1698288/1‑100 (MQ=255)
                                     gTTTCCCTGCTTCACACCCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGG‑CCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGc                                                            <  2:3633290/100‑1 (MQ=255)
                                      tttCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGACACATAGCGG‑CCGCGCCACGGAAGCCGAGCATCGGGTTCTCTGCATCTCGc                                                            >  2:68125/1‑99 (MQ=255)
                                      tttCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGG‑CCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCt                                                           <  1:593759/100‑1 (MQ=255)
                                       ttCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGGCCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCt                                                           <  1:3178249/100‑1 (MQ=255)
                                       tccACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGG‑CCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTc                                                          <  2:2508531/98‑1 (MQ=255)
                                        tCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGG‑CCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCg                                                         >  1:2276356/1‑100 (MQ=255)
                                         cACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGGCCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCg                                                         >  2:3751028/1‑100 (MQ=255)
                                         cACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGG‑CCGCGCCACGTAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGt                                                        >  1:311895/1‑100 (MQ=255)
                                            tGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGG‑CCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAAc                                                     <  1:3606324/100‑1 (MQ=255)
                                            tGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGG‑CCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAAc                                                     <  2:1049188/100‑1 (MQ=255)
                                            tGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGG‑CCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAAc                                                     >  1:1037489/1‑100 (MQ=255)
                                             gCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGGCCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAAc                                                     >  1:2885259/1‑100 (MQ=255)
                                             gCTTCACACTCCAGAGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGG‑CCGCGCCACGGAAGACGAGCATCGGGTTCTCTTCATCTGGCTCGTaa                                                      >  2:3499689/1‑98 (MQ=255)
                                             acTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGG‑CCGCGGCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAAc                                                     >  2:2914153/2‑99 (MQ=255)
                                              cTTCACATTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCTG‑CCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGc                                                   <  2:3607429/100‑1 (MQ=255)
                                                 cacaCTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGGCCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTc                                                 <  2:163685/100‑1 (MQ=255)
                                                 cacaCTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGG‑CCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTca                                                >  2:208183/1‑100 (MQ=255)
                                                 cacaCTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGG‑CCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTca                                                <  1:1531916/100‑1 (MQ=255)
                                                 cacaCTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGG‑CCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTca                                                <  1:738843/100‑1 (MQ=255)
                                                   cacTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGG‑CCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTcacc                                              >  1:2853819/1‑100 (MQ=255)
                                                      tCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGG‑CCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCaccacc                                           <  1:2413137/100‑1 (MQ=255)
                                                      tCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGG‑CCGCGCCACGGAAGCCGAGCAACGGGTTCTCTTCATCTGGCTCGTAAAGATcaccac                                            >  2:3475886/1‑99 (MQ=255)
                                                       ccAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGGCCCGCGCCACGGAAGCCGAGCATCGGGTTGTCTTCATCTGGCTCGTAACGCTCaccacc                                           >  1:203436/1‑100 (MQ=255)
                                                       ccAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGGCCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCaccacc                                           >  1:175010/1‑100 (MQ=255)
                                                       ccAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGGCCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCaccacc                                           <  2:2739815/100‑1 (MQ=255)
                                                       ccAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGG‑CCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCg                                          >  1:490759/1‑100 (MQ=255)
                                                       ccAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGG‑CCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCg                                          >  2:3212408/1‑100 (MQ=255)
                                                         aGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGG‑CCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGAc                                        >  1:276167/1‑100 (MQ=255)
                                                         aGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGG‑CCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGAc                                        >  1:830908/1‑100 (MQ=255)
                                                          gcgcGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGG‑CCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGAcc                                       <  2:3172666/100‑1 (MQ=255)
                                                             cgAAACAGTCGCGGAAGCTGTCGGAAACATAGCGG‑CCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAgg                                    <  1:68534/100‑1 (MQ=255)
                                                              gAAACAGTCGCGGAAGCTGTCGGAAACATAGCGG‑CCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGCCCAGGt                                   >  1:3092794/1‑100 (MQ=255)
                                                                aaCAGTCGCGGAAGCTGTCGGAAACATAGCGGCCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGtt                                  <  1:172722/100‑1 (MQ=255)
                                                                aaCAGTCGCGGAAGCTGTCGGAAACATAGCGG‑CCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTg                                 <  2:3207857/100‑1 (MQ=255)
                                                                 aCAGTCGCGGAAGCTGTCGGAAACATAGCGG‑CCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTgg                                <  1:2333064/100‑1 (MQ=255)
                                                                  cAGTCGCGGAAGCTGTCGGAAACATAGCGG‑CCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGc                               <  2:3807022/100‑1 (MQ=255)
                                                                   aGTCGCGGAAGCTGTCGGAAACATAGCGG‑CCGGGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGCa                              <  2:1003476/100‑1 (MQ=255)
                                                                   aGTCGCGGAAGCTGTCGGAAACATAGCGG‑CCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGCa                              >  1:1424250/1‑100 (MQ=255)
                                                                     tCGCGGAAGCTGTCGGAAACATAGCGGCCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGCat                             <  2:501215/100‑1 (MQ=255)
                                                                        cgGAAGCTGTCGGAAACATAGCGG‑CCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGCATAtt                          <  2:2632619/99‑1 (MQ=255)
                                                                         ggAAGCTGTCGGAAACATAGCGG‑CCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGCATAttcg                        <  1:3684930/100‑1 (MQ=255)
                                                                            aGCTGTCGGAAACATAGCGG‑CCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGCATAttcgttc                     <  1:2423782/100‑1 (MQ=255)
                                                                            aGCTGTCGGAAACATAGCGG‑CCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGCATAttcgttc                     <  2:1779852/100‑1 (MQ=255)
                                                                             gCTGTCGGAAACATAGCGG‑CCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGCATAttcgttcg                    <  1:3070668/100‑1 (MQ=255)
                                                                             gCTGTCGGAAACATAGCGG‑CCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGCATAttcgttcg                    <  1:1454061/100‑1 (MQ=255)
                                                                             gCTGTCGGAAACATAGCGG‑CCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGCATAttcgttcg                    <  2:1991427/100‑1 (MQ=255)
                                                                              cTGTCGGAAACATAGCGGCCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGCATAtgcgttc                     >  1:3720755/1‑99 (MQ=255)
                                                                               tGTCGGAAACATAGCGG‑CCGCGCCACGGAAGCTGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGCATATTCGTTCGAt                  <  2:4089269/100‑1 (MQ=255)
                                                                                gTCGGAAACATAGCGG‑CCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGCATATTCGTTCGAtt                 <  1:1914820/100‑1 (MQ=255)
                                                                                  cGGAAACATAGCGGCCCGCGCCACGGAAGCCGAGCGTCGGGTTCTCTTTATCTGGCTCGTTACGCTCACCACCGACCAGGTTGGCATATTCGTTCGAttt                <  1:3876095/100‑1 (MQ=255)
                                                                                  cGGAAACATAGCGGCCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGCATATTCGTTCGAttt                <  2:2334584/100‑1 (MQ=255)
                                                                                       aCATAGCGG‑CCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGCATATTCGTTCGATTTAAAATc          <  1:3974271/100‑1 (MQ=255)
                                                                                        cATAGCGGCCCGCGCCACGGAAGCCGCGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGACTATTCGTTCGATTTaaaa            >  2:3600104/1‑98 (MQ=255)
                                                                                        cATAGCGG‑CCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGCATATTCGTTCGATTTAAAATCa         <  2:547393/100‑1 (MQ=255)
                                                                                         aTAGCGGCCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGCATATTCGTTCGATTTAAAATCa         <  1:2594229/100‑1 (MQ=255)
                                                                                         aTAGCGG‑CCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGCATATTCGTTCGATTTAAAATCag        >  2:2303963/1‑100 (MQ=255)
                                                                                          tggCGGCCCGCGCCACGCAAGCCGAGCATCGGGTTCTCTTCTTCTGGCTCGTAACGCTCGCAACCGACCAGGTTGGCATATTCGTTCGATTTAAAATCag        <  2:1855024/98‑1 (MQ=255)
                                                                                          tAGCGGCCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGCATATTCGTTCGATTTAAAATCag        >  1:250422/1‑100 (MQ=255)
                                                                                          tAGCGGCCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGCATATTCGTTCGATTTAAAATCag        <  2:2771276/100‑1 (MQ=255)
                                                                                          tAGCGGCCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGCATATTCGTTCGATTTAAAATCag        <  1:185189/100‑1 (MQ=255)
                                                                                           aGCGG‑CCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACTGACCAGGTTGGCATATTCGTTCGATTTAAAATCagag      <  2:2197480/100‑1 (MQ=255)
                                                                                           aGCGG‑CCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGCCTATTCGTTCGATTTAAAATCagag      >  2:478992/1‑100 (MQ=255)
                                                                                           aGCGG‑CCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGCATATTCGTTCGATTTAAAATCagag      <  1:3244715/100‑1 (MQ=255)
                                                                                           aGCGG‑CCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGCATATTCGTTCGATTTAAAATCagag      <  2:44970/100‑1 (MQ=255)
                                                                                           aGCGG‑CCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGCATATTCGTTCGATTTAAAATCagag      <  1:3038370/100‑1 (MQ=255)
                                                                                           aGCGG‑CCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGCATATTCGTTCGATTTAAAATCagag      >  1:947526/1‑100 (MQ=255)
                                                                                            gCGGCCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGCATATTCGTTCGATTTAAAATCagag      >  1:108773/1‑100 (MQ=255)
                                                                                                cccGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGCATATTCGTTCGATTTAAAATCAGAGAGAc  >  1:2020759/1‑100 (MQ=255)
                                                                                                |                                                                                                   
ATCTCAACGTTGGTCAGTCCCATGTCGTTGCGCACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGGCCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGCATATTCGTTCGATTTAAAATCAGAGAGAC  >  NC_000913/1785185‑1785380

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑