Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 2,251,585 A→G 2.6% F94F (TTT→TTC lysP ← lysine transporter

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0077792,251,5850A→G 2.6% 6.7 ‑0.8 ‑0.0246F94F (TTT→TTClysPlysine transporter
Reads supporting (aligned to +/- strand):  new base (2/5):  ref base (107/132):  total (109/137)
Fisher's exact test strand distribution p-value = 4.68e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.93e-01

Only 100 of 254 total aligned reads displayed.

ATCCGGGAACCACCAGCTCATGACCAGCTGAGCTGCAACCAGGTCAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCA  >  NC_007779/2251486‑2251685
                                                                                                   |                                                                                                    
 tCCGGGAACCACCAGCTCATGACCAGCTGAGCTGCAACCAGGTCAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCaaa                                                                                                    <  1:13342521/101‑1 (MQ=255)
 tCCGGGAACCACCAGCTCATGACCAGCTGAGCTGCAACCAGGTCAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCaaa                                                                                                    <  1:14198286/101‑1 (MQ=255)
 tCCGGGAACCACCAGCTCATGACCAGCTGAGCTGCAACCAGGTCAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCaaa                                                                                                    >  1:8705943/1‑101 (MQ=255)
  ccGGGAACCAATAGCGCATCACCAGCGCAGCAGCGACGAGCTCTGCGGCCACAGTAATCGCCCAGTTAAACCAGTAGTTCCAGCCCATGGCAAAGCCGAAg                                                                                                   >  1:5364485/1‑101 (MQ=255)
   cGGGAACCACCAGCTCATGACCAGCTGAGCTGCAACCAGGTCAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGc                                                                                                  <  1:2478342/101‑1 (MQ=255)
      gAACCACCAGCTCATGACCAGCTGAGCTGCAACCAGGTCAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCctt                                                                                               >  1:5396084/1‑101 (MQ=255)
        accaccAGCTCATGACCAGCTGAGCTGCAACCAGGTCAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCcttct                                                                                             <  1:5163181/101‑1 (MQ=255)
         ccaccaGCTCATGACCAGCTGAGCTGCAACCAGGTCAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCcttctt                                                                                            >  1:657325/1‑101 (MQ=255)
           accaGCTCATGACCAGCTGAGCTGCAACCAGGTCAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCa                                                                                          >  1:6508790/1‑101 (MQ=255)
            ccaGCTCATGACCAGCTGAGCTGCAACCAGGTCAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCaa                                                                                         <  1:8800775/101‑1 (MQ=255)
              aGCTCATGACCAGCTGAGCTGCAACCAGGTCAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACa                                                                                       >  1:12909182/1‑101 (MQ=255)
                cTCATGACCAGCTGAGCTGCAACCAGGTCAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATa                                                                                     >  1:3846665/1‑101 (MQ=255)
                  cATCACCAGCGCAGCAGCGACGAGCTCTGCGGCCACAGTAATCGCCCAGTTAAACCAGTAGTTCCAGCCCATGGCAAAGCCGAAGGattctgagacgtatc                                                                                   <  1:9716745/101‑11 (MQ=255)
                   aTGACCAGCTGAGCTGCAACCAGGTCAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGtt                                                                                  >  1:1090930/1‑101 (MQ=255)
                    tGACCAGCTGAGCTGCAACCAGGTCAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTc                                                                                 <  1:2125502/101‑1 (MQ=255)
                    tGACCAGCTGAGCTGCAACCAGGTCAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTc                                                                                 <  1:12015757/101‑1 (MQ=255)
                    tGACCAGCTGAGCTGCAACCAGGTCAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTc                                                                                 <  1:14792309/101‑1 (MQ=255)
                     gACCAGCTGAGCTGCAACCAGGTCAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCt                                                                                >  1:1748992/1‑101 (MQ=255)
                     gACCAGCTGAGCTGCAACCAGGTCAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCt                                                                                >  1:10968708/1‑101 (MQ=255)
                     gACCAGCTGAGCTGCAACCAGGTCAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCt                                                                                <  1:8533692/101‑1 (MQ=255)
                      aCCAGCTGAGCTGCAACCAGGTCAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTg                                                                               <  1:9143156/101‑1 (MQ=255)
                      aCCAGCTGAGCTGCAACCAGGTCAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTg                                                                               <  1:15203019/101‑1 (MQ=255)
                       ccAGCTGAGCTGCAACCAGGTCAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGa                                                                              <  1:5681838/101‑1 (MQ=255)
                        cAGCTGAGCTGCAACCAGGTCAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGAc                                                                             >  1:14604466/1‑101 (MQ=255)
                         agctgagctgCAACCAGGTCAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCCGAcc                                                                            <  1:1796842/101‑1 (MQ=255)
                         agctaagctgCAACCAGGTCAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGAcc                                                                            >  1:12734037/1‑101 (MQ=255)
                           ctgagctgCAACCAGGTCAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGt                                                                          >  1:11834521/1‑101 (MQ=255)
                             gagctgCAACCAGGTCAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTaa                                                                        <  1:13235281/101‑1 (MQ=255)
                                ctgCAACCAGGTCAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTg                                                                     >  1:14109415/1‑101 (MQ=255)
                                   cAACCAGGTCAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCa                                                                  <  1:8862013/101‑1 (MQ=255)
                                    aaCCAGGTCAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCaa                                                                 >  1:1887823/1‑101 (MQ=255)
                                    aaCCAGGTCAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCaa                                                                 <  1:14704315/101‑1 (MQ=255)
                                     aCCAGGTCAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCaaa                                                                >  1:5906650/1‑101 (MQ=255)
                                     aCCAGGTCAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCaaa                                                                <  1:6627479/101‑1 (MQ=255)
                                       cAGGTCAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACg                                                              >  1:13762774/1‑101 (MQ=255)
                                       cAGGTCAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACg                                                              >  1:3482212/1‑101 (MQ=255)
                                         ggTCAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGaa                                                            <  1:7710626/101‑1 (MQ=255)
                                           tCAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAAcc                                                          <  1:7441595/101‑1 (MQ=255)
                                           tCAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAAcc                                                          <  1:11564040/101‑1 (MQ=255)
                                           tCAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGTCCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAAcc                                                          <  1:10614484/101‑1 (MQ=255)
                                             aaCGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCgg                                                        <  1:10251642/101‑1 (MQ=255)
                                             aaCGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCgg                                                        <  1:1605660/101‑1 (MQ=255)
                                             aaCGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCgg                                                        <  1:4868091/101‑1 (MQ=255)
                                             aaCGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCgg                                                        <  1:7902384/101‑1 (MQ=255)
                                              aCGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGa                                                       <  1:6305184/101‑1 (MQ=255)
                                                ggCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGaaa                                                     <  1:8028525/101‑1 (MQ=255)
                                                ggCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGaaa                                                     >  1:10037244/1‑101 (MQ=255)
                                                 gCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAAc                                                    >  1:8116367/1‑101 (MQ=255)
                                                 gCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAAc                                                    <  1:3076133/101‑1 (MQ=255)
                                                  cGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAAcc                                                   <  1:13982889/101‑1 (MQ=255)
                                                  cGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAAcc                                                   >  1:12320127/1‑101 (MQ=255)
                                                  cGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAAcc                                                   <  1:9244952/101‑1 (MQ=255)
                                                  cGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAAcc                                                   <  1:14985938/101‑1 (MQ=255)
                                                    aTAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCgg                                                 >  1:15128976/1‑101 (MQ=255)
                                                     tAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGc                                                <  1:7050879/101‑1 (MQ=255)
                                                      aGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCa                                               >  1:112459/1‑101 (MQ=255)
                                                       gTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCat                                              >  1:574700/1‑101 (MQ=255)
                                                        tCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCata                                             <  1:472818/101‑1 (MQ=255)
                                                         cACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCatat                                            >  1:3346708/1‑101 (MQ=255)
                                                          aCCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCatata                                           >  1:7794715/1‑101 (MQ=255)
                                                          aCCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCatata                                           <  1:1032645/101‑1 (MQ=255)
                                                           ccGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCatatat                                          >  1:5272380/1‑101 (MQ=255)
                                                            cGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATg                                         >  1:9870084/1‑101 (MQ=255)
                                                             gCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGc                                        >  1:4809235/1‑101 (MQ=255)
                                                             gCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGc                                        <  1:4500929/101‑1 (MQ=255)
                                                              cccAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCa                                       >  1:3592528/1‑101 (MQ=255)
                                                               ccAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAg                                      <  1:11378291/101‑1 (MQ=255)
                                                                 aGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGcc                                    <  1:8678457/101‑1 (MQ=255)
                                                                  gTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCa                                   >  1:11105912/1‑101 (MQ=255)
                                                                  gTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCa                                   >  1:12108485/1‑101 (MQ=255)
                                                                    tGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGt                                 >  1:13865567/1‑101 (MQ=255)
                                                                    tGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGt                                 <  1:3142647/101‑1 (MQ=255)
                                                                      tACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTc                               <  1:7150489/101‑1 (MQ=255)
                                                                       aCCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCa                              >  1:14915041/1‑101 (MQ=255)
                                                                        ccAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCACCATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCAc                             <  1:12879613/101‑1 (MQ=255)
                                                                        ccAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCAc                             <  1:4094917/101‑1 (MQ=255)
                                                                        ccAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCACCCAGTTCAc                             <  1:10116798/101‑1 (MQ=255)
                                                                         cAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCAcc                            >  1:10592541/1‑101 (MQ=255)
                                                                           gtagtTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCga                          >  1:11054983/1‑101 (MQ=255)
                                                                            tagtTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCgag                         >  1:14008624/1‑101 (MQ=255)
                                                                             agtTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCgaga                        >  1:9357182/1‑101 (MQ=255)
                                                                                 ccagcccagcGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTgg                    >  1:4283581/1‑101 (MQ=255)
                                                                                     cccagcGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGtcat                >  1:139190/1‑101 (MQ=255)
                                                                                       cagcGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGtcatca              <  1:682279/101‑1 (MQ=255)
                                                                                        agcGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAg             <  1:7613806/101‑1 (MQ=255)
                                                                                        agcGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAg             >  1:11905891/1‑101 (MQ=255)
                                                                                         gcgcGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAgg            <  1:13528873/101‑1 (MQ=255)
                                                                                          cgcgAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGa           <  1:2221075/101‑1 (MQ=255)
                                                                                          cgcgAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGa           <  1:2938405/101‑1 (MQ=255)
                                                                                           gcgAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGaa          >  1:5584590/1‑101 (MQ=255)
                                                                                            cgAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAg         >  1:6258381/1‑101 (MQ=255)
                                                                                            cgAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAg         >  1:2564260/1‑101 (MQ=255)
                                                                                              aaGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTa       >  1:3287985/1‑101 (MQ=255)
                                                                                              aaGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTa       <  1:9573962/101‑1 (MQ=255)
                                                                                               aGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTaa      <  1:1841831/101‑1 (MQ=255)
                                                                                                 ccAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAAc    <  1:7244370/101‑1 (MQ=255)
                                                                                                 ccAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAAc    <  1:9048968/101‑1 (MQ=255)
                                                                                                 ccAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGCCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAAc    >  1:7850327/1‑101 (MQ=255)
                                                                                                  cAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAAcc   <  1:12578015/101‑1 (MQ=255)
                                                                                                   aaaGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCa  >  1:14499352/1‑101 (MQ=255)
                                                                                                   |                                                                                                    
ATCCGGGAACCACCAGCTCATGACCAGCTGAGCTGCAACCAGGTCAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCA  >  NC_007779/2251486‑2251685

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