Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 3,442,207 A→G 3.0% K542K (AAA→AAG thiC → thiamin (pyrimidine moiety) biosynthesis protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0077793,442,2070A→G 3.0% 9.8 0.0 ‑0.1230K542K (AAA→AAGthiCthiamin (pyrimidine moiety) biosynthesis protein
Reads supporting (aligned to +/- strand):  new base (4/3):  ref base (121/102):  total (125/105)
Fisher's exact test strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 8.20e-01

Only 100 of 248 total aligned reads displayed.

ATGTTAAGCAGGGGCTTATCACCTATAAGATTGCTGCCCACGCCGCTGACCTGGCGAAAGGGCATCCGGGCGCGCAAATTCGCGATAACGCCATGTCGAAAGCCCGCTTCGAATTTCGCTGGGAAGACCAGTTTAATCTGGCCCTCGACCCGTTTACCGCCCGCGCTTATCACGATGAAACCCTGCCGCAAGAGTCAGGTA  >  NC_007779/3442107‑3442307
                                                                                                    |                                                                                                    
aTGTTAAGCAGGGGCTTATCACCTATAAGATTGCTGCCCACGCCGCTGACCTGGCGAAAGGGCATCCGGGCGCGCAAATTCGCGATAACGCCATGTCGaaa                                                                                                      >  1:13932359/1‑101 (MQ=255)
  gTTAAGCAGGGGCTTATCACCTATAAGATTGCTGCCCACGCCGCTGACCTGGCGAAAGGGCATCCGGGCGCGCAAATTCGCGATAACGCCATGTCGAAAGc                                                                                                    >  1:12694015/1‑101 (MQ=255)
  gTTAAGCAGGGGCTTATCACCTATAAGATTGCTGCCCACGCCGCTGACCTGGCGAAAGGGCATCCGGGCGCGCAAATTCGCGATAACGCCATGTCGAAAGc                                                                                                    <  1:5152132/101‑1 (MQ=255)
   ttAAGCAGGGGCTTATCACCTATAAGATTGCTGCCCACGCCGCTGACCTGGCGAAAGGGCATCCGGGCGCGCAAATTCGCGATAACGCCATGTCGAAAGcc                                                                                                   >  1:624766/1‑101 (MQ=255)
     aaGCAGGGGCTTATCACCTATAAGATTGCTGCCCACGCCGCTGACCTGGCGAAAGGGCATCCGGGCGCGCAAATTCGCGATAACGCCATGTCGAAAGCCCg                                                                                                 >  1:3425706/1‑101 (MQ=255)
     aaGCAGGGGCTTATCACCTATAAGATTGCTGCCCACGCCGCTGACCTGGCGAAAGGGCATCCGGGCGCGCAAATTCGCGATAACGCCATGTCGAAAGCCCg                                                                                                 >  1:669250/1‑101 (MQ=255)
     aaGCAGGGGCTTATCACCTATAAGATTGCTGCCCACGCCGCTGACCTGGCGAAAGGGCATCCGGGCGCGCAAATTCGCGATAACGCCATGTCGAAAGCCCg                                                                                                 <  1:5329941/101‑1 (MQ=255)
     aaGCAGGGGCTTATCACCTATAAGATTGCTGCCCACGCCGCTGACCTGGCGAAAGGGCATCCGCGCGCGCAAATTCGCGATAACGCCATGTCGAAAGCCCg                                                                                                 >  1:385663/1‑101 (MQ=255)
         aGGGGCTTATCACCTATAAGATTGCTGCCCACGCCGCTGACCTGGCGAAAGGGCATCCGGGCGCGCAAATTCGCGATAACGCCATGTCGAAAGCCCGCTTc                                                                                             <  1:10327064/101‑1 (MQ=255)
            ggCTTATCACCTATAAGATTGCTGCCCACGCCGCTGACCTGGCGAAAGGGCATCCGGGCGCGCAAATTCGCGATAACGCCATGTCGAAAGCCCGCTTCGaa                                                                                          <  1:8387406/101‑1 (MQ=255)
            ggCTTATCACCTATAAGATTGCTGCCCACGCCGCTGACCTGGCGAAAGGGCATCCGGGCGCGCAAATTCGCGATAACGCCATGTCGAAAGCCCGCTTCGaa                                                                                          >  1:10003153/1‑101 (MQ=255)
            ggCTTATCACCTATAAGATTGCTGCCCACGCCGCTGACCTGGCGAAAGGGCATCCGGGCGCGCAAATTCGCGATAACGCCATGTCGAAAGACCGCTTCGaa                                                                                          >  1:660265/1‑101 (MQ=255)
             gCTTATCACCTATAAGATTGCTGCCCACGCCGCTGACCTGGCGAAAGGGCATCCGGGCGCGCAAATTCGCGATAACGCCATGTCGAAAGCCCGCTTCGAAt                                                                                         >  1:4868446/1‑101 (MQ=255)
             gCTTATCACCTATAAGATTGCTGCCCACGCCGCTGACCTGGCGAAAGGGCATCCGGGCGCGCAAATTCGCGATAACGCCATGTCGAAAGCCCGCTTCGAAt                                                                                         <  1:12399882/101‑1 (MQ=255)
              cTTATCACCTATAAGATTGCTGCCCACGCCGCTGACCTGGCGAAAGGGCATCCGGGCGCGCAAATTCGCGATAACGCCATGTCGAAAGCCCGCTTCGAAtt                                                                                        <  1:9323607/101‑1 (MQ=255)
               ttATCACCTATAAGATTGCTGCCCACGCCGCTGACCTGGCGAAAGGGCATCCGGGCGCGCAAATTCGCGATAACGCCATGTCGAAAGCCCGCTTCGAAttt                                                                                       <  1:5664813/101‑1 (MQ=255)
               ttATCACCTATAAGATTGCTGCCCACGCCGCTGACCTGGCGAAAGGGCATCCGGGCGCGCAAATTCGCGATAACGCCATGTCGAAAGCCCGCTTCGAAttt                                                                                       >  1:1549036/1‑101 (MQ=255)
               ttATCACCTATAAGATTGCTGCCCACGCCGCTGACCTGGCGAAAGGGCATCCGGGCGCGCAAATTCGCGATAACGCCATGTCGAAAGCCCGCTTCGAAttt                                                                                       >  1:6286806/1‑101 (MQ=255)
               ttATCACCTATAAGATTGCTGCCCACGCCGCTGACCTGGCGAAAGGGCATCCGGGCGCGCAAATTCGCGATAACGCCATGTCGAAAGCCCGCTTCGAAttt                                                                                       <  1:417594/101‑1 (MQ=255)
                tATCACCTATAAGATTGCTGCCCACTCCGCTGACCTGGCGAAAGGGCATCCGGGCGCGCAAATTCGCGATAACGCCATGTCGAAAGCCCGCTTCGAATTTc                                                                                      >  1:10309996/1‑101 (MQ=255)
                tATCACCTATAAGATTGCTGCCCACGCCGCTGACCTGGCGAAAGGGCATCCGGGCGCGCAAATTCGCGATAACGCCATGTCGAAAGCCCGCTTCGAATTTc                                                                                      <  1:125499/101‑1 (MQ=255)
                tATCACCTATAAGATTGCTGCCCACGCCGCTGACCTGGCGAAAGGGCATCCGGGCGCGCAAATTCGCGATAACGCCATGTCGAAAGCCCGCTTCGAATTTc                                                                                      <  1:375982/101‑1 (MQ=255)
                tATCACCTATAAGATTGCTGCCCACGCCGCTGACCTGGCGAAAGGGCATCCGGGCGCGCAAATTCGCGATAACGCCATGTCGAAAGCCCGCTTCGAATTTc                                                                                      <  1:11482295/101‑1 (MQ=255)
                tATCACCTATAAGATTGCTGCCCACGCCGCTGACCTGGCGAAAGGGCATCCGGGCGCGCAAATTCGCGATAACGCCATGTCGAAAGCCCGCTTCGAATTTc                                                                                      <  1:4460107/101‑1 (MQ=255)
                tATCACCTATAAGATTGCTGCCCACGCCGCTGACCTGGCGAAAGGGCATCCGGGCGCGCAAATTCGCGATAACGCCATGTCGAAAGCCCGCTTCGAATTTc                                                                                      <  1:10642962/101‑1 (MQ=255)
                 aTCACCTATAAGATTGCTGCCCACGCCGCTGACCTGGCGAAAGGGCATCCGGGCGCGCAAATTCGCGATAACGCCATGTCGAAAGCCCGCTTCGAATTTCg                                                                                     >  1:1799543/1‑101 (MQ=255)
                   cACCTATAAGATTGCTGCCCACGCCGCTGACCTGGCGAAAGGGCATCCGGGCGCGCAAATTCGCGATAACGCCATGTCGAAAGCCCGCTTCGAATTTCGCt                                                                                   >  1:10932575/1‑101 (MQ=255)
                    aCCTATAAGATTGCTGCCCACGCCGCTGACCTGGCGAAAGGGCATCCGGGCGCGCAAATTCGCGATAACGCCATGTCGAAAGCCCGCTTCGAATTTCGCTg                                                                                  <  1:3448636/101‑1 (MQ=255)
                     ccTATAAGATTGCTGCCCACGCCGCTGACCTGGCGAAAGGGCATCCGGGCGCGCAAATTCGCGATAACGCCATGTCGAAAGCCCGCTTCGAATTTCGCTgg                                                                                 <  1:8756154/101‑1 (MQ=255)
                      cTATAAGATTGCTGCCCACGCCGCTGACCTGGCGAAAGGGCATCCGGGCGCGCAAATTCGCGATAACGCCATGTCGAAAGCCCGCTTCGAATTTCGCTggg                                                                                >  1:2704308/1‑101 (MQ=255)
                        ataAGATTGCTGCCCACGCCGCTGACCTGGCGAAAGGGCATCCGGGCGCGCAAATTCGCGATAACGCCATGTCGAAAGCCCGCTTCGAATTTCGCTGGGaa                                                                              <  1:10132710/101‑1 (MQ=255)
                        ataAGATTGCTGCCCACGCCGCTGACCTGGCGAAAGGGCATCCGGGCGCGCAAATTCGCGATAACGCCATGTCGAAAGCCCGCTTCGAATTTCGCTGGGaa                                                                              <  1:5062547/101‑1 (MQ=255)
                        ataAGATTGCTGCCCACGCCGCTGACCTGGCGAAAGGGCATCCGGGCGCGCAAATTCGCGATAACGCCATGTCGAAAGCCCGCTTCGAATTTCGCTGGGaa                                                                              <  1:12147646/101‑1 (MQ=255)
                        acaaGATTTCCGCGCACGCCGCTGACCTGGCAAAGGGGCTGCCGCGCGCTCAGGAGCGTGACGACGCGCTGTCGAAGGCCCGTTTCGAGTTCCGTTggcat                                                                              <  1:6073106/99‑4 (MQ=255)
                          aaGATTTCCGCGCACGCCGCTGACCTGGCAAAGGGGCTGCCGCGCGCTCAGGAGCGTGACGACGCGCTGTCGAAGGCCCGTTTCGAGTTCCGTTggcatga                                                                            >  1:5872673/1‑96 (MQ=255)
                             aTTGCTGCCCACGCCGCTGACCTGGCGAAAGGGCATCCGGGCGCGCAAATTCGCGATAACGCCATGTCGAAAGCCCGCTTCGAATTTCGCTGGGAAGACCa                                                                         <  1:3071458/101‑1 (MQ=255)
                              tttccgcgcACGCCGCTGACCTGGCAAAGGGGCTGCCGCGCGCTCAGGAGCGTGACGACGCGCTGTCGAAGGCCCGTTTCGAGTTCCGTTGGCATGACCAg                                                                        <  1:9185394/93‑1 (MQ=255)
                              ttGCTGCCCACGCCGCTGACCTGGCGAAAGGGCATCCGGGCGCGCAAATTCGCGATAACGCCATGTCGAAAGCCCGCTTCGAATTTCGCTGGGAAGACCAg                                                                        <  1:10541257/101‑1 (MQ=255)
                               tgctgcCCACGCCGCTGACCTGGCGAAAGGGCATCCGGGCGCGCAAATTCGCGATAACGCCATGTCGAAAGCCCGCTTCGAATTTCGCTGGGAAGACCAGt                                                                       >  1:1778934/1‑101 (MQ=255)
                                 ctgcCCACGCCGCTGACCTGGCGAAAGGGCATCCGGGCGCGCAAATTCGCGATAACGCCATGTCGAAAGCCCGCTTCGAATTTCGCTGGGAAGACCAGttt                                                                     <  1:14184899/101‑1 (MQ=255)
                                 ccgcgcACGCCGCTGACCTGGCAAAGGGGCTGCCGCGCGCTCAGGAGCGTGACGACGCGCTGTCGAAGGCCCGTTTCGAGTTCCGTTGGCATGACCAGttt                                                                     >  1:6300775/6‑101 (MQ=255)
                                  tgcCCACGCCGCTGACCTGGCGAAAGGGCATCCGGGCGCGCAAATTCGCGATAACGCCATGTCGAAAGCCCGCTTCGAATTTCGCTGGGAAGACCAGTTTa                                                                    >  1:8636240/1‑101 (MQ=255)
                                  tgcCCACGCCGCTGACCTGGCGAAAGGGCATCCGGGCGCGCAAATTCGCGATAACGCCATGTCGAAAGCCCGCTTCGAATTTCGCTGGGAAGACCAGTTTa                                                                    >  1:3302822/1‑101 (MQ=255)
                                    cccACGCCGCTGACCTGGCGAAAGGGCATCCGGGCGCGCAAATTCGCGATAACGCCATGTCGAAAGCCCGCTTCGAATTTCGCTGGGAAGACCAGTTTAAt                                                                  >  1:7424086/1‑101 (MQ=255)
                                      cACGCCGCTGACCTGGCGAAAGGGCATCCGGGCGCGCAAATTCGCGATAACGCCATGTCGAAAGCCCGCTTCGAATTTCGCTGGGAAGACCAGTTTAATCt                                                                <  1:11845346/101‑1 (MQ=255)
                                       aCGCCGCTGACCTGGCGAAAGGGCATCCGGGCGCGCAAATTCGCGATAACGCCATGTCGAAAGCCCGCTTCGAATTTCGCTGGGAAGACCAGTTTAATCTg                                                               <  1:9817250/101‑1 (MQ=255)
                                         gccgcTGACCTGGCGAAAGGGCATCCGGGCGCGCAAATTCGCGATAACGCCATGTCGAAAGCCCGCTTCGAATTTCGCTGGGAAGACCAGTTTAATCTGGc                                                             >  1:7294018/1‑101 (MQ=255)
                                         gccgcTGACCTGGCGAAAGGGCATCCGGGCGCGCAAATTCGCGATAACGCCATGTCGAAAGCCCGCTTCGAATTTCGCTGGGAAGACCAGTTTAATCTGGc                                                             >  1:4323771/1‑101 (MQ=255)
                                         gccgcTGACCTGGCGAAAGGGCATCCGGGCGCGCAAATTCGCGATAACGCCATGTCGAAAGCCCGCTTCGAATTTCGCTGGGAAGACCAGTTTAATCTGGc                                                             >  1:4019823/1‑101 (MQ=255)
                                         gccgcTGACCTGGCGAAAGGGCATCCGGGCGCGCAAATTCGCGATAACGCCATGTCGAAAGCCCGCTTCGAATTTCGCTGGGAAGACCAGTTTAATCTGGc                                                             >  1:1933906/1‑101 (MQ=255)
                                           cgcTGACCTGGCGAAAGGGCATCCGGGCGCGCAAATTCGCGATAACGCCATGTCGAAAGCCCGCTTCGAATTTCGCTGGGAAGACCAGTTTAATCTGGccc                                                           >  1:2339488/1‑101 (MQ=255)
                                             cTGACCTGGCGAAAGGGCATCCGGGCGCGCAAATTCGCGATAACGCCATGTCGAAAGCCCGCTTCGAATTTCGCTGGGAAGACCAGTTTAATCTGGCCCTc                                                         <  1:10486459/101‑1 (MQ=255)
                                              tGACCTGGCGAAAGGGCATCCGGGCGCGCAAATTCGCGATAACGCCATGTCGAAAGCCCGCTTCGAATTTCGCTGGGAAGACCAGTTTAATCTGGCCCTCg                                                        >  1:2196635/1‑101 (MQ=255)
                                                aCCTGGCGAAAGGGCATCCGGGCGCGCAAATTCGCGATAACGCCATGTCGAAAGCCCGCTTCGAATTTCGCTGGGAAGACCAGTTTAATCTGGCCCTCGAc                                                      <  1:11722456/101‑1 (MQ=255)
                                                   tGGCGAAAGGGCATCCGGGCGCGCAAATTCGCGATAACGCCATGTCGAAAGCCCGCTTCGAATTTCGCTGGGAAGACCAGTTTAATCTGGCCCTCGACCCg                                                   >  1:3729348/1‑101 (MQ=255)
                                                            ggCATCCGGGCGCGCAAATTCGCGATAACGCCATGTCGAAAGCCCGCTTCGAATTTCGCTGGGAAGACCAGTTTAATCTGGCCCTCGACCCGTTTAccgcc                                          >  1:1117753/1‑101 (MQ=255)
                                                              cATCCGGGCGCGCAAATTCGCGATAACGCCATGTCGAAAGCCCGCTTCGAATTTCGCTGGGAAGACCAGTTTAATCTGGCCCTCGACCCGTTTAccgcccg                                        >  1:5432572/1‑101 (MQ=255)
                                                               aTCCGGGCGCGCAAATTCGCGATAACGCCATGTCGAAAGCCCGCTTCGAATTTCGCTGGGAAGACCAGTTTAATCTGGCCCTCGACCCGTTTAccgcccgc                                       <  1:11825250/101‑1 (MQ=255)
                                                               aTCCGGGCGCGCAAATTCGCGATAACGCCATGTCGAAAGCCCGCTTCGAATTTCGCTGGGAAGACCAGTTTAATCTGGCCCTCGACCCGTTTAccgcccgc                                       >  1:3667008/1‑101 (MQ=255)
                                                               aTCCGGGCGCGCAAATTCGCGATAACGCCATGTCGAAAGCCCGCTTCGAATTTCGCTGGGAAGACCAGTTTAATCTGGCCCTCGACCCGTTTAccgcccgc                                       >  1:9647336/1‑101 (MQ=255)
                                                                   gggCGCGCAAATTCGCGATAACGCCATGTCGAAAGCCCGCTTCGAATTTCGCTGGGAAGACCAGTTTAATCTGGCCCTCGACCCGTTTACCGCCCGCGCtt                                   >  1:11410720/1‑101 (MQ=255)
                                                                   gggCGCGCAAATTCGCGATAACGCCATGTCGAAAGCCCGCTTCGAATTTCGCTGGGAAGACCAGTTTAATCTGGCCCTCGACCCGTTTACCGCCCGCGCtt                                   <  1:5227467/101‑1 (MQ=255)
                                                                     gcgcgcAAATTCGCGATAACGCCATGTCGAAAGCCCGCTTCGAATTTCGCTGGGAAGACCAGTTTAATCTGGCCCTCGACCCGTTTACCGCCCGCGCTTAt                                 <  1:1310459/101‑1 (MQ=255)
                                                                      cgcgcAAATTCGCGATAACGCCATGTCGAAAGCCCGCTTCGAATTTCGCTGGGAAGACCAGTTTAATCTGGCCCTCGACCCGTTTACCGCCCGCGCTTATc                                >  1:2006819/1‑101 (MQ=255)
                                                                      cgcgcAAATTCGCGATAACGCCATGTCGAAAGCCCGCTTCGAATTTCGCTGGGAAGACCAGTTTAATCTGGCCCTCGACCCGTTTACCGCCCGCGCTTATc                                >  1:5904403/1‑101 (MQ=255)
                                                                      cgcgcAAATTCGCGATAACGCCATGTCGAAAGCCCGCTTCGAATTTCGCTGGGAAGACCAGTTTAATCTGGCCCTCGACCCGTTTACCGCCCGCGCTTATc                                >  1:3272280/1‑101 (MQ=255)
                                                                        cgcAAATTCGCGATAACGCCATGTCGAAAGCCCGCTTCGAATTTCGCTGGGAAGACCAGTTTAATCTGGCCCTCGACCCGTTTACCGCCCGCGCTTATCAc                              <  1:8266665/101‑1 (MQ=255)
                                                                         gcAAATTCGCGATAACGCCATGTCGAAAGCCCGCTTCGAATTTCGCTGGGAAGACCAGTTTAATCTGGCCCTCGACCCGTTTACCGCCCGCGCTTATCACg                             <  1:6508265/101‑1 (MQ=255)
                                                                          cAAATTCGCGATAACGCCATGTCGAAAGCCCGCTTCGAATTTCGCTGGGAAGACCAGTTTAATCTGGCCCTCGACCCGTTTACCGCCCGCGCTTATCACGa                            >  1:7077919/1‑101 (MQ=255)
                                                                           aaaTTCGCGATAACGCCATGTCGAAAGCCCGCTTCGAATTTCGCTGGGAAGACCAGTTTAATCTGGCCCTCGACCCGTTTACCGCCCGCGCTTATCACGAt                           >  1:12811841/1‑101 (MQ=255)
                                                                           aaaTTCGCGATAACGCCATGTCGAAAGCCCGCTTCGAATTTCGCTGGGAAGACCAGTTTAATCTGGCCCTCGACCCGTTTACCGCCCGCGCTTATCACGAt                           <  1:7938633/101‑1 (MQ=255)
                                                                           aaaTTCGCGATAACGCCATGTCGAAAGCCCGCTTCGAATTTCGCTGGGAAGACCAGTTTAATCTGGCCCTCGACCCGTTTACCGCCCGCGCTTATCACGAt                           <  1:1518490/101‑1 (MQ=255)
                                                                                cgcgATAACGCCATGTCGAAAGCCCGCTTCGAATTTCGCTGGGAAGACCAGTTTAATCTGGCCCTCGACCCGTTTACCGCCCGCGCTTATCCCGATAAAAc                      >  1:14025900/1‑101 (MQ=255)
                                                                                cgcgATAACGCCATGTCGAAAGCCCGCTTCGAATTTCGCTGGGAAGACCAGTTTAATCTGGCCCTCGACCCGTTTACCGCCCGCGCTTATCACGATGAAAc                      <  1:10088698/101‑1 (MQ=255)
                                                                                 gcgATAACGCCATGTCGAAAGCCCGCTTCGAATTTCGCTGGGAAGACCAGTTTAATCTGGCCCTCGACCCGTTTACCGCCCGCGCTTATCACGATGAAAcc                     <  1:13544987/101‑1 (MQ=255)
                                                                                  cgATAACGCCATGTCGAAAGCCCGCTTCGAATTTCGCTGGGAAGACCAGTTTAATCTGGCCCTCGACCCGTTTACCGCCCGCGCTTATCACGATGAAAcac                    >  1:1209374/1‑99 (MQ=255)
                                                                                     tAACGCCATGTCGAAAGCCCGCTTCGAATTTCGCTGGGAAGACCAGTTTAATCTGGCCCTCGACCCGTTTACCGCCCGCGCTTATCACGATGAAACCCTGc                 >  1:8191347/1‑101 (MQ=255)
                                                                                      aaCGCCATGTCGAAAGCCCGCTTCGAATTTCGCTGGGAAGACCAGTTTAATCTGGCCCTCGACCCGTTTACCGCCCGCGCTTATCACGATGAAACCCTGcc                >  1:2638909/1‑101 (MQ=255)
                                                                                       aCGCCATGTCGAAAGCCCGCTTCGAATTTCGCTGGGAAGACCAGTTTAATCTGGCCCTCGACCCGTTTACCGCCCGCGCTTATCACGATGAAACCCTGCCg               >  1:1276996/1‑101 (MQ=255)
                                                                                         gCCATGTCGAAAGCCCGCTTCGAATTTCGCTGGGAAGACCAGTTTAATCTGGCCCTCGACCCGTTTACCGCCCGCGCTTATCACGATGAAACCCTGCCGCa             >  1:645996/1‑101 (MQ=255)
                                                                                           caagTCGAAAACCCGCTTCCAATTTCACTGGGAAGACCAGTTTAATCTGGCCCTGGACCCGTTTACCGCCCGCGCTTATCACGATGAAACCCTGCCGCAag           <  1:13378015/98‑1 (MQ=255)
                                                                                           cATGTCGAAAGCCCGCTTCGAATTTCGCTGGGAAGACCAGTTTAATCTGGCCCTCGACCCGTTTACCGCCCGCGCTTATCACGATGAAACCCTGCCGCAag           >  1:4081568/1‑101 (MQ=255)
                                                                                            aTGTCGAAAGCCCGCTTCGAATTTCGCTGGGAAGACCAGTTTAATCTGGCCCTCGACCCGTTTACCGCCCGCGCTTATCACGATGAAACCCTGCCGCAaga          <  1:11068281/101‑1 (MQ=255)
                                                                                            aTGTCGAAAGCCCGCTTCGAATTTCGCTGGGAAGACCAGTTTAATCTGGCCCTCGACCCGTTTACCGCCCGCGCTTATCACGATGAAACCCTGCCGCAaga          <  1:3905920/101‑1 (MQ=255)
                                                                                            aTGTCGAAAGCCCGCTTCGAATTTCGCTGGGAAGACCAGTTTAATCTGGCCCTCGACCCGTTTACCGCCCGCGCTTATCACGATGAAACCCTGCCGCAaga          >  1:4551437/1‑101 (MQ=255)
                                                                                             tGTCGAAAGCCCGCTTCGAATTTCGCTGGGAAGACCAGTTTAATCTGGCCCTCGACCCGTTTACCGCCCGCGCTTATCACGATGAAACCCTGCCGCAagag         <  1:6852889/101‑1 (MQ=255)
                                                                                             tGTCGAAAGCCCGCTTCGAATTTCGCTGGGAAGACCAGTTTAATCTGGCCCTCGACCCGTTTACCGCCCGCGCTTATCACGATGAAACCCTGCCGCAagag         <  1:8911333/101‑1 (MQ=255)
                                                                                             tGTCGAAAGCCCGCTTCGAATTTCGCTGGGAAGACCAGTTTAATCTGGCCCTCGACCCGTTTACCGCCCGCGCTTATCACGATGAAACCCTGCCGCAagag         >  1:9380927/1‑101 (MQ=255)
                                                                                              gTCGAAAGCCCGCTTCGAATTTCGCTGGGAAGACCAGTTTAATCTGGCCCTCGACCCGTTTACCGCCCGCGCTTATCACGATGAAACCCTGCCGCAAGAGt        <  1:564305/101‑1 (MQ=255)
                                                                                                cGAAAGCCCGCTTCGAATTTCGCTGGGAAGACCAGTTTAATCTGGCCCTCGACCCGTTTACCGCCCGCGCTTATCACGATGAAACCCTGCCGCAAGAGTCa      >  1:7925601/1‑101 (MQ=255)
                                                                                                cGAAAGCCCGCTTCGAATTTCGCTGGGAAGACCAGTTTAATCTGGCCCTCGACCCGTTTACCGCCCGCGCTTATCACGATGAAACCCTGCCGCAAGAGTCa      <  1:14089060/101‑1 (MQ=255)
                                                                                                 gAAAGCCCGCTTCGAATTTCGCTGGGAAGACCAGTTTAATCTGGCCCTCGACCCGTTTACCGCCCGCGCTTATCACGATGAAACCCTGCCGCAAGAGTCAg     >  1:3958380/1‑101 (MQ=255)
                                                                                                  aaaGCCCGCTTCGAATTTCGCTGGGAAGACCAGTTTAATCTGGCCCTCGACCCGTTTACCGCCCGCGCTTATCACGATGAAACCCTGCCGCAAGAGTCAgg    <  1:13312314/101‑1 (MQ=255)
                                                                                                  aaaGCCCGCTTCGAATTTCGCTGGGAAGACCAGTTTAATCTGGCCCTCGACCCGTTTACCGCCCGCGCTTATCACGATGAAACCCTGCCGCAAGAGTCAgg    >  1:610989/1‑101 (MQ=255)
                                                                                                   aaGCCCGCTTCGAATTTCGCTGGGAAGACCAGTTTAATCTGGCCCTCGACCCGTTTACCGCCCGCGCTTATCACGATGAAACCCTGCCGCAAGAGTCAGGt   <  1:4664742/101‑1 (MQ=255)
                                                                                                   aaGCCCGCTTCGAATTTCGCTGGGAAGACCAGTTTAATCTGGCCCTCGACCCGTTTACCGCCCGCGCTTATCACGATGAAACCCTGCCGCAAGAGTCAGGt   <  1:11794378/101‑1 (MQ=255)
                                                                                                   aaGCCCGCTTCGAATTTCGCTGGGAAGACCAGTTTAATCTGGCCCTCGACCCGTTTACCGCCCGCGCTTATCACGATGAAACCCTGCCGCAAGAGTCAGGt   >  1:12441342/1‑101 (MQ=255)
                                                                                                   aaGCCCGCTTCGAATTTCGCTGGGAAGACCAGTTTAATCTGGCCCTCGACCCGTTTACCGCCCGCGCTTATCACGATGAAACCCTGCCGCAAGAGTCAGGt   >  1:30375/1‑101 (MQ=255)
                                                                                                   aaGCCCGCTTCGAATTTCACTGGGAAGACCAGTTTAATCTGGCCCTCGACCCGTTTACCGCCCGCGCTTATCACGATGAAACCCTGCCGCAAGAGTCAGGt   >  1:13057560/1‑101 (MQ=255)
                                                                                                    aGCCCGCTTCGAATTTCGCTGGGAAGACCAGTTTAATCTGGCCCTCGACCCGTTTACCGCCCGCGCTTATCACGATGAAACCCTGCCGCAAGAGTCAGGTa  <  1:9257424/101‑1 (MQ=255)
                                                                                                    |                                                                                                    
ATGTTAAGCAGGGGCTTATCACCTATAAGATTGCTGCCCACGCCGCTGACCTGGCGAAAGGGCATCCGGGCGCGCAAATTCGCGATAACGCCATGTCGAAAGCCCGCTTCGAATTTCGCTGGGAAGACCAGTTTAATCTGGCCCTCGACCCGTTTACCGCCCGCGCTTATCACGATGAAACCCTGCCGCAAGAGTCAGGTA  >  NC_007779/3442107‑3442307

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.