Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,885,561 (C)5→4 intergenic (‑319/+40) ygcB ← / ← cysH hypothetical protein/3'‑phosphoadenosine 5'‑phosphosulfate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,885,5570C.99.1% 504.1 / ‑3.3 113intergenic (‑315/+44)ygcB/cysHhypothetical protein/3'‑phosphoadenosine 5'‑phosphosulfate reductase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base . (64/48);  total (65/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.31e-01

Only 100 of 130 total aligned reads displayed.

AGCCTCAAGCACCCGACGCCATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTAAGCCAAAGAAACGCGTT  >  NC_000913/2885465‑2885649
                                                                                            |                                                                                            
aGCCTCAAGCACCCGACGCCATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGc                                                                                              >  2:2197694/1‑93 (MQ=255)
aGCCTCAAGCACCCGACGCCATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGc                                                                                              <  1:577719/93‑1 (MQ=255)
aGCCTCAAGCACCCGACGCCATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGc                                                                                              <  2:3045147/93‑1 (MQ=255)
 gCCTCAAGCACCCGACGCCATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGcc                                                                                             >  1:1486530/1‑93 (MQ=255)
 gCCTCAAGCACCCGACGCCATTCTATACCTGATAATTCTTCCAGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGcc                                                                                             >  1:2511598/1‑93 (MQ=255)
   cTCAAGCACCCGACGCCATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGcccc                                                                                           >  2:1264373/1‑93 (MQ=255)
    tCAAGCACCCGACGCCATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGcccca                                                                                          >  2:2383202/1‑92 (MQ=255)
    tCAAGCACCCGACGCCATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGcccca                                                                                          >  2:708266/1‑92 (MQ=255)
    tCAAGCACCCGACGCCATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGcccca                                                                                          >  2:2975920/1‑92 (MQ=255)
      aaGCACCCGACGCCATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGccccatc                                                                                        >  2:1478405/1‑90 (MQ=255)
      aaGCACCCGACGCCATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTAAAAATCAGGAATGCTGCTCGCAAAATAAAAACGacccatc                                                                                        >  2:1554234/1‑90 (MQ=255)
       aGCACCCGACGCCATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGccccatca                                                                                       >  2:3693176/1‑89 (MQ=255)
        acACCCGACGCCATTCTATTCCTTATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCcatcat                                                                                     <  1:1468276/92‑1 (MQ=255)
          ccgcGACGTCATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCGGCTCGCAAAATAAAAACGCCCCATCAttt                                                                                   <  2:3723747/90‑1 (MQ=255)
             cGACGCCATTCTATAGCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGATCGCAAAATAAAAACGCCCCATCATTTCTg                                                                                <  1:2536150/93‑1 (MQ=255)
             cGACGCCATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTg                                                                                <  2:1940474/93‑1 (MQ=255)
                 gCCATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACag                                                                            >  2:1407964/1‑93 (MQ=255)
                  ccATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACaga                                                                           <  1:2627089/93‑1 (MQ=255)
                  ccATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACaga                                                                           >  2:3482447/1‑93 (MQ=255)
                  ccATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACaga                                                                           >  2:637509/1‑93 (MQ=255)
                  ccATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACaga                                                                           <  1:1468853/93‑1 (MQ=255)
                   cATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACagag                                                                          <  2:3339235/93‑1 (MQ=255)
                   cATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACagag                                                                          <  2:2104620/93‑1 (MQ=255)
                   cATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACagag                                                                          <  1:2113674/93‑1 (MQ=255)
                    aTTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAgg                                                                         >  2:3619940/1‑93 (MQ=255)
                    aTTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAgg                                                                         >  2:3335506/1‑93 (MQ=255)
                      tCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCg                                                                       <  1:1602073/93‑1 (MQ=255)
                          taCCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTa                                                                   >  1:3006201/1‑93 (MQ=255)
                          taCCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTa                                                                   >  1:1548222/1‑93 (MQ=255)
                          gaCCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTa                                                                   <  2:4300788/92‑1 (MQ=255)
                            ccTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAAt                                                                 <  2:1540232/93‑1 (MQ=255)
                              tGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAAttt                                                               >  1:50827/1‑93 (MQ=255)
                              tGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAAttt                                                               >  1:1095098/1‑93 (MQ=255)
                              tGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAAttt                                                               <  2:1146189/93‑1 (MQ=255)
                              tGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAAttt                                                               >  1:3774731/1‑93 (MQ=255)
                              tGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAAttt                                                               >  1:126158/1‑93 (MQ=255)
                               gATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTg                                                              <  2:1466347/93‑1 (MQ=255)
                               gATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTg                                                              >  2:3809912/1‑93 (MQ=255)
                                aTAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGt                                                             >  2:856972/1‑93 (MQ=255)
                                 tAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTc                                                            <  1:3580147/93‑1 (MQ=255)
                                 tAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTc                                                            >  1:4534582/1‑93 (MQ=255)
                                  aaTTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTcc                                                           >  2:1167161/1‑93 (MQ=255)
                                  aaTTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTcc                                                           >  1:2466555/1‑93 (MQ=255)
                                   aTTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCg                                                          <  1:3982769/93‑1 (MQ=255)
                                   aTTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCg                                                          >  1:739630/1‑93 (MQ=255)
                                   aTTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCg                                                          <  2:1884935/93‑1 (MQ=255)
                                   aTTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCg                                                          <  1:3982880/93‑1 (MQ=255)
                                      cttcCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCa                                                       <  1:528149/93‑1 (MQ=255)
                                      cttcCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCa                                                       >  2:3382875/1‑93 (MQ=255)
                                         cccGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAata                                                    >  1:2359969/1‑93 (MQ=255)
                                          ccGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAatat                                                   >  2:2013358/1‑93 (MQ=255)
                                           cgcgGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATAtt                                                  <  1:121894/93‑1 (MQ=255)
                                             cggcgTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTa                                                <  1:2135878/93‑1 (MQ=255)
                                              ggcgTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTAc                                               <  1:1573055/93‑1 (MQ=255)
                                               gcgTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTAcc                                              <  2:1820649/93‑1 (MQ=255)
                                                 gTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCt                                            >  2:872970/1‑93 (MQ=255)
                                                  ttCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCtt                                           >  1:4638268/1‑93 (MQ=255)
                                                  ttCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCtt                                           >  2:3812757/1‑93 (MQ=255)
                                                   tCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTc                                          <  2:915958/93‑1 (MQ=255)
                                                     aCGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCgt                                        >  2:1279650/1‑93 (MQ=255)
                                                     aCGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCgt                                        >  2:3887742/1‑93 (MQ=255)
                                                     aCGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCgt                                        >  1:243216/1‑93 (MQ=255)
                                                      cGTTTAAAATCAGGAATGCTGCTCGCACAATAAAAACGCCCCATCATTTCTGCCAGCGGCGTTTGATTTGTCCGGCAATATTTACCCTTCgtt                                       >  1:4302512/1‑92 (MQ=255)
                                                      cGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCgtg                                       >  1:1234285/1‑93 (MQ=255)
                                                       gTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCgtgt                                      >  1:223814/1‑93 (MQ=255)
                                                        tttAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTa                                     >  1:3677055/1‑93 (MQ=255)
                                                        tttAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTa                                     >  2:2648635/1‑93 (MQ=255)
                                                          tAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAAc                                   <  2:3893363/93‑1 (MQ=255)
                                                            aaaTCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAAccc                                 >  1:580418/1‑93 (MQ=255)
                                                             aaTCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCca                                <  1:2278184/93‑1 (MQ=255)
                                                              aTCAGGAATGGTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCcac                               <  2:2444827/93‑1 (MQ=255)
                                                              aTCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCcac                               >  2:726239/1‑93 (MQ=255)
                                                               tCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCcaca                              >  1:3879608/1‑93 (MQ=255)
                                                                  ggAATGCTGCTCGCAAAATAAAAACGCCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACAtt                            >  2:4397762/1‑93 (MQ=255)
                                                                  ggAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTc                           <  1:3272141/93‑1 (MQ=255)
                                                                  ggAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTc                           >  1:2857284/1‑93 (MQ=255)
                                                                  ggAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTc                           >  1:2656775/1‑93 (MQ=255)
                                                                    aaTGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTccc                         <  2:749430/93‑1 (MQ=255)
                                                                    aaTGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTccc                         <  1:1567770/93‑1 (MQ=255)
                                                                    aaTGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTccc                         >  1:2218934/1‑93 (MQ=255)
                                                                      tgctgcTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCCCtt                       >  2:603982/1‑93 (MQ=255)
                                                                        ctgctCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCCCtttt                     <  1:3630341/93‑1 (MQ=255)
                                                                          gctCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTa                   <  2:640569/93‑1 (MQ=255)
                                                                          gctCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTa                   >  2:279661/1‑93 (MQ=255)
                                                                          gctCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTa                   <  1:4117391/93‑1 (MQ=255)
                                                                           ctCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTaa                  <  1:4661991/93‑1 (MQ=255)
                                                                           ctCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTaa                  <  1:1816056/93‑1 (MQ=255)
                                                                           ctCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTaa                  >  2:620180/1‑93 (MQ=255)
                                                                            tCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTAAg                 >  2:2710141/1‑93 (MQ=255)
                                                                              gCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTAAGcc               >  2:3714489/1‑93 (MQ=255)
                                                                                 aaaTAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTAAGCCaaa            >  1:4171594/1‑93 (MQ=255)
                                                                                  aaTAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTAAGCCAAAg           <  2:2140283/93‑1 (MQ=255)
                                                                                   aTAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTAAGCCAAAGa          >  1:803376/1‑93 (MQ=255)
                                                                                    tAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTAAGCCAAAGaa         <  1:841623/93‑1 (MQ=255)
                                                                                    tAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTAAGCCAAAGaa         >  2:4472259/1‑93 (MQ=255)
                                                                                     aaaaaCGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTAAGCCAAAGaaa        >  1:3631912/1‑93 (MQ=255)
                                                                                      aaaaCGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTAAGCCAAAGAAAc       >  2:4184657/1‑93 (MQ=255)
                                                                                      aaaaCGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTAAGCCAAAGAAAc       >  2:2796896/1‑93 (MQ=255)
                                                                                      aaaaCGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTAAGCCAAAGAAAc       >  1:3759012/1‑93 (MQ=255)
                                                                                       aaaCGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTAAGCCAAAGAAAcg      >  1:381647/1‑93 (MQ=255)
                                                                                            |                                                                                            
AGCCTCAAGCACCCGACGCCATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTAAGCCAAAGAAACGCGTT  >  NC_000913/2885465‑2885649

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: