Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 73,833 A→G I230T (ATC→ACC)  yabK ← putative transport system permease protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091373,8330AG99.5% 628.4 / ‑5.1 200I230T (ATC→ACC) yabKputative transport system permease protein
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (104/95);  total (105/95)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.54e-01

Only 100 of 205 total aligned reads displayed.

AACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCA  >  NC_000913/73743‑73922
                                                                                          |                                                                                         
aaCAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGGtt                                                                                         <  2:3229880/93‑1 (MQ=255)
aaCAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGGtt                                                                                         >  2:3541029/1‑93 (MQ=255)
aaCAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGGtt                                                                                         <  1:3420889/93‑1 (MQ=255)
aaCAGCACCAGCCCGAGGCAGCACACCATCTGGAGAAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGGtt                                                                                         <  1:2063127/93‑1 (MQ=255)
 aCAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGGTTg                                                                                        >  2:2880792/1‑93 (MQ=255)
 aCAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGGTTg                                                                                        >  2:1814672/1‑93 (MQ=255)
 aCAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGGTTg                                                                                        <  1:1601326/93‑1 (MQ=255)
   aGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGGTTGcc                                                                                      <  1:1948693/93‑1 (MQ=255)
    gCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGGTTGCCa                                                                                     >  2:2484810/1‑93 (MQ=255)
    gCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGGTTGCCa                                                                                     >  1:930430/1‑93 (MQ=255)
    gCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGGTTGCCa                                                                                     >  1:3862149/1‑93 (MQ=255)
    gCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGGTTGCCa                                                                                     <  1:34422/93‑1 (MQ=255)
        cAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGGTTGCCAGCTc                                                                                 >  2:2552228/1‑93 (MQ=255)
         aGCCCGAGGCAGCACACCATCTGGAGCGGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGGGTGCCAGCTCg                                                                                >  1:2593406/1‑93 (MQ=255)
          gCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGa                                                                               <  2:2620236/93‑1 (MQ=255)
           cccGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGAt                                                                              <  2:426119/93‑1 (MQ=255)
           cccGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGAt                                                                              <  2:651271/93‑1 (MQ=255)
           cccGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGAt                                                                              >  1:2081343/1‑93 (MQ=255)
            ccGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGATa                                                                             >  2:4272957/1‑93 (MQ=255)
              gAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGt                                                                           <  2:861836/93‑1 (MQ=255)
                ggCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGTgg                                                                         <  2:2174230/93‑1 (MQ=255)
                  cagcaCACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTc                                                                       <  1:119407/93‑1 (MQ=255)
                  cagcaCACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTc                                                                       >  2:1529139/1‑93 (MQ=255)
                        acCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTgg                                                                 >  1:3232341/1‑92 (MQ=255)
                           aTCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGa                                                              <  1:3657950/93‑1 (MQ=255)
                            tCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGGATGCCAGCTCGATAGTGGTCGCCTGCGGAc                                                             >  1:4587968/1‑93 (MQ=255)
                             cTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGAcc                                                            <  2:1078890/93‑1 (MQ=255)
                              tGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCg                                                           >  1:4067645/1‑93 (MQ=255)
                              tGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCg                                                           >  1:2283753/1‑93 (MQ=255)
                                gAGCAGCCCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCCCCTGCGGACCGcc                                                         <  1:399227/93‑1 (MQ=255)
                                     gcgcCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCa                                                    <  2:4143957/93‑1 (MQ=255)
                                     gcgcCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCa                                                    <  2:894087/93‑1 (MQ=255)
                                     gcgcCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCa                                                    >  1:3914941/1‑93 (MQ=255)
                                     gcgcCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCa                                                    <  1:4375192/93‑1 (MQ=255)
                                          aGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGAt                                               >  2:3479539/1‑93 (MQ=255)
                                           gCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATa                                              >  2:2980348/1‑93 (MQ=255)
                                           gCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATa                                              <  1:4511907/93‑1 (MQ=255)
                                             aTTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGc                                            <  1:837623/93‑1 (MQ=255)
                                             aTTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGc                                            >  1:2735220/1‑93 (MQ=255)
                                              ttGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCa                                           >  2:1961548/1‑93 (MQ=255)
                                               tGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCAc                                          <  2:4333436/93‑1 (MQ=255)
                                                   gcgcgGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTg                                      <  2:77891/93‑1 (MQ=255)
                                                   gcgcgGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTg                                      >  2:4467483/1‑93 (MQ=255)
                                                      cgGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCg                                   <  2:1708931/93‑1 (MQ=255)
                                                       gggCAGGATCGAGTCGTAACTCAGCGCCTGATAGGGTGCCAGCTCGATAGTGGTCGCCTGCGGACCCCCCCCCCGCGATAGCACGGGGGCGaa                                 >  2:3942442/1‑93 (MQ=255)
                                                        ggCAGGATCGTAGTCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCgca                                 >  2:362843/1‑91 (MQ=255)
                                                        ggCAGGATCGTAGTCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGaa                                 >  2:4204303/1‑93 (MQ=255)
                                                        ggCAGGATCGTAGTCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGaa                                 <  1:2187454/93‑1 (MQ=255)
                                                        ggCAGGATCGTAGTCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGaa                                 <  2:33928/93‑1 (MQ=255)
                                                         gCAGGATCGTAGTCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAg                                <  2:4107425/93‑1 (MQ=255)
                                                          cAGGATCGTAGTCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGc                               <  2:602647/93‑1 (MQ=255)
                                                          cAGGATCGTAGTCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGc                               <  2:3844361/93‑1 (MQ=255)
                                                           aGGATCGTAGTCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCTGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCt                              >  1:769916/1‑93 (MQ=255)
                                                           aGGATCGTAGTCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCt                              >  2:1350541/1‑93 (MQ=255)
                                                           aGGATCGTAGTCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCt                              >  2:2839164/1‑93 (MQ=255)
                                                           aGGATCGTAGTCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCt                              <  2:1571630/93‑1 (MQ=255)
                                                            ggATCGTAGTCGTATCTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTg                             <  1:1544749/93‑1 (MQ=255)
                                                            ggATCGTAGTCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTg                             >  1:2489535/1‑93 (MQ=255)
                                                            ggATCGTAGTCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTg                             >  1:2593343/1‑93 (MQ=255)
                                                             gATCGTAGTCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTgg                            <  1:172316/93‑1 (MQ=255)
                                                              aTCGTAGTCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGc                           <  2:4296747/93‑1 (MQ=255)
                                                              aTCGTAGTCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGc                           >  1:3174683/1‑93 (MQ=255)
                                                              aTCGTAGTCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGc                           >  1:2822862/1‑93 (MQ=255)
                                                                cGTAGTCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGCGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGa                         >  1:132492/1‑93 (MQ=255)
                                                                 gTAGTCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGaa                        >  1:1571578/1‑93 (MQ=255)
                                                                 gTAGTCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGaa                        >  1:3244857/1‑93 (MQ=255)
                                                                 gTAGTCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGaa                        >  2:2776297/1‑93 (MQ=255)
                                                                  tAGTCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGaga                       >  1:1432230/1‑91 (MQ=255)
                                                                  tAGTCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGaaa                       <  2:1305468/93‑1 (MQ=255)
                                                                  tAGTCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGaaa                       >  1:381308/1‑93 (MQ=255)
                                                                    gTCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACa                     >  1:3560208/1‑93 (MQ=255)
                                                                     tCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACag                    >  2:3450523/1‑93 (MQ=255)
                                                                     tCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACag                    >  1:984782/1‑93 (MQ=255)
                                                                      cGTAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACaga                   <  1:3761675/93‑1 (MQ=255)
                                                                      cGTAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACaga                   <  1:2449266/93‑1 (MQ=255)
                                                                       gTAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACagag                  >  1:3353735/1‑93 (MQ=255)
                                                                       gTAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACagag                  <  1:1807111/93‑1 (MQ=255)
                                                                       gTAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACagag                  >  2:1744802/1‑93 (MQ=255)
                                                                        tAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGc                 <  1:3305942/93‑1 (MQ=255)
                                                                         aaCTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCa                <  2:2432776/93‑1 (MQ=255)
                                                                          aCTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCAt               >  2:47723/1‑93 (MQ=255)
                                                                          aCTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCAt               >  2:967268/1‑93 (MQ=255)
                                                                           cTCAGGGCCTGATAGGTTGCCGGGCCGATAGGGGTGCCCCGCGGACCGCCCCCCAGCAATGGCACGGTGGCGAAGTCGGCGAAACAGTACATa              >  2:760829/1‑93 (MQ=255)
                                                                           cTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATa              >  1:1142747/1‑93 (MQ=255)
                                                                           cTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATa              <  1:1280733/93‑1 (MQ=255)
                                                                              aGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAg           >  2:2109276/1‑93 (MQ=255)
                                                                                cgcCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGAt         <  2:1482042/93‑1 (MQ=255)
                                                                                cgcCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGAt         <  2:3070968/93‑1 (MQ=255)
                                                                                 gcCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATa        >  2:1212115/1‑93 (MQ=255)
                                                                                  ccTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATaa       <  2:1134458/93‑1 (MQ=255)
                                                                                  ccTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATaa       >  2:3803903/1‑93 (MQ=255)
                                                                                  ccTGATAGGTTGCCAGCTCGATAGTGGACGCCTGCGGACCGCGCCCCAGCGACAGCACGGTGGGGAAGCTGGCGAAACAGAGCATAAAGATaa       >  2:2346667/1‑93 (MQ=255)
                                                                                   cTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAg      >  1:645349/1‑93 (MQ=255)
                                                                                   cTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAg      <  2:1790886/93‑1 (MQ=255)
                                                                                    tGATAGGTTGCGAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAgc     <  2:4394242/93‑1 (MQ=255)
                                                                                    tGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAgc     >  1:3793616/1‑93 (MQ=255)
                                                                                    tGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAgc     >  2:2130367/1‑93 (MQ=255)
                                                                                     gATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAgcg    >  1:4226572/1‑93 (MQ=255)
                                                                                     gATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAgag    >  2:2292547/1‑91 (MQ=255)
                                                                                       tAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCgca  >  1:3520123/1‑93 (MQ=255)
                                                                                          |                                                                                         
AACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCA  >  NC_000913/73743‑73922

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