Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 3,326,085 T→A 12.0% M236L (ATG→TTG)  ftsH ← protease, ATP‑dependent zinc‑metallo

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0077793,326,0850T→A12.0% 102.1 0.0 ‑0.4266M236L (ATG→TTG) ftsHprotease, ATP‑dependent zinc‑metallo
Reads supporting (aligned to +/- strand):  new base (16/16):  ref base (115/119):  total (131/135)
Fisher's exact test strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 4.07e-01

Only 100 of 272 total aligned reads displayed.

TCGTGACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACTTTCGCTTCGCCTGCAATCGCTTTC  >  NC_007779/3325985‑3326185
                                                                                                     |                                                                                                    
 cGTGACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATg                                                                                                     <  1:4054564/101‑1 (MQ=255)
 cGTGACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATg                                                                                                     >  1:1646629/1‑101 (MQ=255)
     accaccGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCAc                                                                                                 <  1:8978588/101‑1 (MQ=255)
        accGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGTACATGTCACGaa                                                                                              >  1:8823640/1‑101 (MQ=255)
        accGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGaa                                                                                              <  1:2229197/101‑1 (MQ=255)
         ccGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAac                                                                                             <  1:7119083/101‑1 (MQ=255)
          cGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAaca                                                                                            <  1:5183028/101‑1 (MQ=255)
            cccAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACg                                                                                          <  1:3865080/101‑1 (MQ=255)
              cAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACggg                                                                                        >  1:2020171/1‑101 (MQ=255)
               aGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGa                                                                                       <  1:4369900/101‑1 (MQ=255)
               aGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGa                                                                                       <  1:942836/101‑1 (MQ=255)
                gACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGAt                                                                                      >  1:14857561/1‑101 (MQ=255)
                  ccAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGGTGc                                                                                    <  1:8982299/101‑1 (MQ=255)
                       gcCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGTATGCACCca                                                                               <  1:10121043/101‑1 (MQ=255)
                        ccACGCTGGCGGCCTACGGCGTCGATTTCATCGTTAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCcac                                                                              <  1:8260260/101‑1 (MQ=255)
                        ccACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCcac                                                                              <  1:1943701/101‑1 (MQ=255)
                        ccACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCcac                                                                              <  1:8511224/101‑1 (MQ=255)
                          aCGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCacac                                                                            >  1:9374179/1‑101 (MQ=255)
                             cTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACcga                                                                         <  1:5714375/101‑1 (MQ=255)
                             cTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACcga                                                                         >  1:283375/1‑101 (MQ=255)
                              tGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACcgac                                                                        <  1:9087782/101‑1 (MQ=255)
                                gcggcCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACcgacga                                                                      <  1:1562208/101‑1 (MQ=255)
                                gcggcCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACcgacga                                                                      >  1:14992087/1‑101 (MQ=255)
                                 cggcCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGaa                                                                     >  1:7827197/1‑101 (MQ=255)
                                  ggcCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAAc                                                                    >  1:2706662/1‑101 (MQ=255)
                                  ggcCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAAc                                                                    <  1:11383143/101‑1 (MQ=255)
                                   gcCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACa                                                                   >  1:3786707/1‑101 (MQ=255)
                                    ccTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACAt                                                                  <  1:683027/101‑1 (MQ=255)
                                    ccTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACAt                                                                  >  1:3678872/1‑101 (MQ=255)
                                     cTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACAtt                                                                 >  1:10459382/1‑101 (MQ=255)
                                      tACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACAttt                                                                <  1:5574133/101‑1 (MQ=255)
                                       aCGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTc                                                               >  1:9176575/1‑101 (MQ=255)
                                         ggCGTCAATCTCGTCGATGAAGATAATGCACGGGCTGTTTTCCTTGGCCTGATTAAACAAGTCACGGACACGAGAGGCACCGACACCGACGAACATTTcga                                                             >  1:5819017/1‑99 (MQ=255)
                                          gCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTAc                                                            >  1:5222538/1‑101 (MQ=255)
                                          gCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTAc                                                            <  1:2455865/101‑1 (MQ=255)
                                           cGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACg                                                           >  1:3299817/1‑101 (MQ=255)
                                           cGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACg                                                           >  1:2770722/1‑101 (MQ=255)
                                           cGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACg                                                           >  1:14366126/1‑101 (MQ=255)
                                           cGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACg                                                           >  1:4425780/1‑101 (MQ=255)
                                            gTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGa                                                          >  1:13332871/1‑101 (MQ=255)
                                            gTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGa                                                          <  1:14602175/101‑1 (MQ=255)
                                             tCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGaa                                                         >  1:12629158/1‑101 (MQ=255)
                                             tCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGaa                                                         <  1:6046681/101‑1 (MQ=255)
                                             tCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGaa                                                         >  1:4497830/1‑101 (MQ=255)
                                              caatctcgtCGATGAAGATAATGCACGGGCTGTTTTCCTTGGCCTGATTAAACAAGTCACGGACACGAGAGGCACCGACACCGACGAACATTTCGACGAAg                                                        >  1:10254932/9‑101 (MQ=255)
                                              cGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAg                                                        >  1:1155676/1‑101 (MQ=255)
                                              cGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAg                                                        <  1:6902429/101‑1 (MQ=255)
                                               gATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGt                                                       <  1:66380/101‑1 (MQ=255)
                                                aTTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTc                                                      >  1:11231391/1‑101 (MQ=255)
                                                 tttCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCa                                                     >  1:4640137/1‑101 (MQ=255)
                                                 tctcgtCGATGAAGATAATGCACGGGCTGTTTTCCTTGGCCTGATTAAACAAGTCACGGACACGAGAGGCACCGACACCGACGAACATTTCGACGAAGTcc                                                     >  1:8007029/6‑100 (MQ=255)
                                                  ttCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAg                                                    <  1:14529251/101‑1 (MQ=255)
                                                   tCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGa                                                   <  1:11960045/101‑1 (MQ=255)
                                                    cATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGaa                                                  <  1:2652366/101‑1 (MQ=255)
                                                    cATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGaa                                                  <  1:4687659/101‑1 (MQ=255)
                                                     aTCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAAc                                                 >  1:3335478/1‑101 (MQ=255)
                                                      tCGATGAAGATAATGCACGGGCTGTTTTCCTTGGCCTGATTAAACAAGTCACGGACACGAGAGGCACCGACACCGACGAACATTTCGACGAAGTCCGAAcc                                                >  1:953981/1‑101 (MQ=255)
                                                      tCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAAcc                                                <  1:13756974/101‑1 (MQ=255)
                                                      tCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAAcc                                                >  1:7270895/1‑101 (MQ=255)
                                                      tCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAAcc                                                <  1:15160539/101‑1 (MQ=255)
                                                       cGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCg                                               <  1:11132677/101‑1 (MQ=255)
                                                       cGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCg                                               >  1:3899002/1‑101 (MQ=255)
                                                       cGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCg                                               >  1:13079889/1‑101 (MQ=255)
                                                         atgaaGATAATGCACGGGCTGTTTTCCTTGGCCTGATTAAACAAGTCACGGACACGAGAGGCACCGACACCGACGAACATTTCGACGAAGTCCGAACCCga                                             <  1:14074427/98‑1 (MQ=255)
                                                         atgaaGATAATGCACGGGCTGTTTTCCTTGGCCTGATTAAACAAGTCACGGACACGAGAGGCACCGACACCGACGAACATTTCGACGAAGTCCGAACCCga                                             <  1:7963846/98‑1 (MQ=255)
                                                         aTAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGga                                             >  1:4186581/1‑101 (MQ=255)
                                                         aTAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGga                                             >  1:11066779/1‑101 (MQ=255)
                                                           aaaGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGgaga                                           >  1:1056732/1‑101 (MQ=255)
                                                            aaGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGAt                                          >  1:4422714/1‑101 (MQ=255)
                                                            aaGATAATGCACGGGCTGTTTTCCTTGGCCTGATTAAACAAGTCACGGACACGAGAGGCACCGACACCGACGAACATTTCGACGAAGTCCGAACCCGAGAt                                          <  1:6311117/101‑1 (MQ=255)
                                                               atgatgCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGt                                       <  1:12777580/101‑1 (MQ=255)
                                                               atgatgCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGt                                       <  1:6389673/101‑1 (MQ=255)
                                                                taatgCACGGGCTGTTTTCCTTGGCCTGATTAAACAAGTCACGGACACGAGAGGCACCGACACCGACGAACATTTCGACGAAGTCCGAACCCGAGAtggcg                                      <  1:7500325/99‑5 (MQ=255)
                                                                 gatgCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTaa                                     <  1:7685535/101‑1 (MQ=255)
                                                                 gatgCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTaa                                     >  1:4223769/1‑101 (MQ=255)
                                                                 aatgCACGGGCTGTTTTCCTTGGCCTGATTAAACAAGTCACGGACACGAGAGGCACCGACACCGACGAACATTTCGACGAAGTCCGAACCCGAGAtggcgt                                     <  1:3253299/100‑6 (MQ=255)
                                                                  atgCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTaaa                                    >  1:563598/1‑101 (MQ=255)
                                                                  atgCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTaaa                                    >  1:3456193/1‑101 (MQ=255)
                                                                  atgCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTaaa                                    >  1:3472111/1‑101 (MQ=255)
                                                                       cGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACg                               <  1:13533652/101‑1 (MQ=255)
                                                                          tGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGaa                            <  1:4734786/101‑1 (MQ=255)
                                                                           gCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAAc                           <  1:11711676/101‑1 (MQ=255)
                                                                           gCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAAc                           <  1:12865856/101‑1 (MQ=255)
                                                                           gCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAAc                           <  1:2335757/101‑1 (MQ=255)
                                                                             cGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACtt                         >  1:12180611/1‑101 (MQ=255)
                                                                              gCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACttt                        >  1:8138991/1‑101 (MQ=255)
                                                                              gCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACttt                        >  1:8330256/1‑101 (MQ=255)
                                                                              gCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACttt                        <  1:1311575/101‑1 (MQ=255)
                                                                                tttCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACTttcg                      <  1:2515949/101‑1 (MQ=255)
                                                                                 ttcttcGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACTttcgc                     >  1:9719843/1‑101 (MQ=255)
                                                                                  tcttcGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACTttcgct                    >  1:10039100/1‑101 (MQ=255)
                                                                                    cttcGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAAGGGAACTttcgctt                   <  1:14660242/101‑1 (MQ=255)
                                                                                         ccTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACTTTCGCTTCGCCt              >  1:14154309/1‑101 (MQ=255)
                                                                                          cTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACTTTCGCTTCGCCTg             <  1:5289136/101‑1 (MQ=255)
                                                                                            gTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACTTTCGCTTCGCCTGCa           <  1:3371607/101‑1 (MQ=255)
                                                                                               cGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACTTTCGCTTCGCCTGCAATc        >  1:10892464/1‑101 (MQ=255)
                                                                                                gAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACTTTCGCTTCGCCTGCAATCg       <  1:2562870/101‑1 (MQ=255)
                                                                                                  aCATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACTTTCGCTTCGCCTGCAATCGCt     <  1:5044881/101‑1 (MQ=255)
                                                                                                   cATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACTTTCGCTTCGCCTGCAATCGCtt    <  1:12437578/101‑1 (MQ=255)
                                                                                                   cATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACTTTCGCTTCGCCTGCAATCGCtt    <  1:10302020/101‑1 (MQ=255)
                                                                                                     |                                                                                                    
TCGTGACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACTTTCGCTTCGCCTGCAATCGCTTTC  >  NC_007779/3325985‑3326185

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.