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.

AACCACCTCACCCTGCTCGTTGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCAC  >  NC_000913/3729300‑3729498
                                                                                                    |                                                                                                   
aaCCACCTCACCCTGCTCGTTGAGCAAAATAACTTTTACGCCCGAGGTGCCGAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCgg                                                                                                      <  1:858254/100‑1 (MQ=255)
aaCCACCTCACCCTGCTCGTTGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCgg                                                                                                      <  1:881683/100‑1 (MQ=255)
 accaccTCACCCTGCTCGTTGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑t                                                                                                    >  2:4403338/1‑100 (MQ=255)
 accaccTCACCCTGCGCGTTGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑t                                                                                                    <  2:4196007/100‑1 (MQ=255)
    accTCACCCTGCTTGTTGAGCAAAATAACTTTTACGCCCGAGGTGTCAAGATCTATCCCGATATATATATCGATCGTTCCTTAAAAAAATGCACGG‑TATc                                                                                                 >  1:3917138/1‑100 (MQ=255)
     ccTCACCCTGCTCGTTGAGCAAAATAACTTTTACGCCCTAGGTGCCACGATCTATCCCGATATACATATCGATCGTTCCTTACAAAAATTCCCCG‑TATCg                                                                                                >  1:1366672/1‑100 (MQ=255)
     ccTCACCCTGCTCGTTGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCg                                                                                                <  2:2093537/100‑1 (MQ=255)
      cTCACCCTGCTCGTTGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGc                                                                                               <  2:1699420/100‑1 (MQ=255)
       tCACCCTGCTCGTTGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCt                                                                                              >  1:1032939/1‑100 (MQ=255)
       tCACCCTGCTCGTTGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCt                                                                                              <  1:217070/100‑1 (MQ=255)
        cACCCTGCTCGTTGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTa                                                                                             <  2:3905833/100‑1 (MQ=255)
        cACCCTGCTCGTTGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTa                                                                                             <  2:3924923/100‑1 (MQ=255)
          cccTTCTCGTTGATCAAAATAACTTTTACTCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTAcc                                                                                           <  2:4125508/100‑1 (MQ=255)
          cccTGCTCGTTGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTAcc                                                                                           <  1:1193737/100‑1 (MQ=255)
          cccTGCTCGTTGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTAcc                                                                                           >  1:3380677/1‑100 (MQ=255)
          accTGCTAGTTGAGCAAAATAGCGTTTACGCCCGAGGTGCCAAGATCTATCACGATATACATATCTATCGTTACCTAAACAACTTCACGG‑TATctc                                                                                               >  1:4265669/2‑94 (MQ=255)
             tGCTCGTTGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGAt                                                                                        >  2:4169511/1‑100 (MQ=255)
                tCGTTGAGCAAAATAACTTTTACGACCGAGGTGCCAAGATCTATCCCGATATACATATCGATAAATACTAAAAAAAATGCCCGG‑TATCGCTACCGATa                                                                                       <  1:3931117/98‑1 (MQ=255)
                 cGTTGAGCAAAATAACTTTTACGCAAGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAACAAAATGCCCAG‑TATCGCTACCGATAAcc                                                                                    >  1:2001449/1‑100 (MQ=255)
                    tGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTTATCGGTAGCGAT‑ACCggg                                                                                 <  2:448060/100‑1 (MQ=38)
                    tGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCggg                                                                                 >  1:3727094/1‑100 (MQ=255)
                     gAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGc                                                                                >  2:1328974/1‑100 (MQ=255)
                      aGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGcc                                                                               <  1:1478438/100‑1 (MQ=255)
                         aaaaTAACTTTTACGCCCGAGGTGACAAGATCTATCCCGGTATACATATCGATCGTTCCTTAAAAAAATGCCCGGTTATCGGTAGCGAT‑ACCGGGcattt                                                                            <  1:3801320/100‑5 (MQ=25)
                           aaTAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACgg                                                                          <  1:1275612/100‑1 (MQ=255)
                           aaTAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACgg                                                                          <  2:1202554/100‑1 (MQ=255)
                           aaTAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACgg                                                                          <  2:1155649/100‑1 (MQ=255)
                            attaCTTTTACGCTCGAGGTGCCAAGATCTATGCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTTATCGGTAGCGAT‑ACCGGGcatttttt                                                                         <  1:346447/97‑8 (MQ=21)
                             tAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTTATCGGTAGCGAT‑ACCGGGcattttttt                                                                        <  1:4397985/100‑9 (MQ=25)
                               aCTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACt                                                                       >  1:665857/1‑99 (MQ=255)
                                cTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGc                                                                     >  2:316909/1‑100 (MQ=255)
                                cTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGc                                                                     >  2:3928947/1‑100 (MQ=255)
                                 ttttACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGca                                                                    <  1:189327/100‑1 (MQ=255)
                                  tttACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGcac                                                                   <  1:1839701/100‑1 (MQ=255)
                                      cGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGtt                                                               >  2:4011116/1‑100 (MQ=255)
                                      cGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGtt                                                               >  1:2588937/1‑100 (MQ=255)
                                      cGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGtt                                                               <  1:2296184/100‑1 (MQ=255)
                                       gCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTa                                                              <  1:2804504/100‑1 (MQ=255)
                                        cccGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAg                                                             >  1:559510/1‑100 (MQ=255)
                                        cccGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCATTAAAAAAATGCCCGG‑TATCGCTACAGATAACCGGGCCAACGGAATGCACCGTTAg                                                             >  2:364083/1‑100 (MQ=255)
                                         ccGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGc                                                            <  2:846648/100‑1 (MQ=255)
                                         ccGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGc                                                            >  1:1148236/1‑100 (MQ=255)
                                         ccGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGc                                                            <  1:4067028/100‑1 (MQ=255)
                                          cGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGcc                                                           <  2:1462020/100‑1 (MQ=255)
                                          cGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGc                                                            <  1:3054689/99‑1 (MQ=255)
                                           gAGGTGCTAAAATCTCGCCCGATATACATATCGATCGTGCCTTGAAAAAATGCCCGG‑TGTCGCTACCGATAACCGGGCGAACGGACTGCACAGTTAGCCg                                                          <  1:529375/100‑1 (MQ=255)
                                               tGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGtttt                                                      <  2:4093656/100‑3 (MQ=255)
                                               tGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTa                                                       <  1:3291079/99‑1 (MQ=255)
                                                gCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTAtt                                                     >  2:258506/1‑100 (MQ=255)
                                                gCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTAtt                                                     >  2:704978/1‑100 (MQ=255)
                                                    aGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTc                                                 <  2:2347115/100‑1 (MQ=255)
                                                     gATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCg                                                <  1:3623431/100‑1 (MQ=255)
                                                     gATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCg                                                >  2:1077201/1‑100 (MQ=255)
                                                       tCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGaa                                              <  2:2956115/100‑1 (MQ=255)
                                                        cTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAAc                                             <  1:156067/100‑1 (MQ=255)
                                                             ccGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATa                                        >  2:2536930/1‑100 (MQ=255)
                                                             ccGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCAGTTAATTGTCGacc                                             >  1:3002227/1‑93 (MQ=255)
                                                                 tataCATATCGATCGTTCCTTAAAAAAATGCCCGG‑TAGCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATgg                                    >  1:2354412/1‑100 (MQ=255)
                                                                  ataCATATCGATCGTTCCTTAAAAAAATGCCCGG‑TAGCGCTACCGATATCCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGt                                   <  2:3101833/100‑1 (MQ=255)
                                                                   taCATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGtt                                  >  2:3735699/1‑100 (MQ=255)
                                                                   taCATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGtt                                  >  2:2796441/1‑100 (MQ=255)
                                                                   taCATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGtt                                  >  1:789805/1‑100 (MQ=255)
                                                                   taCATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGATACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGtt                                  >  1:3522131/1‑100 (MQ=255)
                                                                   naCATATCGATCGTTCCTTACAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCATTATTTGTCGAACAGATAATgg                                    >  1:1624053/2‑98 (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)
                                                                      atatCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTAc                               <  2:934309/100‑1 (MQ=255)
                                                                      atatCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTAc                               <  2:311203/100‑1 (MQ=255)
                                                                       tatCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTAcc                              >  1:2213339/1‑100 (MQ=255)
                                                                       tatCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGATTccc                              >  2:2060783/1‑100 (MQ=255)
                                                                          cgagcgTTCCTTAAAAAAATGCCCGG‑TATCGGTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTTGAACAGATAATGGTTTACCAGa                           <  1:2558439/100‑1 (MQ=255)
                                                                              cgTTCCATAAAAAAATGCCCGG‑TATCGCAACCGATAACCGGGCCAACGGACTGCACAGTTATCCGTTATTTGTCGAACAGATAATGGTTTACCAGAtttt                       <  1:1387656/100‑1 (MQ=255)
                                                                               gTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGAtttt                       >  1:676023/1‑99 (MQ=255)
                                                                                  ccTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAg                   <  2:249729/100‑1 (MQ=255)
                                                                                  ccTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAg                   >  1:630957/1‑100 (MQ=255)
                                                                                  ccTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAg                   >  1:1925291/1‑100 (MQ=255)
                                                                                   cTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAgt                  <  2:2643806/100‑1 (MQ=255)
                                                                                   cTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAgt                  >  2:2666931/1‑100 (MQ=255)
                                                                                   cTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAgt                  >  2:3024253/1‑100 (MQ=255)
                                                                                     tAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTTCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAgttg                <  1:2843335/100‑1 (MQ=255)
                                                                                     tAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAgttg                >  2:1873925/1‑100 (MQ=255)
                                                                                      aaaaaaaTGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAgttgt               <  1:3180636/100‑1 (MQ=255)
                                                                                        ccgttggCCCGGTTATCGGTAGCGAT‑ACCGGGCCAACGGACTGCACAGTTAGCCGTTATTCGTCGAACAGCTACTGGTTTACCAGATTTTCCAGTTGTTc             >  1:4338477/7‑100 (MQ=21)
                                                                                        aaaaaTGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTc             >  2:2442521/1‑100 (MQ=255)
                                                                                        aaaaaTGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTc             >  2:3427202/1‑100 (MQ=255)
                                                                                            aTGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTgg         <  1:3560413/100‑1 (MQ=255)
                                                                                              gCCCGGTTATCGGTAGCGAT‑ACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCg       >  2:1483470/1‑100 (MQ=38)
                                                                                              gCCCGGTTATCGGTAGCGAT‑ACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCg       >  2:3502265/1‑100 (MQ=38)
                                                                                              gCCCGGTTATCGGTAGCGAT‑ACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCg       >  2:4293348/1‑100 (MQ=38)
                                                                                              gCCCGGTTATCGGTAGCGAT‑ACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTAGTTCCTGGCg       >  1:2893359/1‑100 (MQ=37)
                                                                                               cccGGTTATCGGTAGCGAT‑ACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGa      >  2:2239613/1‑100 (MQ=38)
                                                                                               cccGGTTATCGGTAGCGAT‑ACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGa      >  2:2758796/1‑100 (MQ=38)
                                                                                               cccGGTTATCGGTACCGAT‑ACCGGGCCAACGGACTGCACAGTTTGCCGTTATTTGTCGAACAGATAATGGCTTACCAGGTTTTCCACTTGTTACTGGCg       >  1:3305594/1‑99 (MQ=255)
                                                                                                ccGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGan     <  1:201319/100‑2 (MQ=255)
                                                                                                ccGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGAc     >  2:110063/1‑100 (MQ=255)
                                                                                                ccGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGAc     <  2:3141431/100‑1 (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)
                                                                                                    |                                                                                                   
AACCACCTCACCCTGCTCGTTGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCAC  >  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: ‑