Read alignment evidence...
	seq id	position		change	freq	score	reads	annotation	genes	product
*	NC_000913	3,508,894	0	G→A	10.5%	11.6 ‑4.9 0.0	96	A307V (GCA→GTA)	yhfW	phosphopentomutase‑related metalloenzyme superfamily protein
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 (9/1):  ref base (38/47):  total (48/48)
Fisher's exact test strand distribution p-value = 7.63e-03
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

GTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCATCACGACCAGGCAATCATCTGGCTGCATCGCCTCAACAAGCCGGGCAAGGTTACGGTCAACGACCTGCAAACGTTCGGCATAACGTGCGACGTCTTCTGCATGACCAGCGAGGTCGGTTTCCTGAATGTTGGTGCAAATAAACGCCGTCGGATGGGTGTTAAATTCGTTGAGGGTGATATCCATAATCCGCTGGCTATCCACCAGATTTTGCCAGCTCACGCCATAAGGATTGTTGACGATATCTGCCACCTTACCCACCAGCACGGTTGGCACGC  >  NC_000913/3508696‑3509083                                                                                                                                                                                                       |                                                                                                                                                                                              gtgCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCATCACGACCAGGCAATCATCTGGCTGCATCGCCTCAACAAGCCGGGCAAGGTTACGGTCAACGACCTGCAAACGTTCGGCATAACGTGc                                                                                                                                                                                              >  1:528378/1‑200 (MQ=255)            gAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCATCACGACCAGGCAATCATCTGGCTGCATCGCCTCAACAAGCCGGGCAAGGTTACGGTCAACGACCTGCAAACGTTCGGCATAACGTGCGACGTCTTCTg                                                                                                                                                                                   >  1:551323/1‑200 (MQ=255)                 gCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCATCACGACCAGGCAATCATCTGGCTGCATCGCCTCAACAAGCCGGGCAAGGTTACGGTCAACGACCTGCAAACGTTCGGCATAACGTGCGACGTCTTCTGCATGa                                                                                                                                                                              >  1:12426/1‑200 (MQ=255)                 gCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCATCACGACCAGGCAATCATCTGGCTGCATCGCCTCAACAAGCCGGGCAAGGTTACGGTCAACGACCTGCAAACGTTCGGCATAACGTGCGACGTCTTCTGCATGa                                                                                                                                                                              >  1:696489/1‑200 (MQ=255)                                 ttGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCATCACGACCAGGCAATCATCTGGCTGCATCGCCTCAACAAGCCGGGCAAGGTTACGGTCAACGACCTGCAAACGTTCGGCATAACGTGCGACGTCTTCTGCATGACCAGCGAGGTCGGttt                                                                                                                                                              <  1:267986/200‑1 (MQ=255)                                                    ggCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCATCACGACCAGGCAATCATCTGGCTGCATCGCCTCAACAAGCCGGGCAAGGTTACGGTCAACGACCTGCAAACGTTCGGCATAACGTGCGACGTCTTCTGCATGACCAGCGAGGTCGGTTTCCTGAATGTTGGTGCAAAt                                                                                                                                           >  1:119770/1‑200 (MQ=255)                                                           cTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCATCACGACCAGGCAATCATCTGGCTGCATCGCCTCAACAAGCCGGGCAAGGTTACGGTCAACGACCTGCAAACGTTCGGCATAACGTGCGACGtctt                                                                                                                                                                                      >  2:346689/1‑150 (MQ=255)                                                               gcgGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCATCACGACCAGGCAATCATCTGGCTGCATCGCCTCAACAAGCCGGGCAAGGTTACGGTCAACGACCTGCAAACGTTCGGCATAACGTGCGACGTCTTCTGc                                                                                                                                                                                  <  2:83412/150‑1 (MQ=255)                                                                  ggTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCATCACGACCAGGCAATCATCTGGCTGCATCGCCTCAACAAGCCGGGCAAGGTTACGGTCAACGACCTGCAAACGTTCGTCATAACGTGCGACGTCTTCTGCATg                                                                                                                                                                               >  2:266026/1‑150 (MQ=255)                                                                  ggTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCATCACGACCAGGCAATCATCTGGCTGCATCGCCTCAACAAGCCGGGCAAGGTTACGGTCAACGACCTGCAAACGTTCGGCATAACGTGCGACGtc                                                                                                                                                                                        >  4:401110/1‑141 (MQ=255)                                                                  ggTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCATCACGACCAGGCAATCATCTGGCTGCATCGCCTCAACAAGCCGGGCAAGGTTACGGTCAACGACCTGCAAACGTTCGGCATAACGTGCGACGTCTTCTGCATGACCAGCGAGGTCGGTTTCCTGAATGTTGGTGCAAATAAACGCCGTCGGAt                                                                                                                             >  1:486462/1‑200 (MQ=255)                                                                    tATGGTGACTGAGACCAAGGGGCGGATCGTGGCCATGAGCCACCGTCACAACCAGGCAAACATCTGGCGGCATCACCGCAACAAGCCGGGCAAGGTTACGGTCAACGACCTGAAAACGTTCGGCATAACGTGCGACGTCTTCTGAAAGAc                                                                                                                                                                             <  2:343708/150‑1 (MQ=255)                                                                      tggtggCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCATCACGACCAGGCAATCATCTGGCTGCATCGCCTCAACAAGCCGGGCAAGGTTACGGTCAACGACCTGCAAACGTTCGGCATCTCGTATGCCGTCTTCTGCTTGAAAAAAAAAATCGGATTCCTGAATGTTGGTTACACTAAACGACGTCGTCTGGgt                                                                                                                         >  1:737277/1‑200 (MQ=255)                                                                      tggtggCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCATCACGACCAGGCAATCATCTGGCTGCATCGCCTCAACAAGCCGGGCAAGGTTACGGTCAACGACCTGCAAACGTTCGGCATCACGTGTGCCGTCTTCTGCATGACCACCGAGGTCGGTTTCCTGAATGTTGGTGCAAATAAACGCCGTCGGATGGgt                                                                                                                         >  1:736995/1‑200 (MQ=255)                                                                      tggtggCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCATCACGACCAGGCAATCATCTGGCTGCATCGCCTCAACAAGCCGGGCAAGGTTACGGTCAACGACCTGCAAACGTTCGGCATAACGTGCGACGTCTGCTGCATGACCa                                                                                                                                                                           >  2:197384/1‑150 (MQ=255)                                                                       ggtggCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCATCACGACCAGGCAATCATCTGGCTGCATCGCCTCAACAAGCCGGGCAAGGTTACGGTCAACGACCTGCAAACGTTCGGCATAACGTGCGACGTCTTCTGCATGACCAGCGAGGTCGGTTTCCTGAATGTTGGTGCAAATAAACGCCGTCGGATGGgtg                                                                                                                        >  1:415901/1‑200 (MQ=255)                                                                           gCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCATCACGACCAGGCAATCATCTGGCTGCATCGCCTCAACAAGCCGGGCAAGGTTACGGTCAACGACCTGCAAACGTTCGGCATCTCGTATGACGTCTTCTGCATGACAAAAAAAATCGGTTTCCTGAATGTTGGTGCAAATAAACGCCGTCGGATGGGTGTTaa                                                                                                                    >  1:813762/1‑200 (MQ=255)                                                                           gCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCATCACGACCAGGCAATCATCTGGCTGCATCGCCTCAACAAGCCGGGCAAGGTTACGGTCAACGACCTGCAAACGTTCGGCATAACGTGCGACGTCTTCTGCATGACCAGCGAGGTCGGTTTCCTGAATGTTGGTGCAAATAAACGCCGTCGGATGGGTGTTaa                                                                                                                    >  1:573707/1‑200 (MQ=255)                                                                             tgtgACCAATGGTCGGATCGTTGCCGTGATCCGCCATCACGACCAGGCAATCATCTGGCTGCATCGCCTCAACAAGCCGGGCAAGGTTACGGTCAACGACCTGCAAACGTTCGGCATAACGTGCGACGTCTTCTGCATGACCAGCGAGGTCGGTTTCCTGAATGTTGGTGCAAATAAACGCCGTCGGATGGGTGTTAAAt                                                                                                                  <  1:346718/200‑1 (MQ=255)                                                                             tgtgACCAATGGTCGGATCGTTGCCGTGATCCGCCATCACGACCAGGCAATCATCTGGCTGCATCGCCTCAACAAGCCGGGCAAGGTTACGGTCAACGACCTGCAAACGTTCGGCATAACGTGCGACGTCTTCTGCATGACCAGCGAGGTCGGTTTCCTGAATGTTGGTGCAAATAAACGCCGTCGGATGGGTGTTAAAt                                                                                                                  <  1:811568/200‑1 (MQ=255)                                                                              gtgACCAATGGTCGGATCGTTGCCGTGATCCGCCATCACGACCAGGCAATCATCTGGCTGCATCGCCTCAACAAGCCGGGCAAGGTTACGGTCAACGACCTGCAAACGTTCGGCATAACGTGCGACGTCTTCTGCATGACCAGCGAGGTCGGTTTCCTGAATGTTGGTGCAAATAAACGCCGTCGGATGGGTGTTAAAtt                                                                                                                 >  1:595722/1‑200 (MQ=255)                                                                                gACCAATGGTCGGATCGTTGCCGTGATCCGCCATCACGACCAGGCAATCATCTGGCTGCATCGCCTCAACAATCCGGGCAAGGTTACGGTCAACGACCTGCAAACGTTCGGCATAACGTGCGACGTCTTCTGCATGCCCACCGAGGTCgg                                                                                                                                                                 >  2:623502/1‑150 (MQ=255)                                                                                   cAATGGTCGGATCGTTGCCGTGATCCGACATCACGACCAGGCAATCATCTGGCTGCATCGCCTCAACAAGCCGGGATAGGTTACAGTCAACGACCTGCAAACGTTCGGCATAACGTGCGACGTCTTCTACATGACCAGCGAGGTCGGttt                                                                                                                                                              <  2:610090/150‑1 (MQ=255)                                                                                     aTGGTCGGATCGTTGCCGTGATCCGCCATCACGACCAGGCAATCATCTGGCTGCATCGCCTCAACAAGCCGGGCAAGGTTACGGTCAACGACCTGCAAACGTTCGGCATAACGTGCGACGTCTTCTGCATGa                                                                                                                                                                              >  3:232325/1‑132 (MQ=255)                                                                                     aTGGTCGGATCGTTGCCGTGATCCGCCATCACGACCAGGCAATCATCTGGCTGCATCGCCTCAACAAGCCGGGCAAGGTTACGGTCAACGACCTGCAAACGTTCGGCATAACGTATGACGTCTTCTGCATGAAAAACAAAGTAGGTTTCCTGAATGTTGGTGCAAATAAACGCCGTCGGATGGGGGTTAAATTCGTGGAg                                                                                                          >  1:801146/1‑200 (MQ=255)                                                                                      tGGTCGGATCGTTGCCGTGATCCGCCATCACGACCAGGCAATCATCTGGCTGCATCGCCTCAACAAGCCGGGCAAGGTTACGGTCAACGACCTGCAAACGTTCGGCATAACGTGCGACGTCTTCTGCATGAc                                                                                                                                                                             >  3:275894/1‑132 (MQ=255)                                                                                      tGGTCGGATCGTTGCCGTGATCCGCCATCACGACCAGGCAATCATCTGGCTGCATCGCCTCAACAAGCCGGGCAAGGTTACGGTCAACGACCTGCAAACGTTCGGCATAACGTGCGACGTCTTCTGCATGAc                                                                                                                                                                             <  3:401743/132‑1 (MQ=255)                                                                                       ggTCGGATCGTTGCCGTGATCCGCCATCACGACCAGGCAATCATCTGGCTGCATCGCCTCAACAAGCCGGGCAAGGTTACGGTCAACGACCTGCAAACGTTCGGCATAACGTGTGACGTCTTCTGCATGACCAGCGAGGTc                                                                                                                                                                   >  4:557433/1‑141 (MQ=255)                                                                                        gTCGGATCGTTGCCGTGATCCGCCATCACGACCAGGCAATCATCTGGCTGCATCGCCTCAACAAGCCGGGCAAGGTTACGGTCAACGACCTGCAAACGTTCGGCATAACGTGCGACGTCTTCTGCATGACCAGCGAGGTCGGTTTCCTGa                                                                                                                                                         <  2:536206/150‑1 (MQ=255)                                                                                           ggATCGTTGCCGTGATCCGCCATCACGACCAGGCAATCATCTGGCTGCATCGCCTCAACAAGCCGGGCAAGGTTACGGTCAACGACCTGCAAACGTTCGGCATAACGTGCGACGTCTTCTGCATGACCAGCg                                                                                                                                                                        <  3:23714/132‑1 (MQ=255)                                                                                            gATCGTTGCCGTGATCCGCCATCACGACCAGGCAATCATCTGGCTGCATCGCCTCAACAAGCCGGGCAAGGTTACGGTCAACGACCTGCAAACGTTCGGCATATCGTATGACGTCTTCTGCATGACCAGCAAGGTCGGttt                                                                                                                                                              >  4:240298/1‑141 (MQ=255)                                                                                              tCGTTGCCGTGATCCGCCATCACGACCAGGCAATCATCTGGCTGCATCGCCTCAACAAGCCGGGCAAGGTTACGGTCAACGACCTGCAAACGTTCGGCATAACGGGCGACGGTTTATGCATGACAAGCGTgg                                                                                                                                                                     >  3:284204/1‑132 (MQ=255)                                                                                                    gcGTGACCCGCCATCACGACCAGGCAATCATCTGGCTGCATCGCCTCAACAAGCCGGGCAAGGTTACGGTCAACGACCTGCAAACGTTCGGCATAACGTGCGACGTCTTCTGCATGACCAGCGAGGTCGGTTTCCTGAATg                                                                                                                                                      <  4:572969/140‑1 (MQ=255)                                                                                                     cGTGATCCGCCATCACGACCAGGCAATCATCTGGCTGCATCGCCTCAACAAGCCGGGCAAGGTTACGGTCAACGACCTGCAAACGTTCGGCATAACGTGCGACGTCTTCTGCATGACCAGCGAGGTCGGTTTCCTGAATGTTGGTGCAAATAAACGCCGTCGGATGGGTGTTAAATTCGTTGAGGGTGATATCCATAATc                                                                                          <  1:595441/200‑1 (MQ=255)                                                                                                      gTGATCCGCCATCACGACCAGGCAATCATCTGGCTGCATCGCCTCAACAAGCCGGGCAAGGTTACGGTCAACGACCTGCAAACGTTCGGCATCTCGTATGCCGTCTTCTGCTTGAAAAAAAAAATCGGTTTACTGAATGTTGGTGCAAATAAACGCCGTCGGATGGGTGTTAAATTCGTTGAGGGTGATATCCATAATcc                                                                                         >  1:818870/1‑200 (MQ=255)                                                                                                       tGATCCGCCATCACGACCAGGCAATCATCTGGCTGCATCGCCTCAACAAGCCGGGCAAGGTTACGGTCAACGACCTGCAAACGTTCGGCATAACGTGCGACGTCTTCTGCATGACCAGCGAGGTCGGTTCCCTGAATGTTGGTGCAAATa                                                                                                                                          >  2:584626/1‑150 (MQ=255)                                                                                                          tCCGCCATCACGACCAGGCAATCATCTGGATGCAGCGCCTCAACAAGCCGGGCAAGGTTACGGTCAACGACCTGCAAACGTTCGGCATAACGTGCGACGTCTTCTGCATGACCAGCGAGGTCGGTTTCCTGa                                                                                                                                                         <  3:566358/132‑1 (MQ=255)                                                                                                           ccGCCATCACGACCAGGCAATCATCTGGCTGCATCGCCTCAACAAGCCGGGCAAGGTTACGGTCAACGACCTGCAAACGTTCGGCATAACGTGCGACGTCTTCTGCATGACCAGCGAGGTCGGTTTCCTGAATGTTGGTGc                                                                                                                                               <  4:245315/141‑1 (MQ=255)                                                                                                             gCCATCACGACCAGGGAATCATCTGGCTGCATCGCCTGAACAAGCCGGGCAAGGTTACGGTCAACGACCTGCAAACGTGCGGCATAACGTGCGACGTCTTCTGCATGACCAGCGAGGTCGGTTTCCTGAATGTTGGTGCAAATAAACGCCGTCGGATGGGTGTTAAATTCGTTGAGGGTGATATCCATAATCCGCTGGCt                                                                                  <  1:181311/200‑1 (MQ=255)                                                                                                             gCCATCACGACCAGGCAATCATCTGGCTGCATCGCCTCAACAAGCCGGGCAAGGTTACGGTCAACGACCTGCAAACGTTCGGCATAACGTGCGACGTCTTCTGCATGACCAGCGAGGTCGGTTTCCTGAATGTTGGTGCAAATAAACGcc                                                                                                                                    >  2:472091/1‑150 (MQ=255)                                                                                                                aTCACGACCAGGCAATCATCTGGCTGCATCGCCTCAACAAGCCGGGCAAGGTTACGGTCAACGACCTGCAAACGTTCGGCATAACGTGCGACGTCTTCTGCATGACCAGCGAGGTCGGTTTCCTGAATGTTGGTGCAAATa                                                                                                                                          >  4:436107/1‑141 (MQ=255)                                                                                                                aTCACGACCAGGCAATCATCTGGCTGCATCGCCTCAACAAGCCGGGCAAGGTTACGGTCAACGACCTGCAAACGTTCGGCATAACGTGCGACGTCTTCTGCATGACCAGCGAGGTCGGGTTCCTGAATGTTGGTGCAAATAAACGCCGTc                                                                                                                                 >  2:585561/1‑150 (MQ=255)                                                                                                                    cGACCAGGCAATCATCTGGCTGCATCGCCTCAACAAGCCGGGCAAGGTTACGGTCAACGACCTGCAAACGTTCGGCATAACGTGCGACGTCTTCTGCATGACCAGCGAGGTCGGTTTCCTGAATGTTGGTGCAAATAAACGCCGTCGGAt                                                                                                                             <  2:317601/150‑1 (MQ=255)                                                                                                                      aCCAGGCAATCATCTGGCTGCATCGCCTCAACAAGCCGGGCAAGGTTACGGTCAACGACCTGCAAACGTTCGGCATAACGTATGCCGTCTTCTGCATGACAAAAAAAATCAGTTTCCTAAATCTTGGTGCAAATAAACGCCATCGGATGGGTGTTAAATTCGTTGAGGGTGATATCCATAATCCGCTGGCTATCCACCAg                                                                         >  1:400776/1‑200 (MQ=255)                                                                                                                         aGGCAATCATCTGGCTGCATCGCCTCAACAAGCCGGGCAAGGTTACGGTCAACGACCTGCAAACGTTCGGCATCACGTATGCCGTCTTCTGCTTGAACAACAAGATCGGTTTCCTGAATGTTGGTGCAAATAAACGCCGTCGGATGGGTGTTAAATTCGTTGAGGGTGATATCCATAATCCGCTGGCTATCCACCAGAtt                                                                      >  1:392627/1‑200 (MQ=255)                                                                                                                         aGGCAATCATCTGGCTGCATCGCCTCAACAAGCCGGGCAAGGTTACGGTCAACGACCTGCAAACGTTCGGCATAACGTGCGACGTCTTCTGCATGACCAGCGAGGTCGGTTTCCTGAATGTTGGTGCAAATAAACGCCGTCGGATGGgtg                                                                                                                        >  2:226445/1‑150 (MQ=255)                                                                                                                              atcagcTGGATGCATCGCCTCAACAAGCCGGGCAAGGTTACGGTCAACGACCTGCAAACGTTCGGCATAACGTGCGACGTCTTCTGCATGACCAGCGAGGTCGGTTTCCTGAATGTTGGTGCAAATAAACGCCGTCGGATGGGTGTTAAATTCGTTGAGGGTGATATCCATAATCCGCTGGCTATCCACCAGATTTTGcc                                                                 <  1:214951/200‑1 (MQ=255)                                                                                                                               tcatcTGGCTGCATCGCCTCAACAAGCCGGGCAAGGTTACGGTCAACGACCTGCAAACGTTCGGCATAACGTGCGACGTCTTCTGCATGACCAGCGAGGTCGGTTTCCTGAATGTTGGTGCAAATAAACGCCGTCGGATGGGTGTTAAATTCGTTGAGGGTGATATCCATAATCCGCTGGCTATCCACCAGATTTTGCCa                                                                >  1:441427/1‑200 (MQ=255)                                                                                                                               tcataTGGCTGCATCGCCGCAACAAGCCGGGCAAGGTAACGGTCAACGCCCTGCAAACGTTCGGCATAACGGGAGGCGTCTTCTGCATGACCAGCGAGGTCGGTTTCCGGAATGTTGGTGCTAATTAACGCCGTCGGATGGGTGTTAAAt                                                                                                                  <  2:119213/150‑1 (MQ=255)                                                                                                                                catcTGGCTGCATCGCCTCAACAAGCCGGGCAAGGTTACGGTCAACGACCTGCAAACGTTCGGCATAACGTCCGTCGTCTTCTGCATGACCAGAGGGGTTGGTTTCCAGAATGTTGGCGCAAAGAAACGCCGTCGGTTggg                                                                                                                          >  4:438125/1‑141 (MQ=255)                                                                                                                                    tGGCTGCATCGCCTCAACAAGCCGGGCAAGGTTACGGTCAACGACCTGCAAACGTTCGGCATAACGTGCGACGTCTTCTGCATGACCAGCGAGGTCGGTGTCCTGAATGTTGGTGCAAATAAACGCCGTCGGATGGGTGTTAAATTCGtt                                                                                                             >  2:712503/1‑150 (MQ=255)                                                                                                                                     ggCTGCATCGCCTCAACAAGCCGGGCAAGGTTACGGTCAACGACCTGCAAACGTTCGGCATAACGTGCGACGTCTTCTGCATGACCAGCGAGGTCGGTTTCCTGAATGTTGGTGCAAATAAACGCCGTCGGATGGGTGTTAAATTCGTTGAGGGTGATATCCATAATCCGCTGGCTATCCACCAGATTTTGCCAGCTCAc                                                          >  1:72639/1‑200 (MQ=255)                                                                                                                                       cTGCATCGCCTCAACAAGCCGGGCAAGGTTACGGTCAACGACCTGCAAACGTTCGGCATAACGTGCGACGTCTTCTGCATGACCAGCGAGGTCGGTTTCCTGAATGTTGGTGCAAATAAACGCCGTCGGATGGGTGTTaaa                                                                                                                   <  4:494800/141‑1 (MQ=255)                                                                                                                                        tGCATCGCCTCAACAAGCCGGGCAAGGTTACGGTCAACGACCTGCAAACGTTCGGCATAACGTGCGACGTCTTCTGCATGACCAGCGAGGTCGGTTTCCTGAATGTTGGTGCAAATAAACGCCGTCGGATGGGTGTTAAAt                                                                                                                  <  4:275905/141‑1 (MQ=255)                                                                                                                                         gCATCGCCTCAACGAGCCGGGCAAGGTTACGGTCAACGACCTGCAATCGTTCGGCATAACGTGCGACGTCTTCTGCATGACCAGCGAGGTCGGTTTCCTGAATGTTGGTGCAAATAAACGCCGTCGGATGGGTGTTAAAtt                                                                                                                 <  4:447710/141‑1 (MQ=255)                                                                                                                                           atggagtCAACAGGCCGGGGAAGGTTACGGGCAACGACCTGCGAAGGTTCGGCGTAACGGGCGACGTCTTCTGGATGACCAGCGAGGTCGGTTTCCTGAATGTTGGTGCAAATAAACGCCGTCGGATGGgtg                                                                                                                        <  3:613967/126‑1 (MQ=255)                                                                                                                                             cGCCTCAACAAGCCGGGCAAGGTTACGGTCAACGACCTGCAAACGTTCGGCATAACGTGCGACGTCTTCTGCATGACCAGCGAGGTCGGTTTCCTGAATGTTGGTGCAAATAAACGCCGTCGGATGGGTGTTAAATTCGTTGAGGGTGATATCCATAATCCGCTGGCTATCCACCAGATTTTGCCAGCTCACGCCATAAg                                                  <  1:142328/200‑1 (MQ=255)                                                                                                                                              gCCTCAACAAGCCGGGCAAGGTTACGGTCAACGACCTGCAAACGTTCGGCATATCGTACGCCGTCTTCTGCATGACCAACAAAGTCGGTTTCCTGAATGTTGGTGCAAATAAACGCCGTCGGATGGGTGTTAAATTCGTTGTGGGTGATATCCATAATCCGCTGTCTATCCACCAGATTTTGCCAGCTCACGCCATAAgg                                                 >  1:383203/1‑200 (MQ=255)                                                                                                                                              gCCTCAACAAGCCGGGCAAGGTTACGGTCAACGACCTGCAAACGTTCGGCATAACGTGCGACGTCTTCTGCATGACCAGCGAGGTCGGTTTCCTGAATGTTGGTGCAAATAAACGCCGTCGGATGGGTGTTAAATTCGTTGAGGGTGATATCCATAATCCGCTGGCTATCCACCAGATTTTGCCAGCTCACGCCATAAgg                                                 >  1:660550/1‑200 (MQ=255)                                                                                                                                               ccgcaacaaGCCGGGCAAGGTTACGGTCAACGGGCTGCAAACGTTCGGCATAGCGGGCGACGTCTTCTGCATGACCAGCGAGGGCGTTTTCCTGAATGTTTGTGCAAATAAACGCCGTCGGATGGGTGTTAAATTCGTTAAGTTTGATATCCATAATCCGCTGGCTATCCACCAGATTTTGCCAGCTCACGCCATAAGGa                                                <  1:440301/197‑1 (MQ=255)                                                                                                                                               ccTCAACAAGCCGGGCAAGGTTACGGTCAACGACCTGCAAACGTTCGGCATAACGTGCGACGTCTTCTGCATGACCAGCGAGGTCGGTTTCCTGAATGTTGGTGCAAATAAACGCCGTCGGATGGGTGTTAAATTCGTTGAGGGTGATATCCATAATCCGCTGGCTATCCACCAGATTTTGCCAGCTCACGCCATAAGGa                                                <  1:216249/200‑1 (MQ=255)                                                                                                                                               ccTCAACAAGCCGGGCAAGGTTACGGTCAACGACCTGCAAACGTTCGGCATAACGTGCGACGTCTTCTGCATGACCAGCGAGGTCGGTTTCCTGAATGTTGGTGCAAATAAACGCCGTCGGATGGGTGTTAAATTCGTTGAGGGTGATATCCATAATCCGCTGGCTATCCACCAGATTTTGCCAGCTCACGCCATAAGGa                                                <  1:545660/200‑1 (MQ=255)                                                                                                                                               ccTCAACAAGCCGGGCAAGGTTACGGTCAACGACCTGCAAACGTTCGGCATAACGTACGCCGTCTTCTGCATGACCAGCGAGGTCGGTTTCCTGAATGTTGGTGCAAATAAACGCCGTCGGATGGGTGTTAAATTCGTTGAGGGTGATATCCATAATCCGCTGGCTATCCACCAGATTTTGCCAGCTCACGCCATAAGGa                                                >  1:307179/1‑200 (MQ=255)                                                                                                                                               ccTCAACAAGCCGGGCAAGGTTACGGTCAACGACCTGCAAACGTGCGGCATAACGTGCGACGTCGTCTGCATGACCAGCGAGGTCGGTTTCCTGAATGTTGGTGCAAATAAACGCCGTCGGATGGGTGTTAAATTCGTTGAGGGTGatat                                                                                                  <  2:372934/150‑1 (MQ=255)                                                                                                                                                cTCAACAAGCCGGGCAAGGTTACGGTCAACGACCTGCAAACGTTCGGCATAACGTGCGACGTCTTCTGCATGACCAGCGAGGTCGGTTTCCTGAATGTTGGTGCAAATAAACGCCGTCGGATGGGTGTTAAATTCGTTGAGGGTGATATc                                                                                                 <  2:528322/150‑1 (MQ=255)                                                                                                                                                 tCAACAGGCCGGGCAAGGATACGGTCAACGACCTGCAAACGGTCGGCATGACGTGCGTAGTATACTGCAAGACCAGCGAGGTCGGTTTCCTGAATGTTGGTGCAAATAAACGCCGTCGGATGGGTGTTAAATTCGTTGAGGGTGATATcc                                                                                                <  2:337349/150‑1 (MQ=255)                                                                                                                                                 tCAACACGCCGGGCAAGGTTACGGTCAACGACCTGCAAACGTTCGGCATAACGTGCGACGTCTTCTGCATGACCAGCGAGGTCGGTTTCCTGAATGTTGGTGCAAATAAACGCCGTCGGATGGGTGTTAAAt                                                                                                                  <  3:69435/132‑1 (MQ=255)                                                                                                                                                 tCAACAAGCCGAGCAAGGTTACGGTCAACGACCTGCAAGCGTTCGGCATAACGTGCGACGTCTTCTGCATGACCAGCGATGTCGGTTTCCTGAATGTTGGTGCAAATAAACGCCGTCGGATGGGTGTTAAAt                                                                                                                  <  3:502249/132‑1 (MQ=255)                                                                                                                                                  caacaaGCCGGGCAAGGTTACGGTCAACGACCTGCAAGCGTGCGGCATAACGTGCGACGTCTTCTGCATGACCGGCGAGGTCGGTTTCCTGAATGTTGGTGCAAATAAACGCCGTCGGATGGGTGTTAAAtt                                                                                                                 <  3:195990/132‑1 (MQ=255)                                                                                                                                                   acgaaGGCGGGCAAGCTTGCGGTAAACGACCTGCAGACGTTCGGCATAACGTGCGACGTCTTCTGCATGACCAGCGAGGTCGGTTTCCTGTATGTTGGTGCAAATGAGCGCCGACGGATGGGTGTTAAATTCGTTGAGGGTGATATCCAt                                                                                              <  2:281605/147‑1 (MQ=255)                                                                                                                                                       aGCCGGGCAAGGTTACGGTCAACGACCTGCAAACGTTCGGCATAACGTGCGACGTCTTCTGCATGACCAGCGAGGTCGGTTTCCTGAATGTTGGTGCAAATAAACGCCGTCGGATGGGTGTTAAATTCGTTGAGGGTGata                                                                                                   >  4:39057/1‑141 (MQ=255)                                                                                                                                                        gCCGGGCAAGGTTACGGTCAACGACCTGCAAACGTTCGGCATAACGTGCGACGTCTTCTGCATGACCAGCGAGGTCGGTTTCCTGAATGTTGGTGCAAATAAACGCCGTCGGATGGGTGTTAAATTCGTTGAGGGTGATATCCATAATcc                                                                                         <  2:509013/150‑1 (MQ=255)                                                                                                                                                            ggCAAGGTTACGGTCAACGACCTGCAAACGTTCGGCATAACGTGCGACGTCTTCTGCATGACCAGCGAGGTCGGTTTCCTGAATGTTGGTGCAAATAAACGCCGTCGGATGGGTGTTAAATTCGTTGAGGGTGATATCCATAATCCGCTGGCTATCCACCAGATTTTGCCAGCTCACGCCATAAGGATTGTTGACGatat                                   >  1:15452/1‑200 (MQ=255)                                                                                                                                                                aGGTTACGGTCAACGACCTGCAAACGTTCGGCATAACGTGCGACGTCTTCTGCATGACCAGCGAGGTCGGTTTCCTGAATGTTGGTGCAAATAAACGCCGTCGGATGGGTGTTAAATTCGTTGAGGGTGATATCCATAATc                                                                                          <  4:243965/141‑1 (MQ=255)                                                                                                                                                                 ggTTACGGTCAACGACCTGCAACCGTTCGGCATAACGTGCGACGTCTTCTGCATGACCAGCGAGGTCGGTTTCATGAATGTTGGTGCAAATAAACGCCGTCGGATGGGTGTTAAATTCGTTGAGGGTGATATCCATAATcc                                                                                         <  4:440877/141‑1 (MQ=255)                                                                                                                                                                 ggTTACGGTCAACGACCTGCAAACGTTCGGCATAACGTGCGACGTCTTCTGCATGACCAGCGAGGTCGGTTTCCTGAATGTTGGTGCTACTAAACGCCGTCGGATGGGTGTTAAATTCGTTGAGGGTGatat                                                                                                  >  3:310723/1‑132 (MQ=255)                                                                                                                                                                   ttACGGTCAACGACCTGCAAACGTTCGGCATAACGTGCGACGTCTTCTGCATGACCAGCGAGGTCGGTTTCCTGAATGTTGGTGCAAATAAACGCCGTCGGATGGGTGTTAAATTCGTTGAGGGTGATATCCATAATCCGCTGGCTATCCACCAGATTTTGCCAGCTCACGCCATAAGGATTGTTGACGATATCTGCCAc                            <  1:266049/200‑1 (MQ=255)                                                                                                                                                                     aCGGTCAACGACCTGCAAACGTTCGGCATAACGTGCGACGTCTTCTGCATGACCAGCGAGGTCGGTTTCCTGAATGTTGGTGCAAATAAACGCCGTCGGATGGGTGTTAAATTCGTTGAGGGTGATATCCAt                                                                                              <  3:115621/132‑1 (MQ=255)                                                                                                                                                                     aCGGGCAACGACATGCAAACGTTCGGGAGAACGTGCGACGACGACAGCATGACCAGCGAGGTCGGTTTCCTGAAGGGTGGTGCAAATAAACGCCGTCGGATGGGTGTTAAATTCGTTGAGGGTGATATCCAt                                                                                              <  3:169369/132‑1 (MQ=255)                                                                                                                                                                        gTCAACGACCTGCAAACGTTCGGCATAACGTGCGACGTCTTCTGCATGACCAGCGAGGTCGGTTTCCTGAATGTTGGTGCAAATAAACGCCGTCGGATGGGTGTTAAATTCGTTGAGGGTGATATCCATAATCCGCTGGCTATCCACCAg                                                                         >  2:661198/1‑150 (MQ=255)                                                                                                                                                                        gTCAACGACCTGCAAACGTTCGGCATAACGTGCGACGTCTTCTGCATGACCAGCGAGGTCGGTTTCCTGAATGTTGGTGCAAATAAACGCCGTCGGATGGGTGTTAAATTCGTTGAGGGTGATATCCATAATCCGCTGGCTATCCACCAGCTTTTGCCAGCTCACGCCATAAGGATTGTTGACGATATCTGCTACCTTAc                       >  1:105158/1‑200 (MQ=255)                                                                                                                                                                        gTCAACGACCTGCAAACGTTCGGCATAACGTGCGACGTCTTCTGCATGACCAGCGAGGTCGGTTTCCTGAATGTTGGTGCAAATAAACGCCGTCGGATGGGTGTTAAATTCGTTGAGGGTGATATCCATAATCCGCTGGCTATCCACCAGATTTTGCCAGCTCACGCCATAAGGa                                                <  1:661264/175‑1 (MQ=255)                                                                                                                                                                        gTCAACGACCTGCAAACGTTCGGCATAACGTGCGACGTCCTCTGCATGACCAGCGAGGTCGGTTTCCTGAATGTTGGCGCAAATAAACGCCGTCGGATGGGTGTTAAATTCGTTGAGGGTGATATCCATAATCCGCTGGCt                                                                                  >  4:201409/1‑141 (MQ=255)                                                                                                                                                                         tCAACGACCTGCAAACGTTCGGCATAACGTGCGACGTCTTCTGCATGACCAGCGAGGTCGGTTTCCTGAATGTTGGTGCAAATAAACGCCGTCGGATGGGTGTTAAATTCGTTGAGGGTGATATCCATAATc                                                                                          <  3:374538/132‑1 (MQ=255)                                                                                                                                                                             cGACCTGCAAACGTTCGGCATAACGTGCGACGTCTTCTGCATGACCAGCGAGGTCGGTTTCCTGAATGTTGGTGCAAATAAACGCCGTCGGATGGGTGTTAAATTCGTTGAGGGTGATATCCATAATCCGCTGGCTATCCACCAGATTTTGCCAGCTCACGCCATAAGGATTGTTGACGATATCTGCCACCTTACccacc                  <  1:578579/200‑1 (MQ=255)                                                                                                                                                                                ccTGCAAACGTTCGGCATAACGTGCGACGTCTTCTGCATGACCAGCGAGGTCGGTTTCCTGAATGTTGGTGCAAATAAACGCCGTCGGATGGGTGTTAAATTCGTTGAGGGTGATATCCATAATCCGCTGGCTATccacca                                                                          >  4:421520/1‑141 (MQ=255)                                                                                                                                                                                ccTGCAAACGTTCGGCATAACGTGCGACGTCTTCTGCATGACCAGCGAGGTCGGTTTCCTGAATGTTGGTGCAAATAAACGCCGTCGGATGGGTGTTAAATTCGTTGAGGGTGATATCCATAATCCGCTGGCTATCCACCAGATTTTGcc                                                                 <  2:595662/150‑1 (MQ=255)                                                                                                                                                                                ccTGCAAACGTTCGGCATAACGTGCGACGTCTTCTGCATGACCAGCGAGGTCGGTTTCCTGAATGTTGGTGCAAATAAACGCCGTCGGATGGGTGTTAAATTCGTTGAGGGTGATATCCATAATCCGCTGGCTATCCACCAGATTTTGCCAGCTCACGCCATAAGGATTGTTGACGATATCTGCCACCTTACCCACCGGc               <  1:206285/200‑1 (MQ=255)                                                                                                                                                                                   gCAAACGTTCGGCATAACGTGCGACGTCTTCTGCATGACCAGCGAGGTCGGTTTCCTGAATGTTGGTGCAAATAAACGCCGTCGGATGGGTGTTAAATTCGTTGAGGGTGATATCCATAATCCGCTGGCTATCCACCAGAt                                                                       <  4:399921/141‑1 (MQ=255)                                                                                                                                                                                         gTTCGGCATAACGTGCGACGTCGTCTGCATGACCAGCGAGGTCGGTTTCCTGAATGTTGGTGCAAATAAACGCCGTCGGATGGGTGTTAAATTCGTTGAGGGTGATATCCATAATCCGCTGGCTATCCACCAGATTTTGcc                                                                 <  4:323899/141‑1 (MQ=255)                                                                                                                                                                                          tgcgggagaaCGGGCGACATCTTCTGCAGGACCCGCGAAGTCGGTTTCCTGGATGTTGGTGCAAATAGACGCCGTGGGATGGGTGTTAAATTCGTTGAGGGTGATATCCATAATCCGCTGGCTATCCACCAg                                                                         <  3:85218/124‑1 (MQ=255)                                                                                                                                                                                          gtCGGAATAACGTAGGACGTCTTCTGCATGACCAGCGGTGTCGGTTTCCTGAATGTTGGTGCAAATAGACGCCGTCGGATGGGTGTTAAATTCGTTGAGGGTGATATGCATAATCCGGTGGCTATCCACCAGATTTTGCCAGCTCACGcc                                                       <  2:690302/149‑1 (MQ=255)                                                                                                                                                                                             ggCATAACGTGCGACGTCTTCTGCATGACCAGCGAGGTCGGTTTCCTGAATGTTGGTGCAAATAAACGCCGTCGGATGGGTGTTAAATTCGTTGAGGGTGATATCCATAATCCGCTGGCTATCCACCAGATTTTGCCAGCTCACGCCATAAGGATTGTTGACGATATCTGCCACCTTACCCACCAGCACGGTTGGCACGc  >  1:161306/1‑200 (MQ=255)                                                                                                                                                                                             ggCATAACGTGCGACGTCTTCTGCATGACCAGCGAGGTCGGTTTCCTGAATATTGGTGCAAATAAACGCCGTCGGATGGGTGTTAAATTCGTTGAGGATGATATCCATAATCCGCTGGCTATCCACCAGAtt                                                                      >  3:462118/1‑132 (MQ=255)                                                                                                                                                                                              gCATAACGTGCGACGTCTTCTGCATGACCAGCGAGGTCGGTTTCCTGAATGTTGGTGCAAATAAACGCCGTCGGATGGGTGTTAAATTCGTTGAGGGTGATATCCATAATCCGCTGGCTATCCACCAGATTTTGCCAGCTc                                                            >  4:460145/1‑141 (MQ=255)                                                                                                                                                                                                aTAACGTGCGACGTCTTCTGCATGACCAGCGAGGTCGGTTTCCTGAATGTTGGTGCAAATAAACGCCGTCGGATGGGTGTTAAATTCGTTGAGGGTGATATCCATAATCCGCTGGCTATCCACCAGATTTTGCCAGCTCACGCCATAAgg                                                 <  2:160365/150‑1 (MQ=255)                                                                                                                                                                                                       gCGACGTCTTCTGCATGACCAGCGAGGTCGGTTTCCTGAATGTTGGTGCAAATAAACGCCGTCGGATGGGTGTTAAATTCGTTGAGGGTGATATCCATAATCCGCTGGCTATCCACCAGATTTTGCCAGCTCACGCCATAAGGAttgttg                                          <  2:573650/150‑1 (MQ=255)                                                                                                                                                                                                       |                                                                                                                                                                                              GTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCATCACGACCAGGCAATCATCTGGCTGCATCGCCTCAACAAGCCGGGCAAGGTTACGGTCAACGACCTGCAAACGTTCGGCATAACGTGCGACGTCTTCTGCATGACCAGCGAGGTCGGTTTCCTGAATGTTGGTGCAAATAAACGCCGTCGGATGGGTGTTAAATTCGTTGAGGGTGATATCCATAATCCGCTGGCTATCCACCAGATTTTGCCAGCTCACGCCATAAGGATTGTTGACGATATCTGCCACCTTACCCACCAGCACGGTTGGCACGC  >  NC_000913/3508696‑3509083

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