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% 1037.5 / 1.6 235intergenic (‑315/+44)ygcB/cysHhypothetical protein/3'‑phosphoadenosine 5'‑phosphosulfate reductase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base . (124/107);  minor base A (2/0);  total (127/108)
Fisher's exact test for biased strand distribution p-value = 5.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 266 total aligned reads displayed.

AGCCTCAAGCACCCGACGCCATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTAAGCCAAAGAAAC  >  NC_000913/2885465‑2885644
                                                                                            |                                                                                       
aGCCTCAAGCACCCGACGCCATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGc                                                                                         >  1:4143442/1‑93 (MQ=255)
 gCCTCAAGCACCCGACGCCATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGcc                                                                                        >  1:7516352/1‑93 (MQ=255)
   cTCAAGCACCCGACGCCATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGcccc                                                                                      <  2:6203721/93‑1 (MQ=255)
    tCAAGCACCCGACGCCATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGcccca                                                                                     >  2:189935/1‑92 (MQ=255)
    tCAAGCACCCGACGCCATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGcccca                                                                                     >  1:7842673/1‑92 (MQ=255)
     cAAGCCCCCGACGCCATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGccccat                                                                                    <  1:142860/93‑3 (MQ=255)
     cAAGCACCCGACGCCATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGccccat                                                                                    >  1:4320665/1‑91 (MQ=255)
      aaGCACCCGACGCCATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAACACGccccatc                                                                                   >  1:1153961/1‑90 (MQ=255)
          aCCCGACGCCATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCAttt                                                                              >  1:3041851/1‑93 (MQ=255)
          aCCCGACGCCATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCAttt                                                                              >  1:385431/1‑93 (MQ=255)
           cccGACGCCATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTc                                                                             <  1:670126/93‑1 (MQ=255)
           cccGACGCCATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTc                                                                             <  1:3316921/93‑1 (MQ=255)
           cccGACGCCATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTc                                                                             <  2:5034043/93‑1 (MQ=255)
            ccGACGCCATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCt                                                                            <  1:971095/93‑1 (MQ=255)
            ccGACGCCATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCt                                                                            >  2:7381064/1‑93 (MQ=255)
             cGACGCCATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTg                                                                           >  1:2408291/1‑93 (MQ=255)
             cGACGCCATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTg                                                                           <  1:4716538/93‑1 (MQ=255)
               aCGCCATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGAc                                                                         >  1:5165945/1‑93 (MQ=255)
                cGCCATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACa                                                                        <  2:2679083/93‑1 (MQ=255)
                 gCCATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACag                                                                       >  1:2344712/1‑93 (MQ=255)
                 gCCATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACACCCCATCATTTCTGACag                                                                       >  1:504849/1‑93 (MQ=255)
                  ccATTCTATACTTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACagg                                                                      >  2:1720470/1‑92 (MQ=255)
                  ccATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACaga                                                                      <  1:482372/93‑1 (MQ=255)
                   cATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACagag                                                                     <  1:7697257/93‑1 (MQ=255)
                    aTTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAACAACGCCCCATCATTTCTGACAGAgg                                                                    >  2:3719655/1‑93 (MQ=255)
                    aTTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAgg                                                                    >  1:7425389/1‑93 (MQ=255)
                     ttCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCTCAAAATAAAAACGCCCCATCATTTCTGACAGAggg                                                                   >  1:2253385/1‑92 (MQ=255)
                     ttCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGc                                                                   >  2:1274658/1‑93 (MQ=255)
                     ttCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGc                                                                   <  1:204430/93‑1 (MQ=255)
                     ttCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGc                                                                   >  1:4571313/1‑93 (MQ=255)
                        tataCCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGtt                                                                <  2:7027695/93‑1 (MQ=255)
                          taCCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTa                                                              <  2:7459863/93‑1 (MQ=255)
                          taCCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTa                                                              <  2:5229820/93‑1 (MQ=255)
                            ccTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAAt                                                            <  1:2291836/93‑1 (MQ=255)
                             cTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAAtt                                                           >  1:947196/1‑93 (MQ=255)
                              tGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAAttt                                                          >  1:6148460/1‑93 (MQ=255)
                               taTAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTg                                                         <  1:6852994/92‑1 (MQ=255)
                                aTAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGt                                                        <  2:836621/93‑1 (MQ=255)
                                 tAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTc                                                       <  2:3203725/93‑1 (MQ=255)
                                 tAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTc                                                       <  1:6511280/93‑1 (MQ=255)
                                 tAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTc                                                       >  2:2960364/1‑93 (MQ=255)
                                  aaTTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTcc                                                      <  2:1093359/93‑1 (MQ=255)
                                   aTTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCg                                                     >  2:4727083/1‑93 (MQ=255)
                                    ttcttcCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGGCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCgg                                                    >  2:4806085/1‑93 (MQ=255)
                                     tcttcCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGc                                                   >  1:1924614/1‑93 (MQ=255)
                                       ttcCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCaa                                                 <  2:3877168/93‑1 (MQ=255)
                                        tcCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAat                                                >  1:5485108/1‑93 (MQ=255)
                                         cccGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAata                                               <  1:6884334/93‑1 (MQ=255)
                                         cccGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAata                                               <  2:5544088/93‑1 (MQ=255)
                                           cgcgGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATAtt                                             <  1:733559/93‑1 (MQ=255)
                                           cgcgGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATAtt                                             >  2:7849276/1‑93 (MQ=255)
                                              ggcgTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTAc                                          <  1:13122/93‑1 (MQ=255)
                                              ggcgTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTAc                                          <  1:1713897/93‑1 (MQ=255)
                                              ggcgTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTAc                                          <  2:5382683/93‑1 (MQ=255)
                                               gcgTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTAcc                                         <  2:441681/93‑1 (MQ=255)
                                                cgTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTAccc                                        <  2:3261059/93‑1 (MQ=255)
                                                 gTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCt                                       <  1:4922128/93‑1 (MQ=255)
                                                 gTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCt                                       <  2:1200496/93‑1 (MQ=255)
                                                  ttCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCtt                                      >  2:6148986/1‑93 (MQ=255)
                                                    cACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCg                                    >  1:4204228/1‑93 (MQ=255)
                                                    cACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCg                                    >  1:5776886/1‑93 (MQ=255)
                                                     aCGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCgt                                   >  1:6381105/1‑93 (MQ=255)
                                                      cGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCgtg                                  <  1:901712/93‑1 (MQ=255)
                                                      cGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCgtg                                  >  2:6868640/1‑93 (MQ=255)
                                                       gTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGAGAGAGGCGTTTAATTTGTCCGGCAATATTTCCCCTTCgtga                                 >  1:5633065/1‑92 (MQ=255)
                                                       gTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCgtga                                 >  2:5462548/1‑92 (MQ=255)
                                                         ttAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTaa                               >  2:3556711/1‑93 (MQ=255)
                                                         ttAAAATCAGGAATGCCGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTaa                               >  1:6676210/1‑93 (MQ=255)
                                                           aaaaTCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAAcc                             >  1:1117397/1‑93 (MQ=255)
                                                             aaTCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCca                           <  1:5670270/93‑1 (MQ=255)
                                                              aTCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCcac                          <  1:3654439/93‑1 (MQ=255)
                                                               tCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCcaca                         <  1:1027005/93‑1 (MQ=255)
                                                               tCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTCCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCcaca                         <  2:6451411/93‑1 (MQ=255)
                                                                cAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACAt                        >  2:7920319/1‑93 (MQ=255)
                                                                 aGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACCtt                       >  1:632757/1‑93 (MQ=255)
                                                                 aGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACAtt                       >  1:1570920/1‑93 (MQ=255)
                                                                  ggAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTc                      >  2:268664/1‑93 (MQ=255)
                                                                      tgttgctCGCAAAATAAAAACTTCCCATCATTTCTGTCAGAGGCGTTTAATTTGTCTGGCAATATTTACCCTTCGTGTAACTCATATTCTCtt                  >  2:7246770/4‑93 (MQ=255)
                                                                      tgctgcTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCCCtt                  >  2:586387/1‑93 (MQ=255)
                                                                      tgctgcTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCCCtt                  >  1:3831587/1‑93 (MQ=255)
                                                                       gctgctCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCCCttt                 >  2:6594431/1‑93 (MQ=255)
                                                                        ctgctCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCCCtttt                >  1:1808418/1‑93 (MQ=255)
                                                                         tgctCGCAAAATAAAAACGACCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCCCttttt               >  2:7708530/1‑93 (MQ=255)
                                                                          gctCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTa              <  1:2199112/93‑1 (MQ=255)
                                                                           ctCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTaa             >  1:3428463/1‑93 (MQ=255)
                                                                            tCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTAAg            >  2:704683/1‑93 (MQ=255)
                                                                            tCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTAAg            <  2:5621618/93‑1 (MQ=255)
                                                                             cGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTAAGc           >  1:5920247/1‑93 (MQ=255)
                                                                             cGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTAAGc           >  2:2037795/1‑93 (MQ=255)
                                                                              gCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACACACATTCCCTTTTTAAGca          >  2:2346657/1‑92 (MQ=255)
                                                                                aaaaTAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTAAGCCaa        >  1:4862000/1‑93 (MQ=255)
                                                                                 aaaTAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTAAGCCaaa       >  1:3102843/1‑93 (MQ=255)
                                                                                 aaaTAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTAAGCCaaa       <  2:5354011/93‑1 (MQ=255)
                                                                                 aaaTAAAAACACCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTAAGCCaaa       >  2:413407/1‑93 (MQ=255)
                                                                                     aaaaaCGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTAAGCCAAAGaaa   <  1:2679827/93‑1 (MQ=255)
                                                                                     aaaaaCGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTAAGCCAAAGaaa   <  1:1052160/93‑1 (MQ=255)
                                                                                     aaaaaCCCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTAAGCCAAAGaaa   <  1:6207249/93‑1 (MQ=255)
                                                                                      aaaaCGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTAAGCCAAAGAAAc  >  2:6500218/1‑93 (MQ=255)
                                                                                      aaaaCGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTAAGCCAAAGAAAc  <  1:5340958/93‑1 (MQ=255)
                                                                                      aaaaCGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTAAGCCAAAGAAAc  >  2:3981794/1‑93 (MQ=255)
                                                                                            |                                                                                       
AGCCTCAAGCACCCGACGCCATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTAAGCCAAAGAAAC  >  NC_000913/2885465‑2885644

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