Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,021,493 C→T A121A (GCC→GCT
A117A (GCC→GCT
tatD →
b4556 →
Mg‑dependent DNase
hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,021,4930CT98.3% 614.6 / ‑2.2 176A121A (GCC→GCT
A117A (GCC→GCT
tatD
b4556
Mg‑dependent DNase
hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (1/2);  new base T (79/94);  total (80/96)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.79e-01

Only 100 of 183 total aligned reads displayed.

GGCGATTGGTGAATGTGGTCTCGACTTTAACCGCAACTTTTCGACGCCGGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCCGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGCCCACGAGCGGTTTATGACATTGCTGGAGCCGTGGCTGGATAAACTGCC  >  NC_000913/4021403‑4021582
                                                                                          |                                                                                         
ggCGATTGGTGAATGTGGTCTCGACTTTAACCGCAACTTTTCGACGCCGGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATtgctgc                                                                                         >  1:2955455/1‑93 (MQ=255)
 gCGATTGGTGAATGTGGTCTCGACTTTAACCGCAACTTTTCGACGCCGGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCa                                                                                        <  2:1984942/93‑1 (MQ=255)
 gCGATTGGTGAATGTGGTCTCGACTTTAACCGCAACTTTTCGACGCCGGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCa                                                                                        <  2:4721086/93‑1 (MQ=255)
  cGATTGGTGAATGTGGTCTCGACTTTAACCGCAACTTTTCGACGCCGGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAg                                                                                       >  2:3346507/1‑93 (MQ=255)
    aTTGGTGAATGTGGTCTCGACTTTAACCGCAACTTTTCGACGCCGGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGAt                                                                                     <  1:1007080/93‑1 (MQ=255)
    aTTGGTGAATGTGGTCTCGACTTTAACCGCAACTTTTCGACACCGGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGAt                                                                                     <  2:3584650/93‑1 (MQ=255)
     ttGGTGAATGTGGTCTCGACTTTAACCGCAACTTTTCGACGCCGGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGAtt                                                                                    >  2:3273712/1‑93 (MQ=255)
     ttGGTGAATGTGGTCTCGACTTTAACCGCAACTTTTCGACGCCGGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGAtt                                                                                    >  2:213901/1‑93 (MQ=255)
     ttGGTGAATGTGGTCTCGACTTTAACCGCAACTTTTCGACGCCGGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGAtt                                                                                    <  2:1670790/93‑1 (MQ=255)
     ttGGTGAATGTGGTCTCGACTTTAACCGCAACTTTGCGACGCCGGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGAtt                                                                                    <  2:2723992/93‑1 (MQ=255)
        gTGAATGTGGTCTCGACTTTAACCGCAACTTTTCGACGCCGGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTaa                                                                                 >  1:1338034/1‑93 (MQ=255)
         tGAATGTGGTCTCGACTTTAACCGCAACTTTTCGACGCCGGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTaaa                                                                                >  2:4363797/1‑93 (MQ=255)
         tGAATGTGGTCTCGACTTTAACCGCAACTTTTCGACGCCGGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTaaa                                                                                >  2:798255/1‑93 (MQ=255)
           aaTGTGGTCTCGACTTTAACCGCAACTTTTCGACGCCGGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACa                                                                              <  1:577616/93‑1 (MQ=255)
            aTGTGGTCTCGACTTTAACCGCAACTTTTCGACGCCGGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACAt                                                                             <  2:1878872/93‑1 (MQ=255)
            aTGTGGTCTCGACTTTAACCGCAACTTTTCGACGCCGGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACAt                                                                             <  1:3477923/93‑1 (MQ=255)
            aTGTGGTCTCGACTTTAACCGCAACTTTTCGACGCCGGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACAt                                                                             >  1:1022086/1‑93 (MQ=255)
              gtgGTCTCGACTTTAACCGCAACTTTTCGACGCCGGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGc                                                                           <  2:4714546/93‑1 (MQ=255)
              gtgGTCTCGACTTTAACCGCAACTTTTCGACGCCGGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGc                                                                           <  1:412985/93‑1 (MQ=255)
               tgGTCTCGACTTTAACCGCAACTTTTCGACGCCGGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGcc                                                                          >  1:2888865/1‑93 (MQ=255)
               tgGTCTCGACTTTAACCGCAACTTTTCGACGCCGGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGcc                                                                          <  2:3944949/93‑1 (MQ=255)
                 gTCTCGACTTTAACCGCAACTTTTCGACGCCGGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCgg                                                                        <  2:399937/93‑1 (MQ=255)
                 gTCTCGACTTTAACCGCAACTTTTCGACGCCGGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCgg                                                                        <  2:1092647/93‑1 (MQ=255)
                    tcGACTTTAACCGCAACTTTTCGACGCCGGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACTTGCCGGTAt                                                                     <  2:3141942/93‑1 (MQ=255)
                    tcGACTTTAACCGCAACTTTTCGACGCCGGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGTAt                                                                     >  2:1722124/1‑93 (MQ=255)
                      gACTTTAACCGCAACTTTTCGACGCCGGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGTAttt                                                                   >  1:1365701/1‑93 (MQ=255)
                        cTTTAACCGCAACTTTTCGACGCCGGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTAt                                                                 <  1:2877346/93‑1 (MQ=255)
                          ttAACCGCAACTTTTCGACGCCGGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGc                                                               >  1:2034478/1‑93 (MQ=255)
                          ttAACCGCAACTTTTCGACGCCGGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGc                                                               <  1:636597/93‑1 (MQ=255)
                          ttAACCGCAACTTTTCGACGCCGGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCCGCAGATTTAAACATGCCGGTATTTATGc                                                               <  1:4500837/93‑1 (MQ=255)
                               cGCAACTTTTCGACGCCGGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGt                                                          <  1:1805046/93‑1 (MQ=255)
                                gCAACTTTTCGACGCCGGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTc                                                         >  1:1058042/1‑93 (MQ=255)
                                 cAACTTTTCGACGCCGGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTcg                                                        >  2:1586775/1‑93 (MQ=255)
                                  caCTTTTCGACGCCGGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTcgc                                                       <  2:1186020/92‑1 (MQ=255)
                                    cTTTTCGACGCCGGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGa                                                     <  1:4027637/93‑1 (MQ=255)
                                    cTTTTCGACGCCGGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGa                                                     >  2:58371/1‑93 (MQ=255)
                                    cTTTTCGACGCCGGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGa                                                     >  2:386482/1‑93 (MQ=255)
                                     ttttCGACGCCGGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGAt                                                    <  2:181206/93‑1 (MQ=255)
                                     ttttCGACGCCGGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGAt                                                    <  1:3438950/93‑1 (MQ=255)
                                      tttCGACGCCGGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATg                                                   <  2:3185614/93‑1 (MQ=255)
                                        tCGACGCCGGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGcc                                                 <  1:4463677/93‑1 (MQ=255)
                                         cGACGCCGGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGccc                                                <  2:1382912/93‑1 (MQ=255)
                                          gACGCCGGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGCCCa                                               >  1:224940/1‑93 (MQ=255)
                                          gACGCCGGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGCCCa                                               <  2:3595433/93‑1 (MQ=255)
                                          gACGCCGGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGCCCa                                               <  1:2511437/93‑1 (MQ=255)
                                           aCGCCGGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGCCCAc                                              <  2:3112760/93‑1 (MQ=255)
                                           aCGCCGGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGCCCAc                                              <  1:4541423/93‑1 (MQ=255)
                                           aCGCCGGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGCCCAc                                              <  1:920734/93‑1 (MQ=255)
                                            cGCCGGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGCCCACg                                             >  1:3352012/1‑93 (MQ=255)
                                             gCCGGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGCCCACGa                                            <  1:4306037/93‑1 (MQ=255)
                                             gCCGGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGCCCACGa                                            >  2:2476643/1‑93 (MQ=255)
                                             gCCAGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGCCCACGa                                            >  1:2529897/1‑93 (MQ=255)
                                               cGGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGCCCACGAGc                                          <  1:2076026/93‑1 (MQ=255)
                                                 gAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGCCCACGAGCgg                                        >  2:1142374/1‑93 (MQ=255)
                                                  aaGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGCCCACGAGCGGt                                       >  2:3627546/1‑93 (MQ=255)
                                                   agagCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGCCCACGAGCggct                                      >  2:4231515/1‑91 (MQ=255)
                                                   agagCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGCCCACGAGCGGtt                                      >  2:2680129/1‑93 (MQ=255)
                                                    gagCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGCCCACGAGCGGttt                                     >  2:3140996/1‑93 (MQ=255)
                                                    aaccAGGAACGCTCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGCCCACGAGCGGttt                                     <  2:2269586/90‑1 (MQ=255)
                                                      gCAGGAACGCGCTTTTGTTGCCCAGATACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGCCCACGAGCGGTTTAt                                   >  1:1452226/1‑93 (MQ=255)
                                                       cAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGCCCACGAGCGGTTTATg                                  <  2:4644080/93‑1 (MQ=255)
                                                        aGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGGTGCCCACGAGCGGTTTATGa                                 >  2:4277698/1‑93 (MQ=255)
                                                        aGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGCCCACGAGCGGTTTATGa                                 >  2:1258011/1‑93 (MQ=255)
                                                        aGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGCCCACGAGCGGTTTATGa                                 >  1:17434/1‑93 (MQ=255)
                                                        aGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGCCCACGAGCGGTTTATGa                                 <  1:326077/93‑1 (MQ=255)
                                                         ggAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGCCCACGAGCGGTTTATGAc                                >  1:1555962/1‑93 (MQ=255)
                                                         ggAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGCCCACGAGCGGTTTATGAc                                <  1:1179153/93‑1 (MQ=255)
                                                         ggAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGCCCACGAGCGGTTTATGAc                                >  1:512189/1‑93 (MQ=255)
                                                         ggAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGCCCACGAGCGGTTTATGAc                                >  2:374196/1‑93 (MQ=255)
                                                            aCGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGCCCACGAGCGGTTTATGACAtt                             <  2:3683575/93‑1 (MQ=255)
                                                             cgcgCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGCCCACGAGCGGTTTATGACATTg                            >  1:4247337/1‑93 (MQ=255)
                                                                 cTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGCCCACGAGCGGTTTATGACATTGCTgg                        >  1:3739279/1‑93 (MQ=255)
                                                                 cTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGCCCACGAGCGGTTTATGACATTGCTgg                        <  2:4253533/93‑1 (MQ=255)
                                                                 cTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGCCCACGAGCGGTTTATGACATTGCTgg                        <  1:2799872/93‑1 (MQ=255)
                                                                 cTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGCCCACGAGCGGTTTATGACATTGCTgg                        >  1:3959956/1‑93 (MQ=255)
                                                                  ttttGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGCCCACGAGCGGTTTATGACATTGCTGGa                       <  2:2952154/93‑1 (MQ=255)
                                                                  ttttGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGCCCACGAGCGGTTTATGACATTGCTGGa                       >  2:1698609/1‑93 (MQ=255)
                                                                  ttttGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGCCCACGAGCGGTTTATGACATTGCTGGa                       >  2:348814/1‑93 (MQ=255)
                                                                    ttgttgCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGCCCACGAGCGGTTTATGACATTGCTGGAGc                     <  2:3453817/93‑1 (MQ=255)
                                                                     tgttgCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGCCCACGAGCGGTTTATGACATTGCTGGAGcc                    <  1:3073817/93‑1 (MQ=255)
                                                                     tgttgCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGCCCACGAGCGGTTTATGACATTGCTGGAGcc                    <  2:2788387/93‑1 (MQ=255)
                                                                       ttgCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGCCCACGAGCGGTTTATGACATTGCTGGAGCCGt                  <  2:4547130/93‑1 (MQ=255)
                                                                       ttgCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGCCCACGAGCGGTTTATGACATTGCTGGAGCCGt                  <  1:1250250/93‑1 (MQ=255)
                                                                        tgCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGCCCACGAGCGGTTTATGACATTGCTGGAGCCGTg                 <  1:2647671/93‑1 (MQ=255)
                                                                        tgCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGCCCACGAGCGGTTTATGACATTGCTGGAGCCGTg                 <  2:2733630/93‑1 (MQ=255)
                                                                          cccAGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGCCCACGAGCGGTTTATGACATTGCTGGAGCCGTGGc               <  1:3113196/93‑1 (MQ=255)
                                                                          cccAGCTACGCATTGCCGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGCCCACGAGCGGTTTATGACATTGCTGGAGCCGTGGc               <  2:3312845/93‑1 (MQ=255)
                                                                            cAGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGCCCACGAGCGGTTTATGACATTGCTGGAGCCGTGGCTg             <  1:2154549/93‑1 (MQ=255)
                                                                             aGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGCCCACGAGCGGTTTATGACATTGCTGGAGCCGTGGCTgg            >  2:4130723/1‑93 (MQ=255)
                                                                             aGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGCCCACGAGCGGTTTATGACATTGCTGGAGCCGTGGCTgg            >  2:4471129/1‑93 (MQ=255)
                                                                             aGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGCCCACGAGCGGTTTATGACATTGCTGGAGCCGTGGCTgg            >  2:1438392/1‑93 (MQ=255)
                                                                              gCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGCCCACGAGCGGTTTATGACATTGCTGGAGCCGTGGCTGGa           <  2:4500008/93‑1 (MQ=255)
                                                                              gCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGCCCACGAGCGGTTTATGACATTGCTGGAGCCGTGGCTGGa           >  1:4617134/1‑93 (MQ=255)
                                                                              gCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGCCCACGAGCGGTTTATGACATTGCTGGAGCCGTGGCTGGa           >  1:2117533/1‑93 (MQ=255)
                                                                                 aCGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGCCCACGAGCGGTTTATGACATTGCTGGAGCCGTGGCTGGAtta        >  2:1955042/1‑91 (MQ=255)
                                                                                     aTTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGCCCACGAGCGGTTTATGACATTGCTGGAGCCGTGGCTGGATAAACTg    <  2:1014325/93‑1 (MQ=255)
                                                                                     aTTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGCCCACGAGCGGTTTATGACATTGCTGGAGCCGTGGCTGGATAAACTg    >  2:1321891/1‑93 (MQ=255)
                                                                                      ttGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGCCCACGAGCGGTTTATGACATTGCTGGAGCCGTGGCTGGATAAACTGc   >  2:1769969/1‑93 (MQ=255)
                                                                                       tgctgcAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGCCCACGAGCGGTTTATGACATTGCTGGAGCCGTGGCTGGATAAACTGcc  <  1:715343/93‑1 (MQ=255)
                                                                                       tgctgcAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGCCCACGAGCGGTTTATGACATTGCTGGAGCCGTGGCTGGATAAACTGcc  <  2:4635735/93‑1 (MQ=255)
                                                                                          |                                                                                         
GGCGATTGGTGAATGTGGTCTCGACTTTAACCGCAACTTTTCGACGCCGGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCCGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGCCCACGAGCGGTTTATGACATTGCTGGAGCCGTGGCTGGATAAACTGCC  >  NC_000913/4021403‑4021582

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