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)
   tcAACGTTGGTCAGTCCCATGTCGTTGCGCACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGGCCGCGcc                                                                                              >  2:1232396/1‑100 (MQ=255)
      aCGTTGGTCAGTCCCATGTCGTTGCGCACACGTTTGACTGCTTCACACTCCAGCGCGAAACAGGCGCGGAAGCTGTCGGAAACATAGCGGCCCGCGCCAc                                                                                            <  2:2593046/100‑1 (MQ=255)
      aCGTTGGTCAGTCCCATGTCGTTGCGCACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGACGCGGAAGCTGTCGGAAACATAGCGGCCGCGCCACg                                                                                           >  1:3937309/1‑100 (MQ=255)
          tGGTCAGTCCCATGTCGTTGCGCACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGGCCGCGCCACGGAAg                                                                                       <  1:2712109/100‑1 (MQ=255)
           ggTCAGTCCCATGTCGTTGCGCACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGGCCGCGCCACGGAAGc                                                                                      >  1:4054354/1‑100 (MQ=255)
           ggTCAGTCCCATGTCGTTGCGCACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGGCCGCGCCACGGAAGc                                                                                      >  1:2845731/1‑100 (MQ=255)
            gTCAGTCCCATGTCGTTGCGCACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGGCCGCGCCACGGAAGcc                                                                                     >  1:3107187/1‑100 (MQ=255)
            gTCAGTCCCATGTCGTTGCGCACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGGCCGCGCCACGGAAGcc                                                                                     >  1:2197493/1‑100 (MQ=255)
            gTCAGTCCCATGTCGTTGCGCACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGGCCGCGCCACGGAAGcc                                                                                     >  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)
                     aTGTCGTTGCGCACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGGCCGCGCCACGGAAGCCGAGCATCgg                                                                            <  1:3120934/100‑1 (MQ=255)
                      tGTCGTTGCGCACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGGCCGCGCCACGGAAGCCGAGCATCggg                                                                           <  2:2724978/100‑1 (MQ=255)
                      tGTCGTTGCGCACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGGCCGCGCCACGGAAGCCGAGCATCggg                                                                           <  1:1981846/100‑1 (MQ=255)
                       gTCGTTGCGCACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGGCCCGCGCCACGGAAGCCGAGCATCggg                                                                           <  2:616126/100‑1 (MQ=255)
                       gTCGTTGCGCACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGGCCGCGCCACGGAAGCCGAGCATCGGGt                                                                          <  2:2149638/100‑1 (MQ=255)
                       gTCGTTGCGCACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGGACGCGCCACGGAAGCCGAGCATCggg                                                                           >  2:4051769/1‑99 (MQ=255)
                        tCGTTGCGCACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGGCCGCGCCACGGAAGCCGAGCATCGGGtt                                                                         <  1:2181645/100‑1 (MQ=255)
                        tCGGTGCGCACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGGCCCGCGCCACGGAAGCCGAGCATCGGGt                                                                          <  2:51494/100‑1 (MQ=255)
                             gcgcACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGGCCGCGCCACGGAAGCCGAGCATCGGGTTCTCtt                                                                    >  2:2320103/1‑100 (MQ=255)
                                cacaCGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGGCCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCa                                                                  >  2:889052/1‑100 (MQ=255)
                                cacaCGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGGCCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCa                                                                  >  1:3100940/1‑100 (MQ=255)
                                cacaCGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGGCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCAt                                                                 >  1:2508545/1‑100 (MQ=255)
                                  cacGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGGCCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTc                                                                   >  1:1548531/1‑97 (MQ=255)
                                  cacGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGGCCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATc                                                                >  1:3554862/1‑100 (MQ=255)
                                    cGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGGCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTgg                                                             >  1:1698288/1‑100 (MQ=255)
                                     gTTTCCCTGCTTCACACCCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGGCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGc                                                            <  2:3633290/100‑1 (MQ=255)
                                      tttCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGACACATAGCGGCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTGCATCTCGc                                                            >  2:68125/1‑99 (MQ=255)
                                      tttCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGGCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCt                                                           <  1:593759/100‑1 (MQ=255)
                                       ttCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGGCCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCt                                                           <  1:3178249/100‑1 (MQ=255)
                                       tccACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGGCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTc                                                          <  2:2508531/98‑1 (MQ=255)
                                        tCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGGCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCg                                                         >  1:2276356/1‑100 (MQ=255)
                                         cACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGGCCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCg                                                         >  2:3751028/1‑100 (MQ=255)
                                         cACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGGCCGCGCCACGTAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGt                                                        >  1:311895/1‑100 (MQ=255)
                                            tGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGGCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAAc                                                     <  1:3606324/100‑1 (MQ=255)
                                            tGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGGCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAAc                                                     <  2:1049188/100‑1 (MQ=255)
                                            tGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGGCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAAc                                                     >  1:1037489/1‑100 (MQ=255)
                                             gCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGGCCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAAc                                                     >  1:2885259/1‑100 (MQ=255)
                                             gCTTCACACTCCAGAGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGGCCGCGCCACGGAAGACGAGCATCGGGTTCTCTTCATCTGGCTCGTaa                                                      >  2:3499689/1‑98 (MQ=255)
                                             acTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGGCCGCGGCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAAc                                                     >  2:2914153/2‑99 (MQ=255)
                                              cTTCACATTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCTGCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGc                                                   <  2:3607429/100‑1 (MQ=255)
                                                 cacaCTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGGCCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTc                                                 <  2:163685/100‑1 (MQ=255)
                                                 cacaCTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGGCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTca                                                >  2:208183/1‑100 (MQ=255)
                                                 cacaCTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGGCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTca                                                <  1:1531916/100‑1 (MQ=255)
                                                 cacaCTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGGCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTca                                                <  1:738843/100‑1 (MQ=255)
                                                   cacTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGGCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTcacc                                              >  1:2853819/1‑100 (MQ=255)
                                                      tCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGGCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCaccacc                                           <  1:2413137/100‑1 (MQ=255)
                                                      tCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGGCCGCGCCACGGAAGCCGAGCAACGGGTTCTCTTCATCTGGCTCGTAAAGATcaccac                                            >  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)
                                                       ccAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGGCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCg                                          >  1:490759/1‑100 (MQ=255)
                                                       ccAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGGCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCg                                          >  2:3212408/1‑100 (MQ=255)
                                                         aGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGGCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGAc                                        >  1:276167/1‑100 (MQ=255)
                                                         aGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGGCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGAc                                        >  1:830908/1‑100 (MQ=255)
                                                          gcgcGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGGCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGAcc                                       <  2:3172666/100‑1 (MQ=255)
                                                             cgAAACAGTCGCGGAAGCTGTCGGAAACATAGCGGCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAgg                                    <  1:68534/100‑1 (MQ=255)
                                                              gAAACAGTCGCGGAAGCTGTCGGAAACATAGCGGCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGCCCAGGt                                   >  1:3092794/1‑100 (MQ=255)
                                                                aaCAGTCGCGGAAGCTGTCGGAAACATAGCGGCCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGtt                                  <  1:172722/100‑1 (MQ=255)
                                                                aaCAGTCGCGGAAGCTGTCGGAAACATAGCGGCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTg                                 <  2:3207857/100‑1 (MQ=255)
                                                                 aCAGTCGCGGAAGCTGTCGGAAACATAGCGGCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTgg                                <  1:2333064/100‑1 (MQ=255)
                                                                  cAGTCGCGGAAGCTGTCGGAAACATAGCGGCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGc                               <  2:3807022/100‑1 (MQ=255)
                                                                   aGTCGCGGAAGCTGTCGGAAACATAGCGGCCGGGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGCa                              <  2:1003476/100‑1 (MQ=255)
                                                                   aGTCGCGGAAGCTGTCGGAAACATAGCGGCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGCa                              >  1:1424250/1‑100 (MQ=255)
                                                                     tCGCGGAAGCTGTCGGAAACATAGCGGCCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGCat                             <  2:501215/100‑1 (MQ=255)
                                                                        cgGAAGCTGTCGGAAACATAGCGGCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGCATAtt                          <  2:2632619/99‑1 (MQ=255)
                                                                         ggAAGCTGTCGGAAACATAGCGGCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGCATAttcg                        <  1:3684930/100‑1 (MQ=255)
                                                                            aGCTGTCGGAAACATAGCGGCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGCATAttcgttc                     <  1:2423782/100‑1 (MQ=255)
                                                                            aGCTGTCGGAAACATAGCGGCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGCATAttcgttc                     <  2:1779852/100‑1 (MQ=255)
                                                                             gCTGTCGGAAACATAGCGGCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGCATAttcgttcg                    <  1:3070668/100‑1 (MQ=255)
                                                                             gCTGTCGGAAACATAGCGGCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGCATAttcgttcg                    <  1:1454061/100‑1 (MQ=255)
                                                                             gCTGTCGGAAACATAGCGGCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGCATAttcgttcg                    <  2:1991427/100‑1 (MQ=255)
                                                                              cTGTCGGAAACATAGCGGCCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGCATAtgcgttc                     >  1:3720755/1‑99 (MQ=255)
                                                                               tGTCGGAAACATAGCGGCCGCGCCACGGAAGCTGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGCATATTCGTTCGAt                  <  2:4089269/100‑1 (MQ=255)
                                                                                gTCGGAAACATAGCGGCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGCATATTCGTTCGAtt                 <  1:1914820/100‑1 (MQ=255)
                                                                                  cGGAAACATAGCGGCCCGCGCCACGGAAGCCGAGCGTCGGGTTCTCTTTATCTGGCTCGTTACGCTCACCACCGACCAGGTTGGCATATTCGTTCGAttt                <  1:3876095/100‑1 (MQ=255)
                                                                                  cGGAAACATAGCGGCCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGCATATTCGTTCGAttt                <  2:2334584/100‑1 (MQ=255)
                                                                                       aCATAGCGGCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGCATATTCGTTCGATTTAAAATc          <  1:3974271/100‑1 (MQ=255)
                                                                                        cATAGCGGCCCGCGCCACGGAAGCCGCGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGACTATTCGTTCGATTTaaaa            >  2:3600104/1‑98 (MQ=255)
                                                                                        cATAGCGGCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGCATATTCGTTCGATTTAAAATCa         <  2:547393/100‑1 (MQ=255)
                                                                                         aTAGCGGCCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGCATATTCGTTCGATTTAAAATCa         <  1:2594229/100‑1 (MQ=255)
                                                                                         aTAGCGGCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGCATATTCGTTCGATTTAAAATCag        >  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)
                                                                                           aGCGGCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACTGACCAGGTTGGCATATTCGTTCGATTTAAAATCagag      <  2:2197480/100‑1 (MQ=255)
                                                                                           aGCGGCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGCCTATTCGTTCGATTTAAAATCagag      >  2:478992/1‑100 (MQ=255)
                                                                                           aGCGGCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGCATATTCGTTCGATTTAAAATCagag      <  1:3244715/100‑1 (MQ=255)
                                                                                           aGCGGCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGCATATTCGTTCGATTTAAAATCagag      <  2:44970/100‑1 (MQ=255)
                                                                                           aGCGGCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGCATATTCGTTCGATTTAAAATCagag      <  1:3038370/100‑1 (MQ=255)
                                                                                           aGCGGCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGCATATTCGTTCGATTTAAAATCagag      >  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: