Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 4,371,467 Δ1 bp 41.1% intergenic (‑202/+104) dcuA ← / ← aspA C4‑dicarboxylate antiporter/aspartate ammonia‑lyase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0077794,371,4670G→.41.1% 0.0 0.0 ‑1.1214intergenic (‑202/+104)dcuA/aspAC4‑dicarboxylate antiporter/aspartate ammonia‑lyase
Reads supporting (aligned to +/- strand):  new base (41/47):  ref base (57/63):  total (101/113)
Fisher's exact test strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 8.43e-02

Only 100 of 217 total aligned reads displayed.

CGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTA  >  NC_007779/4371369‑4371567
                                                                                                  |                                                                                                    
cGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGtt                                                                                                    <  1:8048998/101‑1 (MQ=255)
 gCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTTttt                                                                                                   >  1:499095/1‑101 (MQ=255)
   cTGCAAAACCAATACTTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGttttt                                                                                                 >  1:6642634/1‑101 (MQ=255)
   cTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTTttttt                                                                                                 >  1:8769618/1‑101 (MQ=255)
    tGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTTttttta                                                                                                <  1:2046774/101‑2 (MQ=255)
    tGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGtttttt                                                                                                >  1:3431457/1‑101 (MQ=255)
     gCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTa                                                                                               >  1:3081566/1‑101 (MQ=255)
     gCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTa                                                                                               >  1:7811943/1‑101 (MQ=255)
     gCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTa                                                                                               <  1:1593491/101‑1 (MQ=255)
     gCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTa                                                                                               <  1:5785767/101‑1 (MQ=255)
       aaaaCCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAAc                                                                                             >  1:2236799/1‑101 (MQ=255)
       aaaaCCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTT‑TTTTTTAACa                                                                                            <  1:9123284/101‑1 (MQ=255)
        aaaCCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTT‑TTTTTTAACaa                                                                                           >  1:9035886/1‑101 (MQ=255)
        aaaCCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTT‑TTTTTTAACaa                                                                                           <  1:1058231/101‑1 (MQ=255)
         aaCCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACaa                                                                                           <  1:747744/101‑1 (MQ=255)
           ccAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTT‑TTTTTTAACAAGtt                                                                                        <  1:4673382/101‑1 (MQ=255)
             aaTACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTg                                                                                       <  1:5644759/101‑1 (MQ=255)
             aaTACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTg                                                                                       >  1:5541460/1‑101 (MQ=255)
             aaTACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTT‑TGTTTTGACAAGTTGa                                                                                      >  1:7072887/1‑101 (MQ=255)
               tACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGat                                                                                     >  1:3498753/1‑101 (MQ=255)
                  cTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTT‑TTTTTTAACAAGTTGATATTa                                                                                 >  1:5761485/1‑101 (MQ=255)
                  cTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTT‑TTTTTTAACAAGTTGATATTa                                                                                 <  1:5485798/101‑1 (MQ=255)
                  cTATTCCCCCCAATCTGGCCCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTT‑TTTTTTAACAAGTTGATATTa                                                                                 <  1:3363828/101‑1 (MQ=255)
                  atATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGGATATTAGCCTTCCTTGTTTTTTAACAAGTTGATAtt                                                                                  <  1:2115684/100‑1 (MQ=255)
                    aTTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTT‑TTTTTTAACAAGTTGATATTAGa                                                                               <  1:1006481/101‑1 (MQ=255)
                     ttCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTT‑TTTTTTAACAAGTTGATATTAGAt                                                                              <  1:5418809/101‑1 (MQ=255)
                      tCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTTACAAGTTGATATTAGAt                                                                              <  1:5970943/101‑1 (MQ=255)
                      tCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTT‑TTTTTTAACAAGTTGATATTAGAtt                                                                             >  1:4175616/1‑101 (MQ=255)
                       ccccccAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGAtt                                                                             <  1:3298530/101‑1 (MQ=255)
                       ccccccAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTT‑TTTTTTAACAAGTTGATATTAGATTg                                                                            >  1:2725665/1‑101 (MQ=255)
                          cccAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGtt                                                                          <  1:902265/101‑1 (MQ=255)
                           ccAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGCTa                                                                         <  1:9179117/101‑1 (MQ=255)
                             aaTCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTAtt                                                                       >  1:5483438/1‑101 (MQ=255)
                              aTCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTAttt                                                                      <  1:8375068/101‑1 (MQ=255)
                               tCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTT‑TTTTTTAACAAGTTGATATTAGATTGTTAttttt                                                                    <  1:1302945/101‑1 (MQ=255)
                                cTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTT‑TTTTTTAACAAGTTGATATTAGATTGTTATTTTTa                                                                   <  1:1934697/101‑1 (MQ=255)
                                 tGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTT‑TTTTTTAACAAGTTGATATTAGATTGTTATTTTTaa                                                                  <  1:509836/101‑1 (MQ=255)
                                 tGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTT‑TTTTTTAACAAGTTGATATTAGATTGTTATTTTTaa                                                                  <  1:3258803/101‑1 (MQ=255)
                                  ggCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTaa                                                                  <  1:206300/101‑1 (MQ=255)
                                        cAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACt                                                            <  1:1457369/101‑1 (MQ=255)
                                          agaagaTCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGc                                                          <  1:5872930/101‑1 (MQ=255)
                                          agaagaTCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTT‑TTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCt                                                         <  1:6607630/101‑1 (MQ=255)
                                           gaagaTCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTT‑TTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTc                                                        >  1:5687536/1‑101 (MQ=255)
                                            aagaTCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTT‑TTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTca                                                       >  1:7267171/1‑101 (MQ=255)
                                             agaTCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTT‑TTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTcac                                                      >  1:6673070/1‑101 (MQ=255)
                                                tCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACaa                                                    >  1:7692567/1‑101 (MQ=255)
                                                gcGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTT‑TTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAg                                                   <  1:9585808/100‑1 (MQ=255)
                                                 cGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAg                                                   <  1:9629429/101‑1 (MQ=255)
                                                 cGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTT‑TTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGa                                                  <  1:7865120/101‑1 (MQ=255)
                                                  gCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGa                                                  >  1:2326775/1‑101 (MQ=255)
                                                  gCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTT‑TTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGaa                                                 >  1:6964906/1‑101 (MQ=255)
                                                     aGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGaaaa                                               <  1:8084320/101‑1 (MQ=255)
                                                     aGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGaaaa                                               <  1:7962156/101‑1 (MQ=255)
                                                      gCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTT‑TTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGaaaaaa                                             <  1:1377840/101‑1 (MQ=255)
                                                       cAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTT‑TTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAg                                            <  1:6582770/101‑1 (MQ=255)
                                                        aaaaCTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAg                                            <  1:6094518/101‑1 (MQ=255)
                                                        aaaaCTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAg                                            >  1:6257971/1‑101 (MQ=255)
                                                         aaaCTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTT‑TTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGc                                          >  1:8533789/1‑101 (MQ=255)
                                                          aaCTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGc                                          <  1:7173828/101‑1 (MQ=255)
                                                          aaCTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGc                                          >  1:1651649/1‑101 (MQ=255)
                                                          aaCTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGc                                          <  1:4518684/101‑1 (MQ=255)
                                                          aaCTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGc                                          >  1:8258931/1‑101 (MQ=255)
                                                          aaCTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTT‑TTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCa                                         >  1:8349430/1‑101 (MQ=255)
                                                             tATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACg                                       <  1:6844371/101‑1 (MQ=255)
                                                             tATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACg                                       >  1:7362812/1‑101 (MQ=255)
                                                               tgatgaGTTCTACAACTAGCATATTAGCCTTCCTT‑TTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCa                                    >  1:6954993/1‑101 (MQ=255)
                                                                gatgaGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCa                                    >  1:6421444/1‑101 (MQ=255)
                                                                gatgaGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCa                                    <  1:5949756/101‑1 (MQ=255)
                                                                gatgaGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCa                                    <  1:7434084/101‑1 (MQ=255)
                                                                gatgaGTTCTACAACTAGCATATTAGCCTTCCTT‑TTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCAt                                   >  1:1641049/1‑101 (MQ=255)
                                                                 atgaGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCAt                                   <  1:2536754/101‑1 (MQ=255)
                                                                 atgaGTTCTACAACTAGCATATTAGCCTTCCTT‑TTTTTTAAAAAGTT‑ATATTAGATTTGTATTTTTTAGGTAATTGTTAAAAGAAAAAAGGGACGTCATCt                                 >  1:535677/1‑101 (MQ=255)
                                                                  tgaGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATc                                  <  1:7351767/101‑1 (MQ=255)
                                                                  tgaGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATc                                  >  1:8660928/1‑101 (MQ=255)
                                                                   gaGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCt                                 >  1:8849847/1‑101 (MQ=255)
                                                                   gaGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCt                                 >  1:8019334/1‑101 (MQ=255)
                                                                   gaGTTCTACAACTAGCATATTAGCCTTCCTT‑TTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTTCTGCTCACAAGAAAAAAGGCACGTCATCTg                                <  1:9326433/101‑1 (MQ=255)
                                                                     gTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGAACGTCATCTGa                               >  1:5885385/1‑101 (MQ=255)
                                                                     gTTCTACAACTAGCATATTAGCCTTCCTT‑TTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGAc                              <  1:3645348/101‑1 (MQ=255)
                                                                       tCTACAACTAGCATATTAGCCTTCCTT‑TTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGt                            >  1:6406574/1‑101 (MQ=255)
                                                                        cTACAACTAGCATATTAGCCTTCCTT‑TTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTg                           <  1:1779505/101‑1 (MQ=255)
                                                                           cAACTAGCATATTAGCCTTCCTT‑TTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCt                        >  1:2764825/1‑101 (MQ=255)
                                                                              cTAGCATATTAGCCTTCCTT‑TTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCtttt                     >  1:5368094/1‑101 (MQ=255)
                                                                               tAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTCCTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCtttt                     <  1:4614995/101‑1 (MQ=255)
                                                                               tAGCATATTAGCCTTCCTT‑TTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCttttt                    <  1:4824265/101‑1 (MQ=255)
                                                                                  cATATTAGCCTTCCTT‑TTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTAt                 <  1:3880387/101‑1 (MQ=255)
                                                                                   atatTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTAt                 >  1:3702134/1‑101 (MQ=255)
                                                                                   atatTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTAt                 >  1:9283334/1‑101 (MQ=255)
                                                                                   atatTAGCCTTCCTT‑TTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTACTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTAtt                <  1:5256657/101‑1 (MQ=255)
                                                                                    tatTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTAtt                <  1:4135840/101‑1 (MQ=255)
                                                                                     atTAGCCTTCCTT‑TTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTg              >  1:1333319/1‑101 (MQ=255)
                                                                                      ttAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTg              >  1:2417796/1‑101 (MQ=255)
                                                                                       tAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGt             <  1:6340665/101‑1 (MQ=255)
                                                                                        aGCCTTCCTT‑TTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGtac           <  1:7138027/101‑1 (MQ=255)
                                                                                         gCCTTCCTT‑TTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGtact          >  1:9501260/1‑101 (MQ=255)
                                                                                            ttcctt‑TTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTAcc       >  1:8591898/1‑101 (MQ=255)
                                                                                            ttcctt‑TTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTAcc       >  1:4909646/1‑101 (MQ=255)
                                                                                               cttGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCt     >  1:2673190/1‑101 (MQ=255)
                                                                                                 tGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGt   <  1:2886219/101‑1 (MQ=255)
                                                                                                  gTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTa  <  1:2195440/101‑1 (MQ=255)
                                                                                                  |                                                                                                    
CGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTA  >  NC_007779/4371369‑4371567

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.