Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,599,294 T→C V25A (GTG→GCG)  b2481 → hydrogenase 4 Fe‑S subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,599,2940TC96.2% 609.9 / 4.6 187V25A (GTG→GCG) b2481hydrogenase 4 Fe‑S subunit
Reads supporting (aligned to +/- strand):  ref base T (4/3);  new base C (86/93);  total (91/96)
Fisher's exact test for biased strand distribution p-value = 7.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 198 total aligned reads displayed.

CATTTTCAAGGAGCCGACATGAACCGCTTTGTGGTGGCCGAACCACTGTGGTGTACAGGATGTAATACCTGTCTCGCTGCCTGTTCGGACGTGCATAAAACGCAAGGTTTACAGCAACACCCGCGCCTGGCCCTGGCGAAGACGTCAACAATCACTGCCCCTGTCGTGTGTCATCACTGT  >  NC_000913/2599203‑2599382
                                                                                           |                                                                                        
 aTTTTCAAGGAGCCGACATGAACCGCTTTGTGGTGGCCGAACCACTGTGGTGTACAGGATGTAATACCTGTCTCGCTGCCTGTTCGGACgcgc                                                                                        >  1:1190992/1‑93 (MQ=255)
 aTTTTCAAGGAGCCGACATGAACCGCTTTGTGGTGGCCGAACCACTGTGGTGTACAGGATGTAATACCTGTCTCGCTGCCTGTTCGGACgcgc                                                                                        >  1:1197800/1‑93 (MQ=255)
  ttttCAAGGAGCCGACATGAACCGCTTTGTGGTGGCCGAACCACTGTGGTGTACAGGATGTAATACCTGTCTCGCTGCCTGTTCGGACGCGCa                                                                                       >  1:5179090/1‑93 (MQ=255)
  ttttCAAGGAGCCGACATGAACCGCTTTGTGGTGGCCGAACCACTGTGGTGTACAGGATGTAATACCTGTCTCGCTGCCTGTTCGGACGCGCa                                                                                       >  2:1425416/1‑93 (MQ=255)
    ttCAAGGAGCCGACATGAACCGCTTTGTGGTGGCCGAACCACTGTGGTGTACAGGATGTAATACCTGTCTCGCTGCCTGTTCGGACGCGCATa                                                                                     <  1:3887430/93‑1 (MQ=255)
     tCAAGGAGCCGACATGAACCGCTTTGTGGTGGCCGAACCACTGTGGTGTACAGGATGTAATACCTGTCTCGCTGCCTGTTCGGACGTGCATaa                                                                                    >  1:3943215/1‑93 (MQ=255)
     tCAAGGAGCCGACATGAACCGCTTTGTGGTGGCCGAACCACTGTGGTGTACAGGATGTAATACCTGTCTCGCTGCCTGTTCGGACGTGCATaa                                                                                    >  2:4392904/1‑93 (MQ=255)
     tCAAGGAGCCGACATGAACCGCTTTGTGGTGGCCGAACCACTGTGGTGTACAGGATGTAATACCTGTCTCGCTGCCTGTTCGGACGCGCATaa                                                                                    >  1:3561421/1‑93 (MQ=255)
     tCAAGGAGCCGACATGAACCGCTTTGTGGTGGCCGAACCACTGTGGTGTACAGGATGTAATACCTGTCTCGCTGCCTGTTCGGACGCGCATaa                                                                                    >  1:1770916/1‑93 (MQ=255)
     tCAAGGAGCCGACATGAACCGCTTTGTGGTGGCCGAACCACTGTGGTGTACAGGATGTAATACCTGTCTCGCTGCCTGTTCGGACGCGCATaa                                                                                    >  1:1998724/1‑93 (MQ=255)
      cAAGGAGCCGACATGAACCGCTTTGTGGTGGCCGAACCACTGTGGTGTACAGGATGTAATACCTGTCTCGCTGCCTGTTCGGACGCGCATaaa                                                                                   <  1:1910773/93‑1 (MQ=255)
      cAAGGAGCCGACATGAACCGCTTTGTGGTGGCCGAACCACTGTGGTGTACAGGATGTAATACCTGTCTCGCTGCCTGTTCGGACGCGCATaaa                                                                                   <  1:5589989/93‑1 (MQ=255)
      cAAGGAGCAGACATGAACCGCTTTGTGGTGGCCGAACCACTGTGGTGTACAGGATGTAATACCTGTCTCGCTGCCTGTTCGGACGCGCATaaa                                                                                   <  2:4540691/93‑1 (MQ=255)
       aaGGAGCCGACATGAACCGCTTTGTGGTGGCCGAACCACTGTGGTGTACAGGATGTAATACCTGTCTCGCTGCCTGTTCGGACGCGCATaaaa                                                                                  >  2:968687/1‑93 (MQ=255)
       aaGGAGCCGACATGAACCGCTTTGTGGTGGCCGAACCACTGTGGTGTACAGGATGTAATACCTGTCTCGCTGCCTGTTCGGACGCGCATaaaa                                                                                  <  2:4753839/93‑1 (MQ=255)
        aGGAGCCGACATGAACCGCTTTGTGGTGGCCGAACCACTGTGGTGTACAGGATGTAATACCTGTCTCGCTGCCTGTTCGGACGCGCATAAAAc                                                                                 >  2:2033223/1‑93 (MQ=255)
        aGGAGCCGACATGAACCGCTTTGTGGTGGCCGAACCACTGTGGTGTACAGGATGTAATACCTGTCTCGCTGCCTGTTCGGACGCGCATAAAAc                                                                                 <  1:4759191/93‑1 (MQ=255)
        aGGAGCCGACATGAACCGCTTTGTGGTGGCCGAACCACTGTGGTGTACAGGATGTAATACCTGTCTCGCTGCCTGTTCGGACGCGCATAAAAc                                                                                 <  1:4184928/93‑1 (MQ=255)
         ggAGCCGACATGAACCGCTTTGTGGTGGCCGAACCACTGTGGTGTACAGGATGTAATACCTGTCTCGCGGCCTGTTCGGACGCGCATAAAAcc                                                                                >  1:4877977/1‑92 (MQ=255)
            gCCGACATGAACCGCTTTGTGGTGGCCGAACCACTGTGGTGTACAGGATGTAATACCTGTCTCGCTGCCTGTTCGGACGCGCATAAAACGCaa                                                                             >  1:654861/1‑93 (MQ=255)
            gCCGACATGAACCGCTTTGTGGTGGCCGAACCACTGTGGTGTACAGGATGTAATACCTGTCTCGCTGCCTGTTCGGACGCGCATAAAACGCaa                                                                             >  2:3689255/1‑93 (MQ=255)
            gCCGACATGAACCGCTTTGTGGTGGCCGAACCACTGTGGTGTACAGGATGTAATACCTGTCTCGCTGCCTGTTCGGACGCGCATAAAACGCaa                                                                             >  1:1677869/1‑93 (MQ=255)
            gCCGACATGAACCGCTTTGTGGTGGCCGAACCACTGTGGTGTACAGGATGTAATACCTGTCTCGCTGCCTGTTCGGACGCGCATAAAACGCaa                                                                             <  1:163107/93‑1 (MQ=255)
             ccGACATGAACCGCTTTGTGGTGGCCGAACCACTGTGGTGTACAGGATGTAATACCTGTCTCGCTGCCTGTTCGGACGCGCATAAAACGCAAg                                                                            <  1:652284/93‑1 (MQ=255)
               gACATGAACCGCTTTGTGGTGGCTGAACCACTGTGGTGTACAGGATGTAATACCTGTCTCGCTGCCTGTTCGGACGCGCATAAAACGCAAGGt                                                                          <  1:3499801/93‑1 (MQ=255)
               gACATGAACCGCTTTGTGGTGGCCGAACCACTGTGGTGTACAGGATGTAATACCTGTCTCGCTGCCTGTTCGGACGCGCATAAAACGCAAGGt                                                                          <  2:2548539/93‑1 (MQ=255)
               gACATGAACCGCTTTGTGGTGGCCGAACCACTGTGGTGTACAGGATGTAATACCTGTCTCGCTGCCTGTTCGGACGCGCATAAAACGCAAGGt                                                                          >  1:700720/1‑93 (MQ=255)
               gACATGAACCGCTTTGTGGTGGCCGAACCACTGTGGTGTACAGGATGTAATACCTGTCTCGCTGCCTGTTCGGACGCGCATAAAACGCAAGGt                                                                          >  2:4321234/1‑93 (MQ=255)
                aCATGAACCGCTTTGTGGTGGCCGAACCACTGTGGTGTACAGGATGTAATACCTGTCTCGCTGCCTGTTCGGACGCGCATAAAACGCAAGGtt                                                                         <  2:662807/93‑1 (MQ=255)
                aCATGAACCGCTTTGTGGTGGCCGAACCACTGTGGTGTACAGGATGTAATACCTGTCTCGCTGCCTGTTCGGACGCGCATAAAACGCAAGGtt                                                                         >  2:4498094/1‑93 (MQ=255)
                aCATGAACCGCTTTGTGGTGGCCGAACCACTGTGGTGTACAGGATGTAATACCTGTCTCGCTGCCTGTTCGGACGCGCATAAAACGCAAGGtt                                                                         <  2:4880582/93‑1 (MQ=255)
                 cATGAACCGCTTTGTGGTGGCCGAACCACTGTGGTGTACAGGATGTAATACCTGTCTCGCTGCCTGTTCGGACGCGCATAAAACGCAAGGttt                                                                        >  2:1162587/1‑93 (MQ=255)
                 cATGAACCGCTTTGTGGTGGCCGAACCACTGTGGTGTACAGGATGTAATACCTGTCTCGCTGCCTGTTCGGACGCGCATAAAACGCAAGGttt                                                                        >  1:1088516/1‑93 (MQ=255)
                 cATGAACCGCTTTGTGGTGGCCGAACCACTGTGGTGTACAGGATGTAATACCTGTCTCGCTGCCTGTTCGGACGCGCATAAAACGCAAGGttt                                                                        <  2:4666858/93‑1 (MQ=255)
                   tGAACCGCTTTGTGGTGGCCGAACCACTGTGGTGTACAGGATGTAATACCTGTCTCGCTGCCTGTTCGGACGCGCATAAAACGCAAGGTTTAc                                                                      >  2:1572704/1‑93 (MQ=255)
                   tGAACCGCTTTGTGGTGGCCGAACCACTGTGGTGTACAGGATGTAATACCTGTCTCGCTGCCTGTTCGGACGCGCATAAAACGCAAGGTTTAc                                                                      <  1:2368038/93‑1 (MQ=255)
                    gAACCGCTTTGTGGTGGCCGAACCACTGTGGTGTACAGGATGTAATACCTGTCTCGCTGCCTGTTCGGACGCGCATAAAACGCAAGGTTTACa                                                                     <  1:5310711/93‑1 (MQ=255)
                       ccGCTTTGTGGTGGCCGAACCACTGTGGTGTACAGGATGTAATACCTGTCTCGCTGCCTGTTCGGACGCGCATAAAACGCAAGGTTTACAGCa                                                                  <  2:3232930/93‑1 (MQ=255)
                       ccGCTTTGTGGTGGCCGAACCACTGTGGTGTACAGGATGTAATACCTGTCTCGCTGCCTGTTCGGACGCGCATAAAACGCAAGGTTTACAGCa                                                                  >  2:4594248/1‑93 (MQ=255)
                        cGCTTTGTGGTGGCCGAACCACTGTGGTGTACAGGATGTAATACCTGTCTCGCTGCCTGTTCGGACGCGCATAAAACGCAAGGTTTACAGCaa                                                                 <  1:3258896/93‑1 (MQ=255)
                        cGCTTTGTGGTGGCCGAACCACTGTGGTGTACAGGATGTAATACCTGTCTCGCTGCCTGTTCGGACGCGCATAAAACGCAAGGTTTACAGCaa                                                                 >  2:4030512/1‑93 (MQ=255)
                        cGCTTTGTGGTGGCCGAACCACTGTGGTGTACAGGATGTAATACCTGTCTCGCTGCCTGTTCGGACGCGCATAAAACGCAAGGTTTACAGCaa                                                                 <  2:4618091/93‑1 (MQ=255)
                         gCTTTGTGGTGGCCGAACCACTGTGGTGTACAGGATGTAATACCTGTCTCGCTGCCTGTTCGGACGCGCATAAAACGCAAGGTTTACAGCAac                                                                >  2:622699/1‑93 (MQ=255)
                          cTTTGTGGTGGCCGAACCACTGTGGTGTACAGGATGTAATACCTGTCTCGCTGCCTGTTCGGACGCGCATAAAACGCAAGGTTTACAGCAaca                                                               <  1:854767/93‑1 (MQ=255)
                           tttGTGGTGGCCGAACCACTGTGGTGTACAGGATGTAATACCTGTCTCGCTGCCTGTTCGGACGCGCATAAAACGCAAGGTTTACAGCAacac                                                              <  2:4786249/93‑1 (MQ=255)
                             tgtgGTGGCCGAACCACTGTGGTGTACAGGATGTAATACCTGTCTCGCTGCCTGTTCGGACGCGCATAAAACGCAAGGTTTACAGCAACAccc                                                            <  1:4142173/93‑1 (MQ=255)
                               tggtggCCGAACCACTGTGGTGTACAGGATGTAATACCTGTCTCGCTGCCTGTTCGGACGCGCATAAAACGCAAGGTTTACAGCAACACCcgc                                                          <  1:5631589/93‑1 (MQ=255)
                               tggtggCCGAACCACTGTGGTGTACAGGATGTAATACCTGTCTCGCTGCCTGTTCGGACGCGCATAAAACGCAAGGTTTACAGCAACACCcgc                                                          <  2:357420/93‑1 (MQ=255)
                                ggtggCCGAACCACTGTGGTGTACAGGATGTAATACCTGTCTCGCTGCCTGTTCGGACGCGCATAAAACGCAAGGTTTACAGCAACACCcgcg                                                         >  1:3195023/1‑93 (MQ=255)
                                   ggCCGAACCACTGTGGTGTACAGGATGTATTACCTGTCTCGCTGCCTGTTCGGACGCGCATAAAACGCAAGGTTTACAGCAACACCCGCGcct                                                      <  2:2962998/93‑1 (MQ=255)
                                    gCCGAACCACTGTGGTGTACAGGATGTAATACCTGTCTCGCTGCCTGTTCGGACGCGCATAAAACGCAAGGTTTACAGCAACACCCGCGcctg                                                     >  2:2295672/1‑93 (MQ=255)
                                    gCCGAACCACTGTGGTGTACAGGATGTAATACCTGTCTCGCTGCCTGTTCGGACGCGCATAAAACGCAAGGTTTACAGCAACACCCGCGcctg                                                     <  2:1986858/93‑1 (MQ=255)
                                     ccGAACCACTGTGGTGTACAGGATGTAATACCTGTCTCGCTGCCTGTTCGGACGCGCATAAAACGCAAGGTTTACAGCAACACCCGCGcctgg                                                    <  1:3277033/93‑1 (MQ=255)
                                      cGAACCACTGTGGTGTACAGGATGTAATACCTGTCTCGCTGCCTGTTCGGACGCGCATAAAACGCAAGGTTTACAGCAACACCCGCGcctggc                                                   >  2:3517812/1‑93 (MQ=255)
                                        aaCCACTGTGGTGTACAGGATGTAATACCTGTCTCGCTGCCTGTTCGGACGCGCATAAAACGCAAGGTTTACAGCAACACCCGCGcctggccc                                                 >  1:4540796/1‑93 (MQ=255)
                                        aaCCACTGTGGTGTACAGGATGTAATACCTGTCTCGCTGCCTGTTCGGACGCGCATAAAACGCAAGGTTTACAGCAACACCCGCGcctggccc                                                 >  2:28699/1‑93 (MQ=255)
                                         aCCACTGTGGTGTACAGGATGTAATACCTGTCTCGCTGCCTGTTCGGACGCGCATAAAACGCAAGGTTTACAGCAACACCCGCGcctggccct                                                <  1:794386/93‑1 (MQ=255)
                                         aCCACTGTGGTGTACAGGATGTAATACCTGTCTCGCTGCCTGTTCGGACGCGCATAAAACGCAAGGTTTACAGCAACACCCGCGcctggccct                                                <  1:3102157/93‑1 (MQ=255)
                                          ccACTGTGGTGTACAGGATGTAATACCTGTCTCGCTGCCTGTTCGGACGCGCATAAAACGCAAGGTTTACAGCAACACCCGCGcctggccctg                                               <  1:2976687/93‑1 (MQ=255)
                                          ccACTGTGGTGTACAGGATGTAATACCTGTCTCGCTGCCTGTTCGGACGCGCATAAAACGCAAGGTTTACAGCAACACCCGCGcctggccctg                                               <  1:3401505/93‑1 (MQ=255)
                                          ccACTGTGGTGTACAGGATGTAATACCTGTCTCGCTGCCTGTTCGGACGCGCATAAAACGCAAGGTTTACAGCAACACCCGCGcctggccctg                                               <  1:4980615/93‑1 (MQ=255)
                                          ccACTGTGGTGTACAGGATGTAATACCTGTCTCGCTGCCTGTTCGGACGCGCATAAAACGCAAGGTTTACAGCAACACCCGCGcctggccctg                                               <  2:5739314/93‑1 (MQ=255)
                                           cACTGTGGTGTACAGGATGTAATTTCTGTCTCGCTGCCTGTTCGGACGCGCATAAAACGCAAGGTTTACAGCAACCCCCGCGcctggccctgg                                              <  2:5254111/93‑1 (MQ=255)
                                           cACTGTGGTGTACAGGATGTAATGCCTGTCTCGCTGCCTGTTCGGACGCGCATAAAACGCAAGGTTTACAGCAACACCCGCGcctggccctgg                                              <  1:3332870/93‑1 (MQ=255)
                                            aCTGTGGTGTACAGGATGTAATACCTGTCTCGCTGCCTGTTCGGACGCGCATAAAACGCAAGGTTTACAGCAACACCCGCGcctggccctggc                                             >  1:2460538/1‑93 (MQ=255)
                                            aCTGTGGTGTACAGGATGTAATACCTGTCTCGCTGCCTGTTCGGACGCGCATAAAACGCAAGGTTTACAGCAACACCCGCGcctggccctggc                                             >  1:4340258/1‑93 (MQ=255)
                                                 ggtgTACAGGATTTAATACCTGTCTCGCTGCCTGTTCGGACGCGCATAAAACGCAAGGTTTACAGCAACACCCGCGCCTGGCCCTGGCGAAGa                                        <  2:3121356/93‑1 (MQ=255)
                                                  gtgtACAGGATGTAATACCTGTCTCGCTGCCTGTTCGGACGTGCATAAAACGCAAGGTTTACAGCAACACCCGCGCCTGGCCCTGGCGAAGAc                                       <  2:674554/93‑1 (MQ=255)
                                                     tACAGGATGTAATACCTGTCTCGCTGCCTGTTCGGACGCGCATAAAACGCAAGGTTTACAGCAACACCCGCGCCTGGCCCTGGCGAAGACGta                                    >  1:1176053/1‑92 (MQ=255)
                                                       cAGGATGTAATACCTGTCTCGCTGCCTGTTCGGACGCGCATAAAACGCAAGGTTTACAGCAACACCCGCGCCTGGCCCTGGCGAAGACGTcaa                                  <  2:587284/93‑1 (MQ=255)
                                                       cAGGATGTAATACCTGTCTCGCTGCCTGTTCGGACGCGCATAAAACGCAAGGTTTACAGCAACACCCGCGCCTGGCCCTGGCGAAGACGTcaa                                  <  1:2089191/93‑1 (MQ=255)
                                                         ggATGTAATACCTGTCTCGCTGCCTGTTCGGACGCGCATAAAACGCAAGGTTTACAGCAACACCCGCGCCTGGCCCTGGCGAAGACGTcaaca                                <  1:5080203/93‑1 (MQ=255)
                                                               aaTACCTGCCTCGCTGCCTGTTCGGACGCGCATAAAACGCAAGGTTTACAGCACCACCCGCGCCTGGCCCTGGCGAAGACGTCAACAATCACt                          <  1:5109514/93‑1 (MQ=255)
                                                                    cTGTCTCGCTGCCTGTTCGGACGCGCATAAAACGCAAGGTTTACAGCAACACCCGCGCCTGGCCCTGGCGAAGACGTCAACAATCACTGcccc                     <  2:846444/93‑1 (MQ=255)
                                                                     tGTCTCGCTGCCTGTTCGGACGCGCATAAAACGCAAGGTTTACAGCAACACCCGCGCCTGGCCCTGGCGAAGACGTCAACAATCACTGCCCCt                    >  1:2171100/1‑93 (MQ=255)
                                                                       tctcGCTGCCTGTTCGGACGCGCATAAAACGCAAGGTTTACAGCAACACCCGCGCCTGGCCCTGGCGAAGACGTCAACAATCACTGCCCCtat                  >  2:2830775/1‑91 (MQ=255)
                                                                       tctcGCTGCCTGTTCGGACGCGCATAAAACGCAAGGTTTACAGCAACACCCGCGCCTGGCCCTGGCGAAGACGTCAACAATCACTGCCCCTGt                  <  2:5518585/93‑1 (MQ=255)
                                                                       tcgcGCTGCCTGTTCGGTCGCGCATAAAACGCAAGGTTTACAGCAACACCCGCGCCTGGCCCTGGCGAAGACGTCAACAATCACTGCCCCTGt                  >  1:5483354/4‑93 (MQ=255)
                                                                        ctcGCTGCCTGTTCGGACGCGCATAAAACGCAAGGTTTACAGCAACACCCGCGCCTGGCCCTGGCGAAGACGTCAACAATCACTGCCCCTGTc                 <  1:3732071/93‑1 (MQ=255)
                                                                         tcGCTGCCTGTTCGGACGCGCATAAAACGCAAGGTTTACAGCAACACCCGCGCCTGGCCCTGGCGAAGACGTCAACAATCACTGCCCCTGTCg                >  2:1655867/1‑93 (MQ=255)
                                                                           gCTGCCTGTTCGGACGCGCATAAAACGCAAGGTTTACAGCAACACCCGCGCCTGGCCCTGGCGAAGACGTCAACAATCACTGCCCCTGTCgtg              <  2:5431556/93‑1 (MQ=255)
                                                                            cTGCCTGTTCGGACGCGCATAAAACGCAAGGTTTACAGCAACACCCGCGCCTGGCCCTGGCGAAGACGTCAACAATCACTGCCCCTGTCgtgt             >  2:3831176/1‑93 (MQ=255)
                                                                            cTGCCTGTTCGGACGCGCATAAAACGCAAGGTTTACAGCAACACCCGCGCCTGGCCCTGGCGAAGACGTCAACAATCACTGCCCCTGTCgtgt             <  2:5060334/93‑1 (MQ=255)
                                                                             tGCCTGTTCGGACGCGCATAAAACGCAAGGTTTACAGCAACACCCGCGCCTGGCCCTGGCGAAGACGTCAACAATCACTGCCCCTGTCgtgtg            >  1:385849/1‑93 (MQ=255)
                                                                             tGCCTGTTCGGACGCGCATAAAACGCAAGGTTTACAGCAACACCCGCGCCTGGCCCTGGCGAAGACGTCAACAATCACTGCCCCTGTCgtgtg            <  2:392012/93‑1 (MQ=255)
                                                                              gCCTGTTCGGACGCGCATAAAACGCAAGGTTTACAGCAACACCCGCGCCTGGCCCTGGCGAAGACGTCAACAATCACTGCCCCTGTCgtgtgt           >  1:5766164/1‑93 (MQ=255)
                                                                               ccTGTTCGGACGCGCATAAAACGCAAGGTTTACAGCAACACCCGCGCCTGGCCCTGGCGAAGACGTCAACAATCACTGCCCCTGTCGTGTGtc          >  1:4089111/1‑93 (MQ=255)
                                                                                cTGTTCGGACGCGCATAAAACGCAAGGTTTACAGCAACACCCGCGCCTGGGCCTGGCGAAGACGTCAACAATCACTGCCCCTGTCGTGTGtca         >  1:5727556/1‑93 (MQ=255)
                                                                                cTGTTCGGACGCGCATAAAACGCAAGGTTTACAGCAACACCCGCGCCTGGCCCTGGCGAAGACGTCAACAATCACTGCCCCTGTCGTGTGtca         >  1:3025300/1‑93 (MQ=255)
                                                                                cTGTTCGGACGCGCATAAAACGCAAGGTTTACAGCAACACCCGCGCCTGGCCCTGGCGAAGACGTCAACAATCACTGCCCCTGTCGTGTGtca         <  2:4809873/93‑1 (MQ=255)
                                                                                 tGTTCGGACGCGCATAAAACGCAAGGTTTACAGCAACACCCGCGCCTGGCCCTGGCGAAGACGTCAACAATCACTGCCCCTGTCGTGTGtcat        >  2:1378328/1‑93 (MQ=255)
                                                                                 tGTTCGGACGCGCATAAAACGCAAGGTTTACAGCAACACCCGCGCCTGGCCCTGGCGAAGACGTCAACAATCACTGCCCCTGTCGTGTGtcat        >  2:5008353/1‑93 (MQ=255)
                                                                                    tCGGACGCGCATAAAACGCAAGGTTTACAGCAACACCCGCGCCTGGCCCTGGCGAAGACGTCAACAATCACTGCCCCTGTCGTGTGTCATCAc     >  1:141408/1‑93 (MQ=255)
                                                                                    tCGGACGCGCATAAAACGCAAGGTTTACAGCAACACCCGCGCCTGGCCCTGGCGAAGACGTCAACAATCACTGCCCCTGTCGTGTGTCATCAc     <  2:4205462/93‑1 (MQ=255)
                                                                                     cGGACGCGCATAAAACGCAAGGTTTACAGCAACACCCGCGCCTGGCCCTGGCGAAGACGTCAACAATCACTGCCCCTGTCGTGTGTCATCACt    <  2:3993050/93‑1 (MQ=255)
                                                                                     cGGACGCGCATAAAACGCAAGGTTTACAGCAACACCCGCGCCTGGCCCTGGCGAAGACGTCAACAATCACTGCCCCTGTCGTGTGTCATCACt    <  1:3343556/93‑1 (MQ=255)
                                                                                      ggACGCGCATAAAACGCAAGGTTTACAGCAACACCCGCGCCTGGCCCTGGCGAAGACGTCAACAATCACTGCCCCTGTCGTGTGTCATCACtg   >  1:6048058/1‑93 (MQ=255)
                                                                                       gACGCGCATAAAACGCAAGGTTTACAGCAACCCCCGCGCCTGGCCCTGGCGAAGACGTCAACAATCACTGCCCCTGTCGTGTGTCATCACtgt  <  1:4459400/93‑1 (MQ=255)
                                                                                       gACGCGCATAAAACGCAAGGTTTACAGCAACACCCGCGCCTGGCCCTGGCGAAGACGTCAACAATCACTGCCCCTGTCGTGTGTCATCACtgt  <  1:4518141/93‑1 (MQ=255)
                                                                                       gACGCGCATAAAACGCAAGGTTTACAGCAACACCCGCGCCTGGCCCTGGCGAAGACGTCAACAATCACTGCCCCTGTCGTGTGTCATCACtgt  <  1:3743732/93‑1 (MQ=255)
                                                                                           |                                                                                        
CATTTTCAAGGAGCCGACATGAACCGCTTTGTGGTGGCCGAACCACTGTGGTGTACAGGATGTAATACCTGTCTCGCTGCCTGTTCGGACGTGCATAAAACGCAAGGTTTACAGCAACACCCGCGCCTGGCCCTGGCGAAGACGTCAACAATCACTGCCCCTGTCGTGTGTCATCACTGT  >  NC_000913/2599203‑2599382

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