Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 1,539,898 T→A 3.6% S224C (AGC→TGC)  narY ← nitrate reductase 2 (NRZ) subunit beta

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0077791,539,8980T→A 3.6% 12.9 ‑1.8 ‑0.0246S224C (AGC→TGC) narYnitrate reductase 2 (NRZ) subunit beta
Reads supporting (aligned to +/- strand):  new base (2/7):  ref base (119/117):  total (121/125)
Fisher's exact test strand distribution p-value = 1.72e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.99e-01

Only 100 of 255 total aligned reads displayed.

ACGGTCGGTTGACCGGACTCAATTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAG  >  NC_007779/1539798‑1539998
                                                                                                    |                                                                                                    
acGGTCGGTTGACCGGACTCAATTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCt                                                                                                      >  1:1577527/1‑101 (MQ=255)
acGGTCGGTTGACCGGACTCAATTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCt                                                                                                      >  1:13441945/1‑101 (MQ=255)
 cggtcggtTGACCGGACTCAATTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCtt                                                                                                     >  1:3013210/1‑101 (MQ=255)
  ggtcggtTGACCGGACTCAATTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTa                                                                                                    >  1:8447160/1‑101 (MQ=255)
  ggtcggtTGACCGGACTCAATTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTa                                                                                                    <  1:749781/101‑1 (MQ=255)
   gtcggtTGACCGGACTCAATTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTAt                                                                                                   <  1:7825044/101‑1 (MQ=255)
   gtcggtTGACCGGACTCAATTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTAt                                                                                                   <  1:6556418/101‑1 (MQ=255)
     cggtTGACCGGACTCAATTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGc                                                                                                 <  1:14327844/101‑1 (MQ=255)
        ttGACCGGACTCAATTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACGACCGCTTATGcaca                                                                                              >  1:12479706/1‑101 (MQ=255)
        ttGACCGGACTCAATTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGcaca                                                                                              >  1:11915564/1‑101 (MQ=255)
         tGACCGGACTCAATTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACaa                                                                                             >  1:1734508/1‑101 (MQ=255)
         tGACCGGACTCAATTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACaa                                                                                             <  1:4327354/101‑1 (MQ=255)
          ggcctggcTCAATTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACaaa                                                                                            <  1:3839989/94‑1 (MQ=255)
          gACCGGACTCAATTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACaaa                                                                                            >  1:6063417/1‑101 (MQ=255)
          gACCGGACTCAATTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACaaa                                                                                            >  1:6490046/1‑101 (MQ=255)
           aCCGGACTCAATTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAAc                                                                                           <  1:12905889/101‑1 (MQ=255)
           aCCGGACTCAATTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAAc                                                                                           >  1:6277722/1‑101 (MQ=255)
            ccGGACTCAATTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACg                                                                                          <  1:3228974/101‑1 (MQ=255)
              ggACTCAATTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACgcc                                                                                        >  1:10397753/1‑101 (MQ=255)
                aCTCAATTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACgccag                                                                                      >  1:12798481/1‑101 (MQ=255)
                aCTCAATTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACgccag                                                                                      <  1:2833629/101‑1 (MQ=255)
                aCTCAATTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACgccag                                                                                      <  1:270594/101‑1 (MQ=255)
                    aaTTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCAc                                                                                  >  1:3389380/1‑101 (MQ=255)
                    aaTTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCAc                                                                                  >  1:4997089/1‑101 (MQ=255)
                    aaTTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCAc                                                                                  >  1:10220984/1‑101 (MQ=255)
                     aTTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACg                                                                                 >  1:4179311/1‑101 (MQ=255)
                      ttCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACgg                                                                                >  1:7235817/1‑101 (MQ=255)
                      ttCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACgg                                                                                <  1:10722645/101‑1 (MQ=255)
                      ttCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACgg                                                                                <  1:1382296/101‑1 (MQ=255)
                       tCGTGGGTAACAGAAGATGCATTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCa                                                                              >  1:11628795/1‑101 (MQ=255)
                        cGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCa                                                                              >  1:3514801/1‑101 (MQ=255)
                        cGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCa                                                                              <  1:5535264/101‑1 (MQ=255)
                        cGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCa                                                                              <  1:9763837/101‑1 (MQ=255)
                           gggTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTTCGGACAATCGCTTATGCACAAACGCCAGCCACGGCAttt                                                                           >  1:12094714/1‑101 (MQ=255)
                              tAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATc                                                                        >  1:13641602/1‑101 (MQ=255)
                               aaCAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATcc                                                                       <  1:11102179/101‑1 (MQ=255)
                                aCAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCt                                                                      >  1:2443656/1‑101 (MQ=255)
                                 cAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCTTTTATCCTg                                                                     >  1:7016338/1‑101 (MQ=255)
                                 cAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTg                                                                     >  1:2119939/1‑101 (MQ=255)
                                 cAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTg                                                                     >  1:2550117/1‑101 (MQ=255)
                                   gaagaTGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGat                                                                   >  1:11538834/1‑101 (MQ=255)
                                   gaagaTGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGat                                                                   >  1:7344888/1‑101 (MQ=255)
                                   gaagaTGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGat                                                                   <  1:13770436/101‑1 (MQ=255)
                                     agaTGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGatca                                                                 <  1:11748025/101‑1 (MQ=255)
                                      gaTGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGAACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGatcaa                                                                >  1:8972710/1‑101 (MQ=255)
                                         gCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGatcaatca                                                             <  1:14151333/101‑1 (MQ=255)
                                          cATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAg                                                            <  1:12618704/101‑1 (MQ=255)
                                             ttttCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGcac                                                         >  1:7665387/1‑101 (MQ=255)
                                              tttCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGcaca                                                        >  1:6381007/1‑101 (MQ=255)
                                              tttCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGcaca                                                        >  1:13297577/1‑101 (MQ=255)
                                               ttCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACaa                                                       <  1:5346523/101‑1 (MQ=255)
                                               ttCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACaa                                                       >  1:13928117/1‑101 (MQ=255)
                                                 cTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATg                                                     <  1:4253383/101‑1 (MQ=255)
                                                  tGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGc                                                    <  1:5068123/101‑1 (MQ=255)
                                                     cTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGt                                                 <  1:13802265/101‑1 (MQ=255)
                                                     cTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGt                                                 <  1:1136220/101‑1 (MQ=255)
                                                      ttGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGtc                                                <  1:9479365/101‑1 (MQ=255)
                                                         ccGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGtcttc                                             >  1:6745545/1‑101 (MQ=255)
                                                            cTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCttcttc                                          <  1:9416577/101‑1 (MQ=255)
                                                            cTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCttcttc                                          <  1:2888559/101‑1 (MQ=255)
                                                             ttttCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTcg                                         <  1:11912253/101‑1 (MQ=255)
                                                             ttttCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTcg                                         >  1:1979732/1‑101 (MQ=255)
                                                              tttCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTcgc                                        >  1:3452611/1‑101 (MQ=255)
                                                               ttCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTcgcg                                       <  1:6513696/101‑1 (MQ=255)
                                                                tCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGt                                      >  1:5936822/1‑101 (MQ=255)
                                                                 ccAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGtt                                     <  1:4616483/101‑1 (MQ=255)
                                                                   aGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTg                                   <  1:6826055/101‑1 (MQ=255)
                                                                     ttGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTa                                 >  1:9224065/1‑101 (MQ=255)
                                                                      tGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAg                                >  1:5194211/1‑101 (MQ=255)
                                                                      tGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAg                                >  1:8115910/1‑101 (MQ=255)
                                                                      tGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAg                                <  1:8585698/101‑1 (MQ=255)
                                                                       gAAGTAGACCTTCTTGTAGGGGCAACCGGAGACGCACATGCGCCAACCGCGGCAGCGGTCCTGGTCGACCAGCACAATGCCGTCTTCTGTGCGTTTGTaca                               <  1:5889810/101‑3 (MQ=255)
                                                                        aaGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGAt                              <  1:10085256/101‑1 (MQ=255)
                                                                         aGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATa                             <  1:1506551/101‑1 (MQ=255)
                                                                          gTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAg                            <  1:4660300/101‑1 (MQ=255)
                                                                           tAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAgc                           >  1:31701/1‑101 (MQ=255)
                                                                            aGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAgcg                          <  1:14255788/101‑1 (MQ=255)
                                                                            aGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAgcg                          >  1:8759608/1‑101 (MQ=255)
                                                                            aGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAgcg                          <  1:7471357/101‑1 (MQ=255)
                                                                              aTTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGcgcc                        <  1:5479907/101‑1 (MQ=255)
                                                                                   tttGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTg                   <  1:9174565/101‑1 (MQ=255)
                                                                                   tttGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTg                   >  1:13542449/1‑101 (MQ=255)
                                                                                   tttGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTg                   >  1:9683719/1‑101 (MQ=255)
                                                                                   tttGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTg                   <  1:5121800/101‑1 (MQ=255)
                                                                                    ttGTAGGGGCAACCGGAGACGCACATGCGCCAACCGCGGCAGCGGTCCTGGTCGACCAGCACAATGCCGTCTTCTGTGCGTTTGTACATCGcgcccgatgg                  <  1:6584987/101‑7 (MQ=255)
                                                                                    ttGTAGGGGCAACCGGAGACGCACATGCGCCAACCGCGGCAGCGGTCCTGGTCGACCAGCACAATGCCGTCTTCTGTGCGTTTGTACATCGcgcccgatgg                  <  1:8981800/101‑7 (MQ=255)
                                                                                     tGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTggg                 <  1:12537310/101‑1 (MQ=255)
                                                                                     tGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTggg                 >  1:13120352/1‑101 (MQ=255)
                                                                                        aCGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAg              <  1:168492/101‑1 (MQ=255)
                                                                                         cGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAgg             <  1:3741152/101‑1 (MQ=255)
                                                                                           gACACCCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTc           <  1:10622721/101‑1 (MQ=255)
                                                                                           gACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTc           <  1:4344549/101‑1 (MQ=255)
                                                                                            aCAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCg          <  1:10903644/101‑1 (MQ=255)
                                                                                            aCAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCg          <  1:13601691/101‑1 (MQ=255)
                                                                                              aaCCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGcc        >  1:8065359/1‑101 (MQ=255)
                                                                                                ccGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCAc      <  1:12259979/101‑1 (MQ=255)
                                                                                                ccGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCAc      >  1:11368976/1‑101 (MQ=255)
                                                                                                 cGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACg     <  1:10498978/101‑1 (MQ=255)
                                                                                                 cGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACg     >  1:957351/1‑101 (MQ=255)
                                                                                                   cTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATACCGCCGCTTGGGCAGGTCGCCACGCa   <  1:12734826/101‑1 (MQ=255)
                                                                                                    |                                                                                                    
ACGGTCGGTTGACCGGACTCAATTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAG  >  NC_007779/1539798‑1539998

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ 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.