BRESEQ :: Evidence

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,606,197	0	G	T	37.0%	114.4 / 15.2	80	intergenic (‑32/+3)	leuV/leuP	tRNA‑Leu/tRNA‑Leu
Reads supporting (aligned to +/- strand): ref base G (17/29); new base T (0/27); total (19/61)
Fisher's exact test for biased strand distribution p-value = 1.17e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCAC  >  NC_000913/4606047‑4606344

aaGGAATAAACAAGAACGGAAGTAAAGATAATTGGTGCGCGGGGGGGGCCTGGACCCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCg                                                                                                                                                     <  6:1084990/151‑1 (MQ=21)

aaGGAAGAAACAAGAAAGGAAGTAACTATAATTGGTGCGCGGGGGGGCACTTGACCCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCg                                                                                                                                                     <  7:1545533/151‑1 (MQ=255)

 aGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGt                                                                                                                                                    <  2:838715/151‑1 (MQ=255)

   gaagaaACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGttt                                                                                                                                                  <  3:1508764/151‑1 (MQ=255)

    aaaaaaCAACAACGGAATTAAATATACTTGGGGCCCGGGGGGGCATTTGAACCCCCAGGCCCTTATGCACACTAACACCTCAAGCTAGCGCGTCTACCAATTCCGCCACGTTCGCACAGGCATCTTACTTTTTTTTATATCGCGTCGTTTg                                                                                                                                                 <  2:407101/148‑1 (MQ=16)

     ataaaCAAGAAAGGAAGTAACGATAATTGGTGCGAGGGGGGGTACTTCAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTgg                                                                                                                                                <  4:1724838/149‑1 (MQ=255)

      caaaCAATAAAGGAATTAAAGATAATTGGTGCGAGGGGGGGGACTTTAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGt                                                                                                                                               <  2:1427897/150‑1 (MQ=255)

      caaaCAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGCACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGt                                                                                                                                               <  8:1445746/150‑1 (MQ=255)

       aaaCAAGAAAGTAAGTAAAGATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTg                                                                                                                                              <  6:66555/151‑1 (MQ=255)

       aaaCAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTg                                                                                                                                              <  1:173376/151‑1 (MQ=255)

       aaaCAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTg                                                                                                                                              <  6:74939/151‑1 (MQ=255)

       aaaCAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTg                                                                                                                                              <  6:736343/151‑1 (MQ=255)

        aaCAAGAAAGTAAGTAAAGATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGc                                                                                                                                             <  7:1657327/151‑1 (MQ=255)

        aaCAAGAAAGGAAGTAAAGATAATTGGGGCGCGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGc                                                                                                                                             <  6:188889/151‑1 (MQ=255)

         aCAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGGACTTGCACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCAGCGTACTTTTTTGGATGTGGCGTGGTGTGGGGCg                                                                                                                                            >  8:71958/1‑151 (MQ=18)

          cAAGAAAGGAAGTAAACATAATTGGGGCGAGGGGGGGGACTTGAACCCCCACATCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGa                                                                                                                                           <  5:870677/151‑1 (MQ=255)

          aaaGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGa                                                                                                                                           <  6:499829/150‑1 (MQ=255)

           aaGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAg                                                                                                                                          <  2:919098/151‑1 (MQ=255)

            ataaaGGAATTAAAGATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAgg                                                                                                                                         <  2:874022/149‑1 (MQ=255)

             caaaGGAAGTAAAGATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAggg                                                                                                                                        <  2:194151/150‑1 (MQ=255)

              aaaGGAATTAAACAAAATTGGGGCGCGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAgggg                                                                                                                                       <  2:1005816/151‑1 (MQ=255)

              aaaGGAAGTAAAGATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAgggg                                                                                                                                       <  5:704867/151‑1 (MQ=255)

                agcaaGTAAACAAAATTGGGGCGAGGGGGGGGCCTTCAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAgggggg                                                                                                                                     <  1:1465672/148‑1 (MQ=21)

                 ggAAGTAAACAAAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAggggggg                                                                                                                                    <  8:431918/151‑1 (MQ=255)

                   aactAAAGATAATTGGTGCGAGGGGGGCGAGTTGCACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGa                                                                                                                                  <  5:1536603/148‑1 (MQ=255)

                                 ggTGCGAGGCGGGGGACTTGAACCCCCACGTCCGTAAGGTCACTAACACCTGAAGCTAGCGCGGCAACCAAATCCGCCACGTTCCCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCAc                                                                                                                     <  6:620356/150‑1 (MQ=25)

                                 gcggCGTGGGCGGGTCCTTCAACCCCCACGGCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTATCACAGTCATCTTACTTTTTTTGATATCCCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACg                                                                                                                    <  5:1502703/148‑1 (MQ=25)

                                  gggcctgggggggtCCTTCAACCCCCACGTCCGTAAGGACAGTAACACCTGAAGCTAGCGCGGCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGt                                                                                                                   <  8:1300708/145‑1 (MQ=25)

                                        gggttttgcctacaaCCCCCACGGCCTTAATTACACTAACACCTGACGCTACCCCGTCTACAAATTCCGCCACCTACCCACAGGCATCAAACTTTTTTTGATATCCCCTCTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTa                                                                                                              <  8:946769/137‑1 (MQ=19)

                                        ggggtttgactacaaCCCCCCCGTCCTAAAGTACACTAAAACATGCAGATATCCCTTATACCAATTCCGCCACCAAACCACAGGCAACAAAATTTTTTTGATATCCCCTCTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTaa                                                                                                             <  8:863614/138‑1 (MQ=12)

                                          ggttttcctacaaCCCCCACGGCCGTTATTACAGTAAAACCTGAAGCTATCTCTTAAACCAATTCCGGCACCATCCCACACGCATCAAACTTTTTTTGATCTCCCCTCTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGa                                                                                                           <  1:449425/140‑1 (MQ=16)

                                          gcggggcctacaaCCCCCCCGGCCCTAAATACAGTAACACCTGACGCTCGCCCTGAAACCAATTCCGCCAGCTTCTAACAGGCATGTAACTTTTTTTGATATCGCCTCTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGa                                                                                                           <  6:1281165/140‑1 (MQ=16)

                                              gccctcaaccccccACGGCCGTAAGGACACTAACACCTGAAGCTAGCGGGTCTACAAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCCCCTCTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAacac                                                                                                       <  1:893529/142‑1 (MQ=37)

                                                   acaaCCCCCACTGCCTTATGTACACTAACACATGAAGCTCCCCAGTAAACCAATTCCGCCACCATCCAACAGGCATCTTACTTTTTTTGATATCCCCTCTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAaca                                                                                                  <  8:1417599/149‑1 (MQ=18)

                                                   acaaCCCCCACGGCCTTTATTACATTAACACCTGAAGCTCCATCTGATACAAAATCCGCCACCAACCAACAGGCAACTAACATTTTTTGAATTCCCCACGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAaca                                                                                                  <  1:777904/149‑1 (MQ=9)

                                                    gccgcccctggTCCGTAAGGACCCTAACACCGGATGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACATGCATGTTAATTTTTTGGATCTCGCCTATTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAaca                                                                                                  <  1:1468597/140‑1 (MQ=21)

                                                     accccccACGGCCGTAAGCACACAAAATCCGGCCGCTATCGCGTCAACCCTTTGGCCCACCTTCGCACAGGCAACAAACATTATTTGATATACCCTTCTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACAcc                                                                                                <  3:1390589/149‑1 (MQ=12)

                                                      aCCCCCACGGCCTTAATTACAGTACAACCTGACTATAGCCCTTATACAAATTCCGCCACCTACCAACAGTCATGTAACTTTATTTGATCTCGCCTCTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCt                                                                                               <  8:200804/151‑1 (MQ=12)

                                                       cccccACGGCCGTATGTACACTAAAACCTGAAGATACCACGTAAACCAATTCCGCCACCTTCCCACAGTCATCTTACTTTTTTTGATATCCCCTCTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTg                                                                                              <  1:570411/151‑1 (MQ=25)

                                                        ccccccGTCCGTTACGACAATAAGTGCTGACTCCAGCGCGTATTCCAAATCCGCCACCAAACAACATGCATGTAACTTTTATTGATCTACCCTTCTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGa                                                                                             <  1:668155/151‑1 (MQ=9)

                                                              ggccttatctacacTAACACCTGACGCTAGCGCGTATACCAATTCCGCCACATTCCAACAGGCATCTTACTTTTTTTTATTTCGCCTCTTTTGGGGGGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAg                                                                                       <  8:158439/141‑1 (MQ=21)

                                                                       gACAGTTACACCTGCTGCTTGCTCGTCTACCAATTCCGCCACCTTCGCTCAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGTGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTa                                                                               <  2:1500844/150‑1 (MQ=25)

                                                                        acagTAACACCTGCATGTAGATCGTCTACCAATTCCGCCAACTACTAACAGTCATCTTAGTTTTTTTGAAATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTAc                                                                              <  3:922228/150‑1 (MQ=21)

                                                                        acacTAACACCTGAAGCTCTCTCTTCATTAAATTCCGCCACCAAACCACATTCATCAAACATTTTTTGATATACCCTCCTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTAc                                                                              <  3:246823/150‑1 (MQ=17)

                                                                         agctcacacCTGATGCGCGCGCGTCTACCATTGCCGCCCGCAACGAACAGTAAGCTTAATTTTTTTTAAAGCCCGCCTTTGGGGGCGTGGGGGGGGAGTTCCCCCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCa                                                                            <  8:713946/146‑1 (MQ=9)

                                                                                 ccTGAAGCTCCCCCGGCTACCAATCCCGCCACCAACCCACAGGCATATTACTTTTTTTGAAAGCCCCTGTTTGGGGGCGCGGGGGGGCACTTCACCCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGcc                                                                    <  8:706818/151‑1 (MQ=12)

                                                                                              cgTCTACCAATGCCGCCACCTACCCACAGGCATCTAACTTTTTTTGATATCGCCTCTTTTGGGGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGt                                                       <  1:605354/151‑1 (MQ=18)

                                                                                                tCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTAACAAATCAACCACCtactcaaagga                                                      >  5:1378491/1‑140 (MQ=25)

                                                                                                 cTACCAATTCCGCCACCTTCGCAGACTCAGATTACTTTTAATGATATCGGCGCCCGGGGGCCGAGGGTGGGGGCTTGCACCCGCTCGTCAGTACTTACACTTCCACATCCCGGTAGATCCTCTACCAATTCCGCCAcataccgcaagagcc                                                    >  1:444197/1‑137 (MQ=255)

                                                                                                  tACCAATTCCGCCACCTTCGCAAACAAAGAATACTTTAAAAGATATATCCTAATATCGTGCTAGGGGGTAAACTTGAACCCCAAAGTCGATAAGAACACTAACACCTGAAGATAAAAAAACAACCAATTCCGCCACCTTCGCCCAGGGCg                                                    >  4:141988/1‑150 (MQ=255)

                                                                                                   aCCAATTCCGCCACCTTACAACAGTCATCTTACTTTTTTTGATTTAGCCTATTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGa                                                   <  8:675923/150‑1 (MQ=18)

                                                                                                      aaTTCCGCCACCTTCGCACAGTGATCTTAATGTTATTGATATGGAATCGTTTGGTGCGATGGGGTGGACTAGAACCCCCACGTCCGAAAGAAAAATAACACCTGAAGCTAGCGCGTCTACCAATTCCGGCGCCTTCGCAGAGTGCgggagg                                               >  5:1024781/1‑146 (MQ=17)

                                                                                                         tCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAatat                                            >  8:672351/1‑151 (MQ=255)

                                                                                                         tCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCACTGCGAGGAata                                             >  3:848027/1‑150 (MQ=255)

                                                                                                          ccGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCACCTTCCGTAAGAACACTAACACCTGAAGCTAGAGCGTCTACCAATTCCGCCACCTTCGccgactgggacgtcttt                                            >  6:721451/1‑135 (MQ=21)

                                                                                                           cGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGCCTTGAACCCCGACGTCAGTAAGTACACTACCACCTGAAGCTACAGCGTCTACCAATTCCGCCCCCTTAGCGCAGGGCGAggattttt                                           >  6:1359225/1‑143 (MQ=255)

                                                                                                             gcaacaaataacaTTCATCAAACTTTTTTTGATTTATAGTTTTTTGGTGCGTGGGGGGGGACTTGAACCCCCACGTCCGTAAGTACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAAc                                        <  7:363825/141‑1 (MQ=14)

                                                                                                             ccACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCGAgggtgaagactataaa                                         >  8:621336/1‑135 (MQ=21)

                                                                                                               tccTTCTCACTGTCATCTTACTTTTTTTTATATCGCGTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGcataccatcaat                                                 <  5:647608/139‑12 (MQ=14)

                                                                                                               aCCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGcataccatcaat                                                 >  6:645594/1‑129 (MQ=255)

                                                                                                                 caaaaaacatggaTCTTACTTTAATTGAATTCGAGTCTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTg                                     <  6:25578/138‑1 (MQ=17)

                                                                                                                  aaaaatcattgaagatACTTTTATTGATTTATTCTTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGTACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTgg                                    <  4:1170002/135‑1 (MQ=17)

                                                                                                                   tctatcattcatcaTACTTTTATTGATTTAGAGTTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTgg                                    <  7:530599/141‑1 (MQ=17)

                                                                                                                    atattattgattttACATTTTTTGATTTTTCGTTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGTACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGttt                                 <  7:249629/139‑1 (MQ=17)

                                                                                                                    atatcattcaacatACTTTTATTGATTTATACTTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGttt                                 <  2:250209/138‑1 (MQ=17)

                                                                                                                     aatcactcaagaaaCTTTATTTTATTTCTCCTTTTTTGGTGCTTGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGtttt                                <  7:1642254/138‑1 (MQ=17)

                                                                                                                      atcactcatcaAACTTTTTATGATATAAACTTATTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGttttt                               <  3:1473186/146‑1 (MQ=18)

                                                                                                                        cagccagcaaaCTTTTTTTGATATCCCCTCTTTTGGTGCGTGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTgg                             <  6:1769483/141‑1 (MQ=21)

                                                                                                                         attcatcaAACTTTTTTTGATTTCCCGTATTTTGGTGGGTGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGt                            <  4:70796/149‑1 (MQ=18)

                                                                                                                         acgcATCTAACTTTTATGGATCTCCCGTGATTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGt                            <  1:1607251/148‑1 (MQ=18)

                                                                                                                          gTCATGTTCGTTTTTTTGATTTCCAGTCTTTTGGGGCGCGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTg                           <  8:772882/151‑1 (MQ=18)

                                                                                                                             aTCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCACCCCGAGTAATATAAAAGATGTTTTTCGTGCgc                        >  6:449177/1‑150 (MQ=18)

                                                                                                                              acatACTTTTTTTTATCTCCACTTTTTTGGGGCGTGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAg                       <  4:987878/148‑1 (MQ=18)

                                                                                                                                 tACTTTTTTTGATTTATCCACTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCg                         <  3:321277/146‑1 (MQ=255)

                                                                                                                                 gtcTTTTTTTGATGTGGCCTCTTGTGGGGGGCGGGGGGGGACTTCAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAgggg                    <  8:1560145/149‑1 (MQ=21)

                                                                                                                                  aCTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCACGACCACAAATTTAAACGAGGTTTTTGGTCCGAggggg                   >  1:1633242/1‑151 (MQ=18)

                                                                                                                                   cTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTGTTGGTGCGAgggggg                  >  4:210079/1‑151 (MQ=255)

                                                                                                                                            atatCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTACCACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCGCACTAATAAAAACG‑AGTTTTTGGGCCGAGGGGGGGGACTTGacc        >  1:1503455/1‑149 (MQ=14)

                                                                                                                                             tatCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCTAGGGAGAGCTATATCAAAGAGGTTGTTGGTGCGAggggggggcgtgtccg        >  8:1319797/1‑143 (MQ=18)

                                                                                                                                                cGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGCACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCACTGCCAGCAATATAAACGAGTTTTTTGCTGCGAGGCGGGGGCCtgtccgctc     >  8:682166/1‑143 (MQ=17)

                                                                                                                                                cGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCGAGAGGGAGCAATATCAACGTCGTGATTCGAGCGAGGGGGGGGCCTTGAAccac     >  2:903909/1‑149 (MQ=18)

                                                                                                                                                   cTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCGAGAGACAGGTATATCAAAGAAGTGGTTCGTGGGAggggggggcgaggcgccagcac  >  1:3631/1‑137 (MQ=14)

AAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCAC  >  NC_000913/4606047‑4606344

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 14 ≤` `ATCG/ATCG` `< 15 ≤` `ATCG/ATCG` `< 21 ≤` `ATCG/ATCG` `< 36 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`