Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 3,326,107 C→G 8.8% V228V (GTG→GTC ftsH ← protease, ATP‑dependent zinc‑metallo

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0077793,326,1070C→G 8.8% 59.2 ‑1.3 ‑0.3260V228V (GTG→GTCftsHprotease, ATP‑dependent zinc‑metallo
Reads supporting (aligned to +/- strand):  new base (9/14):  ref base (126/111):  total (135/125)
Fisher's exact test strand distribution p-value = 2.74e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 5.10e-01

Only 100 of 275 total aligned reads displayed.

CGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACTTTCGCTTCGCCTGCAATCGCTTTCGCCAGCAGCGTTTTACCGGTAC  >  NC_007779/3326007‑3326207
                                                                                                     |                                                                                                    
cgcCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCAccc                                                                                                      <  1:7312993/101‑1 (MQ=255)
 gcCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCca                                                                                                     >  1:12411969/1‑101 (MQ=255)
 gcCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCca                                                                                                     >  1:10375653/1‑101 (MQ=255)
  ccACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCcac                                                                                                    <  1:10032269/101‑1 (MQ=255)
   cACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCcaca                                                                                                   >  1:11673391/1‑101 (MQ=255)
     cGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACAcc                                                                                                 >  1:12809736/1‑101 (MQ=255)
     cGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACAcc                                                                                                 <  1:11626583/101‑1 (MQ=255)
     cGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACAcc                                                                                                 <  1:8004793/101‑1 (MQ=255)
      gCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACcg                                                                                                <  1:13485595/101‑1 (MQ=255)
       cTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACcga                                                                                               <  1:4273902/101‑1 (MQ=255)
       cTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACcga                                                                                               >  1:13328187/1‑101 (MQ=255)
       cTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACcga                                                                                               <  1:8437146/101‑1 (MQ=255)
       cTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACcga                                                                                               <  1:11057924/101‑1 (MQ=255)
       cTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACcga                                                                                               <  1:11849932/101‑1 (MQ=255)
        tGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACcgac                                                                                              >  1:1425548/1‑101 (MQ=255)
        tGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACcgac                                                                                              <  1:12610016/101‑1 (MQ=255)
            gggcTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAAc                                                                                          <  1:6751126/98‑1 (MQ=255)
            ggcCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAAc                                                                                          >  1:4977544/1‑101 (MQ=255)
                 aCGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTc                                                                                     >  1:13051216/1‑101 (MQ=255)
                  cGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCt                                                                                    <  1:6629305/101‑1 (MQ=255)
                  cGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCt                                                                                    >  1:6460987/1‑101 (MQ=255)
                  cGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCt                                                                                    <  1:4395246/101‑1 (MQ=255)
                  cGGCGTCAATCTCGTCGATGAAGATAATGCACGGGCTGTTTTCCTTGGCCTGATTAAACAAGTCACGGACACGAGAGGCTCCGACACCGACGAACATTTcg                                                                                    >  1:9845614/1‑100 (MQ=255)
                    gCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTATTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTAc                                                                                  <  1:1128859/101‑1 (MQ=255)
                     cGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACg                                                                                 <  1:11593161/101‑1 (MQ=255)
                     cGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTACGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACg                                                                                 >  1:10290471/1‑101 (MQ=255)
                      gTCAATCTCGTCGATGAAGATAATGCGCGGGCTGTTTTCCTTGGCCTGATTAAACAAGTCACGGACACGAGAGGCACCGACACCGACGAACATTTCGACGa                                                                                >  1:12015129/1‑101 (MQ=255)
                      gTCAATCTCGTCGATGAAGATAATGCACGGGCTGTTTTCCTTGGCCTGATTAAACAAGTCACGGACACGAGAGGCACCGACACCGACGAACATTTCGACGa                                                                                >  1:1503468/1‑101 (MQ=255)
                      gTCAATCTCGTCGATGAAGATAATGCACGGGCTGTTTTCCTTGGCCTGATTAAACAAGTCACGGACACGAGAGGCACCGACACCGACGAACATTTCGACGa                                                                                <  1:13395897/101‑1 (MQ=255)
                       tcaatctcgtCGATGAAGATAATGCACGGGCTGTTTTCCTTGGCCTGATTAAACAAGTCACGGACACGAGAGGCACCGACACCGACGAACATTTCGACGaa                                                                               <  1:12139092/92‑1 (MQ=255)
                       tcaatctcgtCGATGAAGATAATGCACGGGCTGTTTTCCTTGGCCTGATTAAACAAGTCACGGACACGAGAGGCACCGACACCGACGAACATTTCGACGaa                                                                               <  1:195414/92‑1 (MQ=255)
                       tCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGaa                                                                               <  1:4127755/101‑1 (MQ=255)
                       tCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGaa                                                                               <  1:3242331/101‑1 (MQ=255)
                        cGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAg                                                                              >  1:5261449/1‑101 (MQ=255)
                          atctcgtCGATGAAGATAATGCACGGGCTGTTTTCCTTGGCCTGATTAAACAAGTCACGGACACGAGAGGCACCGACACCGACGAACATTTCGACGAAGTc                                                                            >  1:2779320/7‑101 (MQ=255)
                          aTTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTc                                                                            >  1:9507330/1‑101 (MQ=255)
                          aTTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTc                                                                            >  1:6841672/1‑101 (MQ=255)
                           tctcgtCGATGAAGATAATGCACGGGCTGTTTTCCTTGGCCTGATTAAACAAGTCACGGACACGAGAGGCACCGACACCGACGAACATTTCGACGAAGTcc                                                                           <  1:13192193/96‑2 (MQ=255)
                                tCGATGAAGATAATGCACGGGCTGTTTTCCTTGGCCTGATTAAACAAGTCACGGACACGAGAGGCACCGACACCGACGAACATTTCGACGAAGTCCGAAcc                                                                      <  1:1838195/101‑1 (MQ=255)
                                  gATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCgg                                                                    >  1:5328520/1‑101 (MQ=255)
                                   aTAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGga                                                                   >  1:10063475/1‑101 (MQ=255)
                                    tgaaGATAATGCACGGGCTGTTTTCCTTGGCCTGATTAAACAAGTCACGGACACGAGAGGCACCGACACCGACGAACATTTCGACGAAGTCCGAACCCgag                                                                  <  1:2022170/99‑1 (MQ=255)
                                    tAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGgag                                                                  <  1:775500/101‑1 (MQ=255)
                                      aaGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGAt                                                                >  1:1399160/1‑101 (MQ=255)
                                       aGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATa                                                               >  1:6143519/1‑101 (MQ=255)
                                        gatgatGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGTACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAg                                                              <  1:12630227/101‑1 (MQ=255)
                                        gataatGCACGGGCTGTTTTCCTTGGCCTGATTAAACAAGTCACGGACACGAGAGGCACCGACACCGACGAACATTTCGACGAAGTCCGAACCCGAGAtgg                                                              <  1:5924906/101‑3 (MQ=255)
                                         atgatgCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGt                                                             <  1:12753871/101‑1 (MQ=255)
                                          tgatgCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTa                                                            >  1:11557408/1‑101 (MQ=255)
                                          tgatgCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTa                                                            >  1:4553777/1‑101 (MQ=255)
                                          tgatgCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTa                                                            >  1:2571776/1‑101 (MQ=255)
                                            atgCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTaaa                                                          <  1:79236/101‑1 (MQ=255)
                                            atgCACGGGCTGTTTTCCTTGGCCTGATTAAACAAGTCACGGACACGAGAGGCACCGACACCGACGAACATTTCGACGAAGTCCGAACCCGAGAtggcgta                                                          <  1:2962477/101‑7 (MQ=255)
                                               cACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGaa                                                       <  1:5511404/101‑1 (MQ=255)
                                               cACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGaa                                                       >  1:6969373/1‑101 (MQ=255)
                                                aCGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAAc                                                      >  1:6541403/1‑101 (MQ=255)
                                                  ggTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACgg                                                    >  1:702521/1‑101 (MQ=255)
                                                    tGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGaa                                                  >  1:3559758/1‑101 (MQ=255)
                                                     gCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAAc                                                 >  1:8583881/1‑101 (MQ=255)
                                                     gCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAAc                                                 <  1:2883250/101‑1 (MQ=255)
                                                     gCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAAc                                                 >  1:5073447/1‑101 (MQ=255)
                                                       cGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACtt                                               <  1:13648855/101‑1 (MQ=255)
                                                        gCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACttt                                              <  1:1136136/101‑1 (MQ=255)
                                                         cTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACTttc                                             >  1:5389038/1‑101 (MQ=255)
                                                         cTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACTttc                                             <  1:2330502/101‑1 (MQ=255)
                                                           ttcttcGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACTttccc                                           >  1:10893380/1‑99 (MQ=255)
                                                            tcttcGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACTttcgct                                          >  1:1204187/1‑101 (MQ=255)
                                                              cttcGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACTttcgctt                                         <  1:1055175/101‑1 (MQ=255)
                                                               ttcGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACTttcgcttc                                        <  1:12914823/101‑1 (MQ=255)
                                                                tcGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACTttcgcttcg                                       >  1:5768112/1‑101 (MQ=255)
                                                                 cGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACTttcgcttcgc                                      <  1:9927672/101‑1 (MQ=255)
                                                                   ccTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACTTTCGCTTCGCCt                                    >  1:8941062/1‑101 (MQ=255)
                                                                   ccTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACTTTCGCTTCGCCt                                    >  1:8226/1‑101 (MQ=255)
                                                                     tGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACTTTCGCTTCGCCTGc                                  <  1:6100234/101‑1 (MQ=255)
                                                                     tGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACTTTCGCTTCGCCTGc                                  <  1:4018131/101‑1 (MQ=255)
                                                                      gTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACTTTCGCTTCGCCTGCa                                 <  1:10822094/101‑1 (MQ=255)
                                                                       ttCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACTTTCGCTTCGCCTGCaa                                <  1:2239979/101‑1 (MQ=255)
                                                                        tCGAACATGTCACGAACCCGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACTTTCGCTTCGCCTGCAAt                               <  1:10105467/101‑1 (MQ=255)
                                                                           aaCATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACTTTCGCTTCGCCTGCAATCGc                            <  1:12513783/101‑1 (MQ=255)
                                                                             cATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACTTTCGCTTCGCCTGCAATCGCtt                          <  1:5733666/101‑1 (MQ=255)
                                                                               tGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGAAGATAGTAAAGAACGGAACTTTCGCTTCGCCTGCAATCGCTTTc                        <  1:2473342/101‑1 (MQ=255)
                                                                                 tCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACTTTCGCTTCGCCTGCAATCGCTTTCGc                      <  1:808107/101‑1 (MQ=255)
                                                                                  cACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACTTTCGCTTCGCCTGCAATCGCTTTCGcc                     <  1:1070042/101‑1 (MQ=255)
                                                                                  cACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACTTTCGCTTCGCCTGCAATCGCTTTCGcc                     >  1:9473653/1‑101 (MQ=255)
                                                                                    cGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACTTTCGCTTCGCCTGCAATCGCTTTCGCcag                   >  1:3849943/1‑101 (MQ=255)
                                                                                    cGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACTTTCGCTTCGCCTGCAATCGCTTTCGCcag                   >  1:1676148/1‑101 (MQ=255)
                                                                                     gAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACTTTCGCTTCGCCTGCAATCGCTTTCGCcagc                  <  1:6190468/101‑1 (MQ=255)
                                                                                       acacGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACTTTCGCTTCGCCTGCAATCGCTTTCGCcagcag                <  1:7417464/101‑1 (MQ=255)
                                                                                        cacGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACTTTCGCTTCGCCTGCAATCGCTTTCGCCagcagc               >  1:8736343/1‑101 (MQ=255)
                                                                                             gATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACTTTCGCTTCGCCTGCAATCGCTTTCGCCAGCAGCGtttt          >  1:7729572/1‑101 (MQ=255)
                                                                                             gATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACTTTCGCTTCGCCTGCAATCGCTTTCGCCAGCAGCGtttt          <  1:7677566/101‑1 (MQ=255)
                                                                                              aTGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACTTTCGCTTCGCCTGCAATCGCTTTCGCCAGCAGCGTTTTa         <  1:1076287/101‑1 (MQ=255)
                                                                                              aTGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACTTTCGCTTCGCCTGCAATCGCTTTCGCCAGCAGCGTTTTa         >  1:13841917/1‑101 (MQ=255)
                                                                                              aTGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACTTTCGCTTCGCCTGCAATCGCTTTCGCCAGCAGCGTTTTa         >  1:10400496/1‑101 (MQ=255)
                                                                                              aTGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACTTTCGCTTCGCCTGCAATCGCTTTCGCCAGCAGCGTTTTa         <  1:6235975/101‑1 (MQ=255)
                                                                                               tGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACTTTCGCTTCGCCTGCAATCGCTTTCGCCAGCAGCGTTTTAc        <  1:4764463/101‑1 (MQ=255)
                                                                                               tGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACTTTCGCTTCGCCTGCAATCGCTTTCGCCAGCAGCGTTTTAc        >  1:3035856/1‑101 (MQ=255)
                                                                                                gCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACTTTCGCTTCGCCTGCAATCGCTTTCGCCAGCAGCGTTTTAcc       <  1:9137052/101‑1 (MQ=255)
                                                                                                  aCCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACTTTCGCTTCGCCTGCAATCGCTTTCGCCAGCAGCGTTTTACCgg     >  1:2752575/1‑101 (MQ=255)
                                                                                                   cccACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACTTTCGCTTCGCCTGCAATCGCTTTCGCCAGCAGCGTTTTACCGGt    <  1:6309966/101‑1 (MQ=255)
                                                                                                     |                                                                                                    
CGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACTTTCGCTTCGCCTGCAATCGCTTTCGCCAGCAGCGTTTTACCGGTAC  >  NC_007779/3326007‑3326207

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.