Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,729,405 C→G 8.3% intergenic (‑34/+38) xylB ← / ← xylA xylulokinase/D‑xylose isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,729,4050CG 8.3% 596.2 / 37.0 168intergenic (‑34/+38)xylB/xylAxylulokinase/D‑xylose isomerase
Reads supporting (aligned to +/- strand):  ref base C (78/76);  new base G (8/6);  total (86/82)
Fisher's exact test for biased strand distribution p-value = 7.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.48e-01

Only 100 of 170 total aligned reads displayed.

GATGCGGGCGCGAAACGGTCAGCTTTTCCGTTTGCGCAGCAACCACCTCACCCTGCTCGTTGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCGGAGACAAATGATGTTCCTGAGCATATTTGGCTAAATCTGCC  >  NC_000913/3729260‑3729551
                                                                                                                                                  |                                                                                                                                                  
gATGCGGGCGCGAAACGGTCAGCTTTTCCGTTTGCGCAGCAACCACCTCACCCTGCTCGTTGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCCTAAAAAAATGCCCGGTATCCCTa                                                                                                                                                  <  1:956269/148‑1 (MQ=255)
      ggCGCGAAACGGTCAGCTTTTCCGTTTGCGCAGCAACCACCTCACCCTGCTCGTTGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATa                                                                                                                                            >  2:52895/1‑148 (MQ=255)
          cgAAACGGTCAGCTTTTCCGTTTGCGCAGCAACCACCTCACCCTGCTCGTTGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCg                                                                                                                                        <  1:52895/148‑1 (MQ=255)
                     gCTTTTCCGTTTGCGCAGCAACCACCTCACCCTGCTCGTTGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCg                                                                                                                                        <  3:1053698/137‑1 (MQ=255)
                        tttCCGTTTGCGCAGCAACCACCTCACCCTGCTCGTTGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGc                                                                                                                          <  4:384887/148‑1 (MQ=255)
                              tttGCGCAGCAACCACCTCACCCTGCTCGTTGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTTATCGGTAGCGATaccg                                                                                                                                         >  2:336683/1‑125 (MQ=255)
                              tttGCGCAGCAACCACCTCACCCTGCTCGTTGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTTATCGGTAGCGATaccg                                                                                                                                         <  1:336683/128‑4 (MQ=255)
                                         accaccGCACCCTGCTCGTTGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTAttt                                                                                                         <  3:165674/148‑1 (MQ=255)
                                          ccaccTCACCCTGCTCGTTGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTg                                                                                                        >  4:179888/1‑148 (MQ=255)
                                               tCACCCTGCTCGTTGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGc                                                                                                                                     <  3:1074966/114‑1 (MQ=255)
                                                           ttGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTa                                                                                                                                                  >  2:634623/1‑89 (MQ=255)
                                                           ttGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTa                                                                                                                                                  <  1:634623/89‑1 (MQ=255)
                                                           ttGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATaa                                                                                                                                           <  1:711046/96‑1 (MQ=255)
                                                                   aaTAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTa                                                                                                                                                  >  3:893768/1‑81 (MQ=255)
                                                                   aaTAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTa                                                                                                                                                  <  4:893768/81‑1 (MQ=255)
                                                                      aaCTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGca                                                                                                                         >  2:16085/1‑103 (MQ=255)
                                                                          tttACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGa                                                                                                                              >  4:729582/1‑94 (MQ=255)
                                                                          tttACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACt                                                                                                                            <  2:532705/96‑1 (MQ=255)
                                                                          tttACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACt                                                                                                                            >  1:532705/1‑96 (MQ=255)
                                                                                cccGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGcaca                                                                                                                       >  2:1075610/1‑95 (MQ=255)
                                                                                cccGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAg                                                                                                <  1:238479/118‑1 (MQ=255)
                                                                                cccGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGAtttt                                                                            <  4:1040795/138‑1 (MQ=255)
                                                                                          cAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAAcc                                                                                                                                         <  1:16160/67‑1 (MQ=255)
                                                                                          cAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCaa                                                                                                                                  >  2:698043/1‑74 (MQ=255)
                                                                                          cAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTAtt                                                                                                          <  1:1053443/98‑1 (MQ=255)
                                                                                          cAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGa                                                           >  3:708755/1‑145 (MQ=255)
                                                                                          cAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGa                                                           <  4:708755/145‑1 (MQ=255)
                                                                                            aGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGtt                                                                                                             <  3:74609/93‑1 (MQ=255)
                                                                                            aGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGtt                                                                                                             >  4:74609/1‑93 (MQ=255)
                                                                                            aGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGtt                                                                                       >  3:923314/1‑115 (MQ=255)
                                                                                            aGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGtt                                                                                       <  4:923314/115‑1 (MQ=255)
                                                                                                 tATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAAt                                                                                           >  1:294196/1‑106 (MQ=255)
                                                                                                   tCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTAtt                                                                                                          <  1:840158/89‑1 (MQ=255)
                                                                                                    cccGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAg                                                                                                >  3:479328/1‑98 (MQ=255)
                                                                                                    cccGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAg                                                                                                <  4:479328/98‑1 (MQ=255)
                                                                                                    cccGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCa                                              >  3:946471/1‑148 (MQ=255)
                                                                                                    cccGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCa                                              >  1:656428/1‑148 (MQ=255)
                                                                                                          ataCATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCGgag                                        <  2:656428/148‑1 (MQ=255)
                                                                                                            aCATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAg                                                                                                <  4:1101304/90‑1 (MQ=255)
                                                                                                             cATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAAt                                                                                           <  4:545116/94‑1 (MQ=255)
                                                                                                             cATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAAt                                                                                           >  3:545116/1‑94 (MQ=255)
                                                                                                               tatCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTcga                                                                >  1:672282/1‑117 (MQ=255)
                                                                                                               tatCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTcga                                                                <  2:672282/119‑3 (MQ=255)
                                                                                                                 tcgatcgTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGtt                                                                                                             >  3:326779/1‑72 (MQ=255)
                                                                                                                 tcgatcgTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGt                                                                                        >  4:1023182/1‑93 (MQ=255)
                                                                                                                  cgatcgTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTg                                                   >  1:724913/1‑129 (MQ=255)
                                                                                                                  cgatcgTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTg                                                   <  2:724913/129‑1 (MQ=255)
                                                                                                                    atcgTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAg                                                                                 >  2:240467/1‑97 (MQ=255)
                                                                                                                    atcgTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTc                                                                           >  4:197215/1‑103 (MQ=255)
                                                                                                                    atcgTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTc                                                                           <  3:197215/103‑1 (MQ=255)
                                                                                                                     tcgTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGt                                                                                                              >  4:497129/1‑67 (MQ=255)
                                                                                                                     tcgTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTa                                                                                     <  1:407804/92‑1 (MQ=255)
                                                                                                                     tcgTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTa                                                                                     >  2:407804/1‑92 (MQ=255)
                                                                                                                     tcgTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCgg                                          >  1:215593/1‑135 (MQ=255)
                                                                                                                     tcgTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCgg                                          <  2:215593/135‑1 (MQ=255)
                                                                                                                     tcgTTCCTTAAAAAAATGCCAGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGt                                                                                                              <  3:497129/67‑1 (MQ=255)
                                                                                                                      cgTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCGGAGACAAATGATGtt                            >  1:329197/1‑148 (MQ=255)
                                                                                                                        ttCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCGGAGACAAATGATGTTcc                          >  4:432069/1‑148 (MQ=255)
                                                                                                                         tCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGcaca                                                                                                                       >  4:339043/1‑54 (MQ=255)
                                                                                                                          ccTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCGGAGACAAATGATGTTCCTg                        <  2:861965/148‑1 (MQ=255)
                                                                                                                             tAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCGGAGACaa                                    <  3:244146/133‑1 (MQ=255)
                                                                                                                             tAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGCCGAACAGATAATGGTTTACCAGAtttt                                                                            >  2:576396/1‑93 (MQ=255)
                                                                                                                             tAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGCCGAACAGATAATGGTTTACCAGAtttt                                                                            <  1:576396/93‑1 (MQ=255)
                                                                                                                             tAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCAAACTTAGCCGTTATTTGTCTAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCGgaga                                       <  2:808424/130‑1 (MQ=255)
                                                                                                                               aaaaaaTGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCGGAGACAAATGATGTTCCTg                        >  1:462872/1‑143 (MQ=255)
                                                                                                                               aaaaaaTGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCGGAGACAAATGATGTTCCTg                        <  2:462872/143‑1 (MQ=255)
                                                                                                                                aaaaaTGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGt                                                                                                       >  4:494484/1‑63 (MQ=255)
                                                                                                                                aaaaaTGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGt                                                                                                       <  3:494484/63‑1 (MQ=255)
                                                                                                                                aaaaaTGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCAc                                                       >  4:1061620/1‑111 (MQ=255)
                                                                                                                                aaaaaTGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCAc                                                       <  3:1061620/111‑1 (MQ=255)
                                                                                                                                aaaaaTGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCgg                                          >  3:784370/1‑124 (MQ=255)
                                                                                                                                aaaaaTGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCgg                                          <  4:784370/124‑1 (MQ=255)
                                                                                                                                 aaaaTGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTa                                                                                                            <  4:244544/57‑1 (MQ=255)
                                                                                                                                 aaaaTGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTg                                                                                                        <  1:1058969/61‑1 (MQ=255)
                                                                                                                                 aaaaTGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATg                                                                                          >  2:831510/1‑75 (MQ=255)
                                                                                                                                 aaaaTGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATg                                                                                          <  1:831510/75‑1 (MQ=255)
                                                                                                                                  aaaTGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGa                                                           <  2:301274/105‑1 (MQ=255)
                                                                                                                                   aaTGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCGGAGACaaa                                   <  4:1072859/128‑1 (MQ=255)
                                                                                                                                    aTGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCGGAGACAAATGATGTTCCTGAGCATATTTg              <  3:432069/148‑1 (MQ=255)
                                                                                                                                         cggttATCGGTAGCGATACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGa                                                  <  2:835467/103‑1 (MQ=255)
                                                                                                                                         cggttATCGGTAGCGATACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCGGAGACAAATGATGTTCCTGAGCATATTTGGCTaa         >  3:730918/5‑148 (MQ=255)
                                                                                                                                         cggttATCGGTAGCGATACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCGGAGACAAATGATGTTCCTGAGCATATTTGGCTaa         >  3:79387/5‑148 (MQ=255)
                                                                                                                                         cggttATCGGTAGCGATACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCGGAGACAAACGATGTTCCTGAGCATATTTGGCTaa         >  1:808492/5‑148 (MQ=255)
                                                                                                                                         cGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTcc                                                                          >  2:10556/1‑83 (MQ=255)
                                                                                                                                         cGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCActct                                                    <  3:1037514/105‑1 (MQ=255)
                                                                                                                                          ggTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTgg                                                              >  1:304394/1‑94 (MQ=255)
                                                                                                                                          ggTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTgg                                                              <  2:304394/94‑1 (MQ=255)
                                                                                                                                           gttATCGGTAGCGATACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCGGAGACAAATGATGTTCCTGAGCATATTTGGCTAAAt       <  2:170267/146‑1 (MQ=255)
                                                                                                                                           gTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCa                                                        <  1:1124144/99‑1 (MQ=255)
                                                                                                                                           gTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCa                                                        >  2:1124144/1‑99 (MQ=255)
                                                                                                                                           gTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCGGAGACAAATGATGTTCCTGAGCATATTTGGCTAAAt       >  1:771018/1‑148 (MQ=255)
                                                                                                                                             tATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCt                                                                >  1:121462/1‑90 (MQ=255)
                                                                                                                                             tATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGc                                               >  2:1044560/1‑107 (MQ=255)
                                                                                                                                             tATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGc                                               <  1:1044560/107‑1 (MQ=255)
                                                                                                                                             tATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCGGAGACAAATGATGTTCCTGAGCATATTTgg             >  3:1017351/1‑141 (MQ=255)
                                                                                                                                             tATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCGGAGACAAATGATGTTCCTGAGCATATTTGGCTAAATc      >  2:968152/1‑148 (MQ=255)
                                                                                                                                              aTCGGTAGCGATACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCGGAGACAAATGATGTTCCTGAGCATATTTGGCTAAATCTg    <  4:79387/148‑1 (MQ=255)
                                                                                                                                                cGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCGGAGACAAATGATGTTCCTGAGCATAtt                <  3:401084/135‑1 (MQ=255)
                                                                                                                                                 gCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCGGAGACAAATGATGTTCCTGAGCATATTTGGCTAAATCTGcc  >  3:602758/1‑148 (MQ=255)
                                                                                                                                                  cTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAg                                                                                 <  1:945318/68‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
GATGCGGGCGCGAAACGGTCAGCTTTTCCGTTTGCGCAGCAACCACCTCACCCTGCTCGTTGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCGGAGACAAATGATGTTCCTGAGCATATTTGGCTAAATCTGCC  >  NC_000913/3729260‑3729551

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