BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,828,980	A→C	6.8%	Q12H (CAA→CAC)	xanP →	xanthine permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,828,980	0	A	C	6.8%	1490.0 / 8.4	557	Q12H (CAA→CAC)	xanP	xanthine permease
Reads supporting (aligned to +/- strand): ref base A (313/190); new base C (4/50); total (317/240)
Fisher's exact test for biased strand distribution p-value = 1.04e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 580 total aligned reads displayed.

ATATAAAAAAGTGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAA‑CGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTGCCAGGCGCTGGG  >  NC_000913/3828831‑3829129

   taAAAAAGTGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAA‑CGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCg                                                                                                                                                    <  2:12510/150‑1 (MQ=255)

    aaaaaaGTGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAA‑CGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCgg                                                                                                                                                   <  2:3678991/150‑1 (MQ=255)

      aaaaGTGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAA‑CGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGtt                                                                                                                                                 >  2:2608559/1‑150 (MQ=255)

       aaaGTGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAA‑CGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTg                                                                                                                                                <  1:5442424/150‑1 (MQ=255)

         aGTGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAA‑CGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTgcg                                                                                                                                              <  2:1847039/150‑1 (MQ=255)

             cATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAA‑CGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGa                                                                                                                                          <  2:3553710/150‑1 (MQ=255)

             cATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAA‑CGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGa                                                                                                                                          <  2:1011542/150‑1 (MQ=255)

               tttATCTCATTTAAATATTTAGATGAGGATGTGGCAAA‑CGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACt                                                                                                                                        >  2:1893329/1‑150 (MQ=255)

                ttATCGCATTTAAATATTTAGATGAGGATGTGGCAAA‑CTTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACTGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCattttc                                                                                                                                       <  2:2761570/150‑6 (MQ=255)

                 tATCTCATTTAAATATTTAGATGAGGATGTGGCAAA‑CGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCa                                                                                                                                      <  2:2288875/150‑1 (MQ=255)

                  aTCTCATTTAAATATTTAGATGAGGATGTGGCAAA‑CGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCaa                                                                                                                                     <  2:4286007/150‑1 (MQ=255)

                   tctcATTTAAATATTTAGATGAGGATGTGGCAAA‑CGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCaaa                                                                                                                                    >  2:387166/1‑150 (MQ=255)

                     tcATTTAAATATTTAGATGAGGATGTGGCAAA‑CGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCCACCGGTTGCGCAGACTCAAAAc                                                                                                                                  >  2:5050895/1‑150 (MQ=255)

                     tcATTTAAATATTTAGATGAGGATGTGGCAAA‑CGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCCACCGCTTGCCCAGACACAAAAc                                                                                                                                  >  2:3085701/1‑150 (MQ=255)

                     tcATTTAAATATTTAGATGAGGATGTGGCAAA‑CGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAAc                                                                                                                                  >  1:1359538/1‑150 (MQ=255)

                     tcATTTAAATATTTAGATGAGGATGTGGCAAA‑CGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAAc                                                                                                                                  >  1:4178486/1‑150 (MQ=255)

                         ttAAATATTTAGATGAGGATGTGGCAAA‑CGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCg                                                                                                                              >  1:307956/1‑150 (MQ=255)

                          tAAATATTTAGATGAGGATGTGGCAAA‑CGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGa                                                                                                                             >  2:1378076/1‑150 (MQ=255)

                            aaTATTTAGATGAGGATGTGGCAAA‑CGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAAc                                                                                                                           >  2:3194674/1‑150 (MQ=255)

                               atTTAGATGAGGATGTGGCAAA‑CGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGa                                                                                                                        <  2:1477739/150‑1 (MQ=255)

                                 ttAGATGAGGATGTGGCAAA‑CGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGAtt                                                                                                                      >  2:5803476/1‑150 (MQ=255)

                                      tGAGGATGTGGCAAA‑CGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCCACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCg                                                                                                                 >  1:2717497/1‑150 (MQ=255)

                                          gATGTGGAATA‑CGTTTGCTTTTCTCTTTCGGTCAGATAATATGCCCGATTTTATACCAGTGTATTGTTTTCGTTGTCTGTTTCCACCCTCGAGTCAGAAATTGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACGGATTTACCGTCtt                                                                                                             <  2:5495476/150‑1 (MQ=255)

                                              gtgCAAACCGTTTGCTTTTCTCTTTCGGTCAGCTAAAAGGCCCGCTTGGCTACCTGGTTTTTTTTTTCGTTGTCTGTTTCCACCCTCGATTCATAAAAGGCGCAACGGGTGGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGaa                                                                                                          <  1:4072992/148‑1 (MQ=255)

                                               ggCAAA‑CGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGa                                                                                                        >  2:4591057/1‑150 (MQ=255)

                                                gCAAA‑CGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTGTGTTTCCACCCTCGAGTCAGAAAATGCGCAAGCGGTTGGGCAGACTCAAAACAGCGAAGTGTTTTACGGTCTTGTAGAt                                                                                                       >  2:2010908/1‑150 (MQ=255)

                                                  aaa‑CGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCg                                                                                                     >  1:832046/1‑150 (MQ=255)

                                                  aaa‑CGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACACAAAACAGCGAACTGAATTAACGTCTTGAAGAtag                                                                                                     >  1:4522069/1‑148 (MQ=255)

                                                      cGTTTGCTTTTCTCTTTCGGTCAGATAAAATGCCCGCTTTTCTACCACGGTATTGTTTTCGATGTGTTTTTCCACCCTCGAGTCATAAAATGCGCAACGGGTGGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTcc                                                                                                  <  1:3337087/150‑1 (MQ=255)

                                                        tttGCTTTTCTCTTTGGGTCAGCTAAATTGCCCGTTTCGATACCCGGTGTTTTTTTTCGTTGTGTGTTTCCACCCTCGAGTCAGAAACTGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTccgc                                                                                                <  1:3689925/150‑1 (MQ=255)

                                                            cTTTTCTCTTGCGGTCAGCTAAACTGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCACCCGGGTGCGCAGACTCAACACAGCGAACTGAGTTACCGTCTTGACGATCGTCCGCAGCt                                                                                            <  2:690022/150‑1 (MQ=255)

                                                             ttttCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCtt                                                                                           >  2:4395722/1‑150 (MQ=255)

                                                             ttttCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCtt                                                                                           >  2:5138284/1‑150 (MQ=255)

                                                             ttttCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCtt                                                                                           >  2:4731122/1‑150 (MQ=255)

                                                              tttCTCTTTCGGTCAGATAAAATCCCCGCTTTGATACCAGGGTTTTGTTTTCTTTTTCTTTTTCCACCCTCTAGTCATAAAATGCGCAACCGGTGGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTc                                                                                          <  1:5284379/150‑1 (MQ=255)

                                                              tttCTCTTTCGCTCAGAAAAAATGCCCGCTTTGATACCACGGTTTTGTTTTCGATGTCTTTTACAACCCTCCAGACACAAAATGCCCAACCGGTCGCCCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTc                                                                                          <  1:1541666/150‑1 (MQ=255)

                                                                         gTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTg                                                                               >  2:2246246/1‑150 (MQ=255)

                                                                         gTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAAAAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTg                                                                               >  2:3025231/1‑150 (MQ=255)

                                                                          tCAGATAAAATCCCCGCTTTGCTACCACGGTATTGTTTTCGATGTCTGTTTCCACCATCGAGTAACAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGt                                                                              <  2:156568/150‑1 (MQ=255)

                                                                            aGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGttt                                                                            >  2:523209/1‑150 (MQ=255)

                                                                             gCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAAAAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCGTCCTCAAACCCTGTTTg                                                                           >  1:1097864/1‑150 (MQ=255)

                                                                              cTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTgc                                                                          <  2:1151940/150‑1 (MQ=255)

                                                                              cTAAAAGCCCCGCTGGGCTACCACGGGATTTTTTTCGATGTCGGTTTCCACCCTCGCGGCAGAAAATGCGCACCCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTgc                                                                          <  1:3225019/150‑1 (MQ=255)

                                                                               tAAAATGCCCGCTTTGCTACCACGGGTTTGTTTTCGATGTCTGTTCCCACCCTCGAGTCAGAAATTGCGCAACGGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTgcc                                                                         <  1:1028079/150‑1 (MQ=255)

                                                                               tAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAATATCGTCCGCCGCTTCCTCAAACCCTGTTTgcc                                                                         >  1:4892799/1‑150 (MQ=255)

                                                                               tAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAAAAGCGTACTGATTTACCGTGTTGTTGTTGGTGCGCCGCGTCCTAAAACCCTGTTTgcc                                                                         >  1:1241642/1‑150 (MQ=255)

                                                                               tAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAAAAGCGAACGTATTTAACCTCTTGAAGATCGTCCGCCGCGTACACAAACCCGGTTTgcc                                                                         >  1:2038134/1‑150 (MQ=255)

                                                                                aaaaTGCCCGCTGTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGATTCATAAAATGCGCTACCGGTTGCGCAGACTCATAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTgccg                                                                        <  1:3120358/150‑1 (MQ=255)

                                                                                aaaaGGCCCGCTTGGCTACCCCGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTgccg                                                                        <  1:2840823/150‑1 (MQ=255)

                                                                                 aattGCCCGCTTTGCTACCCGGTGTTTTTTTTCGATGTGTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTgccgc                                                                       <  1:912094/147‑1 (MQ=255)

                                                                                   aTGCCCGCTTTGCTACCACGGGTTTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCt                                                                     <  2:4799789/150‑1 (MQ=255)

                                                                                    tGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTTAAGATCGTCCGCCCCTTCCACAAACCCTGTTTgccgccag                                                                    >  1:1172774/1‑148 (MQ=255)

                                                                                      cccGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTc                                                                  <  2:334030/150‑1 (MQ=255)

                                                                                        cGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCCACCGGTTGCCCAGACTCAAAACACCGAACTGATTTAACGTCATGAATATCGTCCGCCGCGACCTCAAACCCGGGTTGGCGCCAGTCAg                                                                >  1:3524716/1‑150 (MQ=255)

                                                                                        cGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGCACTGATTTAACGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCGGTTTGCCGCCTGTCAg                                                                >  1:4627899/1‑150 (MQ=255)

                                                                                         gCTTTGCTACCACGTGTTTGTTTTCGTTGTCTGTTTCCCCCCTCGATTCAGAAATTGCGCAACCGGTGGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGc                                                               <  1:5050895/150‑1 (MQ=255)

                                                                                               cTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGGTGCGCCGACTCAAAACAGCCAACTGATTTACCGTCGTGAAGAACGTCCGCCGCTTCCACAAACCCGGGTTGCCGCCTGTCAGCATctgc                                                         >  2:5376138/1‑150 (MQ=255)

                                                                                                tACCTCGGGATTGTTTTCGATGTCTGTTCCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATctgct                                                        <  1:4783227/150‑1 (MQ=255)

                                                                                                     cGGGATTGTTTTCGTTGTCTGTTTCCACCCGCGAGCCAGAAAATGCGCACCCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGa                                                   <  1:5720605/150‑1 (MQ=255)

                                                                                                       ggATTGTTTTCTATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATg                                                 <  2:2809438/150‑1 (MQ=255)

                                                                                                       ggATGGTTTTGGTTGTCGGTTTCCCCCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGTAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATg                                                 <  2:5682993/150‑1 (MQ=255)

                                                                                                         aTTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAAAAGAGAACTGAATTAACCTCTTGAAGATAGTCCGCCGCGTCCACAAAACCAGTTTGCCGCCGGTCAGCATCTGCTGGCGCTGtt                                               >  1:1914760/1‑150 (MQ=255)

                                                                                                          ttGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTc                                              <  1:2352824/150‑1 (MQ=255)

                                                                                                             ttttCGTTGTCTGTTTCCACCCTCGAGACAGAAAATGCGCACCCGGTGGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGtt                                           <  1:5818465/150‑1 (MQ=255)

                                                                                                             ttttCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGtt                                           >  2:2668745/1‑150 (MQ=255)

                                                                                                              tttCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCCGTTGCGCAGACGCAAAAAAGCGAACTGATTTACCCTCTTGAAGATAGTCCCCCGCATCCACAAACCCAGGTTGCCGCCGGTCAGCAACTGCTGGCGATgttggttg                                          >  1:182297/1‑150 (MQ=255)

                                                                                                               ttCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCCACCGGTTGCGCAGACGCAAAAAACCGAACATAGTTAACGGCGTGAAGATCGTCCGCCGCGTCCTCAAAACCGGGTTTCCGCCTGTCAGCATATGCTGGCGCTGTTCGTTGc                                         >  1:5167569/1‑150 (MQ=255)

                                                                                                               ttCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGc                                         >  2:4652341/1‑150 (MQ=255)

                                                                                                               ttCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGc                                         >  2:948831/1‑150 (MQ=255)

                                                                                                               ttCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAAACCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGc                                         >  2:1752746/1‑150 (MQ=255)

                                                                                                                tCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCg                                        <  1:343354/150‑1 (MQ=255)

                                                                                                                  gATGTCTGTTTCCACCCTCGAGACAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGt                                      <  1:1663966/150‑1 (MQ=255)

                                                                                                                   aTGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCCACCGGGTTCGCAGACACAAAAAAGCCAACAGTTTTAACCTCTTGAAGATCGTCCGCCGCTTCCTCAAAACCTGTTTGGCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTg                                     >  1:3831842/1‑150 (MQ=255)

                                                                                                                     gTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCCGTTGCGCAGACACAAAAAAGCGAACAGATTTAAAGTCTTGTATATCCTCCCCCGCGTCCTCAAACCCAGGTTTCCGCCGGTCAGCAACAGCAGGCGATGTTCGTTGCGGTgtt                                   >  1:2217967/1‑148 (MQ=255)

                                                                                                                       cTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCa                                 >  2:3049864/1‑150 (MQ=255)

                                                                                                                       cTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCa                                 <  2:2393934/150‑1 (MQ=255)

                                                                                                                         gTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACg                               <  1:1610346/150‑1 (MQ=255)

                                                                                                                           ttCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGcc                             >  2:1808464/1‑150 (MQ=255)

                                                                                                                             ccACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAg                           >  2:3458036/1‑150 (MQ=255)

                                                                                                                             ccACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAg                           >  2:638819/1‑150 (MQ=255)

                                                                                                                               aCCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAgcg                         <  2:2145772/150‑1 (MQ=255)

                                                                                                                                cccTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGAATTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAgcgc                        <  2:2861661/150‑1 (MQ=255)

                                                                                                                                 ccTCGAGTCAGAAAATGCGCAACCGGGTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCt                       <  2:510678/150‑1 (MQ=255)

                                                                                                                                  cTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCGTCCTCAAACCCTGTTTGCCGCCTGTCAGCCTCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTa                      >  2:314695/1‑150 (MQ=255)

                                                                                                                                   tCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTAt                     <  2:3778663/150‑1 (MQ=255)

                                                                                                                                        tCAGAAAATGCGCAACCGCTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATc                >  2:5440132/1‑150 (MQ=255)

                                                                                                                                          aGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCTGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTg              >  2:4244165/1‑150 (MQ=255)

                                                                                                                                           gAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTGc             <  2:4836307/150‑1 (MQ=255)

                                                                                                                                           gAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTGc             <  2:5598764/150‑1 (MQ=255)

                                                                                                                                           gAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTGc             <  2:5733832/150‑1 (MQ=255)

                                                                                                                                            aaaaTGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTGcc            >  2:1102764/1‑150 (MQ=255)

                                                                                                                                             aaaTGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTGCCa           <  1:4385553/150‑1 (MQ=255)

                                                                                                                                             aaaTGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACAGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCCGCGCTATTAATCTGCCa           >  2:4407129/1‑150 (MQ=255)

                                                                                                                                              aaTGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTGCCAg          >  2:4899541/1‑150 (MQ=255)

                                                                                                                                                   gcAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTGCCAGGCGCt     >  1:974513/1‑150 (MQ=255)

                                                                                                                                                    cAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTGCCAGGCGCTg    <  2:4513345/150‑1 (MQ=255)

                                                                                                                                                     aaCCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTGCCAGGCGCTgg   >  1:5556916/1‑150 (MQ=255)

                                                                                                                                                     aaCCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTGCCAGGCGCTgg   <  1:4710425/150‑1 (MQ=255)

                                                                                                                                                     aaCCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTGCCAGGCGCTgg   <  2:1296676/150‑1 (MQ=255)

                                                                                                                                                     aaCCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTGCCAGGCGCTgg   >  2:4011288/1‑150 (MQ=255)

ATATAAAAAAGTGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAA‑CGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTGCCAGGCGCTGGG  >  NC_000913/3828831‑3829129

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 11 ≤` `ATCG/ATCG` `< 12 ≤` `ATCG/ATCG` `< 37 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`