Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,597,825 A→G intergenic (‑45/‑101) dapA ← / → b2479 dihydrodipicolinate synthase/GcvR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,597,8250AG98.0% 619.5 / 1.2 197intergenic (‑45/‑101)dapA/b2479dihydrodipicolinate synthase/GcvR
Reads supporting (aligned to +/- strand):  ref base A (1/3);  new base G (95/98);  total (96/101)
Fisher's exact test for biased strand distribution p-value = 6.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.69e-01

Only 100 of 200 total aligned reads displayed.

TTTCATCCATCGGAGTAACAATCGCGACAATACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTATGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTCCTTAAGAACAACTCACCCCTTA  >  NC_000913/2597735‑2597911
                                                                                          |                                                                                      
tttCATCCATCGGAGTAACAATCGCGACAATACTTCCCGGGAACATGGGCCATCCTCTGTGCAAAACAGTGTCTCAATGTGAGGTTTGGtgtg                                                                                      >  2:4231367/1‑93 (MQ=255)
     tccatcGGAGTAACAATCGCGACAATACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCAtt                                                                                 <  1:503946/93‑1 (MQ=255)
      ccatcGGAGTAACAATCGCGACAATACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTa                                                                                >  1:2885881/1‑93 (MQ=255)
       catcGGAGTAACAATCGCGACAATACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTaa                                                                               <  2:2725051/93‑1 (MQ=255)
        atcGGAGTAACAATCGCGACAATACTTCCCGTGAACATGGGCCATCCTGTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTaaa                                                                              >  1:2851329/1‑93 (MQ=255)
        atcGGAGTAACAATCGCGACAATACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTaaa                                                                              >  2:565187/1‑93 (MQ=255)
        atcGGAGTAACAATCGCGACAATACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTaaa                                                                              <  2:3646346/93‑1 (MQ=255)
        atcGGAGTAACAATCGCGACAATACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTaaa                                                                              >  1:3337856/1‑93 (MQ=255)
         tcGGAGTAACAATCGCGACAATACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTaaaa                                                                             <  1:1941099/93‑1 (MQ=255)
         tcGGAGTAACAATCGCGACAATACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTaaaa                                                                             >  1:2057072/1‑93 (MQ=255)
         tcGGAGTAACAATCGCGACAATACTTCCCGTGAACATGGGCCATCCTCTGGGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTaaaa                                                                             <  2:3905725/93‑1 (MQ=255)
            gAGTAACAATCGCGACAATACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAaagca                                                                          <  1:2379117/93‑1 (MQ=255)
             aGTAACAATCGCGACAATACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAaagcaa                                                                         <  1:4038424/93‑1 (MQ=255)
             aGTAACAATCGCGACAATACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAaagcaa                                                                         <  1:1253481/93‑1 (MQ=255)
               tAACAATCGCGACAATACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAaagcaagc                                                                       <  2:4574834/93‑1 (MQ=255)
               tAACAATCGCGACAATACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAaagcaagc                                                                       >  2:1884773/1‑93 (MQ=255)
               tAACAATCGCGACAATACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTATGGCATTAAaagcaagc                                                                       <  2:3862524/93‑1 (MQ=255)
                 aCAATCGCGACAATACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGcag                                                                     <  1:1278463/93‑1 (MQ=255)
                  cAATCGCGACAATACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGcaga                                                                    <  1:3914036/93‑1 (MQ=255)
                   aaTCGCGACAATACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGcagac                                                                   >  1:2822851/1‑93 (MQ=255)
                   aaTCGCGACAATACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGcagac                                                                   <  1:2937382/93‑1 (MQ=255)
                     tCGCGACAATACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGcagacag                                                                 <  2:4008658/93‑1 (MQ=255)
                       gcgACAATACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGaa                                                               >  1:1954186/1‑93 (MQ=255)
                        cgACAATACTTCCCGTGAACATGTGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAAc                                                              <  1:3139230/93‑1 (MQ=255)
                        cgACAATACTTCCCGTGAACATGGGCGATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAAc                                                              <  1:63302/93‑1 (MQ=255)
                        cgACAATACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAAc                                                              <  1:4562834/93‑1 (MQ=255)
                        cgACAATACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAAc                                                              >  1:4594236/1‑93 (MQ=255)
                        cgACAATACGTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAAc                                                              <  1:3139232/93‑1 (MQ=255)
                         gACAATACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAAcc                                                             >  1:4389388/1‑93 (MQ=255)
                         gACAATACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAAcc                                                             <  2:4066688/93‑1 (MQ=255)
                         gACAATACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAAcc                                                             >  2:4718316/1‑93 (MQ=255)
                            aaTACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGtt                                                          >  1:1034764/1‑93 (MQ=255)
                            aaTACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGtt                                                          <  1:1851171/93‑1 (MQ=255)
                            aaTACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGtt                                                          >  2:894369/1‑93 (MQ=255)
                             aTACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTc                                                         <  1:4363996/93‑1 (MQ=255)
                             aTACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTc                                                         >  2:398776/1‑93 (MQ=255)
                              tACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCt                                                        >  1:3051858/1‑93 (MQ=255)
                               tcTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTg                                                       <  2:2859558/92‑1 (MQ=255)
                               aCTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTg                                                       <  2:4016852/93‑1 (MQ=255)
                                cTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTgt                                                      >  2:3335551/1‑92 (MQ=255)
                                cTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGa                                                      <  2:206397/93‑1 (MQ=255)
                                cTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGa                                                      <  1:4467563/93‑1 (MQ=255)
                                 ttCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGAt                                                     >  1:1146465/1‑93 (MQ=255)
                                  tCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGAtt                                                    >  2:1679275/1‑93 (MQ=255)
                                   cccGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGAttg                                                   >  2:3419789/1‑93 (MQ=255)
                                   cccGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGAttg                                                   <  2:2822755/93‑1 (MQ=255)
                                   cccGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGGGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGAttg                                                   >  2:1326396/1‑93 (MQ=255)
                                    ccGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGAttgt                                                  >  1:1361300/1‑93 (MQ=255)
                                    ccGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGAttgt                                                  >  1:4670132/1‑93 (MQ=255)
                                     cGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGAttgtt                                                 >  2:3033279/1‑93 (MQ=255)
                                     cGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGAttgtt                                                 <  2:2966520/93‑1 (MQ=255)
                                     cGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGAttgtt                                                 >  1:2247460/1‑93 (MQ=255)
                                       tGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATtgttgt                                               >  1:3288110/1‑93 (MQ=255)
                                       tGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATtgttgt                                               <  2:3496401/93‑1 (MQ=255)
                                        gAACATGGGCCATCCTCTGTGCAAGCAAGTGTCTCAATGGTACGTATGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCAGATTATTGTa                                              <  1:4215442/93‑1 (MQ=255)
                                        gAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTa                                              <  1:4536580/93‑1 (MQ=255)
                                          aCATGGGCCATCCTCTGTGCAAACCAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATg                                            >  2:708288/1‑93 (MQ=255)
                                            aTGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCa                                          <  2:17179/93‑1 (MQ=255)
                                            aTGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCa                                          <  1:6947/93‑1 (MQ=255)
                                             tGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCAt                                         <  2:3257718/93‑1 (MQ=255)
                                                 ccATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTATGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGttt                                     <  2:2268896/93‑1 (MQ=255)
                                                  cATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGtttt                                    <  1:1894985/93‑1 (MQ=255)
                                                   aTCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGttttt                                   >  2:724612/1‑93 (MQ=255)
                                                   aTCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGAGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGttttt                                   <  2:4201684/93‑1 (MQ=255)
                                                     ccTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGttttttt                                 >  1:4326277/1‑93 (MQ=255)
                                                     ccTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGttttttt                                 >  2:2181023/1‑93 (MQ=255)
                                                      ctctGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGtttttttt                                >  1:1575022/1‑93 (MQ=255)
                                                      ctctGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGtttttttt                                >  1:2705556/1‑93 (MQ=255)
                                                       tctGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTa                               <  2:1523582/93‑1 (MQ=255)
                                                        ctGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTAt                              <  2:1359088/93‑1 (MQ=255)
                                                         tgtgCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATg                             <  1:3554953/93‑1 (MQ=255)
                                                          gtgCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGc                            <  1:87591/93‑1 (MQ=255)
                                                            gCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCtt                          >  2:1352009/1‑93 (MQ=255)
                                                              aaaCAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTc                        >  2:4166693/1‑93 (MQ=255)
                                                               aaCAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTcc                       <  1:3425093/93‑1 (MQ=255)
                                                               aaCAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTcc                       >  2:2653062/1‑93 (MQ=255)
                                                               aaCAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTcc                       <  1:1913742/93‑1 (MQ=255)
                                                                aCAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTCCt                      <  1:684003/93‑1 (MQ=255)
                                                                aCAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTCCt                      >  2:1054125/1‑93 (MQ=255)
                                                                  aaGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTCCTTa                    >  1:1850534/1‑93 (MQ=255)
                                                                  aaGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTCCTTa                    <  2:191597/93‑1 (MQ=255)
                                                                   aGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTCCTTaa                   <  1:4616316/93‑1 (MQ=255)
                                                                   aGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTCCTTaa                   <  2:3164836/93‑1 (MQ=255)
                                                                   aGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTCCTTaa                   <  2:2567084/93‑1 (MQ=255)
                                                                    gtgtCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTCCTTAAg                  >  2:2429055/1‑93 (MQ=255)
                                                                    gtgtCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTCCTTAAg                  <  2:2001492/93‑1 (MQ=255)
                                                                        ctcAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTCCTTAAGaaca              >  2:452609/1‑93 (MQ=255)
                                                                         tcAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCGGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTCCTTAAGaacaa             <  2:3687594/93‑1 (MQ=255)
                                                                         tcAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTCCTTAAGaacaa             >  2:1680081/1‑93 (MQ=255)
                                                                         tcAATGGTACGTTTGGTGTGGCATTAAAAGAAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTCCTTAAGaacaa             >  2:3729021/1‑93 (MQ=255)
                                                                          cAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTCCTTAAGaacaac            <  1:2991123/93‑1 (MQ=255)
                                                                          cAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTCCTTAAGaacaac            <  2:605347/93‑1 (MQ=255)
                                                                          cAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTCTTTATGCTTTCCTTAAGaacaac            >  2:2593971/1‑93 (MQ=255)
                                                                             tGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTCCTTAAGAACAACTCa         >  1:3766247/1‑93 (MQ=255)
                                                                               gTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTCCTTAAGAACAACTCAcc       >  2:1096325/1‑93 (MQ=255)
                                                                                tACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTCCTTAAGAACAACTCAccc      <  1:1507083/93‑1 (MQ=255)
                                                                                tACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTCCTTAAGAACAACTCAccc      >  2:2372900/1‑93 (MQ=255)
                                                                                tACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTCCTTAAGAACAACTCAccc      >  2:1251580/1‑93 (MQ=255)
                                                                                 aCGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTCCTTAAGAACAACTCAcacc     >  1:832335/1‑93 (MQ=255)
                                                                                    tttGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTCCTTAAGAACAACTCACCCCTTa  >  1:3550070/1‑93 (MQ=255)
                                                                                          |                                                                                      
TTTCATCCATCGGAGTAACAATCGCGACAATACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTATGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTCCTTAAGAACAACTCACCCCTTA  >  NC_000913/2597735‑2597911

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