Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,812,035 G→A intergenic (+21/‑86) dut → / → b3641 deoxyuridinetriphosphatase/putative transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,812,0350GA99.8% 1699.5 / ‑5.7 468intergenic (+21/‑86)dut/b3641deoxyuridinetriphosphatase/putative transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (210/256);  total (210/258)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

Only 100 of 477 total aligned reads displayed.

TAATCTGGTGGAAGATTTCGACGCCACCGACCGCGGTGAAGGCGGCTTTGGTCACTCTGGTCGTCAGTAACACATACGCATCCGAATAACGTCATAACATAGCCGCAAACATTTCGTTTGCGGTCATAGCGTGGGTGCCGCCTGGCAAGTGCTTATTTTCAGGGGTATTTTGTAACATGG  >  NC_000913/3811945‑3812124
                                                                                          |                                                                                         
tattCTGGTGGAAGATTTCGACGCCACCGACCGCGGTTAAGGCGGCTTTGGTCACTCTGGTCGTCATTAACACATACGCATCCGAATAACATc                                                                                         >  2:7026331/4‑93 (MQ=255)
tAATCTGGTGGAAGATTTCGACGCCACCGACCGCGGTGAAGGCGGCTTTGGTCACTCTGGTCGTCAGTAACACATACGCATCCGAATAACATc                                                                                         <  1:2954281/93‑1 (MQ=255)
tAATCTGGTGGAAGATTTCGACGCCACCGACCGCGGTGAAGGCGGCTTTGGTCACTCTGGTCGTCAGTAACACATACGCATCCGAATAACATc                                                                                         <  2:6407098/93‑1 (MQ=255)
 aaTCTGGTGGAAGATTTCGACGCCACCGACCGCGGTGAAGGCGGCTTTGGTCACTCTGGTCGTCAGTAACACATACGCATCCGAATAACATCa                                                                                        <  2:4638272/93‑1 (MQ=255)
 aaTCTGGTGGAAGATTTCGACGCCACCGACCGCGGTGAAGGCGGCTTTGGTCACTCTGGTCGTCAGTAACACATACGCATCCGAATAACATCa                                                                                        <  1:5976916/93‑1 (MQ=255)
 aaTCTGGTGGAAGATTTCGACGCCACCGACCGCGGTGAAGGCGGCTTTGGTCACTCTGGTCGTCAGTAACACATACGCATCCGAATAACATCa                                                                                        >  1:1255225/1‑93 (MQ=255)
  aTCTGGTGGAAGATTTCGACGCCACCGACCGCGGTGAAGGCGGCTTTGGTCACTCTGGTCGTCAGTAACACATACGCATCCGAATAAcatcat                                                                                       <  1:1439921/93‑1 (MQ=255)
   tCTGGTGGAAGATTTCGACGCCACCGACCGCGGTGAAGGCGGCTTTGGTCACTCTGGTCGTCAGTAACACATACGCATCCGAATAACATCATa                                                                                      <  1:1028494/93‑1 (MQ=255)
   tCTGGTGGAAGATTTCGACGCCACCGACCGCGGTGAAGGCGGCTTTGGTCACTCTGGTCGTCAGTAACACATACGCATCCGAATAACATCATa                                                                                      <  2:3742637/93‑1 (MQ=255)
    cTGGTGGAAGATTTCGACGCCACCGACCGCGGTGAAGGCGGCTTTGGTCACTCTGGTCGTCAGTAACACATACGCATCCGAATAACATCATaa                                                                                     <  1:793205/93‑1 (MQ=255)
    cTGGTGGAAGATTTCGACGCCACCGACCGCGGTGAAGGCGGCTTTGGTCACTCTGGTCGTCAGTAACACATACGCATCCGAATAACATCATaa                                                                                     <  2:1557298/93‑1 (MQ=255)
    cTGGTGGAAGATTTCGACGCCACCGACCGCGGTGAAGGCGGCTTTGGTCACTCTGGTCGTCAGTAACACATACGCATCCGAATAACATCATaa                                                                                     >  1:3316204/1‑93 (MQ=255)
     tggtggAAGATTTCGACGCCACCGACCGCGGTGAAGGCGGCTTTGGTCACTCTGGTCGTCAGTAACACATACGCATCCGAATAACATCATAAc                                                                                    >  1:3209814/1‑93 (MQ=255)
      ggtggAAGATTTCGACGCCACCGACCGCGGTGAAGGCGGCTTTGGTCACTCTGGTCGTCAGTAACACATACGCATCCGAATAACATCATAACa                                                                                   <  2:5776814/93‑1 (MQ=255)
      ggtggAAGATTTCGACGCCACCGACCGCGGTGAAGGCGGCTTTGGTCACTCTGGTCGTCAGTAACACATACGCATCCGAATAACATCATAACa                                                                                   >  1:3103912/1‑93 (MQ=255)
       gtggAAGATTTCGACGCCACCGACCGCGGTGAAGGCGGCTTTGGTCACTCTGGTCGTCAGTAACACATACGCATCCGAATAACATCATAACAt                                                                                  >  2:6584134/1‑93 (MQ=255)
        tggAAGATTTCGACGCCACCGACCGCGGTGAAGGCGGCTTTGGTCACTCTGGTCGTCAGTAACACATACGCATCCGAATAACATCATAACATa                                                                                 >  1:4868484/1‑93 (MQ=255)
          gAAGATTTCGACGCCACCGACCGCGGTGAAGGCGGCTTTGGTCACTCTGGTCGTCAGTAACACATACGCATCCGAATAACATCATAACATAGc                                                                               <  2:3145809/93‑1 (MQ=255)
          gAAGATTTCGACGCCACCGACCGCGGTGAAGGCGGCTTTGGTCACTCTGGTCGTCAGTAACACATACGCATCCGAATAACATCATAACATAGc                                                                               <  1:6908088/93‑1 (MQ=255)
          gAAGATTTCGACGCCACCGACCGCGGTGAAGGCGGCTTTGGTCACTCTGGTCGTCAGTAACACATACGCATCCGAATAACATCATAACATAGc                                                                               <  2:7555490/93‑1 (MQ=255)
           aaGATTTCGACGCCACCGACCGCGGTGAAGGCGGCTTTGGTCACTCTGGTCGTCAGTAACACATACGCATCCGAATAACATCATAACATAGcc                                                                              <  1:4625119/93‑1 (MQ=255)
            aGATTTCGACGCCACCGACCGCGGTGAAGGCGGCTTTGGTCACTCTGGTCGTCAGTAACACATACGCATCCGAATAACATCATAACATAGCCg                                                                             >  1:5297649/1‑93 (MQ=255)
               tttCGACGCCACCGACCGCGGTGAAGGCGGCTTTGGTCACTCTGGTCGTCAGTAACACATACGCATCCGAATAACATCATAACATAGCCGCaa                                                                          <  2:4924626/93‑1 (MQ=255)
               tttCGACGCCACCGACCGCGGTGAAGGCGGCTTTGGTCACTCTGGTCGTCAGTAACACATACGCATCCGAATAACATCATAACATAGCCGCaa                                                                          >  2:6110465/1‑93 (MQ=255)
                ttCGACGCCACCGACCGCGGTGAAGGCGGCTTTGGTCACTCTGGTCGTCAGTAACACATACGCATCCGAATAACATCATAACATAGCCGCaaa                                                                         <  2:2762334/93‑1 (MQ=255)
                ttCGACGCCACCGACCGCGGTGAAGGCGGCTTTGGTCACTCTGGTCGTCAGTAACACATACGCATCCGAATAACATCATAACATAGCCGCaaa                                                                         <  2:1945911/93‑1 (MQ=255)
                 tCGACGCCACCGACCGCGGTGAAGGCGGCTTTGGTCACTCTGGTCGTCAGTAACACATACGCATCCGAATAACATCATAACATAGCCGCAAAc                                                                        <  2:5559237/93‑1 (MQ=255)
                  cgccGCTACCGACCGCGGTGAAGGCGGCTTTGGTCACGCTGGTCGTCAGTAACACATACGCATCCGAATAACATCATAACATAGCCGCAAACa                                                                       <  2:308032/90‑1 (MQ=255)
                  cGACGCCACCGACCGCGGTGAAGGCGGCTTTGGTCACTCTGGTCGTCAGTAACACATACGCATCCGAATAACATCATAACATAGCCGCAAACa                                                                       <  2:7310968/93‑1 (MQ=255)
                    aCGCCACCGACCGCGGTGAAGGCGGCTTTGGTCACTCTGGTCGTCAGTAACACATACGCATCCGAATAACATCATAACATAGCCGCAAACAtt                                                                     <  1:5521441/93‑1 (MQ=255)
                       ccACCGACCGCGGTGAAGGCGGCTTTGGTCACTCTGGTCGTCAGTAACACATACGCATCCGAATAACATCATAACATAGCCGCAAACATTTCg                                                                  <  2:5351189/93‑1 (MQ=255)
                             accgCGGTGAAGGCGGCTTTGGTCACTCTGGTCGTCAGTAACACATACGCATCCGAATAACATCATAACATAGCCGCAAACATTTCGTTTGCg                                                            >  2:4532075/1‑93 (MQ=255)
                             accgCGGTGAAGGCGGCTTTGGTCACTCTGGTCGTCAGTAACACATACGCATCCGAATAACATCATAACATAGCCGCAAACATTTCGTTTGCg                                                            <  2:1139376/93‑1 (MQ=255)
                               cgcgGTGAAGGCGGCTTTGGTCACTCTGGTCGTCAGTAACACATACGCATCCGAATAACATCATAACATAGCCGCAAACATTTCGTTTGCGGt                                                          >  2:1854314/1‑93 (MQ=255)
                                gcgGTGAAGGCGGCTTTGGTCACTCTGGTCGTCAGTAACACATACGCATCCGAATAACATCATAACATAGCCGCAAACATTTCGTTTGCGGTc                                                         >  2:1227461/1‑93 (MQ=255)
                                 cgGTGAAGGCGGCTTTGGTCACTCTGGTCGTCAGTAACACATACGCATCCGAATAACATCATAACATAGCCGCAAACATTTCGTTTGCGGTCa                                                        >  2:522266/1‑93 (MQ=255)
                                 cgGTGAAGGCGGCTTTGGTCACTCTGGTCGTCAGTAACACATACGCATCCGAATAACATCATAACATAGCCGCAAACATTTCGTTTGCGGTCa                                                        <  2:4777427/93‑1 (MQ=255)
                                 cgGTGAAGGCGGCTTTGGTCACTCTGGTCGTCAGTAACACATACGCATCCGAATAACATCATAACATAGCCGCAAACATTTCGTTTGCCGTCa                                                        >  1:1716353/1‑93 (MQ=255)
                                  ggTGAAGGCGGCTTTGGTCACTCTGGTCGTCAGTAACACATACGCATCCGAATAACATCATAACATAGCCGCAAACATTTCGTTTGCGGTCAt                                                       >  1:3586101/1‑93 (MQ=255)
                                   gTGAAGGCGGCTTTGGTCACTCTGGTCGTCAGTAACACATACGCATCCGAATAACATCATAACATAGCCGCAAACATTTCGTTTGCGGTCATa                                                      <  2:854895/93‑1 (MQ=255)
                                   gTGAAGGCGGCTTTGGTCACTCTGGTCGTCAGTAACACATACGCATCCGAATAACATCATAACATAGCCGCAAACATTTCGTTTGCGGTCATa                                                      >  1:1070547/1‑93 (MQ=255)
                                   gTGAAGGCGGCTTTGGTCACTCTGGTCGTCAGTAACACATACGCATCCGAATAACATCATAACATAGCCGCAAACATTTCGTTTGCGGTCATa                                                      <  1:6096222/93‑1 (MQ=255)
                                    tGAAGGCGGCTTTGGTCACTCTGGTCGTCAGTAACACATACGCATCCGAATAACATCATAACATAGCCGCAAACATTTCGTTTGCGGTCATAg                                                     >  1:3027902/1‑93 (MQ=255)
                                     gAAGGCGGCTTTGGTCACTCTGGTCGTCAGTAACACATACGCATCCGAATAACATCATAACATAGCCGCAAACATTTCGTTTGCGGTCATAGc                                                    <  1:4469384/93‑1 (MQ=255)
                                       aGGCGGCTTTGGTCACTCTGGTCGTCAGTAACACATACGCATCCGAATAACATCATAACATAGCCGCAAACATTTCGTTTGCGGTCATAGCGt                                                  >  2:6781140/1‑93 (MQ=255)
                                        ggcggcTTTGGTCACTCTGGTCGTCAGTAACACATACGCATCCGAATAACATCATAACATAGCCGCAAACATTTCGTTTGCGGTCATAGCGTg                                                 >  1:3856419/1‑93 (MQ=255)
                                        ggcggcTTTGGTCACTCTGGTCGTCAGTAACACATACGCATCCGAATAACATCATAACATAGCCGCAAACATTTCGTTTGCGGTCATAGCGTg                                                 >  2:623352/1‑93 (MQ=255)
                                          cggcTTTGGTCACTCTGGTCGTCAGTAACACATACGCATCCGAATAACATCATAACATAGCCGCAAACATTTCGTTTGCGGTCATAGCGTggg                                               <  2:2222967/93‑1 (MQ=255)
                                              tttGGTCACTCTGGTCGTCAGTAACACATACGCATCCGAATAACATCATAACATAGCCGCAAACATTTCGTTTGCGGTCATAGCGTGGGTgcc                                           >  2:1297994/1‑93 (MQ=255)
                                              tttGGTCACTCTGGTCGTCAGTAACACATACGCATCCGAATAACATCATAACATAGCCGCAAACATTTCGTTTGCGGTCATAGCGTGGGTgcc                                           >  1:7617358/1‑93 (MQ=255)
                                                tGGTCACTCTGGTCGTCAGTAACACATACGCATCCGAATAACATCATAACATAGCCGCAAACATTTCGTTTGCGGTCATAGCGTGGGTgccgc                                         <  1:3793007/93‑1 (MQ=255)
                                                 ggTCACTCTGGTCGTCAGTAACACATACGCATCCGAATAACATCATAACATAGCCGCAAACATTTCGTTTGCGGTCATAGCGTGGGTgccgcc                                        <  2:5971212/93‑1 (MQ=255)
                                                  gTCACTCTGGTCGTCAGTAACACATACGCATCCGAATAACATCATAACATAGCCGCAAACATTTCGTTTGCGGTCATAGCGTGGGTGCCGCCt                                       <  1:6399068/93‑1 (MQ=255)
                                                  gTCACTCTGGTCGTCAGTAACACATACGCATCCGAATAACATCATAACATAGCCGCAAACATTTCGTTTGCGGTCATAGCGTGGGTGCCGCCt                                       >  2:2095161/1‑93 (MQ=255)
                                                   tCACTCTGGTCGTCAGTAACACATACGCATCCGAATAACATCATAACATAGCCGCAAACATTTCGTTTGCGGTCATAGCGTGGGTGCCGCCTg                                      <  2:5075226/93‑1 (MQ=255)
                                                    cgctctGGTCGTCAGTAACACATACGCATCCGAATAACATCATAACATAGCCGCAAACATTTCGTTTGCGGTCATAGCGTGGGTGCCGCCTgg                                     <  2:7686580/91‑1 (MQ=255)
                                                    cACTCTGGTCGTCAGTAACACATACGCATCCGAATAACATCATAACATAGCCGCAAACATTTCGTTTGCGGTCATAGCGTGGGTGCCGCCTgg                                     <  1:7145050/93‑1 (MQ=255)
                                                    cACTCTGGTCGTCAGTAACACATACGCATCCGAATAACATCATAACATAGCCGCAAACATTTCGTTTGCGGTCATAGCGTGGGTGCCGCCTgg                                     >  2:4846726/1‑93 (MQ=255)
                                                     aCTCTGGTCGTCAGTAACACATACGCATCCGAATAACATCATAACATAGCCGCAAACATTTCGTTTGCGGTCATAGCGTGGGTGCCGCCTGGc                                    >  1:6268732/1‑93 (MQ=255)
                                                       tctGGTCGTCAGTAACACATACGCATCCGAATAACATCATAACATAGCCGCAAACATTTCGTTTGCGGTCATAGCGTGGGTGCCGCCTGGCaa                                  <  1:7292574/93‑1 (MQ=255)
                                                       tctGGTCGTCAGTAACACATACGCATCCGAATAACATCATAACATAGCCGCAAACATTTCGTTTGCGGTCATAGCGTGGGTGCCGCCTGGCaa                                  >  2:105815/1‑93 (MQ=255)
                                                       tctGGTCGTCAGTAACACATACGCATCCGAATAACATCATAACATAGCCGCAAACATTTCGTTTGCGGTCATAGCGTGGGTGCCGCCTGGCaa                                  <  2:2825579/93‑1 (MQ=255)
                                                        ctGGTCGTCAGTAACACATACGCATCCGAATAACATCATAACATAGCCGCAAACATTTCGTTTGCGGTCATAGCGTGGGTGCCGCCTGGCAAg                                 <  2:4033257/93‑1 (MQ=255)
                                                         tGGTCGTCTGTAACACATACGCATCCGAATAACATCATAACATAGCCGCAAACATTTCGTTTGCGGTCATAGCGTGGGTGCCGCCTGGCAAGt                                <  2:425880/93‑1 (MQ=255)
                                                         tGGTCGTCAGTAACACATACGCATCCGAATAACATCATAACATAGCCGCAAACATTTCGTTTGCGGTCATAGCGTGGGTGCCGCCTGGCAAGt                                <  2:2141811/93‑1 (MQ=255)
                                                          ggTCGTCAGTAACACATACGCATCCGAATAACATCATAACATAGCCGCAAACATTTCGTTTGCGGTCATAGCGTGGGTGCCGCCTGGCAAGTg                               >  2:2945360/1‑93 (MQ=255)
                                                          ggTCGTCAGTAACACATACGCATCCGAATAACATCATAACATAGCCGCAAACATTTCGTTTGCGGTCATAGCGTGGGTGCCGCCTGGCAAGTg                               >  1:6754866/1‑93 (MQ=255)
                                                           gtcgtcAGTAACACATACGCATCCGAATAGCATCATAACATAGCCGCAAACATTTCGTTTGCGGTCATAGCGTGGGTGCCGCCTGGCAAGTGc                              <  1:2670904/93‑1 (MQ=255)
                                                           gtcgtcAGTAACACATACGCATCCGAATAACATCATAACATAGCCGCAAACATTTCGTTTGCGGTCATAGCGTGGGTGCCGCCTGGCAAGTGc                              >  1:3524642/1‑93 (MQ=255)
                                                           gtcgtcAGTAACACATACGCATCCGAATAACATCATAACATAGCCGCAAACATTTCGTTTGCGGTCATAGCGTGGGTGCCGCCTGGCAAGTGc                              <  2:3539872/93‑1 (MQ=255)
                                                            tcgtcAGTAACACATACGCATCCGAATAACATCATAACATAGCCGCAAACATTTCGTTTGCGGTCATAGCGTGGGTGCCGCCTGGCAAGTGCt                             >  1:1534795/1‑93 (MQ=255)
                                                             cgtcAGTAACACATACGCATCCGAATAACATCATAACATAGCCCAAAACATTTCGTTTGCGGTCATAGCGTGGGTGCCGCCTGGGAAGTGCtt                            >  1:6569641/1‑93 (MQ=255)
                                                               tcAGTAACACATACGCATCCGAATAACATCATAACATAGCCGCAAACATTTCGTTTGCGGTCATAGCGTGGGTGCCGCCTGGCAAGTGCTTAt                          >  2:7888626/1‑93 (MQ=255)
                                                               tcAGTAACACATACGCATCCGAATAACATCATAACATAGCCGCAAACATTTCGTTTGCGGTCATAGCGTGGGTGCCGCCTGGCAAGTGCTTAt                          >  1:7768149/1‑93 (MQ=255)
                                                                cAGTAACACATACGCATCCGAATAACATCATAACATAGCCGCAAACATTTCGTTTGCGGTCATAGCGTGGGTGCCGCCTGGCAAGTGCTTAtt                         >  1:184343/1‑93 (MQ=255)
                                                                cAGTAACACATACGCATCCGAATAACATCATAACATAGCCGCAAACATTTCGTTTGCGGTCATAGCGTGGGTGCCGCCTGGCAAGTGCTTAtt                         <  1:1933521/93‑1 (MQ=255)
                                                                 aGTAACACATACGCATCCGAATAACATCATAACATAGCCGCAAACATTTCGTTTGCGGTCATAGCGTGGGTGCCGCCTGGCAAGTGCTTAttt                        >  2:7470588/1‑93 (MQ=255)
                                                                 aGTAACACATACGCATCCGAATAACATCATAACATAGCCGCAAACATTTCGTTTGCGGTCATAGCGTGGGTGCCGCCTGGCAAGTGCTTAttt                        <  1:2066149/93‑1 (MQ=255)
                                                                 aGTAACACATACGCATCCGAATAACATCATAACATAGCCGCAAACATTTCGTTTGCGGTCATAGCGTGGGTGCCGCCTGGCAAGTGCTTAttt                        <  1:3366267/93‑1 (MQ=255)
                                                                    aaCACATACGCATCCGAATAACATCATAACATAGCCGCAAACATTTCGTTTGCGGTCATAGCGTGGGTGCCGCCTGGCAAGTGCTTATTTTCa                     <  1:5023236/93‑1 (MQ=255)
                                                                        caTACGCATCCGAATAACATCATAACATAGCCGCAAACATTTCGTTTGCGGTCATAGCGTGGGTGCCGCCTGGCAAGTGCTTATTTTCAgggg                 >  1:5948782/1‑93 (MQ=255)
                                                                         aTACGCATCCGAATAACATCATAACATAGCCGCAAACATTTCGTTTGCGGTCATAGCGTGGGTGCCGCCTGGCAAGTGCTTATTTTCAgggga                >  1:788457/1‑92 (MQ=255)
                                                                             gCATCCGAATAACATCATAACATAGCCGCAAACATTTCGTTTGCGGTCATAGCGTGGGTGCCGCCTGGCAAGTGCTTATTTTCAGGGGTAttt            <  2:1447703/93‑1 (MQ=255)
                                                                             gCATCCGAATAACATCATAACATAGCCGCAAACATTTCGTTTGCGGTCATAGCGTGGGTGCCGCCTGGCAAGTGCTTATTTTCAGGGGTAttt            <  1:696983/93‑1 (MQ=255)
                                                                               aTCCGAATAACATCATAACATAGCCGCAAACATTTCGTTTGCGGTCATAGCGTGGGTGCCGCCTGGCAAGTGCTTATTTTCAGGGGTATTTTg          >  1:4099309/1‑93 (MQ=255)
                                                                               aTCCGAATAACATCATAACATAGCCGCAAACATTTCGTTTGCGGTCATAGCGTGGGTGCCGCCTGGCAAGTGCTTATTTTCAGGGGTATTTTg          >  1:5358544/1‑93 (MQ=255)
                                                                               aTCCGAATAACATCATAACATAGCCGCAAACATTTCGTTTGCGGTCATAGCGTGGGTGCCGCCTGGCAAGTGCTTATTTTCAGGGGTATTTTg          >  2:2619074/1‑93 (MQ=255)
                                                                                tCCGAATAACATCATAACATAGCCGCAAACATTTCGTTTGCGGTCATAGCGTGGGTGCCGCCTGGCAAGTGCTTATTTTCAGGGGTATTTTGt         >  1:4047448/1‑93 (MQ=255)
                                                                                tCCGAATAACATCATAACATAGCCGCAAACATTTCGTTTGCGGTCATAGCGTGGGTGCCGCCTGGCAAGTGCTTATTTTCAGGGGTATTTTGt         >  2:7517855/1‑93 (MQ=255)
                                                                                tCCGAATAACATCATAACATAGCCGCAAACATTTCGTTTGCGGTCATAGCGTGGGTGCCGCCTGGCAAGTGCTTATTTTCAGGGGTATTTTGt         >  2:3241279/1‑93 (MQ=255)
                                                                                 ccGAATAACATCATAACATAGCCGCAAACATTTCGTTTGCGGTCATAGCGTGGGTGCCGCCTGGCAAGTGCTTATTTTCAGGGGTATTTTGTa        <  1:748516/93‑1 (MQ=255)
                                                                                   gAATAACATCATAACATAGCCGCAAACATTTCGTTTGCGGTCATAGCGTGGGTGCCGCCTGGCAAGTGCTTATTTTCAGGGGTATTTTGTAAc      >  1:1620595/1‑93 (MQ=255)
                                                                                    aaTAACATCATAACATAGCCGCAAACATTTCGTTTGCGGTCATAGCGTGGGTGCCGCCTGGCAAGTGCTTATTTTCAGGGGTATTTTGTAACa     <  1:4559330/93‑1 (MQ=255)
                                                                                    aaTAACATCATAACATAGCCGCAAACATTTCGTTTGCGGTCATAGCGTGGGTGCCGCCTGGCAAGTGCTTATTTTCAGGGGTATTTTGTAACa     <  2:3223411/93‑1 (MQ=255)
                                                                                    aaTAACATCATAACATAGCCGCAAACATTTCGTTTGCGGTCATAGCGTGGGTGCCGCCTGGCAAGTGCTTATTTTCAGGGGTATTTTGTAACa     >  1:7438339/1‑93 (MQ=255)
                                                                                    aaTAACATCATAACATAGCCGCAAACATTTCGTTTGCGGTCATAGCGTGGGTGCCGCCTGGCAAGTGCTTATTTTCAGGGGTATTTTGTAACa     <  2:3957870/93‑1 (MQ=255)
                                                                                     aTAACATCATAACATAGCCGCAAACATTTCGTTTGCGGTCATAGCGTGGGTGCCGCCTGGCAAGTGCTTATTTTCAGGGGTATTTTGTAACAt    >  1:5669213/1‑93 (MQ=255)
                                                                                      tAACATCATAACATAGCCGCAAACATTTCGTTTGCGGTCATAGCGTGGGTGCCGCCTGGCAAGTGCTTATTTTCAGGGGTATTTTGTAACATg   <  1:2531826/93‑1 (MQ=255)
                                                                                       aaCATCATAACATAGCCGCAAACATTTCGTTTGCGGTCATAGCGTGGGTGCCGCCTGGCAAGTGCTTATTTTCAGGGGTATTTTGTAACATgg  <  2:1686412/93‑1 (MQ=255)
                                                                                       aaCATCATAACATAGCCGCAAACATTTCGTTTGCGGTCATAGCGTGGGTGCCGCCTGGCAAGTGCTTATTTTCAGGGGTATTTTGTAACATgg  <  1:219160/93‑1 (MQ=255)
                                                                                          |                                                                                         
TAATCTGGTGGAAGATTTCGACGCCACCGACCGCGGTGAAGGCGGCTTTGGTCACTCTGGTCGTCAGTAACACATACGCATCCGAATAACGTCATAACATAGCCGCAAACATTTCGTTTGCGGTCATAGCGTGGGTGCCGCCTGGCAAGTGCTTATTTTCAGGGGTATTTTGTAACATGG  >  NC_000913/3811945‑3812124

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