Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,496,451 T→C G329G (GGT→GGC cysG → uroporphyrinogen III methylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,496,4510TC99.7% 1049.4 / ‑5.2 332G329G (GGT→GGCcysGuroporphyrinogen III methylase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (149/182);  total (149/183)
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.59e-01

Only 100 of 347 total aligned reads displayed.

GAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATTCCGTTCTCGGTGGTTCCGGGTATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGCCCAGAGCGTACGCTTAATTACCGGACA  >  NC_000913/3496361‑3496540
                                                                                          |                                                                                         
gAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCAt                                                                                         >  2:1864919/1‑93 (MQ=255)
 aaaGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCAtt                                                                                        <  2:4347460/93‑1 (MQ=255)
   aGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTAc                                                                                      >  2:1099000/1‑93 (MQ=255)
         cGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGc                                                                                <  1:6269268/93‑1 (MQ=255)
         cGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGc                                                                                <  2:7757875/93‑1 (MQ=255)
           aTCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCtt                                                                              >  1:5387961/1‑93 (MQ=255)
            tCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTc                                                                             <  1:6597035/93‑1 (MQ=255)
            tCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTc                                                                             <  2:4887677/93‑1 (MQ=255)
             ccGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCt                                                                            <  1:818684/93‑1 (MQ=255)
               gTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAATCACTGTGCAACGCGGTTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTgg                                                                          <  2:7108639/93‑1 (MQ=255)
                tttATTTTTGGCCGTGGTGGCTAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGt                                                                         >  1:1378772/1‑93 (MQ=255)
                tttATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGt                                                                         <  2:5881902/93‑1 (MQ=255)
                tttATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGt                                                                         >  2:5114797/1‑93 (MQ=255)
                     ttttGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGctc                                                                    <  2:2751164/93‑1 (MQ=255)
                     ttttGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGctc                                                                    <  2:2560824/93‑1 (MQ=255)
                     ttttGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGctc                                                                    <  2:4972455/93‑1 (MQ=255)
                     ttttGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGctc                                                                    >  1:7549166/1‑93 (MQ=255)
                     ttttGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGctc                                                                    <  2:660032/93‑1 (MQ=255)
                      tttGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGctct                                                                   <  1:3912316/93‑1 (MQ=255)
                       ttGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTg                                                                  <  2:5222136/93‑1 (MQ=255)
                           ccGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTa                                                              >  2:4485417/1‑93 (MQ=255)
                             gtggtgGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTAtt                                                            >  2:5339349/1‑93 (MQ=255)
                             gtggtgGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTAtt                                                            <  2:845128/93‑1 (MQ=255)
                             gtggtgGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTAtt                                                            <  1:3548998/93‑1 (MQ=255)
                             gtggtgGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTAtt                                                            <  2:7501970/93‑1 (MQ=255)
                              tggtggCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTc                                                           >  1:1006254/1‑93 (MQ=255)
                                gtggCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCgg                                                         <  2:1985701/93‑1 (MQ=255)
                                 tggCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCggg                                                        >  1:5863248/1‑93 (MQ=255)
                                   gCGAAGAGCTGGAAACCCTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTa                                                      <  2:2183181/93‑1 (MQ=255)
                                   gCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTa                                                      <  2:359539/93‑1 (MQ=255)
                                   gCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTa                                                      >  2:146330/1‑93 (MQ=255)
                                   gCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTa                                                      <  2:42046/93‑1 (MQ=255)
                                   gCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTa                                                      <  1:4404721/93‑1 (MQ=255)
                                    cGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTAt                                                     <  2:469136/93‑1 (MQ=255)
                                     gAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTAtt                                                    <  1:6346377/93‑1 (MQ=255)
                                     gAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTAtt                                                    >  1:6808332/1‑93 (MQ=255)
                                      aaGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTc                                                   >  1:3256571/1‑93 (MQ=255)
                                      aaGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTc                                                   <  2:6663336/93‑1 (MQ=255)
                                      aaGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTc                                                   >  2:1724965/1‑93 (MQ=255)
                                      aaGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTc                                                   <  1:3645193/93‑1 (MQ=255)
                                       agagCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTcc                                                  >  2:6151619/1‑93 (MQ=255)
                                       agagCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTcc                                                  >  2:4719485/1‑93 (MQ=255)
                                       agagCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTcc                                                  >  2:860127/1‑93 (MQ=255)
                                           cTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCTGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTc                                              >  1:6985013/1‑93 (MQ=255)
                                             ggAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCAc                                            >  1:4768415/1‑93 (MQ=255)
                                             ggAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCAc                                            >  1:7793858/1‑93 (MQ=255)
                                              gAAACACTGTGCAACGCGGGTATCTCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACg                                           >  2:3348686/1‑93 (MQ=255)
                                               aaaCACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGc                                          >  1:1698188/1‑93 (MQ=255)
                                                aaCACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCa                                         >  1:2444449/1‑93 (MQ=255)
                                                aaCACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCa                                         >  2:3082415/1‑93 (MQ=255)
                                                aaCACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCa                                         >  1:6935175/1‑93 (MQ=255)
                                                 acacTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCAt                                        >  1:4531509/1‑93 (MQ=255)
                                                  cacTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATc                                       <  2:5764207/93‑1 (MQ=255)
                                                   acTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCATGCATcg                                      >  1:5165833/1‑93 (MQ=255)
                                                   acTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACCCATcg                                      >  1:6424056/1‑93 (MQ=255)
                                                     tgtgCAACGCGGGCATCCCGTTCTCGGTTGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATcgcg                                    >  2:5428367/1‑93 (MQ=255)
                                                        gCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGAtt                                 >  2:1644889/1‑93 (MQ=255)
                                                        gCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGAtt                                 >  1:61292/1‑93 (MQ=255)
                                                         cAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTCCCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTa                                <  2:6325561/93‑1 (MQ=255)
                                                         cAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCCCTCACGCATCGCGATTa                                <  1:2631181/93‑1 (MQ=255)
                                                         cAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTa                                <  2:4216897/93‑1 (MQ=255)
                                                         cAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTa                                <  2:3839649/93‑1 (MQ=255)
                                                            cgcgGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGc                             <  2:4627722/93‑1 (MQ=255)
                                                              cgGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGccc                           <  2:2409757/93‑1 (MQ=255)
                                                              cgGGAATCCCGTTCTCGGTGGTTCCTGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCTACTCACGCATCGCGATTATGccc                           >  2:4002717/1‑93 (MQ=255)
                                                               gggTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGCCCa                          <  1:427924/93‑1 (MQ=255)
                                                                ggTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGCCCag                         >  1:7065201/1‑93 (MQ=255)
                                                                 gTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTACTCCACTCACGCATCGCGATTATGCCCaga                        <  1:6508614/93‑1 (MQ=255)
                                                                  tATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGCCCagag                       <  1:4700695/93‑1 (MQ=255)
                                                                   atcccGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGCCCAGAGc                      <  1:347750/90‑1 (MQ=255)
                                                                   atcccGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGCCCAGAGc                      <  1:4200540/90‑1 (MQ=255)
                                                                    tcccGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGCCCAGAGCg                     >  2:6862526/3‑93 (MQ=255)
                                                                    tcccGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGCCCAGAGCg                     >  2:3044198/3‑93 (MQ=255)
                                                                     cccGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGCCCAGAGCGt                    <  1:6707285/92‑1 (MQ=255)
                                                                      ccGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGCCCAGAGCGTa                   <  2:7609194/93‑1 (MQ=255)
                                                                      ccGTTCTCGGTGGTTCCGGGCAGTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCAGCGAGATTATCCCCAGAGCGTa                   <  2:719662/93‑1 (MQ=255)
                                                                       cGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGCCCAGAGCGTAc                  >  2:2344468/1‑93 (MQ=255)
                                                                       cGTTCTCGGTGGTGCCGCGTATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCCCTCACGCATCGCGATTATGCCCAGAGCGTAc                  <  2:1139579/93‑1 (MQ=255)
                                                                       cGTTCTCGGGGGTTCCGGGCATTCCCGCAGCTTCTGGTTGCTCTGCCTATTCGGCTATTCCACTCACGCATCGCGATTATGCCCAGAGCGTAc                  <  2:375813/93‑1 (MQ=255)
                                                                         ttCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGCCCAGAGCGTACGc                >  2:1899351/1‑93 (MQ=255)
                                                                         ttCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGCCCAGAGCGTACGc                >  1:1040517/1‑93 (MQ=255)
                                                                         ttCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGCCCAGAGCGTACGc                >  2:3301818/1‑93 (MQ=255)
                                                                         ttCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGCCCAGAGCGTACGc                >  1:3174259/1‑93 (MQ=255)
                                                                           ctcGGTGGTTCCGGGCATTCCCGCAGCTTCTGGTTGCTCTGCCTATTAGGGTATTCCACTCACGCATCGCGATTATGCCCAGAGAGTACACtt              <  2:7036725/93‑1 (MQ=255)
                                                                            tcGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGCCCAGAGCGTACGCTTa             >  1:4032794/1‑93 (MQ=255)
                                                                            tcGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGCCCAGAGCGTACGCTTa             >  1:1186685/1‑93 (MQ=255)
                                                                             cGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGCCCAGAGCGTACGCTTaa            <  1:7354902/93‑1 (MQ=255)
                                                                              ggtggtTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGCCCAGAGCGTACGCTTAAt           <  1:1562558/93‑1 (MQ=255)
                                                                                 ggtTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGCCCAGAGCGTACGCTTAATTAc        <  2:1332858/93‑1 (MQ=255)
                                                                                 ggtTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGCCCAGAGCGTACGCTTAATTAc        >  2:7322550/1‑93 (MQ=255)
                                                                                 ggtTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGCCCAGAGCGTACGCTTAATTAc        >  2:5687302/1‑93 (MQ=255)
                                                                                  gtTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGCCCAGAGCGTACGCTTAATTAcc       >  1:2244413/1‑93 (MQ=255)
                                                                                   ttCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGCCCAGAGCGTACGCTTAATTACCg      <  1:2063164/93‑1 (MQ=255)
                                                                                     ccGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGCCCAGAGCGTACGCTTAATTACCGGa    <  1:5494218/93‑1 (MQ=255)
                                                                                      cGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGCCCAGAGCGTACGCTTAATTACCGGac   <  2:6318278/93‑1 (MQ=255)
                                                                                      cGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGCCCAGAGCGTACGCTTAATTACCGGac   >  2:4528126/1‑93 (MQ=255)
                                                                                      cGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGCCCAGAGCGTACGCTTAATTACCGGac   <  2:2707863/93‑1 (MQ=255)
                                                                                       gggCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGCCCAGAGCGTACGCTTAATTACCGGaca  >  2:593082/1‑93 (MQ=255)
                                                                                       gggCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGCCCAGAGCGTACGCTTAATTACCGGaca  >  1:1947946/1‑93 (MQ=255)
                                                                                       gggCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGCCCAGAGCGTACGCTTAATTACCGGaca  >  2:4394781/1‑93 (MQ=255)
                                                                                          |                                                                                         
GAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATTCCGTTCTCGGTGGTTCCGGGTATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGCCCAGAGCGTACGCTTAATTACCGGACA  >  NC_000913/3496361‑3496540

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: