Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 3,441,379 G→C 14.6% T266T (ACG→ACC)  thiC → thiamin (pyrimidine moiety) biosynthesis protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0077793,441,3790G→C14.6% 122.4 ‑3.8 ‑0.1267T266T (ACG→ACC) thiCthiamin (pyrimidine moiety) biosynthesis protein
Reads supporting (aligned to +/- strand):  new base (13/26):  ref base (123/105):  total (136/131)
Fisher's exact test strand distribution p-value = 2.35e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 8.69e-01

Only 100 of 271 total aligned reads displayed.

TGGTAAAAGTTAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCGATCGGTACAGTGCCGATCTACCAGGCGC  >  NC_007779/3441279‑3441479
                                                                                                    |                                                                                                    
tGGTAAAAGTTAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACg                                                                                                      <  1:11125850/101‑1 (MQ=255)
 ggTAAAAGTTAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACgg                                                                                                     <  1:8638510/101‑1 (MQ=255)
   tAAAAGTTAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTg                                                                                                   <  1:13558903/101‑1 (MQ=255)
       aGTTAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATgg                                                                                               >  1:12284735/1‑101 (MQ=255)
       aGTTAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATgg                                                                                               <  1:5352720/101‑1 (MQ=255)
       aGTTAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATgg                                                                                               >  1:6880764/1‑101 (MQ=255)
        gTTAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATggc                                                                                              <  1:14472228/101‑2 (MQ=255)
          tAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGAtc                                                                                            <  1:6683208/101‑1 (MQ=255)
          tAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGAtc                                                                                            <  1:4549898/101‑1 (MQ=255)
           aaCGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGAtct                                                                                           <  1:9260622/101‑1 (MQ=255)
           aaCGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGAtct                                                                                           >  1:2903707/1‑101 (MQ=255)
            aCGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGAtctc                                                                                          >  1:4573758/1‑101 (MQ=255)
             cGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGAtctct                                                                                         <  1:9651994/101‑1 (MQ=255)
             cGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGAtctct                                                                                         >  1:6372279/1‑101 (MQ=255)
               ccAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTcc                                                                                       >  1:2593022/1‑101 (MQ=255)
                cAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCa                                                                                      >  1:14285102/1‑101 (MQ=255)
                 aaTATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCAc                                                                                     >  1:7304500/1‑101 (MQ=255)
                 aaTATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCAc                                                                                     >  1:5710648/1‑101 (MQ=255)
                   tatCGGCAACTCGGCGGTCACCTCTTCCATCGTAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCg                                                                                   <  1:8839886/101‑1 (MQ=255)
                    atCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCgg                                                                                  >  1:1984313/1‑101 (MQ=255)
                    atCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCgg                                                                                  >  1:593296/1‑101 (MQ=255)
                    atCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCgg                                                                                  >  1:486414/1‑101 (MQ=255)
                     tCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGt                                                                                 <  1:5226558/101‑1 (MQ=255)
                     tCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGt                                                                                 >  1:8275148/1‑101 (MQ=255)
                      cGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTc                                                                                >  1:13339685/1‑101 (MQ=255)
                      cGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTc                                                                                <  1:830944/101‑1 (MQ=255)
                       ggCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCg                                                                               >  1:13030555/1‑101 (MQ=255)
                          aaCTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCta                                                                            >  1:13007628/1‑101 (MQ=255)
                           aCTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCtat                                                                           >  1:13281582/1‑101 (MQ=255)
                           aCTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCtat                                                                           >  1:255727/1‑101 (MQ=255)
                            cTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCtata                                                                          >  1:3376171/1‑101 (MQ=255)
                              cggcggTCCCCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATCCGGTGAT‑GATCTCTCCACCGGTCGCTATATTc                                                                       <  1:8959165/101‑1 (MQ=255)
                               ggcggTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTc                                                                       <  1:13447622/101‑1 (MQ=255)
                               ggcggTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTc                                                                       >  1:14840825/1‑101 (MQ=255)
                                  ggTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACg                                                                    >  1:4733432/1‑101 (MQ=255)
                                  ggTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACg                                                                    <  1:12274079/101‑1 (MQ=255)
                                     cACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGaaa                                                                 >  1:5340291/1‑101 (MQ=255)
                                     cACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGaaa                                                                 >  1:13191393/1‑101 (MQ=255)
                                      cccTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAAc                                                                <  1:7417032/100‑1 (MQ=255)
                                      aCCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAAc                                                                >  1:3045128/1‑101 (MQ=255)
                                       ccTCTTCCATCGAAGAAGAAGCGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAAcc                                                               <  1:8373270/101‑1 (MQ=255)
                                          ctTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCcgc                                                            >  1:10430107/1‑101 (MQ=255)
                                          ctTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCcgc                                                            <  1:11108618/101‑1 (MQ=255)
                                           ttCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCcgcg                                                           <  1:7206490/101‑1 (MQ=255)
                                           ttCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCcgcg                                                           <  1:11845316/101‑1 (MQ=255)
                                            tCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGa                                                          >  1:4021795/1‑101 (MQ=255)
                                             ccATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAg                                                         >  1:13655026/1‑101 (MQ=255)
                                              catctccgaagaggTGGAGAAGATGGTGTGGGCGACCCGTTGGGGCGCGGACACCATCATGGACCTGTCCACCGGTAAGGACATTCACGAAACCCGCGAGt                                                        >  1:6488772/8‑101 (MQ=255)
                                              cATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGt                                                        <  1:15039760/101‑1 (MQ=255)
                                              cATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGt                                                        >  1:6792129/1‑101 (MQ=255)
                                                tctccgaagaggTGGAGAAGATGGTGTGGGCGACCCGTTGGGGCGCGGACACCATCATGGACCTGTCCACCGGTAAGGACATTCACGAAACCCGCGAGTgg                                                      <  1:4985198/96‑1 (MQ=255)
                                                tCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTgg                                                      >  1:21321/1‑101 (MQ=255)
                                                  tccgaagaggTGGAGAAGATGGTGTGGGCGACCCGTTGGGGCGCGGACACCATCATGGACCTGTCCACCGGTAAGGACATTCACGAAACCCGCGAGTGGAt                                                    >  1:5551260/4‑101 (MQ=255)
                                                      aagaggTGGAGAAGATGGTGTGGGCGACCCGTTGGGGCGCGGACACCATCATGGACCTGACCACCGGTAAGGACATTCACGAAACCCGCGAGTGGATTCTg                                                >  1:10880419/1‑101 (MQ=255)
                                                      aagaagTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTg                                                >  1:10717848/1‑101 (MQ=255)
                                                      aagaagTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTg                                                >  1:1205864/1‑101 (MQ=255)
                                                       agaagTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGc                                               <  1:161708/101‑1 (MQ=255)
                                                        gaagTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCg                                              >  1:14720665/1‑101 (MQ=255)
                                                         aagTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGt                                             <  1:9244449/101‑1 (MQ=255)
                                                         aagTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGt                                             >  1:14053133/1‑101 (MQ=255)
                                                         aagTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGt                                             >  1:14594355/1‑101 (MQ=255)
                                                          agTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTa                                            <  1:12075842/101‑1 (MQ=255)
                                                          agTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTa                                            <  1:9833529/101‑1 (MQ=255)
                                                           gTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTaa                                           >  1:5473709/1‑101 (MQ=255)
                                                            tGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAAc                                          <  1:13075571/101‑1 (MQ=255)
                                                              gagaaGATGGTGTGGGCGACCCGTTGGGGCGCGGACACCATCATGGACCTGTCCACCGGTAAGGACATTCACGAAACCCGCGAGTGGATTCTGCGTAActc                                        <  1:5798036/98‑3 (MQ=255)
                                                               aaaaGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGc                                       >  1:15115439/1‑101 (MQ=255)
                                                                aaaGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGcc                                      <  1:6995050/101‑1 (MQ=255)
                                                                 aaGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGccc                                     <  1:4101435/101‑1 (MQ=255)
                                                                 aaGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGccc                                     <  1:2385759/101‑1 (MQ=255)
                                                                  agatGGTGTGGGCGACCCGTTGGGGCGCGGACACCATCATGGACCTGTCCACCGGTAAGGACATTCACGAAACCCGCGAGTGGATTCTGCGTAACTCCCCg                                    <  1:613476/98‑1 (MQ=255)
                                                                  aGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCg                                    <  1:7842378/101‑1 (MQ=255)
                                                                    cTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGt                                  >  1:3558040/1‑101 (MQ=255)
                                                                    cTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGt                                  <  1:12428216/101‑1 (MQ=255)
                                                                     tGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTg                                 >  1:2751262/1‑101 (MQ=255)
                                                                      gttATGGTCCCCGCCCTGGGGGCCGGTTCCGTGGTTGGTCCCCCCCCCCGGTCCCTTTA‑TCCCCAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGcc                               <  1:4377175/99‑1 (MQ=255)
                                                                        tATGGTCCAGGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCg                              >  1:10130698/1‑101 (MQ=255)
                                                                        tATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCg                              >  1:11368602/1‑101 (MQ=255)
                                                                          tGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCGAt                            <  1:992744/101‑1 (MQ=255)
                                                                           ggTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCGATc                           <  1:13845646/101‑1 (MQ=255)
                                                                              cgacccgttGGGGCGCGGACACCATCATGGACCTGTCCACCGGTAAGGACATTCACGAAACCCGCGAGTGGATTCTGCGTAACTCCCCGGTGCCGGTTGGt                        <  1:3850763/93‑1 (MQ=255)
                                                                              ccACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCGATCGGt                        >  1:12859614/1‑101 (MQ=255)
                                                                                acccgttGGGGCGCGGACACCATCATGGACCTGTCCACCGGTAAGGACATTCACGAAACCCGCGAGTGGATTCTGCGTAACTCCCCGGTGCCGGTTGGTAc                      <  1:11001462/95‑1 (MQ=255)
                                                                                aCGTGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCGATCGGTAc                      <  1:955583/101‑1 (MQ=255)
                                                                                  ccgttGGGGCGCGGACACCATCATGGACCTGTCCACCGGTAAGGACATTCACGAAACCCGCGAGTGGATTCTGCGTAACTCCCCGGTGCCGGTTGGTAcgg                    <  1:14678439/97‑3 (MQ=255)
                                                                                    gcTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCGATCGGTACAGTg                  <  1:173427/101‑1 (MQ=255)
                                                                                    gcTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCGATCGGTACAGTg                  <  1:13951297/101‑1 (MQ=255)
                                                                                     cTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCGATCGGTACAGTGc                 <  1:7732656/101‑1 (MQ=255)
                                                                                      tGGGGCGCGGACACCATCATGGACCTGTCCACCGGTAAGGACATTCACGAAACCCGCGAGTGGATTCTGCGTAACTCCCCGGTGCCGGTTGGTACGGTGcc                <  1:9108653/101‑1 (MQ=255)
                                                                                           aGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCGATCGGTACAGTGCCGATCt           <  1:4412078/101‑1 (MQ=255)
                                                                                           aGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGACCGGTGCCGATCGGTACAGTGCCGATCt           >  1:9290323/1‑101 (MQ=255)
                                                                                           aGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGACCGGTGCCGATCGGTACAGTGCCGATCt           >  1:6202976/1‑101 (MQ=255)
                                                                                             cGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCGATCGGTACAGTGCCGATCTAc         <  1:12241109/101‑1 (MQ=255)
                                                                                             cGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCGATCGGTACAGTGCCGATCTAc         >  1:1009113/1‑101 (MQ=255)
                                                                                              ggATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCGATCGGTACAGTGCCGATCTAcc        <  1:9397888/101‑1 (MQ=255)
                                                                                              ggATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGACCGGTGCCGATCGGTACAGTGCCGATCTAcc        >  1:3150305/1‑101 (MQ=255)
                                                                                                aTACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCGATCGGTACAGTGCCGATCTACCAg      <  1:2686515/101‑1 (MQ=255)
                                                                                                 tACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCGATCGGTACAGTGCCGATCTACCAgg     >  1:3090999/1‑101 (MQ=255)
                                                                                                  aCGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCGATCGGTACAGTGCCGATCTACCAGgc    >  1:3456959/1‑101 (MQ=255)
                                                                                                   cGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCGATCGGTACAGTGCCGATCTACCAGgcg   >  1:4920291/1‑101 (MQ=255)
                                                                                                    |                                                                                                    
TGGTAAAAGTTAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCGATCGGTACAGTGCCGATCTACCAGGCGC  >  NC_007779/3441279‑3441479

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