BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	953,720	A→T	11.9%	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	11.9%	276.5 / 24.2	101	N249K (AAT→AAA)	focA	formate channel FocA
Reads supporting (aligned to +/- strand): ref base A (49/40); new base T (8/4); total (57/44)
Fisher's exact test for biased strand distribution p-value = 5.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.93e-01

Only 100 of 101 total aligned reads displayed.

GTGGATTTTTTATTTACTGCGTACTTCGACAACCATTAATGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACA  >  NC_000913/953574‑953865

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

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

  ggATTTTTTATTTACTGCGTACTTCGACAACCATTAATGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATtcat                                                                                                                                              <  1:469210/150‑1 (MQ=255)

   gATTTTTTATTTACTGCGTACTTCGACAACCATTAATGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATtcatc                                                                                                                                             <  1:77921/150‑1 (MQ=255)

      tttttATTTACTGCGTACTTCGACAACCATTAATGGTGGTCTTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAATTTCATCAct                                                                                                                                          <  1:629059/150‑2 (MQ=255)

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

         ttatttaCTGCGTACTTCGACAACCATTAATGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAATTTCATCActgat                                                                                                                                       <  2:746597/150‑5 (MQ=255)

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

            tttaCTGCGTACTTCGACAACCATTAATGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAgg                                                                                                                                    <  2:879516/150‑1 (MQ=255)

            tttaCTGCGTACTTCGACAACCATTAATGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAgg                                                                                                                                    <  2:1184580/150‑1 (MQ=255)

              taCTGCGTACTTCGACAACCATTAATGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTATCCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTg                                                                                                                                  <  1:523730/150‑1 (MQ=255)

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

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

                cTGCGTACTTCGACAACCATTAATGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTgag                                                                                                                                >  1:1274589/1‑150 (MQ=255)

                cTGCGTACTTCGACAACCATTAATGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTgag                                                                                                                                >  1:736785/1‑150 (MQ=255)

                  gCGTACTTCGACAACCATTAATGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGaa                                                                                                                              >  1:831809/1‑150 (MQ=255)

                  gCGTACTTCGACAACCATTAATGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGaa                                                                                                                              >  2:1218698/1‑150 (MQ=255)

                  gCGTACTTCGACAAACATTAATGGTGGTCGTTTTCACGCAGTTAAATGACCCAGTATGTAACCCCAACCAACAAACCAACACCGATAATGATGACGATCGTAACCGCAATCAGGTTATAAGTGATGAAATTCATAAAGGTCAGGTGAGaa                                                                                                                              >  2:486020/1‑150 (MQ=255)

                       cTTCGACAACCATTAATGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGTCACCCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAAtt                                                                                                                         >  1:898256/1‑150 (MQ=255)

                        ttCGACTACCATTAATGGTGGTTGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGTTCAGGTGAGAAAAAttt                                                                                                                        <  2:322317/150‑1 (MQ=255)

                         tCGACAACCATTAATGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGATGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAAtttt                                                                                                                       >  1:1521839/1‑150 (MQ=255)

                          cGACAACCATTAATGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTc                                                                                                                      <  1:739370/150‑1 (MQ=255)

                           gACAACCATTAATGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTcc                                                                                                                     >  2:1032306/1‑150 (MQ=255)

                             cAACCATTAATGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCgg                                                                                                                   >  1:1001022/1‑150 (MQ=255)

                                ccATTAATGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAACTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGc                                                                                                                <  1:972104/150‑1 (MQ=255)

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

                                        ggtggtCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGa                                                                                                        >  1:1460052/1‑150 (MQ=255)

                                            gtCGTTTTCACGCAGGTAAATGTCAAAGTATGTAAAGCCAACGAACAAACCCGCACCGATAATTTTGCCGATCGTATCCGGAATGAGTTTTTCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGc                                                                                                    <  2:1303860/150‑1 (MQ=255)

                                            gtCGTTTTAACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGc                                                                                                    >  1:364084/1‑150 (MQ=255)

                                                ttttCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAACCCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTc                                                                                                >  1:1082984/1‑150 (MQ=255)

                                                ttttCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTc                                                                                                <  1:1314285/150‑1 (MQ=255)

                                                ttttCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTc                                                                                                <  1:1310439/150‑1 (MQ=255)

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

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

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

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

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

                                                               aTGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACGGACTGCGGTCCAAAATTCCGGGGAt                                                                                 <  2:1291203/150‑1 (MQ=255)

                                                                  aCCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCg                                                                              <  2:300049/150‑1 (MQ=255)

                                                                           gTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCgg                                                                     <  1:989708/150‑1 (MQ=255)

                                                                           gTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCgg                                                                     <  1:339771/150‑1 (MQ=255)

                                                                           gTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCACTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCgg                                                                     >  2:681789/1‑150 (MQ=255)

                                                                            tCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGa                                                                    >  1:420599/1‑150 (MQ=255)

                                                                            tCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGa                                                                    >  1:777098/1‑150 (MQ=255)

                                                                             caacccaaccAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGAt                                                                   >  1:82574/1‑150 (MQ=255)

                                                                               acccaaccAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACTGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTa                                                                 <  1:148907/150‑1 (MQ=255)

                                                                               acccaaccAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCg                                                                                      >  2:233381/1‑129 (MQ=255)

                                                                               acccaaccAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCg                                                                                      <  1:233381/129‑1 (MQ=255)

                                                                               acccaaccAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCgg                                                                     <  2:206608/146‑1 (MQ=255)

                                                                               acccaaccAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCgg                                                                     >  1:206608/1‑146 (MQ=255)

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

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

                                                                                      caacaaACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATAcc                                                          >  2:1357013/1‑150 (MQ=255)

                                                                                         caaACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCAt                                                       >  2:764452/1‑150 (MQ=255)

                                                                                          aaaCCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATc                                                      <  2:725331/150‑1 (MQ=255)

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

                                                                                                ccaccGATAATGTTGCCGAACGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCACGTGAGAAAAATTTTCCCGTGCAGAACCGACAGCGGTCCAAAATTCCGCGCATCCGAAGTCTCGGATTACAAAACCCATCGGCAta                                                >  2:366458/1‑149 (MQ=255)

                                                                                                 caccGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCa                                               <  2:1001022/150‑1 (MQ=255)

                                                                                                      aTAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCATAAAc                                          <  1:827232/150‑1 (MQ=255)

                                                                                                      aTAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCATAAAc                                          <  2:1470249/150‑1 (MQ=255)

                                                                                                       tAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCATAAACa                                         <  1:1032306/150‑1 (MQ=255)

                                                                                                              gCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCATAAACATGTTTGc                                  >  1:1450881/1‑150 (MQ=255)

                                                                                                              gCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCATAAACATGTTTGc                                  <  1:1357013/150‑1 (MQ=255)

                                                                                                              gCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCATAAACATGTTTGc                                  >  2:1282188/1‑150 (MQ=255)

                                                                                                               ccGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCATAAACATGTTTGCg                                 <  2:1053952/150‑1 (MQ=255)

                                                                                                                     gTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCATAAACATGTTTGCGATACtg                           >  1:917256/1‑150 (MQ=255)

                                                                                                                     gTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCATAAACATGTTTGCGATACtg                           >  1:915173/1‑150 (MQ=255)

                                                                                                                        aCCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGc                        >  1:1213059/1‑150 (MQ=255)

                                                                                                                         ccGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCt                       <  2:831809/150‑1 (MQ=255)

                                                                                                                             aaTCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCaaa                   <  1:1218698/150‑1 (MQ=255)

                                                                                                                             aaTCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCaaa                   >  1:599059/1‑150 (MQ=255)

                                                                                                                               tCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAAc                 <  1:486020/150‑1 (MQ=255)

                                                                                                                                 aGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCg               >  2:658669/1‑150 (MQ=255)

                                                                                                                                 aGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCg               >  2:1217729/1‑150 (MQ=255)

                                                                                                                                 aGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCg               >  1:1083266/1‑150 (MQ=255)

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

                                                                                                                                   cctgaccgTGATGAATTTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCt             >  2:1220545/8‑150 (MQ=255)

                                                                                                                                    ttATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTg            <  1:1200932/150‑1 (MQ=255)

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

                                                                                                                                      aTCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTAAAAACCGCTGGc          >  2:652847/1‑150 (MQ=255)

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

                                                                                                                                       tCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGca         >  1:606807/1‑150 (MQ=255)

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

                                                                                                                                       accgTGATGAATTTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGca         >  2:21453/4‑150 (MQ=255)

                                                                                                                                       accgTGATGAATTTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGca         <  1:558271/147‑1 (MQ=255)

                                                                                                                                       accgTGATGAATTTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGca         >  2:15656/4‑150 (MQ=255)

                                                                                                                                        ccgTGATGAATTTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGcaa        >  2:465699/3‑150 (MQ=255)

                                                                                                                                        ccgTGATGAATTTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGcaa        >  1:1119592/3‑150 (MQ=255)

                                                                                                                                        ccgTGATGAATTTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGcaa        >  1:962569/3‑150 (MQ=255)

                                                                                                                                        cAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGcaa        >  2:683326/1‑150 (MQ=255)

                                                                                                                                         cgTGATGAATTTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGcaac       >  1:1123632/2‑150 (MQ=255)

                                                                                                                                         cgTGATGAATTTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGcaac       >  1:1166414/2‑150 (MQ=255)

                                                                                                                                         aGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACTGCTGGcaac       <  2:674896/150‑1 (MQ=255)

                                                                                                                                         aGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGcaac       >  1:777030/1‑150 (MQ=255)

                                                                                                                                         aGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGcaac       >  1:1069318/1‑150 (MQ=255)

                                                                                                                                            gatgaAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCTCGGATTACAATACCCATCGGGATCATAATCATTTTTGCGATACTGTGCTCAAAACCGCTGGcaacaat    <  2:412704/150‑2 (MQ=255)

                                                                                                                                            gatgaAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCaacaaa    >  2:149755/1‑150 (MQ=255)

                                                                                                                                            gatgaAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCaacaaa    >  1:476782/1‑150 (MQ=255)

                                                                                                                                             atgaAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCaacaaac   <  1:1341788/150‑1 (MQ=255)

                                                                                                                                              tgaAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCaacaaaca  >  1:470490/1‑150 (MQ=255)

GTGGATTTTTTATTTACTGCGTACTTCGACAACCATTAATGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACA  >  NC_000913/953574‑953865

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