Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,729,408 C→G 15.8% intergenic (‑37/+35) 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,4080CG15.8% 494.4 / 79.2 177intergenic (‑37/+35)xylB/xylAxylulokinase/D‑xylose isomerase
Reads supporting (aligned to +/- strand):  ref base C (80/69);  new base G (19/9);  total (99/78)
Fisher's exact test for biased strand distribution p-value = 2.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.16e-02

Only 100 of 186 total aligned reads displayed.

CCCTGCTCGTTGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCA  >  NC_000913/3729310‑3729507
                                                                                                   |                                                                                                   
cccTGCTCGTTGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTAcc                                                                                                    >  1:3380677/1‑100 (MQ=255)
   tGCTCGTTGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGAt                                                                                                 >  2:4169511/1‑100 (MQ=255)
       cGTTGAGCAAAATAACTTTTACGCAAGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAACAAAATGCCCAG‑TATCGCTACCGATAAcc                                                                                             >  1:2001449/1‑100 (MQ=255)
        gTTGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCg                                                                                            >  1:869572/1‑100 (MQ=255)
          tGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCggg                                                                                          >  1:3727094/1‑100 (MQ=255)
          tGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCggg                                                                                          <  2:4357134/100‑1 (MQ=255)
           gAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGc                                                                                         >  2:939765/1‑100 (MQ=255)
            aGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGcc                                                                                        <  1:1478438/100‑1 (MQ=255)
            aGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGcc                                                                                        >  1:3208697/1‑100 (MQ=255)
             gCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTTATCGGTAGCGAT‑ACCGGGcat                                                                                       <  2:3685349/100‑3 (MQ=37)
               aaaaTAACTTTTACGCCCGAGGTGACAAGATCTATCCCGGTATACATATCGATCGTTCCTTAAAAAAATGCCCGGTTATCGGTAGCGAT‑ACCGGGcattt                                                                                     <  1:3801320/100‑5 (MQ=25)
                 aaTAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACgg                                                                                   <  2:1202554/100‑1 (MQ=255)
                   tAACTTTTACGCCCGAGGTGACAAGATCTATCCCGATATACATATCGATCGTTCCTTACAAAAATGCCCGGTTATCGGTAGCGAT‑ACCGGGcattttttt                                                                                 <  1:1793325/100‑9 (MQ=21)
                     aCTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTg                                                                               >  1:665802/1‑100 (MQ=255)
                      cTTTTACGCCCGAGGTGCCAAGATCTATCCCGCTATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTg                                                                               <  2:1913970/99‑1 (MQ=255)
                      cTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGc                                                                              >  2:3928947/1‑100 (MQ=255)
                      cTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGc                                                                              >  2:316909/1‑100 (MQ=255)
                       ttttACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGca                                                                             <  1:189327/100‑1 (MQ=255)
                            cGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGtt                                                                        >  1:2588937/1‑100 (MQ=255)
                            cGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGtt                                                                        >  2:4011116/1‑100 (MQ=255)
                            cGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGtt                                                                        <  1:2296184/100‑1 (MQ=255)
                            cGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGtt                                                                        >  2:1143102/1‑100 (MQ=255)
                             gCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTa                                                                       <  1:2804504/100‑1 (MQ=255)
                               ccGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAg                                                                      >  2:3604996/1‑99 (MQ=255)
                               ccGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGc                                                                     >  1:1148236/1‑100 (MQ=255)
                               ccGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGc                                                                     <  1:4067028/100‑1 (MQ=255)
                               ccGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGc                                                                     <  1:4314844/100‑1 (MQ=255)
                                cGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGcc                                                                    <  2:1462020/100‑1 (MQ=255)
                                cGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGc                                                                     <  1:3054689/99‑1 (MQ=255)
                                   ggTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGtt                                                                 <  1:1616026/100‑1 (MQ=255)
                                     tGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGtttt                                                               <  2:4093656/100‑3 (MQ=255)
                                      gCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTAtt                                                              <  2:3135526/100‑1 (MQ=255)
                                      gCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTAtt                                                              >  2:704978/1‑100 (MQ=255)
                                      gCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTAtt                                                              >  2:258506/1‑100 (MQ=255)
                                      gCCAAGATCTATCCCGATATACATATCGATCGTTACTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTAt                                                               <  1:49239/99‑1 (MQ=255)
                                        cAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTg                                                            >  2:2408958/1‑100 (MQ=255)
                                          aGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTc                                                          <  2:2347115/100‑1 (MQ=255)
                                          aGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTc                                                          >  2:4162946/1‑100 (MQ=255)
                                           gATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCg                                                         <  1:3623431/100‑1 (MQ=255)
                                           gATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCg                                                         >  2:1077201/1‑100 (MQ=255)
                                            aTCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTc                                                          >  2:3027967/1‑98 (MQ=255)
                                             tCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGaa                                                       <  2:4207955/100‑1 (MQ=255)
                                               tATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACa                                                     >  2:1418744/1‑100 (MQ=255)
                                                   ccGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATa                                                 >  2:2536930/1‑100 (MQ=255)
                                                   ccGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCAGTTAATTGTCGacc                                                      >  1:3002227/1‑93 (MQ=255)
                                                    cggtataCAGATCGAGCGTTCCTTAAAAAGATGCCCGG‑TATCGCTACCGATAACCGGGACAACGGACTGCACCGTTAGCCGTTATTTGTCGAAc                                                      <  1:779940/91‑1 (MQ=255)
                                                     gATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTTCACAGTTAGCCGTTATTTGTCGAACAGATAAt                                               >  2:2857498/1‑100 (MQ=255)
                                                       tataCATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATgg                                             <  1:1273126/100‑1 (MQ=255)
                                                       tataCATATCGATCGTTCCTTAAAAAAATGCCCGG‑TAGCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATgg                                             >  1:2354412/1‑100 (MQ=255)
                                                       tataCATATCGATCCTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATgn                                             <  1:3556749/100‑2 (MQ=255)
                                                        ataCATATCGATCGTTCCTTAAAAAAATGCCCGG‑TAGCGCTACCGATATCCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGt                                            <  2:3101833/100‑1 (MQ=255)
                                                         taCATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGtt                                           >  1:789805/1‑100 (MQ=255)
                                                          aCATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGttt                                          >  1:1162824/1‑100 (MQ=255)
                                                           cATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTa                                         >  1:2505580/1‑100 (MQ=255)
                                                           cATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTa                                         >  1:3855714/1‑100 (MQ=255)
                                                           cATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACAGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTa                                         >  2:2711342/1‑100 (MQ=255)
                                                             tatCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTAcc                                       >  1:2213339/1‑100 (MQ=255)
                                                             tatCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGATTccc                                       >  2:2060783/1‑100 (MQ=255)
                                                              atcgatcgTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCa                                      <  1:612262/100‑1 (MQ=255)
                                                                cgagcgTTCCTTAAAAAAATGCCCGG‑TATCGGTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTTGAACAGATAATGGTTTACCAGa                                    <  1:2558439/100‑1 (MQ=255)
                                                                   tcgTTCCTTAAAAAAATGCCAGG‑TATCGATACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGAttt                                 >  1:529394/1‑100 (MQ=255)
                                                                     gTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGAtttt                                >  1:676023/1‑99 (MQ=255)
                                                                     gTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTc                               >  2:3437896/1‑100 (MQ=255)
                                                                      ttCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTcc                              <  1:3501053/100‑1 (MQ=255)
                                                                        ccTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAATGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAg                            <  2:1614736/100‑1 (MQ=255)
                                                                        ccTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAg                            >  1:1925291/1‑100 (MQ=255)
                                                                         cTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATACTGGTTTACCAGATTTTCCAgt                           >  1:1285917/1‑100 (MQ=255)
                                                                         cTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAgt                           <  2:2643806/100‑1 (MQ=255)
                                                                          ttAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAgtt                          >  1:3402665/1‑100 (MQ=255)
                                                                           tAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTTCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAgttg                         <  1:2843335/100‑1 (MQ=255)
                                                                            aaaaaaaTGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAgttgt                        <  1:3164076/100‑1 (MQ=255)
                                                                              ccgttggCCCGGTTATCGGTAGCGAT‑ACCGGGCCAACGGACTGCACAGTTAGCCGTTATTCGTCGAACAGCTACTGGTTTACCAGATTTTCCAGTTGTTc                      >  1:4338477/7‑100 (MQ=21)
                                                                                    gCCCGGTTATCGGTAGCGAT‑ACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCg                >  1:3865094/1‑100 (MQ=38)
                                                                                    gCCCGGTTATCGGTAGCGAT‑ACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTAGTTCCTGGCg                >  1:2893359/1‑100 (MQ=37)
                                                                                     cccGGTTATCGGTAGCGAT‑ACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGa               >  2:2762663/1‑100 (MQ=38)
                                                                                     cccGGTTATCGGTAGCGAT‑ACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGa               >  2:898526/1‑100 (MQ=38)
                                                                                     cccGGTTATCGGTAGCGAT‑ACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGa               >  2:2239613/1‑100 (MQ=38)
                                                                                     cccGGTTATCGGTAGCGAT‑ACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGa               >  2:2758796/1‑100 (MQ=38)
                                                                                     cccGGTTATCGGTAGCGAT‑ACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGa               >  2:826655/1‑100 (MQ=38)
                                                                                     cccGGTTATCGGTACCGAT‑ACCGGGCCAACGGACTGCACAGTTTGCCGTTATTTGTCGAACAGATAATGGCTTACCAGGTTTTCCACTTGTTACTGGCg                >  1:3305594/1‑99 (MQ=255)
                                                                                      ccggttATCGGTAGCGAT‑ACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGAc              >  2:2804526/6‑100 (MQ=255)
                                                                                      ccGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGan              <  1:201319/100‑2 (MQ=255)
                                                                                      ccGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGAc              <  2:3141431/100‑1 (MQ=255)
                                                                                      cagg‑TATCGCTACAGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGAc              >  2:2966454/3‑100 (MQ=255)
                                                                                       cGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGAcc             >  2:3278728/1‑100 (MQ=255)
                                                                                       cGG‑TATCGCTACCGATAACCGGGCCAACGGACTACACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGAcc             >  2:3582684/1‑100 (MQ=255)
                                                                                         g‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCAc           <  1:2390101/100‑1 (MQ=255)
                                                                                              cGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCGGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTg       >  2:2448543/1‑100 (MQ=255)
                                                                                              cGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCActc         >  1:3982408/1‑98 (MQ=255)
                                                                                               gCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGa      >  2:3548153/1‑100 (MQ=255)
                                                                                                cTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGAt     <  1:3819874/100‑1 (MQ=255)
                                                                                                 tACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATg    <  2:3738597/100‑1 (MQ=255)
                                                                                                 tACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATg    <  1:2463129/100‑1 (MQ=255)
                                                                                                 tACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATg    <  1:1841668/100‑1 (MQ=255)
                                                                                                  aCCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGc   <  2:3313401/100‑1 (MQ=255)
                                                                                                  aCCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGc   <  1:3455741/100‑1 (MQ=255)
                                                                                                  aCCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGc   <  1:3811380/100‑1 (MQ=255)
                                                                                                  aCCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGc   >  2:1864017/1‑100 (MQ=255)
                                                                                                   ccGATAACCGGGCCAACGGACTGCACAGTTAGTCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCa  >  2:1269170/1‑100 (MQ=255)
                                                                                                   ccGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCa  <  1:4168117/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
CCCTGCTCGTTGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCA  >  NC_000913/3729310‑3729507

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