Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,862,488 G→A A471A (GCC→GCT yidE ← putative transport protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,862,4880GA97.8% 487.6 / ‑1.9 139A471A (GCC→GCTyidEputative transport protein
Reads supporting (aligned to +/- strand):  ref base G (2/1);  new base A (70/66);  total (72/67)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.84e-01

Only 100 of 142 total aligned reads displayed.

TGCCAGCATCCCGCACATGGTCAGGTAATTCATTTTGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGGGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAATCCCCACCAGATTTCAGACCAACGACCGAGAGGAAC  >  NC_000913/3862398‑3862578
                                                                                          |                                                                                          
tGCCAGCATCCCGCACATGGTCAGGTAATTCATTTTGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGc                                                                                          >  1:1242694/1‑93 (MQ=255)
 gCCAGCATCCCGCACATGGTCAGGTAATTCATTTTGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCa                                                                                         >  1:902865/1‑93 (MQ=255)
    aGCATCCCGCACATGGTCAGGTAATTCATTTTGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCa                                                                                      <  1:1363101/93‑1 (MQ=255)
     gCATCCCGCACATGGTCAGGTAATTCATTTTGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCAt                                                                                     <  2:4032230/93‑1 (MQ=255)
       aTCCCGCACATGGTCAGGTAATTCATTTTGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATaa                                                                                   <  2:3626769/93‑1 (MQ=255)
        tCCCGCACATGGTCAGGTAATTCATTTTGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAAc                                                                                  >  1:2177071/1‑93 (MQ=255)
          ccGCACAAGGTCAGGTAATTCATTTTGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCa                                                                                >  1:1874780/1‑93 (MQ=255)
           cGCACATGGTCAGGTAATTCATTTTGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCaa                                                                               <  2:1621488/93‑1 (MQ=255)
            gCACATGGTCAGGTAATTCATTTTGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAAt                                                                              >  1:501039/1‑93 (MQ=255)
               caTGGTCAGGTAATTCATTTTGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCa                                                                           <  2:3370919/93‑1 (MQ=255)
                aTGGTCAGGTAATTCATTTTGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAg                                                                          >  2:947701/1‑93 (MQ=255)
                aTGGTCAGGTAATTCATTTTGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAg                                                                          <  2:3888696/93‑1 (MQ=255)
                  ggTCAGGTAATTCATTTTGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCt                                                                        >  2:1981974/1‑93 (MQ=255)
                  ggTCAGGTAATTCATTTTGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCt                                                                        <  1:3737618/93‑1 (MQ=255)
                   gTCAGGTAATTCATTTTGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCtt                                                                       <  1:814660/93‑1 (MQ=255)
                   gTCAGGTAATTCATTTTGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCtt                                                                       <  2:2167693/93‑1 (MQ=255)
                    tCAGGTAATTCATTTTGGCTAACATCCGCGCCAGAATGCAAACAGTAATTAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTa                                                                      <  1:3059664/93‑1 (MQ=255)
                    gcAGGTAATTCATTTTGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTa                                                                      <  2:1739177/92‑1 (MQ=255)
                       ggTAATTCATTTTGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGcc                                                                   <  2:3793829/93‑1 (MQ=255)
                       ggTAATTCATTTTGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGcc                                                                   >  1:2080076/1‑93 (MQ=255)
                       ggTAATTCATTTTGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGcc                                                                   >  1:3325177/1‑93 (MQ=255)
                        gTAATTCATTTTGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGccc                                                                  <  2:486052/93‑1 (MQ=255)
                        gTAATTCATTTTGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGccc                                                                  >  2:2961063/1‑93 (MQ=255)
                          aaTTCATTTTGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCtt                                                                <  2:421538/93‑1 (MQ=255)
                           aTTCATTTTGGCTAACATCCGCGCCAGAATGCCAGCAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTc                                                               <  1:2905412/93‑1 (MQ=255)
                           aTTCATTTTGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTc                                                               <  1:3972655/93‑1 (MQ=255)
                             tCATTTTGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGc                                                             <  1:2971333/93‑1 (MQ=255)
                               aTTTTGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGcct                                                           >  2:1605669/1‑92 (MQ=255)
                               aTTTTGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCa                                                           >  2:3690555/1‑93 (MQ=255)
                                ttttGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCAt                                                          <  1:4636022/93‑1 (MQ=255)
                                ttttGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCAt                                                          >  2:2285038/1‑93 (MQ=255)
                                ttttGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCAt                                                          >  2:1631500/1‑93 (MQ=255)
                                ttttGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCAt                                                          <  2:1041174/93‑1 (MQ=255)
                                ttttGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCAt                                                          <  1:2594882/93‑1 (MQ=255)
                                ttttGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCAt                                                          >  1:2325030/1‑93 (MQ=255)
                                ttttGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCAt                                                          >  2:4476798/1‑93 (MQ=255)
                                 tttGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCAtt                                                         <  1:4420164/93‑1 (MQ=255)
                                  ttGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGGGGTGAGCAGAGCGCCATAACCAAGCCAGCTTAGCCCTTCGCCATTg                                                        <  2:402517/93‑1 (MQ=255)
                                  ttGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTg                                                        <  2:402518/93‑1 (MQ=255)
                                   tGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGa                                                       >  1:3475302/1‑93 (MQ=255)
                                    ggCTAACATCCGCGCCAGAGTGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGAc                                                      <  2:2658097/93‑1 (MQ=255)
                                    ggCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGAc                                                      <  1:1795471/93‑1 (MQ=255)
                                    ggCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGAc                                                      >  1:1871958/1‑93 (MQ=255)
                                     gCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGAcc                                                     >  1:4478515/1‑93 (MQ=255)
                                      cTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCGGCTTAGCCCTTCGCCATTGACCa                                                    >  1:2426302/1‑93 (MQ=255)
                                      cTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCa                                                    <  2:28773/93‑1 (MQ=255)
                                       tAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAg                                                   >  2:3645773/1‑93 (MQ=255)
                                       tAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAg                                                   <  1:609703/93‑1 (MQ=255)
                                       tAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAg                                                   <  1:3806200/93‑1 (MQ=255)
                                        aaCATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAgg                                                  >  1:2797532/1‑93 (MQ=255)
                                         aCATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAggg                                                 <  1:4690608/93‑1 (MQ=255)
                                         aCATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAggg                                                 <  1:167656/93‑1 (MQ=255)
                                          cATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAgggg                                                >  1:141424/1‑92 (MQ=255)
                                           aTCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAACCCTTCGCCATTGACCAGGGTa                                               >  1:120816/1‑93 (MQ=255)
                                            tCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTAt                                              <  2:3979393/93‑1 (MQ=255)
                                               gcgcCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTca                                           >  1:1168965/1‑93 (MQ=255)
                                               gcgcCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTca                                           >  1:3941595/1‑93 (MQ=255)
                                                cgcCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTcac                                          <  1:1907408/93‑1 (MQ=255)
                                                 gcCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTcaca                                         >  1:488025/1‑93 (MQ=255)
                                                  ccAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACaa                                        >  2:4588142/1‑93 (MQ=255)
                                                   cAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACaaa                                       >  2:359309/1‑93 (MQ=255)
                                                   cAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACaaa                                       <  1:4569864/93‑1 (MQ=255)
                                                    aGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACTAGGGTATTCACaaaa                                      >  2:3262584/1‑93 (MQ=255)
                                                    aGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACaaaa                                      >  2:3096763/1‑93 (MQ=255)
                                                    aGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACaaaa                                      >  1:3762417/1‑93 (MQ=255)
                                                     gAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAAt                                     <  1:2148750/93‑1 (MQ=255)
                                                      aaTGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAATc                                    <  2:1196722/93‑1 (MQ=255)
                                                       aTGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAATcc                                   <  2:2296289/93‑1 (MQ=255)
                                                        tGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAATccc                                  <  2:2479239/93‑1 (MQ=255)
                                                        tGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAATccc                                  >  1:4053903/1‑93 (MQ=255)
                                                         gCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAATcccc                                 >  2:1727119/1‑93 (MQ=255)
                                                          ccAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCCCAAAATccccc                                >  2:410345/1‑92 (MQ=255)
                                                           cAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAATCCccac                               >  2:2151755/1‑93 (MQ=255)
                                                           cAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAATCCccac                               >  1:3358435/1‑93 (MQ=255)
                                                           cAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAATCCccac                               <  1:291926/93‑1 (MQ=255)
                                                            aaCAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAATCCccacc                              >  1:2791488/1‑93 (MQ=255)
                                                             aCAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAATCCccacca                             >  1:367470/1‑93 (MQ=255)
                                                             aCAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAATCCccacca                             >  2:2467850/1‑93 (MQ=255)
                                                             aCAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAATCCccacca                             <  1:388557/93‑1 (MQ=255)
                                                               aGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAATCCCCACCAGa                           <  2:3043838/93‑1 (MQ=255)
                                                                gTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAATCCCCACCAGAt                          >  1:3636829/1‑93 (MQ=255)
                                                                 tAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAATCCCCACCAGAtt                         <  2:613065/93‑1 (MQ=255)
                                                                  aaTCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAATCCCCACCAGAttt                        >  2:2768140/1‑93 (MQ=255)
                                                                  aaTCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAATCCCCACCAGAttt                        >  2:553665/1‑93 (MQ=255)
                                                                  aaTCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAATCCCCACCAGAttt                        >  2:4728055/1‑93 (MQ=255)
                                                                    acAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAATCCCCACCAGATTTCa                      <  2:1872302/92‑1 (MQ=255)
                                                                       gCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAATCCCCACCAGATTTCAGAc                   >  1:4664127/1‑93 (MQ=255)
                                                                          gAACGGCGGTGATCAGAGCGCCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAATCCCCACCAGATTTCAGACCaa                <  1:935574/93‑1 (MQ=255)
                                                                          gAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAATCCCCACCAGATTTCAGACCaa                <  1:1637913/93‑1 (MQ=255)
                                                                          gAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAATCCCCACCAGATTTCAGACCaa                <  1:4476798/93‑1 (MQ=255)
                                                                              ggcggTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAAACCCCACCAGATTTCAGACCAACGAc            <  2:2115386/93‑1 (MQ=255)
                                                                               gcggTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAATCCCCACCAGATTTCAGACCAACGAcc           <  2:2934055/93‑1 (MQ=255)
                                                                                 ggTGATCAGAGCGCCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAATCCCCACCAGATTTCAGACCAACGACCga         <  2:318150/93‑1 (MQ=255)
                                                                                 ggTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAATCCCCACCAGATTTCAGATCAACGACCga         <  2:445620/93‑1 (MQ=255)
                                                                                  gTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAATCCCCACCAGATTTCAGACCAACGACCgag        <  1:3183393/93‑1 (MQ=255)
                                                                                  gTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAATCCCCACCAGATTTCAGACCAACGACCgag        >  1:435572/1‑93 (MQ=255)
                                                                                    gATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAATCCCCACCAGATTTCAGACCAACGACCGagag      >  1:3678191/1‑93 (MQ=255)
                                                                                      tCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAATCCCCACCAGATTTCAGACCAACGACCGAGAGGa    >  2:947136/1‑93 (MQ=255)
                                                                                      acAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAATCCCCACCAGATTTCAGACCAACGACCGAGAGGa    <  1:2117541/92‑1 (MQ=255)
                                                                                        aGGGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAATCCCCACCAGATTTCAGACCAACGACCGAGAGGAAc  >  1:4043008/1‑93 (MQ=255)
                                                                                          |                                                                                          
TGCCAGCATCCCGCACATGGTCAGGTAATTCATTTTGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGGGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAATCCCCACCAGATTTCAGACCAACGACCGAGAGGAAC  >  NC_000913/3862398‑3862578

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