Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 613,409 A→G A10A (GCA→GCG)  entF → ATP‑dependent serine activating enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913613,4090AG99.4% 526.8 / ‑4.8 169A10A (GCA→GCG) entFATP‑dependent serine activating enzyme
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (72/96);  total (73/96)
Fisher's exact test for biased strand distribution p-value = 4.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.13e-01

Only 100 of 176 total aligned reads displayed.

GGCTGGAAGCCCACTGGCGTACTCTGACACCGACGAATTTTACCCAGTTGCAGGAGGCACAATGAGCCAGCATTTACCTTTGGTCGCCGCACAGCCCGGCATCTGGATGGCAGAAAAACTGTCAGAATTACCCTCCGCCTGGAGCGTGGCGCATTACGTTGAGTTAACCGGAGAGGTTGA  >  NC_000913/613319‑613498
                                                                                          |                                                                                         
 gCTGGAAGCCCACTGGCGTACTCTGACACCGACGAATTTTACCCAGTTGCAGGAGGCACAATGAGCCAGCATTTACCTTTGGTCGCCGCGCAg                                                                                        <  1:261142/93‑1 (MQ=255)
   tGGAAGCCCACTGGCGTACTCTGACACCGACGAATTTTACCCAGTTGCAGGAGGCACAATGAGCCAGCATTTACCTTTGGTCGCCGCGCAGcc                                                                                      <  1:3331202/93‑1 (MQ=255)
     gAAGCCCACTGGCGTACTCTGACACCGACGAATTTTACCCAGTTGCAGGAGGCACAATGAGCCAGCATTTACCTTTGGTCGCCGCGCAGCCCg                                                                                    >  2:357350/1‑93 (MQ=255)
     gAAGCCCACGCGAGTACTCTGACACAGACGAATTGTACCCAGTTGCAGGAGGCACAATGAGCCAGCATTTACCTTTGGTCGCCGCGCAGCCCg                                                                                    <  2:1603715/93‑1 (MQ=255)
       aGCCCACTGGCGTACTCTGACACCGACGAATTTTACCCAGTTGCAGGAGGCACAATGAGCCAGCATTTACCTTTGGTCGCCGCGCAGCCCGGc                                                                                  <  1:1081818/93‑1 (MQ=255)
       aGCCCACTGGCGTACTCTGACACCGACGAATTTTACCCAGTTGCAGGAGGCACAATGAGCCAGCATTTACCTTTGGTCGCCGCGCAGCCCGGc                                                                                  <  1:1377562/93‑1 (MQ=255)
         cccACTGGCGTACTCTGACACCGACGAATTTTACCCAGTTGCAGGAGGCACAATGAGCCAGCATTTACCTTTGGTCGCCGCGCAGCCCGGCAt                                                                                >  2:4072224/1‑93 (MQ=255)
          ccACTGGCGTACTCTGACACCGACGAATTTTACCCAGTTGCAGGAGGCACAATGAGCCAGCATTTACCTTTGGTCGCCGCGCAGCCCGGCATc                                                                               >  1:3449728/1‑93 (MQ=255)
            aCTGGCGTACTCTGACACCGACGAATTTTACCCAGTTGCAGGAGGCACAATGAGCCAGCATTTACCTTTGGTCGCCGCGCAGCCCGGCATCTg                                                                             >  1:3355505/1‑93 (MQ=255)
            aCTGGCGTACTCTGACACCGACGAATTTTACCCAGTTGCAGGAGGCACAATGAGCCAGCATTTACCTTTGGTCGCCGCGCAGCCCGGCATCTg                                                                             <  1:4660861/93‑1 (MQ=255)
            aCTGGCGTACTCTGACACCGACGAATTTTACCCAGTTGCAGGAGGCACAATGAGCCAGCATTTACCTTTGGTCGCCGCGCAGCCCGGCATCTg                                                                             <  1:104743/93‑1 (MQ=255)
              tGGCGTACTCTGACACCGACGAATTTTACCCAGTTGCAGGAGGCACAATGAGCCAGCATTTACCTTTGGTCGCCGCGCAGCCCGGCATCTGGa                                                                           <  1:475045/93‑1 (MQ=255)
              tGGCGTACTCTGACACCGACGAATTTTACCCAGTTGCAGGAGGCACAATGAGCCAGCATTTACCTTTGGTCGCCGCGCAGCCCGGCATCTGGa                                                                           <  1:2717614/93‑1 (MQ=255)
              tGGCGTACTCTGACACCGACGAATTTTACCCAGTTGCAGGAGGCACAATGAGCCAGCATTTACCTTTGGTCGCCGCGCAGCCCGGCATCTGGa                                                                           >  2:3811692/1‑93 (MQ=255)
              tGGCGTACTCTGACACCGACGAATTTTACCCAGTTGCAGGAGGCACAATGAGCCAGCATTTACCTTTGGTCGCCGCGCAGCCCGGCATCTGGa                                                                           <  2:78202/93‑1 (MQ=255)
               ggCGTACTCTGACACCGACGAATTTTACCCAGTTGCAGGAGGCACAATGAGCCAGCATTTACCTTTGGTCGCCGCGCAGCCCGGCATCTGGAt                                                                          <  1:4392144/93‑1 (MQ=255)
                gCGTACTCTGACACCGACGAATTTTACCCAGTTGCAGGAGGCACAATGAGCCAGCATTTACCTTTGGTCGCCGCGCAGCCCGGCATCTGGATg                                                                         <  1:2298743/93‑1 (MQ=255)
                  gTACTCTGACACCGACGAATTTTACCCAGTTGCAGGAGGCACAATGAGCCAGCATTTACCTTTGGTCGCCGCGCAGCCCGGCATCTGGATGGc                                                                       <  1:835504/93‑1 (MQ=255)
                  gTACTCTGACACCGACGAATTTTACCCAGTTGCAGGAGGCACAATGAGCCAGCATTTACCTTTGGTCGCCGCGCAGCCCGGCATCTGGATGGc                                                                       >  1:1914413/1‑93 (MQ=255)
                  gTACTCTGACACCGACGAATTTTACCCAGTTGCAGGAGGCACAATGAGCCAGCATTTACCTTTGGTCGCCGCGCAGCCCGGCATCTGGATGGc                                                                       <  2:3304672/93‑1 (MQ=255)
                  gTACTCTGACACCGACGAATTTTACCCAGTTGCAGGAGGCACAATGAGCCAGCATTTACCTTTGGTCGCCGCGCAGCCCGGCATCTGGATGGc                                                                       >  1:392244/1‑93 (MQ=255)
                    aCTCTGACACCGACGAATTTTACCCAGTTGCAGGAGGCACAATGAGCCAGCATTTACCTTTGGTCGCCGCGCAGCCCGGCATCTGGATGGCAg                                                                     >  1:1599302/1‑93 (MQ=255)
                     ctctGACACCGACGAATTTTACCCAGTTGCAGGAGGCACAATGAGCCAGCATTTACCTTTGGTCGCCGCGCAGCCCGGCATCTGGATGGCAGa                                                                    >  2:762064/1‑93 (MQ=255)
                      tctGACACCGACGAATTTTACCCAGTTGCAGGAGGCACAATGAGCCAGCATTTACCTTTGGTCGCCGCGCAGCCCGGCATCTGGATGGCAGaa                                                                   >  1:4393640/1‑93 (MQ=255)
                        tGACACCGACGAATTTTACCCAGTTGCAGGAGGCACAATGAGCCAGCATTTACCTTTGGTCGCCGCGCAGCCCGGCATCTGGATGGCAGaaaa                                                                 >  1:3912825/1‑93 (MQ=255)
                        tGACACCGACGAATTTTACCCAGTTGCAGGAGGCACAATGAGCCAGCATTTACCTTTGGTCGCCGCGCAGCCCGGCATCTGGATGGCAGaaaa                                                                 >  1:992080/1‑93 (MQ=255)
                         gACACCGACGAATTTTACCCAGTTGCAGGAGGCACAATGAGCCAGCATTTACCTTTGGTCGCCGCGCAGCCCGGCATCTGGATGGCAGaaaaa                                                                <  1:4301600/93‑1 (MQ=255)
                          acacCGACGAATTTTACCCAGTTGCAGGAGGCACAATGAGCCAGCATTTACCTTTGGTCGCCGCGCAGCCCGGCATCTGGATGGCAGAAAAAc                                                               <  2:4351523/93‑1 (MQ=255)
                          acacCGACGAATTTTACCCAGTTGCAGGAGGCACAATGAGCCAGCATTTACCTTTGGTCGCCGCGCAGCCCGGCATCTGGATGGCAGAAAAAc                                                               <  2:214185/93‑1 (MQ=255)
                            acCGACGAATTTTACCCAGTTGCAGGAGGCACAATGAGCCAGCATTTACCTTTGGTCGCCGCGCAGCCCGGCATCTGGATGGCAGAAAAACTg                                                             <  2:2482569/93‑1 (MQ=255)
                             ccGACGAATTTTACCCAGTTGCAGGAGGCACAATGAGCCAGCATTTACCTTTGGTCGCCGCGCAGCCCGGCATCTGGATGGCAGAAAAACTGt                                                            <  1:845679/93‑1 (MQ=255)
                             ccGACGAATTTTACCCAGTTGCAGGAGGCACAATGAGCCAGCATTTACCTTTGGTCGCCGCGCAGCCCGGCATCTGGATGGCAGAAAAACTGt                                                            <  1:57497/93‑1 (MQ=255)
                             ccGACGAATTTTACCCAGTTGCAGGAGGCACAATGAGCCAGCATTTACCTTTGGTCGCCGCGCAGCCCGGCATCTGGATGGCAGAAAAACTGt                                                            <  2:3871619/93‑1 (MQ=255)
                              cgacgaATTTTACCCAGTTGCAGGAGGCACAATGAGCCAGCATTTACCTTTGGTCGCCGCGCAGCCCGGCATCTGGATGGCAGAAAAACTGTc                                                           <  2:386252/93‑1 (MQ=255)
                               gacgaATTTTACCCAGTTGCAGGAGGCACAATGAGCCAGCATTTACCTTTGGTCGCCGCGCAGCCCGGCATCTGGATGGCAGAAAAACTGTCa                                                          <  1:1226039/93‑1 (MQ=255)
                                  gaATTTTACCCAGTTGCAGGAGGCACAATGAGCCAGCATTTACCTTTGGTCGCCGCGCAGCCCGGCATCTGGATGGCAGAAAAACTGTCAGaa                                                       >  2:4194308/1‑93 (MQ=255)
                                  gaATTTTACCCAGTTGCAGGAGGCACAATGAGCCAGCATTTACCTTTGGTCGCCGCGCAGCCCGGCATCTGGATGGCAGAAAAACTGTCAGaa                                                       <  2:4144155/93‑1 (MQ=255)
                                   aaTTTTACCCAGTTGCAGGAGGCACAATGAGCCAGCATTTACCTTTGGTCGCCGCGCAGCCCGGCATCTGGATGGCAGAAAAACTGTCAGAAt                                                      <  2:3364451/93‑1 (MQ=255)
                                   aaTTTTACCCAGTTGCAGGAGGCACAATGAGCCAGCATTTACCTTTGGTCGCCGCGCAGCCCGGCATCTGGATGGCAGAAAAACTGTCAGAAt                                                      <  1:3613477/93‑1 (MQ=255)
                                    aTTTTACCCAGTTGCAGGAGGCACAATGAGCCAGCATTTACCTTTGGTCGCCGCGCAGCCCGGCATCTGGATGGCAGAAAAACTGTCAGAAtt                                                     <  1:2942308/93‑1 (MQ=255)
                                     ttttACCCAGTTGCAGGAGGCACAATGAGCCAGCATTTACCTTTGGTCGCCGCGCAGCCCGGCATCTGGATGGCAGAAAAACTGTCAGAATTa                                                    >  1:4028013/1‑93 (MQ=255)
                                     ttttACCCAGTTGCAGGAGGCACAATGAGCCAGCATTTACCTTTGGTCGCCGCGCAGCCCGGCATCTGGATGGCAGAAAAACTGTCAGAATTa                                                    >  1:1412533/1‑93 (MQ=255)
                                     ttttACCCAGTTGCAGGAGGCACAATGACCCAGCATTTACCTTTGGTCGCCGCGCAGCCCGGCATCTGGATGGCAGAAAAACTGTCAGAATTa                                                    >  2:2278804/1‑93 (MQ=255)
                                      tttACCCAGTTGCAGGAGGCACAATGAGCCAGCATTTACCTTTGGTCGCCGCGCAGCCCGGCATCTGGATGGCAGAAAAACTGTCAGAATTAc                                                   <  1:2804613/93‑1 (MQ=255)
                                        tACCCAGTTGCAGGAGGCACAATGAGCCAGCATTTACCTTTGGTCGCCGCGCAGCCCGGCATCTGGATGGCAGAAAAACTGTCAGAATTAccc                                                 <  2:2936479/93‑1 (MQ=255)
                                        tACCCAGTTGCAGGAGGCACAATGAGCCAGCATTTACCTTTGGTCGCCGCGCAGCCCGGCATCTGGATGGCAGAAAAACTGTCAGAATTAccc                                                 <  2:845800/93‑1 (MQ=255)
                                        tACCCAGTTGCAGGAGGCACAATGAGCCAGCATTTACCTTTGGTCGCCGCGCAGCCCGGCATCTGGATGGCAGAAAAACTGTCAGAATTAccc                                                 <  2:1216472/93‑1 (MQ=255)
                                        tACCCAGTTGCAGGAGGCACAATGAGCCAGCATTTACCTTTGGTCGCCGCGCAGCCCGGCATCTGGATGGCAGAAAAACTGTCAGAATTAccc                                                 <  2:3517742/93‑1 (MQ=255)
                                          cccAGTTGCAGGAGGCACAATGAGCCAGCATTTACCTTTGGTCGCCGCGCAGCCCGGCATCTGGATGGCAGAAAAACTGTCAGAATTACCCTc                                               <  1:1727515/93‑1 (MQ=255)
                                            cAGTTGCAGGAGGCACAATGAGCCAGCATTTACCTTTGGTCGCCGCGCAGCCCGGCATCTGGATGGCAGAAAAACTGTCAGAATTACCCTCCg                                             <  2:1219527/93‑1 (MQ=255)
                                            cAGTTGCAGGAGGCACAATGAGCCAGCATTTACCTTTGGTCGCCGCGCAGCCCGGCATCTGGATGGCAGAAAAACTGTCAGAATTACCCTCCg                                             <  2:1024166/93‑1 (MQ=255)
                                              gTTGCAGGAGGCACAATGAGCCAGCATTTACCTTTGGTCGCCGCGCAGCCCGGCATCTGGATGGCAGAAAAACTGTCAGAATTCCCCTCCGcc                                           <  1:248761/93‑1 (MQ=255)
                                              gTTGCAGGAGGCACAATGAGCCAGCATTTACCTTTGGTCGCCGCGCAGCCCGGCATCTGGATGGCAGAAAAACTGTCAGAATTACCCTCCGcc                                           <  2:1622003/93‑1 (MQ=255)
                                               ttGCAGGAGGCACAATGAGCCAGCATTTACCTTTGGTCGCCGCGCAGCCCGGCATCTGGATGGCAGAAAAACTGTCAGAATTACCCTCCGCCt                                          >  1:4378703/1‑93 (MQ=255)
                                                 gCAGGAGGCACAATGAGCCAGCATTTACCTTTGGTCGCCGCGCAGCCCGGCATCTGGATGGCAGAAAAACTGTCAGAATTACCCTCCGCCTgg                                        <  2:2599631/93‑1 (MQ=255)
                                                  cAGGAGGCACAATGAGCCAGCATTTACCTTTGGTCGCGGCGCAGCCCGGCATCTGGATGGCAGAAAAACTGTCAGAATTACCCTCCGCCTGGa                                       <  2:3635844/93‑1 (MQ=255)
                                                  cAGGAGGCACAATGAGCCAGCATTTACCTTTGGTCGCCGCGCAGCCCGGCATCTGGATGGCAGAAAAACTGTCAGAATTACCCTCCGCCTGGa                                       >  1:1239843/1‑93 (MQ=255)
                                                   aggaggCACAATGAGCCAGCATTTACCTTTGGTCGCCGCGCAGCCCGGCATCTGGATGGCAGAAAAACTGTCAGAATTACCCTCCGCCTGGAg                                      <  1:523559/93‑1 (MQ=255)
                                                    ggaggCACAATGAGCCAGCATTTACCTTTGGTCGCCGCGCAGCCCGGCATCTGGATGGCAGAAAAACTGTCAGAATTACCCTCCGCCTGGAGc                                     <  2:2746052/93‑1 (MQ=255)
                                                      aggCACAATGAGCCAGCATTTACCTTTGGTCGCCGCGCAGCCCGGCATCTGGATGGCAGAAAAACTGTCAGAATTACCCTCCGCCTGGAGCGt                                   <  2:3993162/93‑1 (MQ=255)
                                                       ggCACAATGAGCCAGCATTTACCTTTGGTCGCCGCGCAGCCCGGCATCTGGATGGCAGAAAAACTGTCAGAATTACCCTCCGCCTGGAGCGTg                                  >  2:3043283/1‑93 (MQ=255)
                                                       ggCACAATGAGCCAGCATTTACCTTTGGTCGCCGCGCAGCCCGGCATCTGGATGGCAGAAAAACTGTCAGAATTACCCTCCGCCTGGAGCGTg                                  <  2:4411335/93‑1 (MQ=255)
                                                       gccacaATGAGCCAGCATTTACCTTTGGTCGCCGCGCAGCCCGGCATCTGGATGGCAGAAAACCTGTCAGAATTACCCTCCGCCTGGAGCGTg                                  <  2:424545/91‑1 (MQ=255)
                                                         cacaATGAGCCAGCATTTACCTTTGGTCGCCGCGCAGCCCGGCATCTGGATGGCAGAAAAACTGTCAGAATTACCCTCCGCCTGGAGCGTGgc                                <  2:747484/93‑1 (MQ=255)
                                                         cacaATGAGCCAGCATTTACCTTTGGTCGCCGCGCAGCCCGGCATCTGGATGGCAGAAAAACTGTCAGAATTACCCTCCGCCTGGAGCGTGgc                                <  2:2973560/93‑1 (MQ=255)
                                                          acaATGAGCCAGCATTTACCTTTGGTCGCCGCGCAGCCCGGCATCTGGATGGCAGAAAAACTGTCAGAATTACCCTCCGCCTGGAGCGTGgcg                               >  1:3145735/1‑93 (MQ=255)
                                                           caATGAGCCAGCATTTACCTTTGGTCGCCGCGCAGCCCGGCATCTGGATGGCAGAAAAACTGTCAGAATTACCCTCCGCCTGGAGCGTGgcgc                              >  1:2210491/1‑93 (MQ=255)
                                                            aaTGAGCCAGCATTTACCTTTGGTCGCCGCGCAGCCCGGCATCTGGATGGCAGAAAAACTGTCAGAATTACCCTCCGCCTGGAGCGTGGCGCa                             >  2:2059426/1‑93 (MQ=255)
                                                            aaTGAGCCAGCATTTACCTTTGGTCGCCGCGCAGCCCGGCATCTGGATGGCAGAAAAACTGTCAGAATTACCCTCCGCCTGGAGCGTGGCGCa                             >  2:2171874/1‑93 (MQ=255)
                                                            aaTGAGCCAGCATTTACCTTTGGTCGCCGCGCAGCCCGGCATCTGGATGGCAGAAAAACTGTCAGAATTACCCTCCGCCTGGAGCGTGGCGCa                             >  2:1857772/1‑93 (MQ=255)
                                                             aTGAGCCAGCATTTACCTTTGGTCGCCGCGCAGCCCGGCATCTGGATGGCAGAAAAACTGTCAGAATTACCCTCCGCCTGGAGCGTGGCGCAt                            >  2:3408859/1‑93 (MQ=255)
                                                               gAGCCAGCATTTACCTTTGGTCGCCGCGCAGCCCGGCATCTGGATGGCAGAAAAACTGTCAGAATTACCCTCCGCCTGGAGCGTGGCGCATTa                          <  1:2078287/93‑1 (MQ=255)
                                                                  ccAGCATTTACCTTTGGTCGCCGCGCAGCCCGGCATCTGGATGGCAGAAAAACTGTCAGAATTACCCTCCGCCTGGAGCGTGGCGCATTACGt                       <  2:1108726/93‑1 (MQ=255)
                                                                     gCATTTACCTTTGGTCGCCGCGCAGCCCGGCATCTGGATGGCAGAAAAACTGTCAGAATTACCCTCCGCCTGGAGCGTGGCGCATTACGTTGa                    <  1:422911/93‑1 (MQ=255)
                                                                     gCATTTACCTTTGGTCGCCGCGCAGCCCGGCATCTGGATGGCAGAAAAACTGTCAGAATTACCCTCCGCCTGGAGCGTGGCGCATTACGTTGa                    <  2:4671008/93‑1 (MQ=255)
                                                                     gCATTTACCTTTGGTCGCCGCGCAGCCCGGCATCTGGATGGCAGAAAAACTGTCAGAATTACCCTCCGCCTGGAGCGTGGCGCATTACGTTGa                    >  1:3143202/1‑93 (MQ=255)
                                                                      cATTTACCTTTGGTCGCCGCGCAGCCCGGCATCTGGATGGCAGAAAAACTGTCAGAATTACCCTCCGCCTGGAGCGTGGCGTATTACGTTGAg                   >  2:2008272/1‑93 (MQ=255)
                                                                      cATTTACCTTTGGTCGCCGCGCAGCCCGGCATCTGGATGGCAGAAAAACTGTCAGAATTACCCTCCGCCTGGAGCGTGGCGCATTACGTTGAg                   >  2:4592485/1‑93 (MQ=255)
                                                                       aTTTACCTTTGGTCGCCGCGCAGCCCGGCATCTGGATGGCAGAAAAACTGTCAGAATTACCCTCCGCCTGGAGCGTGGCGCATTACGTTGAGt                  >  1:256212/1‑93 (MQ=255)
                                                                        tttACCTTTGGTCGCCGCGCAGCCCGGCATCTGGATGGCAGAAAAACTGTCAGAATTACCCTCCGCCTGGAGCGTGGCGCATTACGTTGAGtt                 <  2:2175582/93‑1 (MQ=255)
                                                                        tttACCTTTGGTCGCCGCGCAGCCCGGCATCTGGATGGCAGAAAAACTGTCAGAATTACCCACCGCCTGGAGCGTGGCGCATTACGTTGAGtt                 <  2:3756521/93‑1 (MQ=255)
                                                                          tACCTTTGGTCGCCGCGCAGCCCGGCATCTGGATGGCAGAAAAACTGTCAGAATTACCCTCCGCCTGGAGCGTGGCGCATTACGTTGAGTTaa               <  1:2468428/93‑1 (MQ=255)
                                                                            ccTTTGGTCGCCGCGCAGCCCGGCATCTTGATGGGAGAAAAACTGTCAGAATTACACTCCGCCTGGAGCGCGGCGCATTACGTTGAGTTAAcc             >  2:4437573/1‑93 (MQ=255)
                                                                            ccTTTGGTCGCCGCGCAGCCCGGCATCTGGATGGCAGAAAAACTGTCAGAATTACCCTCCGCCTGGAGCGTGGCGCATTACGTTGAGTTAAcc             <  1:1542315/93‑1 (MQ=255)
                                                                            ccTTTGGTCGCCGCACAGCCCGGCATCTGGATGGCAGAAAAACTGTCAGAATTACCCTCCGCCTGGAGCGTGGCGCATTACGTTGAGTTAAcc             >  1:3509534/1‑93 (MQ=255)
                                                                            ccTTTGGTAGCCGCGCAGACCGGAATCTGGATGGAAGAAAAACTGTCAGAATTACGCTCCGCCTGGAGCGTGGGGCGTTACGTTGAGTTAAcg             >  2:1449208/1‑92 (MQ=255)
                                                                             cTTTGGTCGCCGCGCAGCCCGGCATCTGGATGGCAGAAAAACTGTCAGAATTACCCTCCGCCTGGAGCGTGGCGCATTACGTTGAGTTAACCg            <  1:3816916/93‑1 (MQ=255)
                                                                             cTTTGGTCGCCGCGCAGCCCGGCATCTGGATGGCAGAAAAACTGTCAGAATTACCCTCCGCCTGGAGCGTGGCGCATTACGTTGAGTTAACCg            >  2:801313/1‑93 (MQ=255)
                                                                             cTTTGGTCGCCGCGCAGCCCGGCATCTGGATGGCAGAAAAACTGTCAGAATTACCCTCCGCCTGGAGCGTGGCGCATTACGTTGAGTTAACCg            >  1:4004471/1‑93 (MQ=255)
                                                                             cTTTGGTCGCCGCGCAGCCCGGCATCTGGATGGCAGAAAAACTGTCAGAATTACCCTCCGCCTGGAGCGTGGCGCATTACGTTGAGTTAACCg            >  2:2044461/1‑93 (MQ=255)
                                                                             cTTTGGTCGCCGCGCAGCCCGGCATCTGGATGGCAGAAAAACTGTCAGAATTACCCTCCGCCTGGAGCGTGGCGCATTACGTTGAGTTAACCg            >  1:637890/1‑93 (MQ=255)
                                                                              tttGGTCGCCGCGCAGCCCGGCATCTGGATGGCAGAAAAACTGTCAGAATTACCCTCCGCCTGGAGCGTGGCGCATTACGTTGAGTTAACCgg           >  2:1749322/1‑93 (MQ=255)
                                                                              tttGGTCGCCGCGCAGCCCGGCATCTGGATGGCAGAAAAACTGTCAGAATTACCCTCCGCCTGGAGCGTGGCGCATTACGTTGAGTTAACCgg           <  1:4060929/93‑1 (MQ=255)
                                                                               ttGGTCGCCGCGCAGCCCGGCATCTGGATGGCAGAAAAACTGTCAGAATTACCCTCCGCCTGGAGCGTGGCGCATTACGTTGAGTTAACCGga          <  2:1660593/93‑1 (MQ=255)
                                                                               ttGGTCGCCGCGCAGCCCGGCATCTGGATGGCAGAAAAACTGTCAGAATTACCCTCCGCCTGGAGCGTGGCGCATTACGTTGAGTTAACCGga          <  2:1848598/93‑1 (MQ=255)
                                                                                 ggTCGCCGCGCAGCCCGGCATCTGGATGGCAGAAAAACTGTCAGAATTACCCTCCGCCTGGAGCGTGGCGCATTACGTTGAGTTAACCGgaga        >  1:453430/1‑93 (MQ=255)
                                                                                    cgccgcGCAGCCCGGCATCTGGATGGCAGAAAAACTGTCAGAATTACCCTCCGCCTGGAGCGTGGCGCATTACGTTGAGTTAACCGGAGAGGt     <  2:2075945/93‑1 (MQ=255)
                                                                                    cgccgcGCAGCCCGGCATCTGGATGGCAGAAAAACTGTCAGAATTACCCTCCGCCTGGAGCGTGGCGCATTACGTTGAGTTAACCGGAGAGGt     <  2:277866/93‑1 (MQ=255)
                                                                                       cgcgCAGCCCGGCATCTGGATGGCAGAAAAACTGTCAGAATTACCCTCCGCCTGGAGCGTGGCGCATTACGTTGAGTTAACCGGAGAGGTTGa  >  2:1580350/1‑93 (MQ=255)
                                                                                       cgcgCAGCCCGGCATCTGGATGGCAGAAAAACTGTCAGAATTACCCTCCGCCTGGAGCGTGGCGCATTACGTTGAGTTAACCGGAGAGGTTGa  <  1:3192872/93‑1 (MQ=255)
                                                                                          |                                                                                         
GGCTGGAAGCCCACTGGCGTACTCTGACACCGACGAATTTTACCCAGTTGCAGGAGGCACAATGAGCCAGCATTTACCTTTGGTCGCCGCACAGCCCGGCATCTGGATGGCAGAAAAACTGTCAGAATTACCCTCCGCCTGGAGCGTGGCGCATTACGTTGAGTTAACCGGAGAGGTTGA  >  NC_000913/613319‑613498

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