Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 3,441,967 G→C 7.2% L462L (CTG→CTC thiC → thiamin (pyrimidine moiety) biosynthesis protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0077793,441,9670G→C 7.2% 43.8 ‑1.6 ‑1.8223L462L (CTG→CTCthiCthiamin (pyrimidine moiety) biosynthesis protein
Reads supporting (aligned to +/- strand):  new base (11/5):  ref base (101/106):  total (112/111)
Fisher's exact test strand distribution p-value = 1.93e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.59e-02

Only 100 of 241 total aligned reads displayed.

AGGTGATGATTGAAGGCCCAGGCCACGTGCCGATGCAGATGATCCGCCGCAATATGACCGAGGAGTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACG  >  NC_007779/3441867‑3442067
                                                                                                    |                                                                                                    
aGGTGATGATTGAAGGCCCAGGCCACGTGCCGATGCAGATGATCCGCCGCAATATGACCGAGGAGTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTg                                                                                                      >  1:3828526/1‑101 (MQ=255)
aGGTGATGATTGAAGGCCCAGGCCACGTGCCGATGCAGATGATCCGCCGCAATATGACCGAGGAGTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTg                                                                                                      >  1:2500787/1‑101 (MQ=255)
aGGTGATGATTGAAGGCCCAGGCCACGTGCCGATGCAGATGATCCGCCGCAATATGACCGAGGAGTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTg                                                                                                      <  1:8169107/101‑1 (MQ=255)
 ggTGATGATTGAAGGCCCAGGCCACGTGCCGATGCAGATGATCCGCCGCAATATGACCGAGGAGTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTgg                                                                                                     <  1:3502335/101‑1 (MQ=255)
    gatgatTGAAGGCCCAGGCCACGTGCCGATGCAGATGATCCGCCGCAATATGACCGAGGAGTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGc                                                                                                  >  1:7180593/1‑101 (MQ=255)
     atgatTGAAGGCCCAGGCCACGTGCCGATGCAGATGATCCGCCGCAATATGACCGAGGAGTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGcc                                                                                                 <  1:6060159/101‑1 (MQ=255)
     atgatTGAAGGCCCAGGCCACGTGCCGATGCAGATGATCCGCCGCAATATGACCGAGGAGTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGcc                                                                                                 <  1:4427446/101‑1 (MQ=255)
         ttGAAGGCCCAGGCCACGTGCCGATGCAGATGATCCGCCGCAATATGACCGAGGAGTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTa                                                                                             <  1:5808719/101‑1 (MQ=255)
         ttGAAGGCCCAGGCCACGTGCCGATGCAGATGATCCGCCGCAATATGACCGAGGAGTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTa                                                                                             <  1:10361288/101‑1 (MQ=255)
         ttGAAGGCCCAGGCCACGTGCCGATGCAGATGATCCGCCGCAATATGACCGAGGAGTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTa                                                                                             <  1:4029106/101‑1 (MQ=255)
         ttGAAGGCCCAGGCCACGTGCCGATGCAGATGATCCGCCGCAATATGACCGAGGAGTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTa                                                                                             <  1:4913298/101‑1 (MQ=255)
         ttGAAGGCCCAGGCCACGTGCCGATGCAGATGATCCGCCGCAATATGACCGAGGAGTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTa                                                                                             <  1:789169/101‑1 (MQ=255)
           gAAGGCCCAGGCCACGTGCCGATGCAGATGATCCGCCGCAATATGACCGAGGAGTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAAc                                                                                           >  1:125499/1‑101 (MQ=255)
            aaGGCCCAGGCCACGTGCCGATGCAGATGATCCGCCGCAATATGACCGAGGAGTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACt                                                                                          >  1:5662039/1‑101 (MQ=255)
             aGGCCCAGGCCACGTGCCGATGCAGATGATCCGCCGCAATATGACCGAGGAGTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTa                                                                                         <  1:4134955/101‑1 (MQ=255)
             aGGCCCAGGCCACGTGCCGATGCAGATGATCCGCCGCAATATGACCGAGGAGTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTa                                                                                         <  1:7389096/101‑1 (MQ=255)
             aGGCCCAGGCCACGTGCCGATGCAGATGATCCGCCGCAATATGACCGAGGAGTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTa                                                                                         <  1:5908197/101‑1 (MQ=255)
               gCCCAGGCCACGTGCCGATGCAGATGATCCGCCGCAATATGACCGAGGAGTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTAcc                                                                                       >  1:11695370/1‑101 (MQ=255)
                cccAGGCCACGTGCCGATGCAGATGATCCGCCGCAATATGACCGAGGAGTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCg                                                                                      <  1:626289/101‑1 (MQ=255)
                    ggCCACGTGCCGATGCAGATGATCCGCCGCAATATGACCGAGGAGTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGatat                                                                                  <  1:11844424/101‑1 (MQ=255)
                    ggCCACGTGCCGATGCAGATGATCCGCCGCAATATGACCGAGGAGTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGatat                                                                                  >  1:5427219/1‑101 (MQ=255)
                     gCCACGTGCCGATGCAGATGATCCTCCGCAATATGACCGAGGAGTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGATGCCGCTAACTACCGATAtt                                                                                 >  1:13374261/1‑101 (MQ=255)
                     gCCACGTGCCGATGCAGATGATCCGCCGCAATATGACCGAGGAGTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATAtt                                                                                 >  1:375862/1‑101 (MQ=255)
                     gCCACGTGCCGATGCAGATGATCCGCCGCAATATGACCGAGGAATTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATAtt                                                                                 <  1:3204386/101‑1 (MQ=255)
                      ccACGTGCCGATGCAGATGATCCGCCGCAATATGACCGAGGAGTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTg                                                                                <  1:4522979/101‑1 (MQ=255)
                       cACGTGCCGATGCAGATGATCCGCCGCAATATGACCGAGGAGTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTgc                                                                               >  1:14155768/1‑101 (MQ=255)
                       cACGTGCCGATGCAGATGATCCGCCGCAATATGACCGAGGAGTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTgc                                                                               >  1:13915940/1‑101 (MQ=255)
                       cACGTGCCGATGCAGATGATCCGCCGCAATATGACCGAGGAGTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTgc                                                                               >  1:6419834/1‑101 (MQ=255)
                            gCCGATGCAGATGATCCGCCGCAATATGACCGAGGAGTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCgg                                                                          >  1:6598314/1‑101 (MQ=255)
                            gCCGATGCAGATGATCCGCCGCAATATGACCGAGGAGTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCgg                                                                          <  1:7717827/101‑1 (MQ=255)
                              cGATGCAGATGATCCGCCGCAATATGACCGAGGAGTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGc                                                                        >  1:8264628/1‑101 (MQ=255)
                              cGATGCAGATGATCCGCCGCAATATGACCGAGGAGTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGc                                                                        <  1:12338647/101‑1 (MQ=255)
                              cGATGCAGATGATCCGCCGCAATATGACCGAGGAGTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGc                                                                        <  1:6986944/101‑1 (MQ=255)
                                gtgctctaactCCGCCGCAATATGACCGAGGAGTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTa                                                                      >  1:11238656/11‑101 (MQ=255)
                                aTGCAGATGATCCGCCGCAATATGACCGAGGAGTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTa                                                                      >  1:10559342/1‑101 (MQ=255)
                                aTGCAGATGATCCGCCGCAATATGACCGAGGAGTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTa                                                                      >  1:9009488/1‑101 (MQ=255)
                                  gCAGATGATCCGCCGCAATATGACCGAGGAGTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATg                                                                    >  1:2403646/1‑101 (MQ=255)
                                   cAGATGATCCGCCGCAATATGACCGAGGAGTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGa                                                                   <  1:8384865/101‑1 (MQ=255)
                                   cAGATGATCCGCCGCAATATGACCGAGGAGTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGa                                                                   >  1:8044305/1‑101 (MQ=255)
                                    aGATGATCCGCCGCAATATGACCGAGGAGTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGAc                                                                  <  1:4568068/101‑1 (MQ=255)
                                    aGATGATCCGCCGCAATATGACCGAGGAGTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGAc                                                                  >  1:7467821/1‑101 (MQ=255)
                                      atgatCCGCCGCAATATGACCGAGGAGTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCa                                                                >  1:14085103/1‑101 (MQ=255)
                                      atgatCCGCCGCAATATGACCGAGGAGTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCa                                                                <  1:307976/101‑1 (MQ=255)
                                        gatCCGCCGCAATATGACCGAGGAGTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACt                                                              >  1:10854722/1‑101 (MQ=255)
                                             gccgcAATATGACCGAGGAGTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACg                                                         >  1:5150908/1‑101 (MQ=255)
                                              ccgcAATATGACCGAGGAGTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGt                                                        <  1:13989340/101‑1 (MQ=255)
                                               cgcAATATGACCGAGGAGTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTc                                                       <  1:13257584/101‑1 (MQ=255)
                                                gcAATATGACCGAGGAGTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCg                                                      >  1:10371043/1‑101 (MQ=255)
                                                  aaTATGACCGAGGAGTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCggg                                                    <  1:4843432/101‑1 (MQ=255)
                                                  aaTATGACCGAGGAGTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCggg                                                    <  1:13430808/101‑1 (MQ=255)
                                                     atGACCGAGGAGTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGAt                                                 >  1:6243427/1‑101 (MQ=255)
                                                      tGACCGAGGAGTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGAtt                                                <  1:520305/101‑1 (MQ=255)
                                                         ccGAGGAGTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGt                                             <  1:7744808/101‑1 (MQ=255)
                                                          cGAGGAGTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTg                                            >  1:10130077/1‑101 (MQ=255)
                                                          cGAGGAGTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTg                                            >  1:1669206/1‑101 (MQ=255)
                                                             ggagGTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTGCCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTgcgg                                         >  1:12945847/1‑101 (MQ=255)
                                                               agTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTgcggcg                                       >  1:2734374/1‑101 (MQ=255)
                                                               agTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTgcggcg                                       <  1:9530957/101‑1 (MQ=255)
                                                                gTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCga                                      >  1:467872/1‑101 (MQ=255)
                                                                gTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCga                                      <  1:5078680/101‑1 (MQ=255)
                                                                  tAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCgatg                                    <  1:13805272/101‑1 (MQ=255)
                                                                   agagCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCgatga                                   >  1:13253561/1‑101 (MQ=255)
                                                                   agagCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCgatga                                   >  1:11993120/1‑101 (MQ=255)
                                                                   agagCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCgatga                                   >  1:10972058/1‑101 (MQ=255)
                                                                      gCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTg                                >  1:1407081/1‑101 (MQ=255)
                                                                      gCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTg                                <  1:11645228/101‑1 (MQ=255)
                                                                       cACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTgg                               >  1:1309990/1‑101 (MQ=255)
                                                                        ggtgtgaggAAGCTCCGTTCTACACCCTCGGCCCGCTGGCAACCGACATCGCGCCGGGCTATGACCACATCACGTCGGCTATCGGCGCGGCAATCATTgcg                              <  1:3019799/93‑3 (MQ=255)
                                                                        aCTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGc                              >  1:11314013/1‑101 (MQ=255)
                                                                          tgtgaggAAGCTCCGTTCTACACCCTCGGCCCGCTGGCAACCGACATCGCGCCGGGCTATGACCACATCACGTCGGCTATCGGCGCGGCAATCATTgcgca                            <  1:11374226/95‑5 (MQ=255)
                                                                           gtgaggAAGCTCCGTTCTACACCCTCGGCCCGCTGGCAACCGACATCGCGCCGGGCTATGACCACATCACGTCGGCTATCGGCGCGGCAATCATTgcgcag                           >  1:9639816/6‑96 (MQ=255)
                                                                           gtgaggAAGCTCCGTTCTACACCCTCGGCCCGCTGGCAACCGACATCGCGCCGGGCTATGACCACATCACGTCGGCTATCGGCGCGGCAATCATTgcgcag                           >  1:573771/6‑96 (MQ=255)
                                                                           gCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTgg                           <  1:8626153/101‑1 (MQ=255)
                                                                            tgaggAAGCTCCGTTCTACACCCTCGGCCCTCTGGCAACCGACATCGCGCCGGGCTATGACCACATCACGTCGGCTATCGGCGCGGCAATCATTgcccagg                          <  1:8332312/97‑7 (MQ=255)
                                                                            tgaggAAGCTCCGTTCTACACCCTCGGCCCGCTGGCAACCGACATCGCGCCGGGCTATGACCACATCACGTCGGCTATCGGCGCGGCAATCATTgcgcagg                          >  1:553835/5‑95 (MQ=255)
                                                                            tgaggAAGCTCCGTTCTACACCCTCGGCCCGCTGGCAACCGACATCGCGCCGGGCTATGACCACATCACGTCGGCTATCGGCGCGGCAATCATTgcgcagg                          >  1:11026797/5‑95 (MQ=255)
                                                                            gcACGAAGAGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACAAATTCACGTCGGGGATTGGTGCGGCGGTGTATGGATGGt                          >  1:4012016/2‑101 (MQ=255)
                                                                            ccACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGt                          >  1:1517994/1‑101 (MQ=255)
                                                                             cACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGtt                         <  1:6453886/101‑1 (MQ=255)
                                                                             cACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGtt                         >  1:890446/1‑101 (MQ=255)
                                                                               ggAAGCTCCGTTCTACACCCTCGGCCCGCTGGCAACCGACATCGCGCCGGGCTATGACCACATCACGTCGGCTATCGGCGCGGCAATCATTgcgcaggctg                       >  1:7270025/2‑92 (MQ=255)
                                                                                gAAGCTCCGTTCTACGCCCTCGGCCCGCTGGCAGCCGTCATCGCGCCGGGCTATGACCACATCACGTCGGCTATCGGCGCGGCAATCATTgcgcaggctgg                      <  1:8505637/101‑11 (MQ=255)
                                                                                gAAGCTCCGTTCTACACCCTCGGCCCGCTGGCAACCGACATCGCGCCGGGCTATGACCACATCACGTCGGCTATCGGCGCGGCAATCATTgcgcaggctgg                      >  1:13708744/1‑91 (MQ=255)
                                                                                   gcgcCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGTGATGATTGGCTGGTTTGGCTg                   <  1:7649784/101‑1 (MQ=255)
                                                                                   gcgcCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTg                   <  1:12116639/101‑1 (MQ=255)
                                                                                    cgcCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTgc                  <  1:7117127/101‑1 (MQ=255)
                                                                                    cgcCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTgc                  >  1:6749588/1‑101 (MQ=255)
                                                                                      ccGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTgcgc                >  1:10194910/1‑101 (MQ=255)
                                                                                       cGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGcgcg               >  1:3307960/1‑101 (MQ=255)
                                                                                       cGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGcgcg               <  1:10712738/101‑1 (MQ=255)
                                                                                        gTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGGTTGGCTGCGCGa              >  1:505797/1‑101 (MQ=255)
                                                                                         ttttACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGAt             <  1:5173937/101‑1 (MQ=255)
                                                                                         ttttACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGGTTGGCTGCGCGAt             >  1:1358037/1‑101 (MQ=255)
                                                                                          tttACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATg            <  1:4583599/101‑1 (MQ=255)
                                                                                          tttACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATg            >  1:12870077/1‑101 (MQ=255)
                                                                                            tACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGct          <  1:232252/101‑1 (MQ=255)
                                                                                             acacTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGctc         <  1:5687325/101‑1 (MQ=255)
                                                                                             acacTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGctc         <  1:10032789/101‑1 (MQ=255)
                                                                                                atctGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGt      <  1:3686834/100‑1 (MQ=255)
                                                                                                    ggggCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACg  >  1:9255151/1‑101 (MQ=255)
                                                                                                    |                                                                                                    
AGGTGATGATTGAAGGCCCAGGCCACGTGCCGATGCAGATGATCCGCCGCAATATGACCGAGGAGTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACG  >  NC_007779/3441867‑3442067

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