Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 1,485,434 T→C 7.7% R201R (CGT→CGC)  hrpA → ATP‑dependent helicase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0077791,485,4340T→C 7.7% 40.4 0.0 ‑0.3232R201R (CGT→CGC) hrpAATP‑dependent helicase
Reads supporting (aligned to +/- strand):  new base (8/10):  ref base (100/114):  total (108/124)
Fisher's exact test strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 5.11e-01

Only 100 of 243 total aligned reads displayed.

TGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCGCGGCGTCCTGACCTAAAAATCATTATCACTTCCGCGACTATCGACCCGGAACGCTTTTCGCGCCACTTTAATAATGCGCCGATTATTGAAGTCTCCGGTCGGAC  >  NC_007779/1485334‑1485532
                                                                                                    |                                                                                                  
 gctgATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCGCGGCGTc                                                                                                   >  1:4858426/1‑101 (MQ=255)
   tgATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCGCGGCGTCCt                                                                                                 <  1:1297600/101‑1 (MQ=255)
     aTGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCGCGGCGTCCTGa                                                                                               <  1:4027613/101‑1 (MQ=255)
      tGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCGCGGCGTCCTGAc                                                                                              >  1:12246987/1‑101 (MQ=255)
       tgaGTACGACACCATCATCATCGATGAGGCGCACGAACGCAGCCTCAACATCGACTTCCTGCTGGGATATTTAAAGCATATTCTGCCGCAGCGCCCAGatc                                                                                             >  1:1778138/3‑99 (MQ=255)
       gCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCGCGGCGTCCTGAcc                                                                                             >  1:8267182/1‑101 (MQ=255)
        cAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCGCGGCGTCCTGACCt                                                                                            <  1:13250012/101‑1 (MQ=255)
         aGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCGCGGCGTCCTGACCTa                                                                                           <  1:5003838/101‑1 (MQ=255)
         aGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCGCGGCGTCCTGACCTa                                                                                           <  1:5380702/101‑1 (MQ=255)
          gTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCGCGGCGTCCTGACCTaa                                                                                          >  1:8352613/1‑101 (MQ=255)
            aCGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCGCGGCGTCCTGACCTaaaa                                                                                        >  1:12996960/1‑101 (MQ=255)
              gACACCATCATCATCGATGAGGCGCACGAACGCAGCCTCAACATCGACTTCCTGCTGGGATATTTAAAGCATATTCTGCCGCAGCGCCCAGATCTGAAGAt                                                                                      >  1:3854819/1‑101 (MQ=255)
               acacTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCGCGGCGTCCTGACCTAAAAATc                                                                                     <  1:6156072/101‑1 (MQ=255)
               acacTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCGCGGCGTCCTGACCTAAAAATc                                                                                     >  1:12007162/1‑101 (MQ=255)
               acacTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCGCGGCGTCCTGACCTAAAAATc                                                                                     >  1:7400935/1‑101 (MQ=255)
                cacTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCGCGGCGTCCTGACCTAAAAATCa                                                                                    <  1:11410267/101‑1 (MQ=255)
                 acTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCGCGGCGTCCTGACCTAAAAATCAt                                                                                   <  1:7732474/101‑1 (MQ=255)
                  gtATCATTAATGACTAAGCGCGCCATGACATAATTAATCTCGGCTTCTTGCCAGGCTATTTTAAAGAGTTTCTTCCGCGGCGTCCTGACCGAAAAATCAtt                                                                                  <  1:6538454/100‑1 (MQ=255)
                   tATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCGCGGCGTCCTGACCTAAAAATCATTa                                                                                 <  1:10398723/101‑1 (MQ=255)
                     tCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCGCGGCGTCCTGACCTAAAAATCATTATc                                                                               <  1:769259/101‑1 (MQ=255)
                      cATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCGCGGCGTCCTGACCTAAAAATCATTATCa                                                                              >  1:4691625/1‑101 (MQ=255)
                         tattGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCGCGGCGTCCTGACCTAAAAATCATTATCACtt                                                                           >  1:6404124/1‑101 (MQ=255)
                          attGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCGCGGCGTCCTGACCTAAAAATCATTATCACTTc                                                                          <  1:6144027/101‑1 (MQ=255)
                          attGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCGCGGCGTACTGACCTAAAAATCATTATCACTTc                                                                          >  1:6320237/1‑101 (MQ=255)
                           ttGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCGCGGCGTCCTGACCTAAAAATCATTATCACTTcc                                                                         <  1:8839584/101‑1 (MQ=255)
                            cgatgaggcgcACGAACGCAGCCTCAACATCGACTTCCTGCTGGGATATTTAAAGCATATTCTGCCGCAGCGCCCAGATCTGAAGATCATTATCACTTcgg                                                                        <  1:4587160/94‑3 (MQ=255)
                              atggggcgcACGAACGCAGCCTCAACATCGACTTCCTGCTGGGATATTTAAAGCATATTCTGCCGCAGCGCCCAGATCTGAAGATCATTATCACTTCggcg                                                                      <  1:5251474/96‑1 (MQ=255)
                              aCGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCGCGGCGTCCTGACCTAAAAATCATTATCACTTCcgcg                                                                      <  1:226990/101‑1 (MQ=255)
                              aCGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCGCGGCGTCCTGACCTAAAAATCATTATCACTTCcgcg                                                                      <  1:9688801/101‑1 (MQ=255)
                              aCGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCGCGGCGTCCTGACCTAAAAATCATTATCACTTCcgcg                                                                      <  1:10166025/101‑1 (MQ=255)
                               cGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCGCGGCGTCCTGACCTAAAAATCATTATCACTTCCGCGa                                                                     >  1:11295689/1‑101 (MQ=255)
                                 aggcgcACGAACGCAGCCTCAACATCGACTTCCTGCTGGGATATTTAAAGCATATTCTGCCGCAGCGCCCAGATCTGAAGATCATTATCACTTCGGCGAcg                                                                   <  1:3379391/99‑2 (MQ=255)
                                 aggcgcACGAACGCAGCCTCAACATCGACTTCCTGCTGGGATATTTAAAGCATATTCTGCCGCAGCGCCCAGATCTGAAGATCATTATCACTTCGGCGAcg                                                                   <  1:5563052/99‑2 (MQ=255)
                                 aggcgcACGAACGCAGCCTCAACATCGACTTCCTGCTGGGATATTTAAAGCATATTCTGCCGCAGCGCCCAGATCTGAAGATCATTATCACTTCGGCGAcg                                                                   <  1:4607799/99‑2 (MQ=255)
                                 aaGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCGCGGCGTCCTGACCTAAAAATCATTATCACTTCCGCGACt                                                                   <  1:1906045/101‑1 (MQ=255)
                                  aGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCGCGGCGTCCTGACCTAAAAATCATTATCACTTCCGCGACTa                                                                  <  1:13925156/101‑1 (MQ=255)
                                  aGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCGCGGCGTCCTGACCTAAAAATCATTATCACTTCCGCGACTa                                                                  <  1:10451273/101‑1 (MQ=255)
                                   gcgcACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCGCGGCGTCCTGACCTAAAAATCATTATCCCTTCCGCGACTAt                                                                 <  1:2502390/101‑1 (MQ=255)
                                   gcgcACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCGCGGCGTCCTGACCTAAAAATCATTATCACTTCCGCGACTAt                                                                 <  1:13689600/101‑1 (MQ=255)
                                   gcgcACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCGCGGCGTCCTGACCTAAAAATCATTATCACTTCCGCGACTAt                                                                 <  1:14284363/101‑1 (MQ=255)
                                   gcgcACGAACGCAGCCTCAACATCGACTTCCTGCTGGGATATTTAAAGCATATTCTGCCGCAGCGCCCAGATCTGAAGATCATTATCACTTCGGCGACGAt                                                                 <  1:3158456/101‑1 (MQ=255)
                                    ggctcGAACGCATTCTGAATATCGATTTGTTGCTCGGGTATTTGAAAGAGTTGCTGCGGCGGCGGCGTGACCTAAAAAGCATTATCACTGCCTCGACTATc                                                                <  1:8322062/97‑1 (MQ=255)
                                    cgcACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCGCGGCGTCCTGACCTAAAAATCATTATCACTTCCGCGACTATc                                                                >  1:5958377/1‑101 (MQ=255)
                                     gcACGAACGCAGCCTCAACATCGACTTCCTGCTGGGATATTTAAAGCATATTCTGCCGCAGCGCCCAGATCTGAAGATCATTATCACTTCGGCGACGAttg                                                               >  1:10767205/1‑99 (MQ=255)
                                      cACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCGCGGCGTCCTGACCTAAAAATCATTATCACTTCCGCGACTATCGa                                                              >  1:13125526/1‑101 (MQ=255)
                                      cACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCGCGGCGTCCTGACCTAAAAATCATTATCACTTCCGCGACTATCGa                                                              >  1:11663715/1‑101 (MQ=255)
                                       aCGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCGCGGCGTCCTGACCTAAAAATCATTATCACTTCCGCGACTATCGAc                                                             >  1:5908903/1‑101 (MQ=255)
                                       aCGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCGCGGCGTCCTGACCTAAAAATCATTATCACTTCCGCGACTATCGAc                                                             >  1:4133374/1‑101 (MQ=255)
                                       aCGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCGCGGCGTCCTGACCTAAAAATCATTATCACTTCCGCGACTATCGAc                                                             >  1:11081177/1‑101 (MQ=255)
                                          aaCGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCGCGGCGTCCTGACCTAAAAATCATTATCACTTCCGCGACTATCGACCCg                                                          >  1:8235378/1‑101 (MQ=255)
                                          aaCGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCGCGGCGTCCTGACCTAAAAATCATTATCACTTCCGCGACTATCGACCCg                                                          <  1:2370246/101‑1 (MQ=255)
                                          aaCGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCGCGGCGTCCTGACCTAAAAATCATTATCACTTCCGCGACTATCGACCCg                                                          <  1:3827884/101‑1 (MQ=255)
                                           aCGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCGCGGCGTCCTGACCTAAAAATCATTATCACTTCCGCGACTATCGACCCgg                                                         <  1:6763396/101‑1 (MQ=255)
                                           aCGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCGCGGCGTCCTGACCTAAAAATCATTATCACTTCCGCGACTATCGACCCgg                                                         <  1:2191070/101‑1 (MQ=255)
                                              cAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCGCGGCGTCCTGACCTAAAAATCATTATCACTTCCGCGACTATCGACCCGGAAc                                                      <  1:1350946/101‑1 (MQ=255)
                                              cAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCGCGGCGTCCTGACCTAAAAATCATTATCACTTCCGCGACTATCGACCCGGAAc                                                      >  1:9435370/1‑101 (MQ=255)
                                               aGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCGCGGCGTCCTGACCTAAAAATCATTATCACTTCCGCGACTATCGACCCGGAACg                                                     >  1:2847210/1‑101 (MQ=255)
                                               aGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCGCGGCGTCCTGACCTAAAAATCATTATCACTTCCGCGACTATCGACCCGGAACg                                                     <  1:7857152/101‑1 (MQ=255)
                                               aGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCGCGGCGTCCTGACCTAAAAATCATTATCACTTCCGCGACTATCGACCCGGAACg                                                     <  1:10351066/101‑1 (MQ=255)
                                                gCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCGCGGCGTCCTGACCTAAAAATCATTATCACTTCCGCGACTATCGACCCGGAACGc                                                    >  1:8697590/1‑101 (MQ=255)
                                                 ccTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCGCGGCGTCCTGACCTAAAAATCATTATCACTTCCGCGACTATCGACCCGGAACGCt                                                   <  1:7553924/101‑1 (MQ=255)
                                                  cTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCGCGGCGTCCTGACCTAAAAATCATTATCACTTCCGCGACTATCGACCCGGAACGCtt                                                  <  1:10720689/101‑1 (MQ=255)
                                                   tGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCGCGGCGTCCTGACCTAAAAATCATTATCCCTTCCGCGACTATCGACCCGGAACGCttt                                                 <  1:1666208/101‑1 (MQ=255)
                                                    gAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCGCGGCGTCCTGACCTAAAAATCATTATCACTTCCGCGACTATCGACCCGGAACGCtttt                                                <  1:7182221/101‑1 (MQ=255)
                                                    gAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCGCGGCGTCCTGACCTAAAAATCATTATCACTTCCGCGACTATCGACCCGGAACGCtttt                                                >  1:13238620/1‑101 (MQ=255)
                                                        atCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCGCGGCGTCCTGACCTAAAAATCATTATCACTTCCGCGACTATCGACCCGGAACGCTTTTcgcg                                            >  1:12318335/1‑101 (MQ=255)
                                                        atCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCGCGGCGTCCTGACCTAAAAATCATTATCACTTCCGCGACTATCGACCCGGAACGCTTTTcgcg                                            >  1:11722854/1‑101 (MQ=255)
                                                         tCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCGCGGCGTCCTGACCTAAAAATCATTATCACTTCCGCGACTATCGACCCGGAACGCTTTTCgcgc                                           >  1:14174249/1‑101 (MQ=255)
                                                          cGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCGCGGCGTCCTGACCTAAAAATCATTATCACTTCCGCGACTATCGACCCGGAACGCTTTTCGCGcc                                          >  1:3065220/1‑101 (MQ=255)
                                                            aTTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCGCGGCGTCCTGACCTAAAAATCATTATCACTTCCGCGACTATCGACCCGGAACGCTTTTCGCGCCAc                                        <  1:259200/101‑1 (MQ=255)
                                                            aTTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCGCGGCGTCCTGACCTAAAAATCATTATCACTTCCGCGACTATCGACCCGGAACGCTTTTCGCGCCAc                                        >  1:1694727/1‑101 (MQ=255)
                                                            aTTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCGCGGCGTCCTGACCTAAAAATCATTATCACTTCCGCGACTATCGACCCGGAACGCTTTTCGCGCCAc                                        <  1:1402093/101‑1 (MQ=255)
                                                             ttttttGCTCGGCTATTTGAAAGAGTTGCTGCCGCGGCGTCCTGACCTAAAAATCATTATCACTTCCGCGACTATCGACCCGGAACGCTTTTCGCGCCACt                                       >  1:10973188/1‑101 (MQ=255)
                                                                tttGCTCGCCTATTTGAAAGAGTGGCGGCCGCGGCGTCCTGACCTAAAAATCATTATCCCTTCCGCGATTA‑CCACCCGGAACGCTTTTCGCGCCACTTtaa                                   <  1:12392007/101‑1 (MQ=255)
                                                                 ttGCTCGGCTATTTGAAAGAGTTGCTGCCGCGGCGTCCTGACCTAAAAATCATTATCACTTCCGCGACTATCGACCCGGAACGCTTTTCGCGCCACTTtaa                                   <  1:9482048/101‑1 (MQ=255)
                                                                  tGCTCGGCTATTTGAAAGAGTTGCTGCCGCGGCGTCCTGACCTAAAAATCATTATCACTTCCGCGACTATCGACCCGGAACGCTTTTCGCGCCACTTtaat                                  <  1:1723812/101‑1 (MQ=255)
                                                                  tGCTCGGCTATTTGAAAGAGTTGCTGCCGCGGCGTCCTGACCTAAAAATCATTATCACTTCCGCGACTATCGACCCGGAACGCTTTTCGCGCCACTTtaat                                  >  1:12461068/1‑101 (MQ=255)
                                                                   gCTCGGCTATTTGAAAGAGTTGCTGCCGCGGCGTCCTGACCTAAAAATCATTATCACTTCCGCGACTATCGACCCGGAACGCTTTTCGCGCCACTTtaata                                 >  1:1961279/1‑101 (MQ=255)
                                                                   gCTCGGCTATTTGAAAGAGTTGCTGCCGCGGCGTCCTGACCTAAAAATCATTATCACTTCCGCGACTATCGACCCGGAACGCTTTTCGCGCCACTTtaata                                 <  1:6937323/101‑1 (MQ=255)
                                                                     tCGGCTATTTGAAAGAGTTGCTGCCGCGGCGTCCTGACCTAAAAATCATTATCACTTCCGCGACTATCGACCCGGAACGCTTTTCGCGCCACTTTaataat                               >  1:9291700/1‑101 (MQ=255)
                                                                      cGGCTATTTGAAAGAGTTGCTGCCGCGGCGTCCTGACCTAAAAATCATTATCACTTCCGCGACTATCGACCCGGAACGCTTTTCGCGCCACTTTAATAATg                              <  1:5880585/101‑1 (MQ=255)
                                                                            tttGAAAGAGTTGCTGCCGCGGCGTCCTGACCTAAAAATCATTATCACTTCCGCGACTATCGACCCGGAACGCTTTTCGCGCCACTTTAATAATGCGCCGa                        >  1:12732486/1‑101 (MQ=255)
                                                                            tttGAAAGAGTTGCTGCCGCGGCGTCCTGACCTAAAAATCATTATCACTTCCGCGACTATCGACCCGGAACGCTTTTCGCGCCACTTTAATAATGCGCCGa                        <  1:8977106/101‑1 (MQ=255)
                                                                               gAAAGAGTTGCTGCCGCGGCGTCCTGACCTAAAAATCATTATCACTTCCGCGACTATCGACCCGGAACGCTTTTCGCGCCACTTTAATAATGCGCCGatta                     <  1:2161207/101‑1 (MQ=255)
                                                                               gAAAGAGTTGCTGCCGCGGCGTCCTGACCTAAAAATCATTATCACTTCCGCGACTATCGACCCGGAACGCTTTTCGCGCCACTTTAATAATGCGCCGatta                     >  1:12803911/1‑101 (MQ=255)
                                                                                  agagTTGCTGCCGCGGCGTCCTGACCTAAAAATCATTATCACTTCCGCGACTATCGACCCGGAACGCTTTTCGCGCCACTTTAATAATGCGCCGATTATTg                  <  1:7947492/101‑1 (MQ=255)
                                                                                  agagTTGCTGCCGCGGCGTCCTGACCTAAAAATCATTATCACTTCCGCGACTATCGACCCGGAACGCTTTTCGCGCCACTTTAATAATGCGCCGATTATTg                  <  1:2084737/101‑1 (MQ=255)
                                                                                    agTTGCTGCCGCGGCGTCCTGACCTAAAAATCATTATCACTTCCGCGACTATCGACCCGGAACGCTTTTCGCGCCACTTTAATAATGCGCCGATTATTGaa                <  1:5213028/101‑1 (MQ=255)
                                                                                    agTTGCTGCCGCGGCGTCCTGACCTAAAAATCATTATCACTTCCGCGACTATCGACCCGGAACGCTTTTCGCGCCACTTTAATAATGCGCCGATTATTGaa                <  1:4642196/101‑1 (MQ=255)
                                                                                      ttGCTGCCGCGGCGTCCTGACCTAAAAATCATTATCACTTCCGCGACTATCGACCCGGAACGCTTTTCGCGCCACTTTAATAATGCGCCGATTATTGAAGt              <  1:11053721/101‑1 (MQ=255)
                                                                                       tgctgcCGCGGCGTCCTGACCTAAAAATCATTATCACTTCCGCGACTATCGACCCGGAACGCTTTTCGCGCCACTTTAATAATGCGCCGATTATTGAAGtc             <  1:8156923/101‑1 (MQ=255)
                                                                                           gcCGCGGCGTCCTGACCTAAAAATCATTAACACTTCCGCGACTATCGACCCGGAACGCTTTTCGCGCCACTTTAATAATGCGCCGATTATTGAAGTCTCCg         >  1:4278567/1‑101 (MQ=255)
                                                                                             cgcgGCGTCCTGACCTAAAAATCATTATCACTTCCGCGACTATCGACCCGGAACGCTTTTCGCGCCACTTTAATAATGCGCCGATTATTGAAGTCTCCGGt       >  1:9725706/1‑101 (MQ=255)
                                                                                              gcggcgTCCTGACCTAAAAATCATTATCACTTCCGCGACTATCGACCCGGAACGCTTTTCGCGCCACTTTAATAATGCGCCGATTATTGAAGTCTCCGGTc      <  1:935672/101‑1 (MQ=255)
                                                                                              gcggcgTCCTGACCTAAAAATCATTATCACTTCCGCGACTATCGACCCGGAACGCTTTTCGCGCCACTTTAATAATGCGCCGATTATTGAAGTCTCCGGTc      <  1:2143524/101‑1 (MQ=255)
                                                                                               cggcgTCCTGACCTAAAAATCATTATCACTTCCGCGACTATCGACCCGGAACGCTTTTCGCGCCACTTTAATAATGCGCCGATTATTGAAGTCTCCGGTCg     >  1:706044/1‑101 (MQ=255)
                                                                                               cggcgTCCTGACCTAAAAATCATTATCACTTCCGCGACTATCGACCCGGAACGCTTTTCGCGCCACTTTAATAATGCGCCGATTATTGAAGTCTCCGGTCg     >  1:1482983/1‑101 (MQ=255)
                                                                                               cggcgTCCTGACCTAAAAATCATTATCACTTCCGCGACTATCGACCCGGAACGCTTTTCGCGCCACTTTAATAATGCGCCGATTATTGAAGTCTCCGGTCg     <  1:100778/101‑1 (MQ=255)
                                                                                                ggcgTCCTGACCTAAAAATCATTATCACTTCCGCGACTATCGACCCGGAACGCTTTTCGCGCCACTTTAATAATGCGCCGATTATTGAAGTCTCCGGTCgg    <  1:4917160/101‑1 (MQ=255)
                                                                                                ggcgTCCTGACCTAAAAATCATTATCACTTCCGCGACTATCGACCCGGAACGCTTTTCGCGCCACTTTAATAATGCGCCGATTATTGAAGTCTCCGGTCgg    >  1:3361400/1‑101 (MQ=255)
                                                                                                    |                                                                                                  
TGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCGCGGCGTCCTGACCTAAAAATCATTATCACTTCCGCGACTATCGACCCGGAACGCTTTTCGCGCCACTTTAATAATGCGCCGATTATTGAAGTCTCCGGTCGGAC  >  NC_007779/1485334‑1485532

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