BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,720,793	G→T	5.2%	Q38H (CAG→CAT)	insJ →	IS150 transposase A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,720,793	0	G	T	5.2%	481.4 / 11.8	135	Q38H (CAG→CAT)	insJ	IS150 transposase A
Reads supporting (aligned to +/- strand): ref base G (59/69); new base T (2/5); total (61/74)
Fisher's exact test for biased strand distribution p-value = 4.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.44e-01

Only 100 of 141 total aligned reads displayed.

ATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGGCTGCTGCTC  >  NC_000913/3720645‑3720941

aTTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAg                                                                                                                                                      <  3:420308/149‑1 (MQ=255)

 ttttGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAgg                                                                                                                                                     >  3:72244/1‑149 (MQ=255)

 ttttGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTAGTGTCCCCCGAACCCAgg                                                                                                                                                     >  3:167811/1‑149 (MQ=255)

       tgGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCACGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGa                                                                                                                                               <  4:875418/149‑1 (MQ=255)

          aCGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACAt                                                                                                                                            <  2:45450/149‑1 (MQ=255)

          aCGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACAt                                                                                                                                            >  2:627785/1‑149 (MQ=255)

            gATGAAAGAGAATAGCCCCTAATATGTCAAAGCCAATATACCCTATTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGATACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATgg                                                                                                                                          <  2:237342/149‑1 (MQ=255)

               gAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGtt                                                                                                                                       <  1:146558/149‑1 (MQ=255)

               gAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGtt                                                                                                                                       <  1:55297/149‑1 (MQ=255)

                aaaTGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCATGTCCTGACCTGGGttc                                                                                                                                      <  2:201677/149‑2 (MQ=255)

                 aaTGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGc                                                                                                                                     >  4:171907/1‑149 (MQ=255)

                  aTGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGcc                                                                                                                                    >  1:435222/1‑149 (MQ=255)

                    ggAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCATGTCCTGACCTGGGttcgggg                                                                                                                                  <  3:168504/149‑6 (MQ=255)

                        tAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTAt                                                                                                                              >  1:615635/1‑149 (MQ=255)

                            cccTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGaaa                                                                                                                          >  2:520017/1‑149 (MQ=255)

                                aaTATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACa                                                                                                                      >  4:362863/1‑149 (MQ=255)

                                  tatGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATg                                                                                                                    <  1:627785/149‑1 (MQ=255)

                                  tatGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCACCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATg                                                                                                                    <  4:167811/149‑1 (MQ=255)

                                   atGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATgg                                                                                                                   >  2:483095/1‑149 (MQ=255)

                                        aaaGCCAAAATTCCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGaaa                                                                                                              >  3:597775/1‑149 (MQ=255)

                                        aaaGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGaaa                                                                                                              >  3:916943/1‑149 (MQ=255)

                                           gCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCATGTCCTGACCTGGGt                                                                                                                                        >  2:372441/1‑120 (MQ=255)

                                            ccAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAgg                                                                                                          <  3:341438/149‑1 (MQ=255)

                                              aaaaTACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGtt                                                                                                        >  4:307398/1‑149 (MQ=255)

                                              aaaaTACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGtt                                                                                                        <  1:36460/149‑1 (MQ=255)

                                              aaaaTACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGtt                                                                                                        <  4:903028/149‑1 (MQ=255)

                                               aaaTACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGttt                                                                                                       <  1:470756/149‑1 (MQ=255)

                                               aaaTACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGttt                                                                                                       >  2:425629/1‑149 (MQ=255)

                                                 aTACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATgg                                                                                                                                          <  3:546134/112‑1 (MQ=255)

                                                 aTACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATgg                                                                                                                                          >  4:546134/1‑112 (MQ=255)

                                                 aTACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACGTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTaa                                                                                                     <  2:92290/149‑1 (MQ=255)

                                                     ccTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCTGGTCAGGAGATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTcc                                                                                                 >  3:797542/1‑149 (MQ=255)

                                                       ttttGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCa                                                                                               >  2:914460/1‑149 (MQ=255)

                                                                    cTTGATGTCGTGAATCACTACTTCACAACTGATGATGGGTACAGGATCATATCGGTACGGTTTGGTGTCCCCCGAACACAGGTCAGGACATGGGTTGCCCTCTCTGAAAAACTTGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGt                                                                                  >  2:390651/1‑149 (MQ=255)

                                                                    cTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGt                                                                                  >  1:253833/1‑149 (MQ=255)

                                                                    cTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGt                                                                                  <  1:434914/149‑1 (MQ=255)

                                                                      tGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTa                                                                                <  3:482465/149‑1 (MQ=255)

                                                                       gAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAg                                                                               >  4:465234/1‑149 (MQ=255)

                                                                           tCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCt                                                                           >  4:224328/1‑149 (MQ=255)

                                                                           tCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCt                                                                           <  2:435222/149‑1 (MQ=255)

                                                                             gTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGa                                                                         <  1:725278/149‑1 (MQ=255)

                                                                              tGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGAt                                                                        <  2:392822/149‑1 (MQ=255)

                                                                              tGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGAt                                                                        <  2:716450/149‑1 (MQ=255)

                                                                                aaTCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATcc                                                                      <  2:574327/149‑1 (MQ=255)

                                                                                 aTCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAgg                                                                                                                                                     <  2:759577/69‑1 (MQ=255)

                                                                                  tCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAAc                                                                                                                       <  2:156973/98‑1 (MQ=255)

                                                                                  tCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCag                                                                    >  3:131514/1‑149 (MQ=255)

                                                                                   cACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCaga                                                                   <  4:879527/149‑1 (MQ=255)

                                                                                      tactTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGtt                                                                >  2:458056/1‑149 (MQ=255)

                                                                                      tactTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGtt                                                                >  4:439114/1‑149 (MQ=255)

                                                                                         ttCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCg                                                             >  2:184486/1‑149 (MQ=255)

                                                                                         ttCACAACTGAGGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCg                                                             >  1:694627/1‑149 (MQ=255)

                                                                                            acaACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTAt                                                          <  2:315977/149‑1 (MQ=255)

                                                                                            acaACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTAt                                                          >  2:274930/1‑149 (MQ=255)

                                                                                             caACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTAtt                                                         <  2:192191/149‑1 (MQ=255)

                                                                                              aaCTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGt                                                            >  3:263763/1‑145 (MQ=255)

                                                                                              aaCTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTa                                                        >  3:728692/1‑149 (MQ=255)

                                                                                                   atgatgGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTc                                                   <  1:204986/149‑1 (MQ=255)

                                                                                                     gatgGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGt                                                 <  3:153873/149‑1 (MQ=255)

                                                                                                     gatgGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGt                                                 <  1:640835/149‑1 (MQ=255)

                                                                                                      atgGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGc                                                                                    >  3:187899/1‑113 (MQ=255)

                                                                                                          ttACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAg                                            <  1:731557/149‑1 (MQ=255)

                                                                                                          ttACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAg                                            <  1:158877/149‑1 (MQ=255)

                                                                                                          ttACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAg                                            <  4:220754/149‑1 (MQ=255)

                                                                                                            aCAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCt                                          <  3:815961/149‑1 (MQ=255)

                                                                                                             cAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGt                                                    <  2:105754/138‑1 (MQ=255)

                                                                                                             cAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGt                                                    >  1:105754/1‑138 (MQ=255)

                                                                                                              aGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCaga                                                                   >  2:156300/1‑122 (MQ=255)

                                                                                                              aGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCaga                                                                   <  1:156300/122‑1 (MQ=255)

                                                                                                              aGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCtgt                                        <  3:106040/149‑1 (MQ=255)

                                                                                                                gATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGa                                      >  3:458399/1‑149 (MQ=255)

                                                                                                                 atcatcTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGAt                                     <  3:36267/149‑1 (MQ=255)

                                                                                                                 atcatcTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGAt                                     <  4:458399/149‑1 (MQ=255)

                                                                                                                  tcatcTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATc                                    <  1:160622/149‑1 (MQ=255)

                                                                                                                    atcTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAGAAAGATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTGTTAAGGTCGTGAAAGCTGTGATCGa                                  >  3:904274/1‑149 (MQ=255)

                                                                                                                     tctcGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGa                                  <  1:790279/148‑1 (MQ=255)

                                                                                                                     tctcGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGa                                  >  2:790279/1‑148 (MQ=255)

                                                                                                                       tcGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGagca                               >  3:210356/1‑149 (MQ=255)

                                                                                                                             cGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACAt                         <  1:520017/149‑1 (MQ=255)

                                                                                                                               ttttGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGt                       >  1:71696/1‑149 (MQ=255)

                                                                                                                               ttttGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGt                       <  4:210356/149‑1 (MQ=255)

                                                                                                                                       tCCCCCGGACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAAt               >  2:617899/1‑149 (MQ=255)

                                                                                                                                       tCCCCCGAACCCAGGTCAGGACATGGGTTGCCctct                                                                                                                                >  1:451859/1‑36 (MQ=255)

                                                                                                                                        cccccGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATc              <  1:352949/149‑1 (MQ=255)

                                                                                                                                         ccccGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCa             >  4:730200/1‑149 (MQ=255)

                                                                                                                                         ccccGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCa             >  4:915936/1‑149 (MQ=255)

                                                                                                                                         ccccGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCa             >  2:224631/1‑149 (MQ=255)

                                                                                                                                          cccGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAg                                                                               >  2:699331/1‑82 (MQ=255)

                                                                                                                                          cccGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAg            >  2:731677/1‑149 (MQ=255)

                                                                                                                                           tcGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAgg           <  4:797542/148‑1 (MQ=255)

                                                                                                                                           ccGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAgg           >  1:515899/1‑149 (MQ=255)

                                                                                                                                             gAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAAGCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGgct         <  4:523554/149‑1 (MQ=255)

                                                                                                                                              aaCCCATGTCCTGACCTGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTc                                                   <  1:800785/106‑1 (MQ=255)

                                                                                                                                              aaCCCATGTCCTGACCTGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGATTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTc                                                   >  2:800785/1‑106 (MQ=255)

                                                                                                                                              aaCCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGgctg        <  1:483095/149‑1 (MQ=255)

                                                                                                                                               aCCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGgctgc       >  2:753325/1‑149 (MQ=255)

                                                                                                                                               aCCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGgctgc       >  2:698244/1‑149 (MQ=255)

                                                                                                                                                 ccAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGgctgctg     >  1:256278/1‑149 (MQ=255)

                                                                                                                                                  cAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGgctgctgc    <  3:895325/149‑1 (MQ=255)

                                                                                                                                                    ggTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGGCTGCTGCTc  <  4:72244/149‑1 (MQ=255)

ATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGGCTGCTGCTC  >  NC_000913/3720645‑3720941

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 16 ≤` `ATCG/ATCG` `< 34 ≤` `ATCG/ATCG` `< 39 ≤` `ATCG/ATCG` `< 40 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`