Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,729,400:1 +T 15.2% intergenic (‑29/+43) 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,3991.T15.2% 563.7 / 107.6 171intergenic (‑28/+44)xylB/xylAxylulokinase/D‑xylose isomerase
Reads supporting (aligned to +/- strand):  ref base . (78/67);  new base T (17/9);  total (95/76)
Fisher's exact test for biased strand distribution p-value = 2.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.53e-01

Only 100 of 179 total aligned reads displayed.

AACCACCTCACCCTGCTCGTTGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCAC  >  NC_000913/3729300‑3729498
                                                                                                    |                                                                                                   
aaCCACCTCACCCTGCTCGTTGAGCAAAATAACTTTTACGCCCGAGGTGCCGAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCgg                                                                                                      <  1:858254/100‑1 (MQ=255)
aaCCACCTCACCCTGCTCGTTGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCgg                                                                                                      <  1:881683/100‑1 (MQ=255)
 accaccTCACCCTGCTCGTTGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGt                                                                                                    >  2:4403338/1‑100 (MQ=255)
 accaccTCACCCTGCGCGTTGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGt                                                                                                    <  2:4196007/100‑1 (MQ=255)
    accTCACCCTGCTTGTTGAGCAAAATAACTTTTACGCCCGAGGTGTCAAGATCTATCCCGATATATATATCGATCGTTCCTTAAAAAAATGCACGGTATc                                                                                                 >  1:3917138/1‑100 (MQ=255)
     ccTCACCCTGCTCGTTGAGCAAAATAACTTTTACGCCCTAGGTGCCACGATCTATCCCGATATACATATCGATCGTTCCTTACAAAAATTCCCCGTATCg                                                                                                >  1:1366672/1‑100 (MQ=255)
     ccTCACCCTGCTCGTTGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCg                                                                                                <  2:2093537/100‑1 (MQ=255)
      cTCACCCTGCTCGTTGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGc                                                                                               <  2:1699420/100‑1 (MQ=255)
       tCACCCTGCTCGTTGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCt                                                                                              >  1:1032939/1‑100 (MQ=255)
       tCACCCTGCTCGTTGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCt                                                                                              <  1:217070/100‑1 (MQ=255)
        cACCCTGCTCGTTGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTa                                                                                             <  2:3905833/100‑1 (MQ=255)
        cACCCTGCTCGTTGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTa                                                                                             <  2:3924923/100‑1 (MQ=255)
          cccTTCTCGTTGATCAAAATAACTTTTACTCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTAcc                                                                                           <  2:4125508/100‑1 (MQ=255)
          cccTGCTCGTTGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTAcc                                                                                           <  1:1193737/100‑1 (MQ=255)
          cccTGCTCGTTGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTAcc                                                                                           >  1:3380677/1‑100 (MQ=255)
          accTGCTAGTTGAGCAAAATAGCGTTTACGCCCGAGGTGCCAAGATCTATCACGATATACATATCTATCGTTACCTAAACAACTTCACGGTATctc                                                                                               >  1:4265669/2‑94 (MQ=255)
             tGCTCGTTGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGAt                                                                                        >  2:4169511/1‑100 (MQ=255)
                tCGTTGAGCAAAATAACTTTTACGACCGAGGTGCCAAGATCTATCCCGATATACATATCGATAAATACTAAAAAAAATGCCCGGTATCGCTACCGATa                                                                                       <  1:3931117/98‑1 (MQ=255)
                 cGTTGAGCAAAATAACTTTTACGCAAGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAACAAAATGCCCAGTATCGCTACCGATAAcc                                                                                    >  1:2001449/1‑100 (MQ=255)
                    tGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTTATCGGTAGCGATACCggg                                                                                 <  2:448060/100‑1 (MQ=38)
                    tGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCggg                                                                                 >  1:3727094/1‑100 (MQ=255)
                     gAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGc                                                                                >  2:1328974/1‑100 (MQ=255)
                      aGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGcc                                                                               <  1:1478438/100‑1 (MQ=255)
                         aaaaTAACTTTTACGCCCGAGGTGACAAGATCTATCCCGGTATACATATCGATCGTTCCTTAAAAAAATGCCCGGTTATCGGTAGCGATACCGGGcattt                                                                            <  1:3801320/100‑5 (MQ=25)
                           aaTAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACgg                                                                          <  1:1275612/100‑1 (MQ=255)
                           aaTAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACgg                                                                          <  2:1202554/100‑1 (MQ=255)
                           aaTAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACgg                                                                          <  2:1155649/100‑1 (MQ=255)
                            attaCTTTTACGCTCGAGGTGCCAAGATCTATGCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTTATCGGTAGCGATACCGGGcatttttt                                                                         <  1:346447/97‑8 (MQ=21)
                             tAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTTATCGGTAGCGATACCGGGcattttttt                                                                        <  1:4397985/100‑9 (MQ=25)
                               aCTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACt                                                                       >  1:665857/1‑99 (MQ=255)
                                cTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGc                                                                     >  2:316909/1‑100 (MQ=255)
                                cTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGc                                                                     >  2:3928947/1‑100 (MQ=255)
                                 ttttACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGca                                                                    <  1:189327/100‑1 (MQ=255)
                                  tttACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGcac                                                                   <  1:1839701/100‑1 (MQ=255)
                                      cGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGtt                                                               >  2:4011116/1‑100 (MQ=255)
                                      cGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGtt                                                               >  1:2588937/1‑100 (MQ=255)
                                      cGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGtt                                                               <  1:2296184/100‑1 (MQ=255)
                                       gCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTa                                                              <  1:2804504/100‑1 (MQ=255)
                                        cccGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAg                                                             >  1:559510/1‑100 (MQ=255)
                                        cccGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCATTAAAAAAATGCCCGGTATCGCTACAGATAACCGGGCCAACGGAATGCACCGTTAg                                                             >  2:364083/1‑100 (MQ=255)
                                         ccGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGc                                                            <  2:846648/100‑1 (MQ=255)
                                         ccGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGc                                                            >  1:1148236/1‑100 (MQ=255)
                                         ccGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGc                                                            <  1:4067028/100‑1 (MQ=255)
                                          cGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGcc                                                           <  2:1462020/100‑1 (MQ=255)
                                          cGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGc                                                            <  1:3054689/99‑1 (MQ=255)
                                           gAGGTGCTAAAATCTCGCCCGATATACATATCGATCGTGCCTTGAAAAAATGCCCGGTGTCGCTACCGATAACCGGGCGAACGGACTGCACAGTTAGCCg                                                          <  1:529375/100‑1 (MQ=255)
                                               tGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGtttt                                                      <  2:4093656/100‑3 (MQ=255)
                                               tGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTa                                                       <  1:3291079/99‑1 (MQ=255)
                                                gCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTAtt                                                     >  2:258506/1‑100 (MQ=255)
                                                gCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTAtt                                                     >  2:704978/1‑100 (MQ=255)
                                                    aGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTc                                                 <  2:2347115/100‑1 (MQ=255)
                                                     gATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCg                                                <  1:3623431/100‑1 (MQ=255)
                                                     gATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCg                                                >  2:1077201/1‑100 (MQ=255)
                                                       tCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGaa                                              <  2:2956115/100‑1 (MQ=255)
                                                        cTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAAc                                             <  1:156067/100‑1 (MQ=255)
                                                             ccGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATa                                        >  2:2536930/1‑100 (MQ=255)
                                                             ccGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCAGTTAATTGTCGacc                                             >  1:3002227/1‑93 (MQ=255)
                                                                 tataCATATCGATCGTTCCTTAAAAAAATGCCCGGTAGCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATgg                                    >  1:2354412/1‑100 (MQ=255)
                                                                  ataCATATCGATCGTTCCTTAAAAAAATGCCCGGTAGCGCTACCGATATCCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGt                                   <  2:3101833/100‑1 (MQ=255)
                                                                   taCATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGtt                                  >  2:3735699/1‑100 (MQ=255)
                                                                   taCATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGtt                                  >  2:2796441/1‑100 (MQ=255)
                                                                   taCATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGtt                                  >  1:789805/1‑100 (MQ=255)
                                                                   taCATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGATACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGtt                                  >  1:3522131/1‑100 (MQ=255)
                                                                   naCATATCGATCGTTCCTTACAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCATTATTTGTCGAACAGATAATgg                                    >  1:1624053/2‑98 (MQ=255)
                                                                     cATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTa                                >  1:2505580/1‑100 (MQ=255)
                                                                     cATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTa                                >  1:3855714/1‑100 (MQ=255)
                                                                     cATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACAGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTa                                >  2:2711342/1‑100 (MQ=255)
                                                                      atatCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTAc                               <  2:934309/100‑1 (MQ=255)
                                                                      atatCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTAc                               <  2:311203/100‑1 (MQ=255)
                                                                       tatCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTAcc                              >  1:2213339/1‑100 (MQ=255)
                                                                       tatCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGATTccc                              >  2:2060783/1‑100 (MQ=255)
                                                                          cgagcgTTCCTTAAAAAAATGCCCGGTATCGGTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTTGAACAGATAATGGTTTACCAGa                           <  1:2558439/100‑1 (MQ=255)
                                                                              cgTTCCATAAAAAAATGCCCGGTATCGCAACCGATAACCGGGCCAACGGACTGCACAGTTATCCGTTATTTGTCGAACAGATAATGGTTTACCAGAtttt                       <  1:1387656/100‑1 (MQ=255)
                                                                               gTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGAtttt                       >  1:676023/1‑99 (MQ=255)
                                                                                  ccTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAg                   <  2:249729/100‑1 (MQ=255)
                                                                                  ccTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAg                   >  1:630957/1‑100 (MQ=255)
                                                                                  ccTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAg                   >  1:1925291/1‑100 (MQ=255)
                                                                                   cTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAgt                  <  2:2643806/100‑1 (MQ=255)
                                                                                   cTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAgt                  >  2:2666931/1‑100 (MQ=255)
                                                                                   cTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAgt                  >  2:3024253/1‑100 (MQ=255)
                                                                                     tAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTTCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAgttg                <  1:2843335/100‑1 (MQ=255)
                                                                                     tAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAgttg                >  2:1873925/1‑100 (MQ=255)
                                                                                      aaaaaaaTGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAgttgt               <  1:3180636/100‑1 (MQ=255)
                                                                                        ccgttggCCCGGTTATCGGTAGCGATACCGGGCCAACGGACTGCACAGTTAGCCGTTATTCGTCGAACAGCTACTGGTTTACCAGATTTTCCAGTTGTTc             >  1:4338477/7‑100 (MQ=21)
                                                                                        aaaaaTGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTc             >  2:2442521/1‑100 (MQ=255)
                                                                                        aaaaaTGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTc             >  2:3427202/1‑100 (MQ=255)
                                                                                            aTGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTgg         <  1:3560413/100‑1 (MQ=255)
                                                                                              gCCCGGTTATCGGTAGCGATACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCg       >  2:1483470/1‑100 (MQ=38)
                                                                                              gCCCGGTTATCGGTAGCGATACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCg       >  2:3502265/1‑100 (MQ=38)
                                                                                              gCCCGGTTATCGGTAGCGATACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCg       >  2:4293348/1‑100 (MQ=38)
                                                                                              gCCCGGTTATCGGTAGCGATACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTAGTTCCTGGCg       >  1:2893359/1‑100 (MQ=37)
                                                                                               cccGGTTATCGGTAGCGATACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGa      >  2:2239613/1‑100 (MQ=38)
                                                                                               cccGGTTATCGGTAGCGATACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGa      >  2:2758796/1‑100 (MQ=38)
                                                                                               cccGGTTATCGGTACCGATACCGGGCCAACGGACTGCACAGTTTGCCGTTATTTGTCGAACAGATAATGGCTTACCAGGTTTTCCACTTGTTACTGGCg       >  1:3305594/1‑99 (MQ=255)
                                                                                                ccGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGan     <  1:201319/100‑2 (MQ=255)
                                                                                                ccGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGAc     >  2:110063/1‑100 (MQ=255)
                                                                                                ccGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGAc     <  2:3141431/100‑1 (MQ=255)
                                                                                                 cGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGAcc    >  2:3278728/1‑100 (MQ=255)
                                                                                                 cGGTATCGCTACCGATAACCGGGCCAACGGACTACACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGAcc    >  2:3582684/1‑100 (MQ=255)
                                                                                                   gTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCAc  <  1:2390101/100‑1 (MQ=255)
                                                                                                    |                                                                                                   
AACCACCTCACCCTGCTCGTTGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCAC  >  NC_000913/3729300‑3729498

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: