BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,396,961	C→A	5.3%	A279S (GCG→TCG)	rng ←	RNase G

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,396,961	0	C	A	5.3%	634.1 / 11.7	186	A279S (GCG→TCG)	rng	RNase G
Reads supporting (aligned to +/- strand): ref base C (83/91); new base A (9/2); total (92/94)
Fisher's exact test for biased strand distribution p-value = 3.26e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 195 total aligned reads displayed.

CGTCGCTTCAATATTGGTATTGAAAATGGTGTCGTCCAGATTGCGATGACCGACAAACGCTCCGGTATTGATGTCCACGGTGGTCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGAATGTACTCCGAGGTGAACTCAAGTAACGCTTCGTAAGTCAGGCGTGAGTCAACGCGAATGCGGT  >  NC_000913/3396812‑3397110

cGTCGCTTCAATATTGGTATTGAAAATGGTGTCGTCCAGATTGCGATGACCGACAAACGCTCCGGTATTGATGTCCACGGTGGTCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAgcgc                                                                                                                                                       >  2:301717/1‑150 (MQ=255)

  tCGCTTCAATATTGGTATTGAAAATGGTGTCGTCCAGATTGCGATGACCGACAAACGCTCCGGTATTGATGTCCACGGTGGTCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGcgctc                                                                                                                                                     <  1:12527/150‑1 (MQ=255)

    gCTTCAATATTGGTATTGAAAATGGTGTCGTCCAGATTGCGATGACCGACAAACGCTCCGGTATTGATGTCCACGGTGGTCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGcgctcgc                                                                                                                                                   >  1:985933/1‑150 (MQ=255)

    gCTTCAATATTGGTATTGAAAATGGTGTCGTCCAGATTGCGATGACCGACAAACGCTCCGGTATTGATGTCCACGGTGGTCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGcgctcgc                                                                                                                                                   <  1:1293914/150‑1 (MQ=255)

     cTTAAATATTGGTATTGAAAATGGTGTAGTCCAGATTGCGATGACCGACAAACGCTCAGCTATTGATGTCCACGGTGGTAATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGcgctcgct                                                                                                                                                  >  2:521835/1‑150 (MQ=255)

      ttCAATATTGGTATTGAAAATGGTGTCGTCCAGATTGCGATGACCGACAAACGCTCCGGTATTGATGTCCACGGTGGTCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTg                                                                                                                                                 >  2:940673/1‑150 (MQ=255)

       tCAATATTGGTATTGAAAATGGTGTCGTCCAGATTGCGATGACCGACAAACGCTCCGGTATTGATGTCCACGGTGGTCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTgg                                                                                                                                                >  1:1832625/1‑150 (MQ=255)

        cAATATTGGTATTGAAAATGGTGTCGTCCAGATTGCGATGACCGACAAACGCTCCGGTATTGATGTCCACGGTGGTCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGa                                                                                                                                               >  1:1474652/1‑150 (MQ=255)

         aaTATTGGTATTGAAAATGGTGTCGTCCAGATTGCGATGACCGACAAACGCTCCGGTATTGATGTCCACGGTGGTCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGAt                                                                                                                                              >  1:2211309/1‑150 (MQ=255)

          atatTGGTATTGAAAATGGTGTCGTCCAGATTGCGATGACCGACAAACGCTCCGGTATTGATGTCCACGGTGGTCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGAtt                                                                                                                                             >  2:37879/1‑150 (MQ=255)

           tatTGGTATTGAAAATGGTGTCGTCCAGATTGCGATGACCGACAAACGCTCCGGTATTGATGTCCACGGTGGTCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGAGCGCTGGaacg                                                                                                                                            <  2:161604/150‑4 (MQ=255)

           aatTGGTATTGAAAATGGTGTCGTCCAGATTGCGATGACCGACAAACGCTCCGGTATTGATGTCCACGGTGGTCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGAGCGCTGGaacg                                                                                                                                            <  1:1077402/149‑4 (MQ=255)

             ttGGTATTGAAAATGGTGTCGTCCAGATTGCGATGACCGACAAACGCTCCGGTATTGATGTCCACGGTGGTCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCg                                                                                                                                          <  1:771567/150‑1 (MQ=255)

              tGGTATTGAAAATGGTGTCGTCCAGATTGCGATGACCGACAAACGCTCCGGTATTGATGTCCACGGTGGTCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGt                                                                                                                                         <  2:210323/150‑1 (MQ=255)

               ggTATTGAAAATGGTGTCGTCCAGATTGCGATGACCGACAAACGCTCCGGTATTGATGTCCACGGTGGTCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGtt                                                                                                                                        >  1:1170762/1‑150 (MQ=255)

                gTATTGAAAATGGTGTCGTCCAGATTGCGATGACCGACAAACGCTCCGGTATTGATGTCCACGGTGGTCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGttt                                                                                                                                       >  1:25914/1‑150 (MQ=255)

                    tGAAAATGGTGTCGTCCAGATTGCGATGACCGACAAACGCTCCGGTATTGATGTCCACGGTGGTCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGa                                                                                                                                   >  1:1224709/1‑150 (MQ=255)

                      aaaaTGGTGTCGTCCAGATTGCGATGACCGACAAACTCTCCGGTATTGATGTCCACGGTGGTCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACa                                                                                                                                 <  2:753357/150‑1 (MQ=255)

                             tgacgtcCAGATTTCGATGACCGACAAACGCTCCGGTATTGATGTCCACGGTGGTCATCGCTTCGGTCTTGTCGATAATGAGATATCCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCTACATCAAaga                                                                                                                          <  2:1003468/147‑1 (MQ=255)

                                 gtcCAGATTGCGATGACCGACAAACGCTCCGGTATTGATGTCCACGGTGGTCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATc                                                                                                                      <  2:1474652/150‑1 (MQ=255)

                                  tcCAGATTGCGATGACCGACAAACGCTCCGGTATTGATGTCCACGGTGGTCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCg                                                                                                                     <  1:1035816/150‑1 (MQ=255)

                                    cAGATTGCGATGACCGACAAACGCTCCGGTATTGATGTCCACGGTGGTCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGaa                                                                                                                   <  2:2414380/150‑1 (MQ=255)

                                     aGATTGCGATGACCGACAAACGCTCCGGTATTGATGTCCACGGTGGTCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGaaa                                                                                                                  >  1:476016/1‑150 (MQ=255)

                                     aGATTGCGATGACCGACAAACGCTCCGGTATTGATGTCCACGGTGGTCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTATACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGaaa                                                                                                                  >  2:214975/1‑150 (MQ=255)

                                          gCGATGACCGACAAACGCTCCGGTATTGATGTCCACGGTGGTCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCgg                                                                                                             >  1:2113842/1‑150 (MQ=255)

                                          gCGATGACCGACAAACGCTCCGGTATTGATGTCCACGGTGGTCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCgg                                                                                                             >  1:762533/1‑150 (MQ=255)

                                            gATGACCGACAAACGCTCCGGTATTGATGTCCACGGTGGTCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCt                                                                                                           <  1:1238982/150‑1 (MQ=255)

                                              tGACCGACAAACGCTCCGGTATTGATGTCCACGGTGGTCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTgg                                                                                                         >  1:873229/1‑150 (MQ=255)

                                              tGACCGACAAACGCTCCGGTATTGATGTCCACGGTGGTCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTgg                                                                                                         <  1:794719/150‑1 (MQ=255)

                                               gACCGACAAACGCTCCGGTATTGATGTCCACGGTGGTCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGc                                                                                                        >  2:848582/1‑150 (MQ=255)

                                                 ccGACAAACGCTCCGGTATTGATGTCCACGGTGGTCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGt                                                                                                      >  1:2447436/1‑150 (MQ=255)

                                                 ccGACAAACGCTCCGGTATTGATGTCCACGGTGGTCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGt                                                                                                      >  2:1790938/1‑150 (MQ=255)

                                                   gACAAACGCTCCGGTATTGATGTCCACGGTGGTCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTcc                                                                                                    >  2:600141/1‑150 (MQ=255)

                                                     cAAACGCTCCGGTATTGATGTCCACGGTGGTCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCtg                                                                                                  <  1:37879/150‑1 (MQ=255)

                                                     cAAACGCTCCGGTATTGATGTCCACGGTGGTCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCtg                                                                                                  <  1:1463261/150‑1 (MQ=255)

                                                        aCGCTCCGGTATTGATGTCCACGGTGGTCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTgtgt                                                                                               <  1:491748/150‑1 (MQ=255)

                                                         cGCTCCGGTATTGATGTCCACGGTGGTCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTa                                                                                              <  2:1832625/150‑1 (MQ=255)

                                                               ggTATTGATGTCCACGGTGGTCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCTGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTc                                                                                        <  2:1038243/150‑1 (MQ=255)

                                                                 tATTGATGTCCACGGTGGTCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCa                                                                                      <  2:1460011/150‑1 (MQ=255)

                                                                  aTTGATGTCCACGGTGGTCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAg                                                                                     <  1:592/150‑1 (MQ=255)

                                                                  aTTGATGTCCACGGTGGTCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAg                                                                                     <  2:1176465/150‑1 (MQ=255)

                                                                  aTTGATGTCCACGGTGGTCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAg                                                                                     <  1:940673/150‑1 (MQ=255)

                                                                    tGATGTCCACGGTGGTCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTTGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAgct                                                                                   >  1:1343946/1‑150 (MQ=255)

                                                                    tGATGTCCACGGTGGTCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAgct                                                                                   >  2:119896/1‑150 (MQ=255)

                                                                     gATGTCCACGGTGGTCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGTCATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAgctt                                                                                  >  2:248913/1‑150 (MQ=255)

                                                                       tGTCTACGGTGGTCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAgcttgc                                                                                <  2:1352985/150‑1 (MQ=255)

                                                                       tGTCCACGGTGGTCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCATTTCTACTTTGCTTTCCATCGCTCGCTTTATTTCGTTTTCGACTTCAAAGAGATCTAAAATCGTCTTGGGTCCTGTGTAATTTTCCAgattgc                                                                                <  2:725913/150‑1 (MQ=255)

                                                                       tGTCCACGGTGGTCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTCCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAgcttgc                                                                                <  1:22346/150‑1 (MQ=255)

                                                                          ccACGGTGGTCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGcttgcttg                                                                             <  2:2048572/150‑1 (MQ=255)

                                                                          ccACGGTGGTCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGcttgcttg                                                                             >  2:2368590/1‑150 (MQ=255)

                                                                          ccACGGTGGTCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGcttgcttg                                                                             <  1:1110053/150‑1 (MQ=255)

                                                                           cACGGTGGTCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGt                                                                            >  1:120855/1‑150 (MQ=255)

                                                                             cGGTGGTCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCa                                                                          <  2:25914/150‑1 (MQ=255)

                                                                             cGGTGGTCAACGATTCAGTCTAGTAGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCAGTGTAATGCTCCAGCTTGCTTGTCa                                                                          >  2:1169538/1‑150 (MQ=255)

                                                                                  gtCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTCTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCg                                                                     >  2:1562524/1‑150 (MQ=255)

                                                                                    cATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGGTCGCTGGATTTCGTTTTCGACATCATAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATGTCggg                                                                   <  1:1595048/150‑1 (MQ=255)

                                                                                    cATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCggg                                                                   >  2:399918/1‑150 (MQ=255)

                                                                                      tCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGaa                                                                 >  1:2045964/1‑150 (MQ=255)

                                                                                       cGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGAAt                                                                >  2:2322494/1‑150 (MQ=255)

                                                                                         cTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGAATGt                                                              <  2:985933/150‑1 (MQ=255)

                                                                                         cTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGAATGt                                                              >  2:1226376/1‑150 (MQ=255)

                                                                                             ggTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGAATGTACTc                                                          >  2:1638275/1‑150 (MQ=255)

                                                                                              gTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGAATGTACTcc                                                         <  2:2213631/150‑1 (MQ=255)

                                                                                                cTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGAATGTACTCCGa                                                       >  2:620782/1‑150 (MQ=255)

                                                                                                cTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGAATGTACTCCGa                                                       >  1:320807/1‑150 (MQ=255)

                                                                                                  ggTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGAATGTACTCCGAgg                                                     <  1:521835/150‑1 (MQ=255)

                                                                                                  agacGATAAAGAGATAACCACCGGAATTCAGATATACTTTGCGTTCCAGCGCTCGCTGTATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGAATGTACTCCGAgg                                                     >  2:2469334/4‑150 (MQ=255)

                                                                                                    tCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGAATGTACTCCGAGGTg                                                   >  1:2127460/1‑150 (MQ=255)

                                                                                                     cGATAATGAGATAACCACTGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGAATGTACTCCGAGGTGa                                                  >  1:154393/1‑150 (MQ=255)

                                                                                                       aTAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTAGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCACGGGAATGTACTCCGAGGTGAAc                                                >  2:1911886/1‑150 (MQ=255)

                                                                                                            tagaTAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGAATGTACTCCGAGGTGAACTCAAg                                           <  1:1041379/149‑1 (MQ=255)

                                                                                                            gagaTAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGAATGTACTCCGAGGTGAACTCAAg                                           <  1:399918/150‑1 (MQ=255)

                                                                                                                tAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGAATGTACTCCGAGGTGAACTCAAGTAAc                                       >  2:582471/1‑150 (MQ=255)

                                                                                                                tAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGAATGTACTCCGAGGTGAACTCAAGTAAc                                       <  2:1954794/150‑1 (MQ=255)

                                                                                                                  accaccGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGAATGTACTCCGAGGTGAACTCAAGTAACGc                                     <  2:1773980/150‑1 (MQ=255)

                                                                                                                      ccGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGAATGTACTCCGAGGTGAACTCAAGTAACGCTTCg                                 <  1:1194739/150‑1 (MQ=255)

                                                                                                                      ccGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGAATGTACTCCGAGGTGAACTCAAGTAACGCTTCg                                 <  2:338744/150‑1 (MQ=255)

                                                                                                                           tttCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGAATGTACTCCGAGGTGAACTCAAGTAACGCTTCGTAAGt                            >  1:2187775/1‑150 (MQ=255)

                                                                                                                            ttCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGAATGTACTCCGAGGTGAACTCAAGTAACGCTTCGTAAGTc                           >  2:219235/1‑150 (MQ=255)

                                                                                                                              cAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGAATGTACTCCGAGGTGAACTCAAGTAACGCTTCGTAAGTCAg                         <  2:476016/150‑1 (MQ=255)

                                                                                                                                 ttCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGAATGTACTCCGAGGTGAACTCAAGTAACGCTTCGTAAGTCAGGCg                      <  2:388864/150‑1 (MQ=255)

                                                                                                                                  tCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGAATGTACTCCGAGGTGAACTCAAGTAACGCTTCGTAAGTCAGGCGt                     >  2:14284/1‑150 (MQ=255)

                                                                                                                                   cTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGAATGTACTCCGAGGTGAACTCAAGTAACGCTTCGTAAGTCAGGCGTg                    <  1:2175734/150‑1 (MQ=255)

                                                                                                                                       tttGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGAATGTACTCCGAGGTGAACTCAAGTAACGCTTCGTAAGTCAGGCGTGAGTc                <  1:1927441/150‑1 (MQ=255)

                                                                                                                                       aacgaaatCCAGCGAGCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGAATGTACTCCGAGGTGAACTCAAGTAACGCTTCGTAAGTCAGGCGTGAGTc                >  2:2324169/8‑150 (MQ=255)

                                                                                                                                        acgaaatCCAGCGAGCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGAATGTACTCCGAGGTGAACTCAAGTAACGCTTCGTAAGTCAGGCGTGAGTCa               >  1:1065333/7‑150 (MQ=255)

                                                                                                                                         cgaaatCCAGCGAGCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGAATGTACTCCGAGGTGAACTCAAGTAACGCTTCGTAAGTCAGGCGTGAGTCaa              >  1:1739322/6‑150 (MQ=255)

                                                                                                                                         cgaaatCCAGCGAGCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGAATGTACTCCGAGGTGAACTCAAGTAACGCTTCGTAAGTCAGGCGTGAGTCaa              >  2:1395383/6‑150 (MQ=255)

                                                                                                                                           cGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGAATGTACTCCGAGGTGAACTCAAGTAACGCTTCGTAAGTCAGGCGTGAGTCAAcg            >  2:1853984/1‑150 (MQ=255)

                                                                                                                                           cGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGAATGTACTCCGAGGTGAACTCAAGTAACGCTTCGTAAGTCAGGCGTGAGTCAAcg            <  1:1856353/150‑1 (MQ=255)

                                                                                                                                           aaatCCAGCGAGCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGAATGTACTCCGAGGTGAACTCAAGTAACGCTTCGTAAGTCAGGCGTGAGTCAAcg            >  2:1107692/4‑150 (MQ=255)

                                                                                                                                            gTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGAATGTACTCCGAGGTGAACTCAAGTAACGCTTCGTAAGTCAGGCGTGAGTCAAcgc           >  1:664584/1‑150 (MQ=255)

                                                                                                                                            aatCCAGCGAGCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGAATGTACTCCGAGGTGAACTCAAGTAACGCTTCGTAAGTCAGGCGTGAGTCAAcgc           >  2:1257048/3‑150 (MQ=255)

                                                                                                                                             ttCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGAATGTACTCCGAGGTGAACTCAAGTAACGCTTCGTAAGTCAGGCGTGAGTCAAcgcg          <  2:1531192/150‑1 (MQ=255)

                                                                                                                                                 aGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGAATGTACTCCGAGGTGAACTCAAGTAACGCTTCGTAAGTCAGGCGTGAGTCAACGCGAATg      >  2:677277/1‑150 (MQ=255)

                                                                                                                                                 aGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGAATGTACTCCGAGGTGAACTCAAGTAACGCTTCGTAAGTCAGGCGTGAGTCAACGCGAATg      >  1:136994/1‑150 (MQ=255)

                                                                                                                                                 aGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGAATGTACTCCGAGGTGAACTCAAGTAACGCTTCGTAAGTCAGGCGTGAGTCAACGCGAATg      <  2:498397/150‑1 (MQ=255)

                                                                                                                                                 aGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGAATGTACTCCGAGGTGAACTCAAGTAACGCTTCGTAAGTCAGGCGTGAGTCAACGCGAATg      <  2:154393/150‑1 (MQ=255)

                                                                                                                                                  gcgcTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGAATGTATTCCGAGGTGAACTCAAGTAACGCTTCGTAAGTCAGGCGTGAGTCAACGCGAATGc     <  2:873229/150‑1 (MQ=255)

                                                                                                                                                    gctcgctGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGAATGTACTCCGAGGTGAACTCAAGTAACGCTTCGTAAGTCAGGCGTGAGTCAACGCGAATGCgg   <  1:2368590/150‑1 (MQ=255)

CGTCGCTTCAATATTGGTATTGAAAATGGTGTCGTCCAGATTGCGATGACCGACAAACGCTCCGGTATTGATGTCCACGGTGGTCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGAATGTACTCCGAGGTGAACTCAAGTAACGCTTCGTAAGTCAGGCGTGAGTCAACGCGAATGCGGT  >  NC_000913/3396812‑3397110

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