Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0009133,875,5210T→C20.5% 10.5 ‑4.2 ‑0.040A28A (GCT→GCCyidXputative lipoprotein
Rejected: Polymorphism prediction not supported by minimum number of reads on both strands.
Rejected: Prediction has biased read strand distribution supporting polymorphism.
Reads supporting (aligned to +/- strand):  new base (8/0):  ref base (16/15):  total (24/16)
Fisher's exact test strand distribution p-value = 1.46e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.87e-01

AAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATATTCAGGTT  >  NC_000913/3875313‑3875728
                                                                                                                                                                                                                |                                                                                                                                                                                                               
aaaGAACAATGCATAAAAAAATGACATGCGTCGGGAAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGctc                                                                                                                                                                                                                <  2:373741/210‑1 (MQ=255)
  aGAACAAAGCATAAAAAAAAGAGATGCGTCGGGAAGAAATCTGAAAAGGGATATCAGGAGATAAACAGTAGGGCAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTg                                                                                                                                                                                                              <  2:282514/210‑1 (MQ=255)
        ttccgatctaaaaTGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGAGCTAAACAGGAGGGAAAGAGGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAAAGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTcc                                                                                                                                                                                                        <  1:167624/201‑1 (MQ=255)
              acaaaaaTGACATGCGTCGGGCAGTAAGCTGAAAACGGTTATCAGTCGCTAAACAGGAGGGCAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATGAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATAcc                                                                                                                                                                                                  <  2:445718/208‑1 (MQ=255)
               aaaaaaTGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCg                                                                                                                                                                                                 >  2:353731/1‑210 (MQ=255)
                 aaaaTGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTccagatcggaa                                                                                                                                                                                               >  2:167617/1‑201 (MQ=255)
                          tGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGt                                                                                                                                                                                      >  2:174921/1‑210 (MQ=255)
                                 ggCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTCGTCTCCCATACCGCAAAGAGTAGTTCAGGaa                                                                                                                                                                               >  1:411463/1‑210 (MQ=255)
                                   cAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAAGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAg                                                                                                                                                                             <  1:254662/210‑1 (MQ=255)
                                       aaTTTGAAAAGGGAAATGAGGAGATAAAAAGAAAGGAAAGAAGAATATGATTTCAACGACATGGCTAATCGTTTAAAGAATTAAGAGGGAAGGTGAATTTTAAGGGATTGTTTAAGGCTGCCGGTTTATTCCCACGGGAGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTTTc                                                                                                                                                                         <  1:70243/210‑1 (MQ=255)
                                       aaTCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATc                                                                                                                                                                         <  1:174928/210‑1 (MQ=255)
                                       aaTCAGAAAAGAAAAATCAGAAGCTAAACAAAAAGGAAAGAAGTGTGTCAGAACACCGGAATAGCAAATCGTTTAAAAGATTAATCATGAAGTTGAATTTTCAGTCATTTTTTTAGGCTGCCGGTTAATTCCCAATGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATc                                                                                                                                                                         <  1:312911/210‑1 (MQ=255)
                                          cTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAAt                                                                                                                                                                      >  1:247059/1‑210 (MQ=255)
                                          cTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACAGCAAAGGGTAGTTCAGGAAAGTATCAAt                                                                                                                                                                      >  1:202951/1‑210 (MQ=255)
                                                  ggATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCGTCGTTTCCCTGAAAAAAAAAAGTAGTTCAGGAAAGTATCAATCTCAGTCa                                                                                                                                                              >  1:461735/1‑210 (MQ=255)
                                                         aGGCGCGAAACAGGGAGGAAAGAAGAGTATGCTTTCAATGGCTTAGCTAATCGTTTAAAGCATTAATCATGAGGTTGAATTTTAAGGTATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACAcct                                                                                                                                                       <  2:159023/210‑2 (MQ=255)
                                                               tAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTg                                                                                                                                                 >  2:288723/1‑210 (MQ=255)
                                                                 aaCAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCGTCTTCTGCTTTAAAAAAAAAAGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCAGCACTggg                                                                                                                                               >  1:157838/1‑210 (MQ=255)
                                                                          gAAAGAAGAATATGCTTTCAACGGATTAAATAATCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAg                                                                                                                                      <  1:353754/210‑1 (MQ=255)
                                                                                      tGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCGTCTTCTGCTTTAAAAAAAAAAGTAGTTCAGGAAAGTATCAATCTCATTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGt                                                                                                                          >  1:105707/1‑210 (MQ=255)
                                                                                        cTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCGTCTTTTGCTTTACAAAAAAAAGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGAAGGCAAACGATAGTaa                                                                                                                        >  1:377225/1‑210 (MQ=255)
                                                                                                         aCTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATg                                                                                                       >  1:13156/1‑210 (MQ=255)
                                                                                                            cGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGttt                                                                                                    >  1:375087/1‑210 (MQ=255)
                                                                                                             gTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGATATGtttt                                                                                                   >  1:326009/1‑210 (MQ=255)
                                                                                                             gTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTCTGCCTTACAGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGTGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGtttt                                                                                                   >  1:221858/1‑210 (MQ=255)
                                                                                                                gaaaGGAGTAATCATGAAGTTGAATTTTAAGAGATTGTTTAAGTCTGCCGGTTTAGTCCAACGGGCGCTGATGCTTTCAGGGTGTATCTCGTATTCTATGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGt                                                                                                <  2:221468/209‑1 (MQ=255)
                                                                                                                    ggATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAggg                                                                                            >  2:31270/1‑210 (MQ=255)
                                                                                                                           tCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGc                                                                                     <  1:324382/210‑1 (MQ=255)
                                                                                                                            cATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCt                                                                                    >  2:101430/1‑210 (MQ=255)
                                                                                                                              tGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTgg                                                                                  <  2:325987/210‑1 (MQ=255)
                                                                                                                                   ttGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTAc                                                                             >  2:236137/1‑210 (MQ=255)
                                                                                                                                          ttAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGATATGTTTTTGTAGGGGAATCGCGGGATTACCTTATCa                                                                      >  1:15893/1‑210 (MQ=255)
                                                                                                                                             aGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTg                                                                   >  1:388992/1‑210 (MQ=255)
                                                                                                                                               ggATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAATAAGTAATGGAAAAAAAGGCTTTGTTTTTGTAGGGGAATGGCTGGATTACTTTATACCTGAt                                                                 >  1:233952/1‑210 (MQ=255)
                                                                                                                                                          ggCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCGTCTTCTGCTTGAAAAAAAAAAGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATAAAGGATGGTGCCGATGa                                                      >  1:26168/1‑210 (MQ=255)
                                                                                                                                                               ccGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCg                                                 <  2:114483/210‑1 (MQ=255)
                                                                                                                                                                      aTTACCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTATCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGAt                                          <  1:384014/210‑1 (MQ=255)
                                                                                                                                                                           cACTGGCACGTATGCTTTCCGGCTGTATCTCGTATGCTCTGGTGTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCa                                     <  2:375064/210‑1 (MQ=255)
                                                                                                                                                                                     gATGCTTTCAGGCTGTATCTCGTATGCCGTCTTCTGCTTGAAAAAAAAAAGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGATTTTGTAGGGGAATCGCTGGATTACCTTAAAACTGATGGTGCTGATGACATCGTGAAGATGCTCAATTATACAga                           >  1:451231/1‑209 (MQ=255)
                                                                                                                                                                                                           tATGCACTGGTTTCCCATACCGAAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGAAACCATCACTGGGCTATCGAAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATATTCAg     <  2:55939/210‑1 (MQ=255)
                                                                                                                                                                                                              gCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATATTCAGGtt  >  1:147602/1‑210 (MQ=255)
                                                                                                                                                                                                                |                                                                                                                                                                                                               
AAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATATTCAGGTT  >  NC_000913/3875313‑3875728

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