BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,597,364	T→C	pseudogene (624/3861 nt)	yneO ←	pseudogene, AidA homolog

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,597,364	0	T	C	96.0%	1143.9 / 18.4	335	pseudogene (624/3861 nt)	yneO	pseudogene, AidA homolog
Reads supporting (aligned to +/- strand): ref base T (7/6); new base C (137/178); total (149/186)
Fisher's exact test for biased strand distribution p-value = 5.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.17e-02

Only 100 of 335 total aligned reads displayed.

GTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGCCTAAAATGGCGACGATTGAACTCGTGACGTAACCCTTATCCGTAATATTTAATGAACCTGTGCCATAGCT  >  NC_000913/1597216‑1597510

  aCCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATCGAg                                                                                                                                                 <  6:1010684/150‑1 (MQ=255)

   ccGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATCGAggt                                                                                                                                               <  7:1345960/151‑2 (MQ=255)

   ccGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATCGAGGa                                                                                                                                               <  8:1437319/151‑1 (MQ=255)

    cGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTACCCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTGGCTTACCAATTTGAAATTCAATCGAGGa                                                                                                                                               >  5:1368916/1‑150 (MQ=255)

     gTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATCGAGGAAt                                                                                                                                             >  1:33693/1‑151 (MQ=255)

       gatgaCAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATCGAGGAATc                                                                                                                                            <  1:1622452/150‑1 (MQ=255)

        atgaCAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCAGCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATCGAGGAATCa                                                                                                                                           <  4:944547/150‑1 (MQ=255)

         tgaCAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATCGAGGAATCatt                                                                                                                                         >  6:1489821/1‑151 (MQ=255)

          gaCAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCGACCGATAATAGTATTTTAAGCCGTAACGAACCCCAAATCGCGAATAGTCGCCTACCCAGTTACTTTATTACCAAATTGAAATTCAATCGAGGTATCatt                                                                                                                                         >  8:1517923/1‑150 (MQ=255)

           aCAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGACACCGATAATCGTATTTTCAGCCGTAACAAACCCCCCCTCGCGAATAGTCGCCTACGCAGATCCGTGATTAGAAATTTGAAATTAAATCGAGGAATCatta                                                                                                                                        >  2:983360/1‑150 (MQ=255)

            cAGAGTCTTGATCCTGCACATTCAGTGTTCCGATACCTGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATCGAGGAATCattatt                                                                                                                                      <  3:979386/151‑1 (MQ=255)

             agagTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATCGAGGAATCATTAttt                                                                                                                                     >  3:902862/1‑151 (MQ=255)

              gagTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATCGAGGAATCATTAtttt                                                                                                                                    <  1:1393369/151‑1 (MQ=255)

                    tGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCATCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTGCTTGATTACCAATTTGAAATTCAATGGAGGAATCATTTTTTTttattt                                                                                                                              >  1:1155807/1‑150 (MQ=255)

                     gATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTACCCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATCGAGGAATCATTATTTTttatt                                                                                                                               >  3:1333382/1‑149 (MQ=255)

                     gATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTACCCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATCGAGGAATCATTATTTTTTATTAg                                                                                                                             <  3:27609/151‑1 (MQ=255)

                     gATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATCGAGGAATCATTATTTTttatta                                                                                                                              <  1:1261369/150‑1 (MQ=255)

                       tCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTGCCCAGTTCCTTGATTACCAATTTGAAATTCAATCGAGCAATCATTATTTTTTATTAGcc                                                                                                                           >  3:71233/1‑151 (MQ=255)

                        ccTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCCACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAATTTCAATCGAGGAATCATTATTTTTTATTAGcc                                                                                                                           >  3:225782/1‑150 (MQ=255)

                        ccTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATCGAGGAATCATTATTTTTTATTAGCCa                                                                                                                          <  2:1289881/151‑1 (MQ=255)

                           gCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATCGAGGAATCATTATTTTTTATTAGCCACt                                                                                                                        <  1:1347162/150‑1 (MQ=255)

                               aTTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCACCCCCCCCTCGCGAAGAGTCGCCTCCCCAGGTGCTTGATTAGCAATTTGAAATTCAATCGCGGCATCATTATTTTTTATTAGCCACTCGccc                                                                                                                   >  7:787252/1‑150 (MQ=255)

                                  cAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATCGAGGAATCATTATTTTTTATTAGCCACTCGCCAccc                                                                                                                <  3:359871/151‑1 (MQ=255)

                                        tCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCCCCCCCCCCTCGCGACTAGTCGGCGCCCCCGTTGGGGGATGACCAATTTGATCTTCCATCGCGGCATGAGTATTTTTTAGTAGCCCCTCGGCCCCCGTTTcc                                                                                                          >  8:846470/1‑150 (MQ=255)

                                         ccGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCa                                                                                                          >  2:763957/1‑150 (MQ=255)

                                           gATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACg                                                                                                       <  2:611844/151‑1 (MQ=255)

                                           gATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATCGAGGAATCATTATTTTTTATTAGCGAGTCGCGACCCTTTTCAACg                                                                                                       >  6:846161/1‑151 (MQ=255)

                                           gATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATCGAGGAATCATTATTTTTTAGTAGCCACTCGCCACCCTTTTAAACg                                                                                                       >  4:1076888/1‑151 (MQ=255)

                                            aTACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATCGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGa                                                                                                      <  4:1362593/151‑1 (MQ=255)

                                              aCCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCACCCCCCCCTCGCGAATACTCGCCTGCCGAGTGCCGTGATTACCAAGTTGAAATTCGATCGAGCAATCATGATTTTTTAATAGCCCCTCGCGACGCTTTTCAACGAc                                                                                                     >  5:1088551/1‑150 (MQ=255)

                                               ccGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATCGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACaa                                                                                                   <  1:379589/151‑1 (MQ=255)

                                                cGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACGCCCCCTCGCGATTAGTGGCGTCCCGACGTACGTGAATTCCAATTTGATAATCAATCGCGGGATGATTATTTTTTATTATCCCCTAGCGCACCTTATCAACGACCAc                                                                                                  >  8:564641/1‑151 (MQ=255)

                                                cGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGCATAGTCGCCTCCCCAGTTGCTTGCTTACCAAGTTGAAATTCAATCGCGGCCTCATTATTTTTTTTTTGCCCCTGGCCCCCCTTTTccccgccccc                                                                                                  >  5:431367/1‑142 (MQ=255)

                                                 ggTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATCGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAAcc                                                                                                 <  3:626938/151‑1 (MQ=255)

                                                 ggTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATCGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCt                                                                                                <  4:1219015/151‑1 (MQ=255)

                                                    ggCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATCGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTgg                                                                                              <  8:930836/151‑1 (MQ=255)

                                                       aTTGCCACCGATAATCGTATTTTCAGCCGTAACCACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGcc                                                                                            >  8:1099162/1‑150 (MQ=255)

                                                       aTTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATCGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGcc                                                                                            >  1:474244/1‑150 (MQ=255)

                                                         tGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATCGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACACCCTGGCCAtt                                                                                         >  2:135281/1‑151 (MQ=255)

                                                          gCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATCGCGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTa                                                                                        <  6:578822/151‑1 (MQ=255)

                                                            cACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATCGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACt                                                                                      <  6:531680/151‑1 (MQ=255)

                                                                 aTAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCGTTGATTACCAATTTGAAATTCAATCGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGTCCATTACTGCCCg                                                                                 >  6:1304006/1‑151 (MQ=255)

                                                                   aaTCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATCGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGcc                                                                               <  3:468447/151‑1 (MQ=255)

                                                                    aTCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATCGAGGAATCATTATTTTTTATTAGCCCCTCGCCACCCTTTTCAACGCCAACCTGGCCATTACTGCCCGCCt                                                                              >  1:968821/1‑151 (MQ=255)

                                                                      cGTATTTTCAGCCGTAACCACCCCCCCCTCGCGAATAGTCGCCTCCCCACTTCCATGATTACCAATTTGAAATTCAATCCCCGAATCATTATTTTTTAATAGCCCCTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCtc                                                                             >  1:1477405/1‑149 (MQ=255)

                                                                      cGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATCGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCAGCCTg                                                                             >  4:587341/1‑150 (MQ=255)

                                                                       gTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATCGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGAt                                                                           <  8:478197/151‑1 (MQ=255)

                                                                               aGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGAATACCAATTTGAAATTCAATCCAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGCCTaaa                                                                   <  5:278998/151‑1 (MQ=255)

                                                                                   gTAACCAACCCCCACTCGCTAATAGTCGCCTACCCAGTTCCTTGATTACCAATTTGAAATTCAATCGAGGAATCATTATTTTTTATTAGCCACTAGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGCCTAAAATgg                                                               <  8:1587809/151‑1 (MQ=255)

                                                                                    tAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCGTGATTACCAAATTGAAATTAAATCGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCGTTTAAACGCCCACCTGGCCATTACTGCCCGCCTGATAGCCTAAAATGGc                                                              >  5:354281/1‑151 (MQ=255)

                                                                                     aaccaaccGCCCCTCTCGAATAGTCTCCTCCCCAGTTGCTTGATTACCAATTTGAAATTCAATCGAGGAATGATTATTTTTTATAAGCCACTCGCCAGCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGCCTAAAATGGcg                                                             <  1:164375/151‑1 (MQ=255)

                                                                                        caaccCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATCCAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTATTCAACGACAACCTGGCCATTACTGCCCGCCTCATAGCCTAAAATGGcgacg                                                          >  7:953213/1‑151 (MQ=255)

                                                                                         aaccCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATCGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGCCTAAAATGGcgacg                                                          >  2:1104313/1‑150 (MQ=255)

                                                                                          accCCCCCTCGCGAATAGTCGCCTCACCAGTTCCTTGATTACCAATTTGAAATTCAATCGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGCCTAAAATGGCGACGAt                                                        <  8:731791/151‑1 (MQ=255)

                                                                                              ccccTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATCGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGCCTAAAATGGCGACGATTGa                                                     >  7:247676/1‑150 (MQ=255)

                                                                                               cccTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATCGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGCCTAAAATGGCGACGATTGAAc                                                   <  2:499985/151‑1 (MQ=255)

                                                                                                 cTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATCGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGCCTAAAATGGCGACGATTGAACTc                                                 >  6:1194958/1‑151 (MQ=255)

                                                                                                   cgcgAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATCGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGCCTAAAATGGCGACGATTGAACTc                                                 <  6:1599186/149‑1 (MQ=255)

                                                                                                   cgcgAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATCGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGCCTAAAATGGCGACGATTGAACTCGt                                               >  4:1668456/1‑151 (MQ=255)

                                                                                                    gcgAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATCGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGCCTAAAATGGCGACGATTGAACTCGTg                                              >  7:1225788/1‑151 (MQ=255)

                                                                                                     cgAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATCGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCGGATAGCCTAGAATGGCGACGATTGAACTCGTg                                              >  3:1516466/1‑150 (MQ=255)

                                                                                                      gAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATCGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGCCTAAAATGGCGACGATTGAACTCGTGAc                                            <  3:1107806/151‑1 (MQ=255)

                                                                                                       aaTAGTCGCCTCCCCAGTTCCTTTATTACCAATTTGAAATTCAATCGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGCCTAAAATGGCGACGATTGAACTCGTGACg                                           <  8:1483993/151‑1 (MQ=255)

                                                                                                       aaTAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATCGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGCCTAAAATGGCGACGATTGAACTCGAGACg                                           >  6:212868/1‑151 (MQ=255)

                                                                                                         tAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATCGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGCCTAAAATGGCGACGATTGAACTCGTGACGt                                          >  3:1065444/1‑150 (MQ=255)

                                                                                                           gTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATCGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGCCTAAAATGGCGACGATTGAACTCGTGACGTAAc                                       <  4:472872/151‑1 (MQ=255)

                                                                                                             cGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATCGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGCCTAAAATGGCGACGATTGAACTCGTGACGTAAccc                                     >  7:1283305/1‑151 (MQ=255)

                                                                                                              gCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATCGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGCCTAAAATGGCGACGATTGAACTCGTGACGTAACCCt                                    >  5:635615/1‑151 (MQ=255)

                                                                                                              gCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATCGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGCCTAAAATGGCGACGATTGAACTCGTGACGTAACCCt                                    >  2:1529196/1‑151 (MQ=255)

                                                                                                               ccTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATCGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGCCTAAAATGGCGACGATTGAACTCGTGACGTAACCCtt                                   >  8:1705440/1‑151 (MQ=255)

                                                                                                               ccTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATCGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGCCTAAAATGGCGACGATTGAACTCGTGACGTAACCCtt                                   >  5:1658093/1‑151 (MQ=255)

                                                                                                                cTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATCGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGCCTAAAATGGCGACGATTGAACTCGTGACGTAACCCtt                                   >  7:1109361/1‑150 (MQ=255)

                                                                                                                 cccccAGTTCCTTGATTACCAATTTGAAATTCAATCGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGCCTAAAATGGCGACGATTGAACTCGTGACGTAACCCTTAt                                 <  1:895955/150‑1 (MQ=255)

                                                                                                                  ccccAGTTCCTTGATTACCAATTTGAAATTCAATCGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGCCTAAAATGGCGACGATTGAACTCGTGACGTAACCCTTAt                                 <  2:1023411/150‑1 (MQ=255)

                                                                                                                  ccccAGTTCCTTGATTACCAATTTGAAATTCAATCGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGCCTAAAATGGCGACGATTGAACTCGTGACGTAACCCTTAt                                 <  5:1252599/150‑1 (MQ=255)

                                                                                                                  ccccAGTTCCTTGATTACCAATTTGAAATTCAATCGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGCCTAAAATGGCGACGATTGAACTCGTGACGTAACCCTTATc                                >  4:255549/1‑151 (MQ=255)

                                                                                                                       gTTCCTTGATTACCAATTTGAAATTCAATCGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGCCTAAAATGGCGACGATTGAACTCGTGACGTAACCCTTATCCGTaa                           <  7:1432793/151‑1 (MQ=255)

                                                                                                                        ttCCTTGATTACCAATTTGAAATTCAATCGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGCCTAAAATGGCGACGATTGAACTCGTGACGTAACCCTTATCCGTAat                          >  6:66244/1‑151 (MQ=255)

                                                                                                                        ttCCTTGATTACCAATTTGAAATTCAATCGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGCCTAAAATGGCGACGATTGAACTCGTGACGTAACCCTTATCCGTAat                          >  7:457912/1‑151 (MQ=255)

                                                                                                                          ccTTGATTACCAATTTGAAATTCAATCGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGCCTAAAATGGCGACGATTGAACTCGTGACGTAACCCTTATCCGTAatat                        >  7:842653/1‑151 (MQ=255)

                                                                                                                             tGATTACCAATTTGAAATTCAATCGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGCCTAAAATGGCGACGATTGAACTCGTGACGTAACCCTTATCCGTAATAttt                      <  5:874142/150‑1 (MQ=255)

                                                                                                                               aTTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGCCTAAAATGGCGACGATTGAACTCGTGACGTAACCCTTATCCGTAATATTTAAt                   >  7:1581533/1‑151 (MQ=255)

                                                                                                                               aTTACCAATTTGAAATTCAATCGAGGAATCATTATTTTTTCTTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGCCTAAAATGGCGACGATTGAACTCGTGACGTAACCCTTATCCGTAATATTTAAt                   <  7:565109/151‑1 (MQ=255)

                                                                                                                               aTTACCAATTTGAAATTCAATCGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGCCTAAAATGGCGACGATTGAACTCGTGACGTAACCCTTATCCGTAATATTTAAt                   <  4:112271/151‑1 (MQ=255)

                                                                                                                                ttAGCAATTTGATATTCTATCGAGGAATCATTATTTTTTATTAGCCTCTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGCCTAAAATGGCGACGATTGAACTCGTGACGTAACCCTTATCCGTAATATTTAAt                   <  6:274581/150‑1 (MQ=255)

                                                                                                                                  aCCAATTTGAAATTCAATCGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGCCTAAAATGGCGACGATTGAACTCGTGACGTAACCCTTATCCGTAATATTTAATGaa                <  8:1389251/151‑1 (MQ=255)

                                                                                                                                    cAATTTGAAATTCAATCGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGCCTAAAATGGCGACGATTGAACTCGTGACGTAACGCTTATCCGTAATATTTAATGAAcc              >  6:1420088/1‑151 (MQ=255)

                                                                                                                                     aaTTTGAAATTCAATCGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGCCTAAAATGGCGACGATTGAACTCGTGACGTAACCCTTATCCGTAATATTTAATGAACCt             <  8:1197384/151‑1 (MQ=255)

                                                                                                                                     aaTTTGAAATTCAATCGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGCCTAAAATGGCGACGATTGAACTCGTGACGTAACCCTTATCCGTAATATTTAATGAACCt             >  4:999473/1‑151 (MQ=255)

                                                                                                                                      aTTTGAAATTCAATCGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGCCTAAAATGGCGACGATTGAACTCGTGACGTAACCCTTATCCGTAATATTTAATGAACCtg            >  6:1652153/1‑151 (MQ=255)

                                                                                                                                       tttGAAATTCAATCGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGCCTAAAATGGCGACGATTGAACTCGTGACGTAACCCTTATCCGTAATATTTAATGAACCtgt           >  8:318628/1‑151 (MQ=255)

                                                                                                                                       tttGAAATTCAATCGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGCCTAAAATGGCGACGATTGAACTCGTGACGTAACCCTTATCCGTAATATTTAATGAACCtgt           >  2:221227/1‑151 (MQ=255)

                                                                                                                                         tGAAATTCAATCGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGCCTAAAATGGCGACGATTGAACTCGTGACGTAACCCTTATCCGTAATATTTAATGAACCtgtg          <  1:585388/150‑1 (MQ=255)

                                                                                                                                          gAAATTCAATCGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGCTTAAAATGGCGACGATTGAACTCGTGACGTAACCCTTATCCGTAATATTTAATGAACCTGTGcc        <  4:771619/151‑1 (MQ=255)

                                                                                                                                           aaaTTCAATCGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGCCTAAAATGGCGACGATTGAACTCGTGACGTAACCCTTATCCGTAATATTTAATGAACCTGTGca        >  3:567332/1‑149 (MQ=255)

                                                                                                                                            aaTTCAATCGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGCCTAAAATGGCGACGATTGAACTCGTGACGTAACCCTTATCCGTAATATTTAATGAACCTGTGCCAt      <  7:1018287/151‑1 (MQ=255)

                                                                                                                                              ttCAATCGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGCCTAAAATGGCGACGATTGAACTCGTGACGTAACCCTTATCCGTAATATTTAATGAACCTGTGCCAt      >  4:832563/1‑149 (MQ=255)

                                                                                                                                               tCAATCGAGGAATCATTATTTTTTATTAGCCACTCGCCACCATTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGCCTAAAATGGCGACGATTGAACTCGTGACGTAACCCTTATCCGTAATATTTAATGAACCTGTGCCATAg    <  4:628360/150‑1 (MQ=255)

                                                                                                                                                cAATCGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCGATTACTGCCCGCCTGATAGCCTAAAATGGCGACGATTGAACTCGTGACGTAACCCTTATCCGTAATATTTAATGAACCTGTGCCATAGc   <  1:1022709/150‑1 (MQ=255)

                                                                                                                                                cAATCGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGCCTAAAATGGCGACGATTGAACTCGTGACGTAACCCTTATCCGTAATATTTAATGAACCTGTGCCATAGCt  >  4:1103740/1‑151 (MQ=255)

GTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGCCTAAAATGGCGACGATTGAACTCGTGACGTAACCCTTATCCGTAATATTTAATGAACCTGTGCCATAGCT  >  NC_000913/1597216‑1597510

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