BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	953,720	A→T	7.4%	N249K (AAT→AAA)	focA ←	formate channel FocA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	953,720	0	A	T	7.4%	314.1 / 10.9	103	N249K (AAT→AAA)	focA	formate channel FocA
Reads supporting (aligned to +/- strand): ref base A (47/48); new base T (3/5); total (50/53)
Fisher's exact test for biased strand distribution p-value = 7.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

Only 100 of 103 total aligned reads displayed.

GTGGATTTTTTATTTACTGCGTACTTCGACAACCATTAATGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACAT  >  NC_000913/953574‑953866

gTGGATTTTTTATTTACTGCGTACTTCGACAACCATTAATGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATtc                                                                                                                                                 >  1:497583/1‑150 (MQ=255)

gTGGATTTTTTATTTACTGCGTACTTCGACAACCATTAATGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATtc                                                                                                                                                 <  2:1294954/150‑1 (MQ=255)

 tGGATTTTTTATTTACTGCGTACTTCGACAACCATTAATGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATtca                                                                                                                                                >  2:1012487/1‑150 (MQ=255)

   taTTTTTTATTTACTGCGTACTTCGACAACCATTAATGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATtcatc                                                                                                                                              <  1:712247/149‑1 (MQ=255)

     ttttttATTTACTGCGTACTTCGACAACCATTAATGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCAc                                                                                                                                            <  2:894244/150‑1 (MQ=255)

      tttttATTTACTGCGTACTTCGACAACCATTAATGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAATTTCATCAct                                                                                                                                           <  1:325168/150‑2 (MQ=255)

      tttttATTTACTGCGTACTTCGACAACCATTAATGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACg                                                                                                                                           <  2:478480/150‑1 (MQ=255)

       ttttATTTACTGCGTACTTCGACAACCATTAATGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAATTTCATCActg                                                                                                                                          <  2:1470531/150‑3 (MQ=255)

          tatttaCTGCGTACTTCGACAACCATTAATGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAATTTCATCActgata                                                                                                                                       <  1:93275/150‑6 (MQ=255)

          tatttaCTGCGTACTTCGACAACCATTAATGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCa                                                                                                                                       <  2:102318/150‑1 (MQ=255)

              taCTGCGTACTTCGACAACCATTAATGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTg                                                                                                                                   <  1:1568454/150‑1 (MQ=255)

               aCTGCGTACTTCGACAACCATTAATGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTga                                                                                                                                  <  1:965607/150‑1 (MQ=255)

               aCTGCGTACTTCGACAACCATTAATGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTga                                                                                                                                  <  1:959890/150‑1 (MQ=255)

                cTGCGTACTTCGACAACCATTAATGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTgag                                                                                                                                 <  1:940022/150‑1 (MQ=255)

                 tGCGTACTTCGACAACCATTAATGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTgaga                                                                                                                                >  2:864602/1‑150 (MQ=255)

                     tACTTCGACAACCATTAATGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGaaaaa                                                                                                                            >  1:128101/1‑150 (MQ=255)

                     tACTTCGACAACCATTAATGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGaaaaa                                                                                                                            >  1:632121/1‑150 (MQ=255)

                       cTTCGACAACCATTAATGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAGTCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAAtt                                                                                                                          >  2:1311532/1‑150 (MQ=255)

                           gACAACCATTAATGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTcc                                                                                                                      <  1:1254323/150‑1 (MQ=255)

                                ccATTAATGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAATTTCATCActg                                                                                                                                          >  2:450400/1‑123 (MQ=255)

                                ccATTAATGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAATTTCATCActg                                                                                                                                          <  1:450400/125‑3 (MQ=255)

                                ccATTAATGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGc                                                                                                                 >  1:999488/1‑150 (MQ=255)

                                 cATTAATGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCa                                                                                                                >  1:557548/1‑150 (MQ=255)

                                 cATTAATGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCa                                                                                                                <  1:311355/150‑1 (MQ=255)

                                 cATTAATGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCa                                                                                                                <  2:142290/150‑1 (MQ=255)

                                     aaTGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAAc                                                                                                            >  2:1371109/1‑150 (MQ=255)

                                      aTGGTGGTCGTTTTCAAGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAAcc                                                                                                           <  2:727146/150‑1 (MQ=255)

                                          tggtCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGCGAGAAAAATTTTCCGGTGCAGAACCGACt                                                                                                       >  2:970179/1‑150 (MQ=255)

                                           ggtCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCa                                                                                                                                       <  1:67605/117‑1 (MQ=255)

                                           ggtCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTg                                                                                                      >  1:1007276/1‑150 (MQ=255)

                                           ggtCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACAGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCa                                                                                                                                       >  2:67605/1‑117 (MQ=255)

                                            gtCGTTTTCTCGCAGGTAAATGACCCAGTATGTCAAGCCTTCCAACAAACCACGACGGATAATGTTGCCGATCGTTACCGGAATCAGGTTTTCAGTGATGAATTTCATCACGGTCAGGTGAGAAAAATTTTCTGGTGCAGAACCGTCTGc                                                                                                     <  2:102664/150‑1 (MQ=255)

                                            gtCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGc                                                                                                     <  1:1371109/150‑1 (MQ=255)

                                             tCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCg                                                                                                    >  2:628341/1‑150 (MQ=255)

                                             tCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCg                                                                                                    <  1:1500148/150‑1 (MQ=255)

                                              cGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATTGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCgg                                                                                                   <  1:727709/150‑1 (MQ=255)

                                               gTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGt                                                                                                  <  1:619137/150‑1 (MQ=255)

                                               cgaagtaCGCA‑GTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTc                                                                                                 >  1:1305329/7‑150 (MQ=255)

                                                  ttCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCa                                                                                               >  1:650186/1‑150 (MQ=255)

                                                   tCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCaa                                                                                              <  1:1012487/150‑1 (MQ=255)

                                                    cACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCaaa                                                                                             >  2:1357371/1‑150 (MQ=255)

                                                      cGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAAt                                                                                           >  1:1192843/1‑150 (MQ=255)

                                                      cGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAAt                                                                                           >  1:766322/1‑150 (MQ=255)

                                                      cGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATAGTAACCGTAATCCGGTTATCCGTGAAGAAAGTGATCACGGACAGGTGAGAACAATTTGCCTGTGCCGACCCCACTGCGCTCCAAAAt                                                                                           >  2:1290935/1‑150 (MQ=255)

                                                           gTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCgg                                                                                      <  1:100610/150‑1 (MQ=255)

                                                            tAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCggg                                                                                     <  1:545459/150‑1 (MQ=255)

                                                              aaTGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCTCGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGa                                                                                   >  2:1332158/1‑150 (MQ=255)

                                                                tGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATg                                                                                 >  2:568481/1‑150 (MQ=255)

                                                                 gACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGc                                                                                <  2:1007276/150‑1 (MQ=255)

                                                                  aCCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCg                                                                               >  2:653279/1‑150 (MQ=255)

                                                                   cccAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAGCCGACTGCGGTCCAAAATTCCGGGGATGCGa                                                                              <  1:1309346/150‑1 (MQ=255)

                                                                   cccAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGa                                                                              <  1:653279/150‑1 (MQ=255)

                                                                         aTGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGTGGTCCAAAATTCCGGGGATGCGAAGTcgc                                                                        >  2:933331/1‑150 (MQ=255)

                                                                          tGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTcgcg                                                                       <  1:1404314/150‑1 (MQ=255)

                                                                           gTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCgg                                                                      >  2:895569/1‑150 (MQ=255)

                                                                              aacccaaccAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGAtt                                                                   >  2:421997/1‑150 (MQ=255)

                                                                               acccaaccAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTa                                                                  >  2:862307/1‑150 (MQ=255)

                                                                                cccaaccAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTAc                                                                 <  1:421997/150‑1 (MQ=255)

                                                                                 ccaaccaaCAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACa                                                                >  1:1076327/1‑150 (MQ=255)

                                                                                 ccaaccaaCAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACa                                                                <  2:1425590/150‑1 (MQ=255)

                                                                                 ccaaccaaCAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACa                                                                >  2:684981/1‑150 (MQ=255)

                                                                                  caaccaacAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACaa                                                               <  2:513387/150‑1 (MQ=255)

                                                                                  caaccaacAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACaa                                                               <  1:864602/150‑1 (MQ=255)

                                                                                  caaccaacAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGA‑TACAAt                                                              <  2:1008339/150‑1 (MQ=255)

                                                                                  caaccaacAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGCGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACaa                                                               <  1:970179/150‑1 (MQ=255)

                                                                                   aaccaacAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAAt                                                              <  1:577425/150‑1 (MQ=255)

                                                                                   aaccaacAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAAt                                                              <  1:1529865/150‑1 (MQ=255)

                                                                                     ccaacAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATAc                                                            <  1:159054/150‑1 (MQ=255)

                                                                                     ccaacAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGACGCGGATTACAATAc                                                            >  2:560697/1‑150 (MQ=255)

                                                                                      caacaaACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATAcc                                                           <  2:632121/150‑1 (MQ=255)

                                                                                      caacaaACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATAcc                                                           <  2:422946/150‑1 (MQ=255)

                                                                                          aaaCCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATc                                                       >  2:450385/1‑150 (MQ=255)

                                                                                           aaCCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAATTTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCg                                                      >  1:560489/1‑150 (MQ=255)

                                                                                               accaccGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGAt                                                  >  1:52711/1‑150 (MQ=255)

                                                                                                   ccGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCATa                                              >  1:122674/1‑150 (MQ=255)

                                                                                                         aTGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCATAAACATg                                        >  1:855545/1‑150 (MQ=255)

                                                                                                          tGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCATAAACATGt                                       >  1:355200/1‑150 (MQ=255)

                                                                                                               ccGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCATAAACATGTTTGCg                                  >  2:1372696/1‑150 (MQ=255)

                                                                                                               ccGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCATAAACATGTTTGCg                                  >  1:1471404/1‑150 (MQ=255)

                                                                                                                   tCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCATAAACATGTTTGCGATAc                              <  1:652603/150‑1 (MQ=255)

                                                                                                                     gTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCATAAACATGTTTGCGATACtg                            <  2:128101/150‑1 (MQ=255)

                                                                                                                       aaCCGGAATCAGGTTATCAGTGATGAAATTCATCTCGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCATAAACATGTTTGCGATACtgtg                          <  1:1332158/150‑1 (MQ=255)

                                                                                                                        aCCGGAATCAGTTTATCATTGATGAAATTCATCACGGTCAGGTGATAAAAATTTTCCGGTGCAGAAACGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCTGGATCATAAACATGTTTGCGATACTGTGc                         <  2:557548/150‑1 (MQ=255)

                                                                                                                        aCCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGAGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGc                         >  2:1309507/1‑150 (MQ=255)

                                                                                                                              aTCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCaaaa                   >  1:994896/1‑150 (MQ=255)

                                                                                                                              aTCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCaaaa                   >  1:331961/1‑150 (MQ=255)

                                                                                                                                cAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAAcc                 <  2:331961/150‑1 (MQ=255)

                                                                                                                                   gTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCt              >  1:736356/1‑150 (MQ=255)

                                                                                                                                     tATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTgg            <  1:684981/150‑1 (MQ=255)

                                                                                                                                     tATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTgg            <  1:54136/150‑1 (MQ=255)

                                                                                                                                      gaccgTGATGAATTTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCCTCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGc           >  1:1317059/5‑150 (MQ=255)

                                                                                                                                      aTCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGc           <  2:766322/150‑1 (MQ=255)

                                                                                                                                       tCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGca          <  1:568481/150‑1 (MQ=255)

                                                                                                                                       tCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGca          <  1:471143/150‑1 (MQ=255)

                                                                                                                                        ccgTGATGAATTTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAACCCGCTGGcaa         >  1:861331/3‑150 (MQ=255)

                                                                                                                                          gTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGcaaca       >  2:146250/1‑150 (MQ=255)

                                                                                                                                          gTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGcaaca       >  1:5995/1‑150 (MQ=255)

                                                                                                                                           tgatgaAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGTGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGcaacaa      <  1:933331/150‑1 (MQ=255)

                                                                                                                                            gatgaAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCaacaaa     <  1:1311532/150‑1 (MQ=255)

                                                                                                                                               gaAATTAATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCTACTGAGGTCCAAAATTCCGGGGATCCGAAGTCGCGGATTACAATACCCATCGGGATCATAAACATGTTTTCGATACAGTGCTCAAAACCGCTGGCAACAAACAt  >  1:191341/1‑150 (MQ=255)

GTGGATTTTTTATTTACTGCGTACTTCGACAACCATTAATGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACAT  >  NC_000913/953574‑953866

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: `‑`