Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,313,090 T→G intergenic (‑99/‑260) alsR ← / → rpiB d‑allose‑inducible als operon transcriptional repressor; autorepressor; repressor of rpiR/ribose 5‑phosphate isomerase B/allose 6‑phosphate isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,313,0900TG97.0% 937.4 / 7.5 273intergenic (‑99/‑260)alsR/rpiBd‑allose‑inducible als operon transcriptional repressor; autorepressor; repressor of rpiR/ribose 5‑phosphate isomerase B/allose 6‑phosphate isomerase
Reads supporting (aligned to +/- strand):  ref base T (4/4);  new base G (137/126);  total (142/131)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 291 total aligned reads displayed.

CGTTCTCTTGTCCTTGAATTGAAGACGCCACTATATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTTTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTAATTAGTTTAAAAACATGATGTTATTTGTTTTTCTGGTTTTTTAGATG  >  NC_000913/4312994‑4313185
                                                                                                |                                                                                               
    ctctTGTCCTTGAATTGAAGACGCCACTATATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGtgtatttta                                                                                          >  1:4922163/1‑100 (MQ=255)
    ctctTGTCCTTGAATTGAAGACGCCACTATATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGtgtatttta                                                                                          >  2:4872814/1‑100 (MQ=255)
     tctTGTCCTTGAATTGAAGACGCCACTATATCACCACAGCAAGAACAGAAAATATGCTGGCGACTTCCCAAAACAATACAATTTATAAAGtgtattttat                                                                                         >  2:3910034/1‑100 (MQ=255)
      ctTGTCCTTGAATTGAAGACGCCACTATATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTgtattttatt                                                                                        >  1:1064027/1‑100 (MQ=255)
      ctTGTCCTTGAATTGAAGACGCCACTATATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTgtattttatt                                                                                        <  1:1817111/100‑1 (MQ=255)
       ttGTCCTTGAATTGAAGGCTCGACTTTATCACCACATCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTCTATTTTATtg                                                                                       <  2:5202001/100‑1 (MQ=255)
       ttGTCCTTGAATTGAAGACGCCACTATATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTGTATTTTATtg                                                                                       <  1:875918/100‑1 (MQ=255)
       ttGTCCTTGAATTGAAGACGCCACTATATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTGTATTTTATtg                                                                                       <  2:3522382/100‑1 (MQ=255)
        tGTCCTTGAATTGAAGACGCCACTATATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTGTATTTTATtgt                                                                                      <  1:5488462/100‑1 (MQ=255)
         gTCTTTCAATTGAAGACGCCACTATATCACAACATCAAGCACAGAAAAGATGCTGGCCACTTCACAACAAAATAAAATATATAAAGTGTATTTTATtgtg                                                                                     >  2:3456522/1‑100 (MQ=255)
         gTCCTTGAATTGAAGACGCCACTATATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTGTATTTTATtgtg                                                                                     >  2:5465184/1‑100 (MQ=255)
         gTCCTTGAATTGAAGACGCCACTATATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTGTATTTTATtgtg                                                                                     >  1:1976600/1‑100 (MQ=255)
         gTCCTTGAATTGAAGACGCCACTATATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTGTATTTTATtgtg                                                                                     <  1:999820/100‑1 (MQ=255)
            cTTGAATTGAAGACGCCACTATATCGCCACAGCAAGAACAGAAAAGATTCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTGTATTTTATTGTGAAc                                                                                  <  1:3927505/100‑1 (MQ=255)
            cTTGAATTGAAGACGCCACTATATCACCACAGCAAGAACAGAAAAGATGCTGGCTACTTCACAAAAAAATAAAATTTGTAAAGTGTATTGTATTGTGAAc                                                                                  >  2:1260929/1‑100 (MQ=255)
            cTTGAATTGAAGACGCCACTATATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTGTATTTTATTGTGAAc                                                                                  >  2:373243/1‑100 (MQ=255)
             ttgaattgaaGACGCCACTATATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTGTATTTTATTGTGAACa                                                                                 <  1:2196619/100‑1 (MQ=255)
             ttgaattgaaGACGCCACTATATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTGTATTTTATTGTGAACa                                                                                 <  1:4376901/100‑1 (MQ=255)
              tgaattgaaGACGCAACTATATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTGTATTTTATTGTGAACaa                                                                                <  2:4695787/100‑1 (MQ=255)
                aattgaaGACGCCACTATATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAAGTTATAAAGTGTATTTTATTGTGAACAATg                                                                              >  1:5093830/1‑100 (MQ=255)
                 attgaaGACGCCACTATATCACCACAGCAAGAACAGAAAAGCTGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTGTATTTTTTTTTTAACAAttg                                                                             >  2:5232421/1‑98 (MQ=255)
                 attgaaGACGCCACTATATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTGTATTTTATTGTGAACAATgg                                                                             >  2:1632700/1‑100 (MQ=255)
                  ttgaaGACGCCACTATATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTGTATTTTATTGTGAACAATGGc                                                                            <  2:4826022/100‑1 (MQ=255)
                    gaaGACGCCACTATATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTGTATTTTATTGTGAACAATGGCGa                                                                          >  2:336055/1‑100 (MQ=255)
                    gaaGACGCCACTATATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTGTATTTTATTGTGAACAATGGCGa                                                                          >  1:3492247/1‑100 (MQ=255)
                    gaaGACGCCACTATATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTGTATTTTATTGTGAACAATGGCGa                                                                          >  2:793201/1‑100 (MQ=255)
                     aaGACGCCACTATATCACCACAGCAAGAACAGAAAAGATGCTAGCGACTTCACACAAAAATAAAATTTATAAAGTGTATTTTATCGTGAACAATGGCGAg                                                                         <  1:4844849/100‑1 (MQ=255)
                       gACGCCACTATATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTGTATTTTATTGTGAACAATGGCGAGTg                                                                       <  1:3260121/100‑1 (MQ=255)
                       gACGCCACTATATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTGTATTTTATTGTGAACAATGGCGAGTg                                                                       >  1:2901991/1‑100 (MQ=255)
                       gACGCCACTATATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTGTATTTTATTGTGAACAATGGCGAGTg                                                                       >  2:5014746/1‑100 (MQ=255)
                       gACGCCACTATATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTGTATTTTATTGTGAACAATGGCGAGTg                                                                       <  2:2701773/100‑1 (MQ=255)
                         cGCCACTATATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTGTATTTTATTGTGAACAATGGCGAGTGGc                                                                     >  2:1731988/1‑100 (MQ=255)
                          tccgcTATATCACGACAGCAAGAACAGAAAAGAAGCTGGCGACTTCACAAAAAAATAAAATTTAGAAAGTGTATTTTATTGTGAACAATGGCGAGGGGCt                                                                    <  2:2476500/96‑1 (MQ=255)
                           ccACTATATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTGTATTTTATTGTGAACAATGGCGAGTGGCtt                                                                   >  2:3294415/1‑100 (MQ=255)
                               tataTCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTGTATTTTATTGTGAACAATGGCGAGTGGCTTaaaa                                                               <  2:1049938/100‑1 (MQ=255)
                                atatCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTGTATTTTATTGTGAACAATGGCGAGTGGCTTAAAAt                                                              >  1:5020464/1‑100 (MQ=255)
                                atatCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTGTATTTTATTGTGAACAATGGCGAGTGGCTTAAAAt                                                              >  1:4320858/1‑100 (MQ=255)
                                atatCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTGTATTTTATTGTGAACAATGGCGAGTGGCTTAAAAt                                                              <  1:2747657/100‑1 (MQ=255)
                                 tatCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTGTATTTTATTGTGAACAATGGCGAGTGGCTTAAAAtt                                                             <  2:2819106/100‑1 (MQ=255)
                                 tatCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTGTATTTTATTGTGAACAATGGCGAGTGGCTTAAAAtt                                                             >  1:3634099/1‑100 (MQ=255)
                                   tCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTGTATTTTATTGTGAACAATGGCGAGTGGCTTAAAAtttt                                                           <  2:2807505/100‑1 (MQ=255)
                                     accacAGCAAGCCCGGAAAACATGATGGCGACTTCACAACAAAATAACATTTATAACGTGTATTTTCTTGTGAACAATGGCGAGTGGCTTACAATTTTAt                                                         <  2:253351/100‑1 (MQ=255)
                                     accacAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTGTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTAt                                                         <  2:4818747/100‑1 (MQ=255)
                                     accacAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTGTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTAt                                                         <  2:3945642/100‑1 (MQ=255)
                                      ccacAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTGTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTAtt                                                        >  1:2401939/1‑100 (MQ=255)
                                       cacaGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTGTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTAttt                                                       <  2:180683/100‑1 (MQ=255)
                                       cacaGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTGTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTAttt                                                       <  1:5407032/100‑1 (MQ=255)
                                        acaGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTGTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTc                                                      >  2:3570626/1‑100 (MQ=255)
                                         caGCAAGTACAGAAAAGATGCTGGCGACGGCCCAAAAAACTATAATCTATAAAGTCTCGTTGATTGTGGAAAATTGCAAAAGCCATAAAATTTTATTTCa                                                     >  2:653498/1‑100 (MQ=255)
                                           gCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTGTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCAtt                                                   >  2:4850712/1‑100 (MQ=255)
                                           gCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTGTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCAtt                                                   >  2:2280442/1‑100 (MQ=255)
                                            cAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTGTATTTTATTGTGAACAATGGCGTGTGGCTTAAAATTTTATTTCAttt                                                  >  2:1922360/1‑100 (MQ=255)
                                            cAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTGTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCAttt                                                  >  1:2545986/1‑100 (MQ=255)
                                             aaGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTGTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTa                                                 >  1:954354/1‑100 (MQ=255)
                                               gAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTTTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTAAt                                               >  2:1118464/1‑100 (MQ=255)
                                                   aGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTTTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTAATTAGt                                           <  1:1133542/100‑1 (MQ=255)
                                                    gAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTGTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTAATTAGtt                                          <  2:2602092/100‑1 (MQ=255)
                                                    gAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTGTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTAATTAGtt                                          <  1:4011797/100‑1 (MQ=255)
                                                    gAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTGTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTAATTAGtt                                          >  2:4804869/1‑100 (MQ=255)
                                                     aaaaGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTGTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTAATTAGttt                                         >  1:4068984/1‑100 (MQ=255)
                                                      aaaGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTGTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTAATTAGTTTa                                        >  1:657276/1‑100 (MQ=255)
                                                          aTGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTGTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTAATTAGTTTaaaaa                                    <  1:2223045/100‑1 (MQ=255)
                                                           tGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTGTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTAATTAGTTTAAAAAc                                   >  1:4638878/1‑100 (MQ=255)
                                                           tGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTGTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTAATTAGTTTAAAAAc                                   >  1:155485/1‑100 (MQ=255)
                                                             cTGGCGACTTCACAAAAAAATAAAATTTATAAAGTGTATTTTATTGTTAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTAATTAGTTTAAAAACat                                 >  2:2157526/1‑100 (MQ=255)
                                                              tGGCGACTTCACAAAAAAATAAAATTTATAAAGTGTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTAATTAGTTTAAAAACatg                                <  1:265064/100‑1 (MQ=255)
                                                              tGGCGACTTCACAAAAAAATAAAATTTATAAAGTGTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTAATTAGTTTAAAAACatg                                <  1:418514/100‑1 (MQ=255)
                                                               ggCGACTTCACAAAAAAATAAAATTTATAAAGTGTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTCATTAGTTTAAAAACatgc                               >  2:4553764/1‑99 (MQ=255)
                                                                      tCACAAAAAAATAAAATTTATAAAGTGTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTAATTAGTTTAAAAACATGATGTTAtt                        <  1:1965290/100‑1 (MQ=255)
                                                                       cacaAAAAAATAAAATTTATAAAGTGTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTAATTAGTTTAAAAACATGATGTTAttt                       >  2:702339/1‑100 (MQ=255)
                                                                        acaAAAAAATAAAATTTATAAAGTGTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTAATTAGTTTAAAAACATGATGTTATTTg                      <  2:153974/100‑1 (MQ=255)
                                                                        acaAAAAAATAAAATTTATAAAGTGTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTAATTAGTTTAAAAACATGATGTTATTTg                      >  1:1409273/1‑100 (MQ=255)
                                                                         caAAAAAATAAAATTTATAAAGTGTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTAATTAGTTTAAAAACATGATGTTATTTGt                     >  1:1645296/1‑100 (MQ=255)
                                                                         caAAAAAATAAAATTTATAAAGTGTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTAATTAGTTTAAAAACATGATGTTATTTGt                     >  1:2007679/1‑100 (MQ=255)
                                                                         caAAAAAATAAAATTTATAAAGTGTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTAATTAGTTTAAAAACATGATGTTATTTGt                     >  2:5118186/1‑100 (MQ=255)
                                                                         aaaaaaaaTAAAATTTATAAAGTGTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTAATTAGTTTAAAAACATGATGTTATTTGt                     >  2:4686741/2‑100 (MQ=255)
                                                                           aaaaaaTAAAATTTATAAAGTGTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTAATTAGTTTAAAAACATGATGTTATTTGttt                   <  1:120325/100‑1 (MQ=255)
                                                                           aaaaaaTAAAATTTATAAAGTGTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTAATTAGTTTAAAAACATGATGTTATTTGttt                   <  2:4383182/100‑1 (MQ=255)
                                                                           aaaaaaTAAAATTTATAAAGTGTATTTTATTGTGAACAATGGCGAATGGCTTAAAATTTTATTTTATTTAATTAGTTTAAAAACATGATTTTAttttttc                   >  2:4187852/1‑99 (MQ=255)
                                                                            aaaaaTAAAATTTATAAAGTGTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTAATTAGTTTAAAAACATGATGTTATTTGtttt                  <  1:214398/100‑1 (MQ=255)
                                                                            aaaaaTAAAATTTATAAAGTGTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTAATTAGTTTAAAAACATGATGTTATTTGtttt                  <  2:1482805/100‑1 (MQ=255)
                                                                               aataaaatTTATAAAGTGTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTAATTAGTTTAAAAACATGATGTTATTTGTTTTTCt               <  2:2378203/100‑1 (MQ=255)
                                                                               aataaaatTTATAAAGTGTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTAATTAGTTTAAAAACATGATGTTATTTGTTTTTCt               >  2:4920743/1‑100 (MQ=255)
                                                                                ataaagtTTATCAAGAGTACTGTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTCCATTTAATTAGTTTAAAAACATGATGTTACTTGTTTTTCTg              <  2:3972956/100‑1 (MQ=255)
                                                                                ataaaatTTATAAAGTGTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTAATTAGTTTAAAAACATGATGTTATTTGTTTTTCTg              <  2:2764974/100‑1 (MQ=255)
                                                                                ataaaatTTATAAAGTGTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTAATTAGTTTAAAAACATGATGTTATTTGTTTTTCTg              <  1:526437/100‑1 (MQ=255)
                                                                                ataaaatTTATAAAGTGTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTAATTAGTTTAAAAACATGATGTTATTTGTTTTTCTg              <  2:3650844/100‑1 (MQ=255)
                                                                                    aatTTATAAAGTGTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTAATTAGTTTAAAAACATGATGTTATTTGTTTTTCTGGttt          >  1:4695810/1‑100 (MQ=255)
                                                                                    aatTTATAAAGTGCATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTAATTAGTTTAAAAACATGATGTTATTTGTTTTTCTGGttt          >  1:4898224/1‑100 (MQ=255)
                                                                                     atTTATAAAGTGTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTAATTAGTTTCAAAACATGATGTTATTTCTTTTTCTGGtttt         >  2:4305555/1‑100 (MQ=255)
                                                                                     atTTATAAAGTGTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTAATTAGTTTAAAAACATGATGTTATTTGTTTTTCTGGtttt         <  1:1861576/100‑1 (MQ=255)
                                                                                      tttATAAAGTGTATTTTATTGTGAACAATGGCGAGTGGGTTAAAATTTTATTTCATTTAATTAGTTTAAAAACATGATGTTATTTGTTTTTCTGGttttt        >  2:2058675/1‑100 (MQ=255)
                                                                                      tttATAAAGTGTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTAATTAGTTTAAAAACATGATGTTATTTGTTTTTCTGGttttt        >  2:891002/1‑100 (MQ=255)
                                                                                       ttATAAAGTGTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTAATTAGTTTAAAAACATGATGTTATTTGTTTTTCTGGtttttt       <  1:5189306/100‑1 (MQ=255)
                                                                                        tataAAGTGTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTAATTAGTTTAAAAACATGATGTTATTTGTTTTTCTGGTTTTTTa      <  1:3353030/100‑1 (MQ=255)
                                                                                        tataAAGTGTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTAATTAGTTTAAAAACATGATGTTATTTGTTTTTCTGGTTTTTTa      <  1:711092/100‑1 (MQ=255)
                                                                                          taAAGTGTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTAATCAGTTTAAAAACATGATGTTATTTGTTTTTCTGGTTTTTTAga    <  2:3093410/100‑1 (MQ=255)
                                                                                           aaaGTTTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTAATTAGTTTAAAAACATGATGTTATTTGTTTTTCTGGTTTTTTAgat   >  1:2857637/1‑100 (MQ=255)
                                                                                           aaaGTGTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTAATTAGTTTAAAAACATGATGTTATTTGTTTTTCTGGTTTTTTAgat   <  2:1898980/100‑1 (MQ=255)
                                                                                           aaaGTGTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTAATTAGTTTAAAAACATGATGTTATTTGTTTTTCTGGTTTTTTAgat   <  2:1959065/100‑1 (MQ=255)
                                                                                                |                                                                                               
CGTTCTCTTGTCCTTGAATTGAAGACGCCACTATATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTTTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTAATTAGTTTAAAAACATGATGTTATTTGTTTTTCTGGTTTTTTAGATG  >  NC_000913/4312994‑4313185

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