Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,519,159 A→T 70.3% T45S (ACA→TCA)  curA → curcumin/dihydrocurcumin reductase, NADPH‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,519,1590AT70.3% 139.6 / 86.5 103T45S (ACA→TCA) curAcurcumin/dihydrocurcumin reductase, NADPH‑dependent
Reads supporting (aligned to +/- strand):  ref base A (20/10);  new base T (35/37);  total (56/47)
Fisher's exact test for biased strand distribution p-value = 1.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.24e-02

Only 100 of 103 total aligned reads displayed.

CTGGCCTCGCGTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCACAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCACCGCCTGTTGATATTGGCGGCGTGATGGTCGGC  >  NC_000913/1519063‑1519254
                                                                                                |                                                                                               
cTGGCCTCGCGTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAg                                                                                              >  1:2149492/1‑100 (MQ=255)
cTGGCCTCGCGTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAg                                                                                              <  2:2216271/100‑1 (MQ=255)
 tGGCCTCGCGTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGCTGTGGCCACACCGGGTGAAGGACAGGTGTTACTGCGCCCAg                                                                                              >  1:3092816/1‑99 (MQ=255)
  ggCCTCGCGTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTAAAGGACAGGAATTACTGCGCACAGtt                                                                                            >  1:2945755/1‑100 (MQ=255)
   gCCTCGCGTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGttt                                                                                           >  1:2065256/1‑100 (MQ=255)
   gCCTCGCGTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCACAGttt                                                                                           >  2:1321197/1‑100 (MQ=255)
   gCCTCGCGTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCACAGttt                                                                                           <  1:1443081/100‑1 (MQ=255)
     cTCGCGTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCACAGTTTAt                                                                                         <  2:3856211/100‑1 (MQ=255)
      tCGCGTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTAtt                                                                                        <  1:2290254/100‑1 (MQ=255)
      tCGCGTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTAtt                                                                                        <  1:1222483/100‑1 (MQ=255)
       cgcgTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTAttt                                                                                       >  1:1649165/1‑100 (MQ=255)
        gcgTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTg                                                                                      >  2:3612197/1‑100 (MQ=255)
          gTCCACATGGCGCACCTGTTCCGGAGAATTTTCGTCTCGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTc                                                                                    <  2:2782817/100‑1 (MQ=255)
            ccACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCACAGTTTATTTGTccc                                                                                  >  2:645885/1‑100 (MQ=255)
             cacaTGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCt                                                                                 <  2:1873725/100‑1 (MQ=255)
               caTGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCACAGTTTATTTGTCCCTgg                                                                               >  1:4079000/1‑100 (MQ=255)
                aTGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCACAGTTTATTTGTCCCTGGa                                                                              <  1:1190714/100‑1 (MQ=255)
                 tGGCGCACCTGTTCCGGGGAATTTCCGTCTTGAAGAAGATGCTGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGAc                                                                             <  1:1571902/100‑1 (MQ=255)
                  ggCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGAcc                                                                            <  2:1649153/100‑1 (MQ=255)
                  ggCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCACAGTTTATTTGTCCCTGGAcc                                                                            <  2:2018126/100‑1 (MQ=255)
                    cgcACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCg                                                                          <  1:1226537/100‑1 (MQ=255)
                       aCCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGtat                                                                       <  2:1667756/100‑1 (MQ=255)
                        ccTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGGAGGACAGGTGTTACTGCGGTCAGTTTATTTGTCCCTGGACCCGtata                                                                      >  1:627603/1‑100 (MQ=255)
                          tGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACCCCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTatat                                                                     >  2:1426282/1‑99 (MQ=255)
                            ttCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCg                                                                  <  1:3612220/100‑1 (MQ=255)
                              ccGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTg                                                                <  2:2279085/100‑1 (MQ=255)
                              ccGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTg                                                                >  1:3183485/1‑100 (MQ=255)
                              ccGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTg                                                                >  1:1553406/1‑100 (MQ=255)
                               cGGAGTGTTTCCGGCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGGCAGGTGATACTGCGATCAGGTTAGTTGTCCGTGGACCCGGATATGCg                                                                  >  2:3684365/1‑97 (MQ=255)
                               cGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTgg                                                               >  2:285791/1‑100 (MQ=255)
                               cGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCACAGTTTATTTGTCCCTGGACCCGTATATGCGTgg                                                               >  2:2665913/1‑100 (MQ=255)
                               cGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCACAGTTTATTTGTCCCTGGACCCGTATATGCGTgg                                                               >  1:4096543/1‑100 (MQ=255)
                                ggAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCACAGTTTATTTGTCCCTGGACCCGTATATGCGTGGa                                                              >  2:2774250/1‑100 (MQ=255)
                                ggAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCACAGTTTATTTGTCCCTGGACCCGTATATGCGTGGa                                                              >  1:3626364/1‑100 (MQ=255)
                                ggAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCACAGTTTATTTGTCCCTGGACCCGTATATGCGTGGa                                                              <  1:1488234/100‑1 (MQ=255)
                                  agaATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACg                                                            <  2:993788/100‑1 (MQ=255)
                                     aTTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTa                                                          <  2:985005/99‑1 (MQ=255)
                                     aTTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTAt                                                         >  2:3543606/1‑100 (MQ=255)
                                       ttCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGa                                                       >  2:1100775/1‑100 (MQ=255)
                                        tCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAg                                                      <  2:2778125/100‑1 (MQ=255)
                                         ccGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGCTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGc                                                     >  2:2473039/1‑100 (MQ=255)
                                           gTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGa                                                   >  1:301772/1‑100 (MQ=255)
                                              ttGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGa                                                >  1:1572530/1‑100 (MQ=255)
                                                gaagaaGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGc                                              <  1:3629037/100‑1 (MQ=255)
                                                gaagaaGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGc                                              <  2:3761510/100‑1 (MQ=255)
                                                 aagaagATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGcc                                             >  2:3427672/1‑100 (MQ=255)
                                                   gaagaTGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCAt                                           <  1:2756290/100‑1 (MQ=255)
                                                   gaagaTGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCAt                                           <  1:1292308/100‑1 (MQ=255)
                                                      gatgatGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCtt                                        <  1:1426791/100‑1 (MQ=255)
                                                        tgatgTCGCCACGCCGGGTGAATGAGAGGTGTTACTGCGCTCATTTTATTGGTCCCGGGACCGGTGTATGTGTGGACTTATGAGCGATGAGCCATCTTAt                                      <  2:801721/100‑1 (MQ=255)
                                                        tgatgTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTAt                                      >  1:1123872/1‑100 (MQ=255)
                                                        tgatgTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTAt                                      <  1:1272452/100‑1 (MQ=255)
                                                          atgTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCACAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTc                                    <  1:2560830/100‑1 (MQ=255)
                                                           tgTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCGCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCa                                   <  2:3349038/100‑1 (MQ=255)
                                                           tgTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCa                                   >  1:2043715/1‑100 (MQ=255)
                                                           tgTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCa                                   <  1:614194/100‑1 (MQ=255)
                                                           tgTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCACAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCa                                   >  2:1127301/1‑100 (MQ=255)
                                                             tCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCAcc                                 <  1:3555918/100‑1 (MQ=255)
                                                              cGCCACACCGGGTGAAGGACAGGTGTTACTGCGCACAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCACCg                                >  1:3820472/1‑100 (MQ=255)
                                                              cGCCACACCGGGTGAAGGACAGGTGTTACTGCGCACAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCACCg                                >  1:2159136/1‑100 (MQ=255)
                                                               gCCACACCGGGTGAAGGACAGGTGTTACTGCGCACAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCACCGc                               >  1:1389781/1‑100 (MQ=255)
                                                               gCCACACCGGGTGAAGGACAGGTGTTACTGCGCACAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCACCGc                               >  1:971546/1‑100 (MQ=255)
                                                               gCCACACCGGGTGAAGGACAGGTGTTACTGCGCACAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCACCGc                               >  2:1408173/1‑100 (MQ=255)
                                                                ccACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCACCGcc                              <  2:2373426/100‑1 (MQ=255)
                                                                  acacCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCACCGCCTg                            >  2:175557/1‑100 (MQ=255)
                                                                   ccccGGGTGAAGGACAGGTGTTACTGCGCACAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCACCGCCTGt                           <  2:2376458/98‑1 (MQ=255)
                                                                   cacCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCACCGCCTGt                           >  1:443426/1‑100 (MQ=255)
                                                                     ccGGGTGAAGGACAGGTGTTACTGCGCACAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCACCGCCTGTTg                         <  1:4050190/100‑1 (MQ=255)
                                                                      cGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGACCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCACCGCCTGTTGa                        <  2:3846729/100‑1 (MQ=255)
                                                                      cGGGTGAAGGACAGGGGTTACTGCGCTCAGTTTATTAGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCACCGCCTGTTg                         <  2:252334/99‑1 (MQ=255)
                                                                       gggTGAAGGACAGGTGTTACTGCGCACAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCACCGCCTGTTGat                       >  1:2115271/1‑100 (MQ=255)
                                                                          tGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCACCGCCTGTTGATAtt                    >  2:3646824/1‑100 (MQ=255)
                                                                          tGAAGGACAGGTGTTACTGCGCACAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCACCGCCTGTTGATAtt                    >  1:1064353/1‑100 (MQ=255)
                                                                           gAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCACCGCCTGTTGATATTg                   <  2:2065243/100‑1 (MQ=255)
                                                                           gAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCACCGCCTGTTGATATTg                   <  1:1801069/100‑1 (MQ=255)
                                                                           gAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCACCGCCTGTTGATATTg                   >  1:910857/1‑100 (MQ=255)
                                                                            aaGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCACCGCCTGTTGATATTgg                  <  1:937333/100‑1 (MQ=255)
                                                                             aGGAGATGTGTTATTGCGCTCAGTTTATTTGTCCCTTGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCACCGCCTGTTGATATTggc                 <  2:3626125/100‑1 (MQ=255)
                                                                              ngACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCACCGCCTGTTGATATTggcg                >  2:2057407/2‑100 (MQ=255)
                                                                              ggACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCACCGCCTGTTGATATTggcg                >  1:3344391/1‑100 (MQ=255)
                                                                               gACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCACCGCCTGTTGATATTggcgg               >  1:1889399/1‑100 (MQ=255)
                                                                                aCAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCACCGCCTGTTGATATTggcggc              <  1:2175340/100‑1 (MQ=255)
                                                                                 gatctgtTACTGCGCACAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCACCGCCTGTTGATATTGgcggcg             <  2:596565/96‑1 (MQ=255)
                                                                                 cAGGTGTTACTGCGCACAGTTTATTTGTCCCTGCACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCACCGCCTGTTGATATTGgcggcg             <  2:3906852/100‑1 (MQ=255)
                                                                                  aGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCACCGCCTGTTGATATTggcggc              <  2:2772469/98‑1 (MQ=255)
                                                                                   ggTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCACCGCCTGTTGATATTGGCGGCGTg           >  2:839949/1‑100 (MQ=255)
                                                                                   ggTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCACCGCCTGTTGATATTGGCGGCGTg           >  2:420142/1‑100 (MQ=255)
                                                                                   ggTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCACCGCCTGTTGATATTGGCGGCGTg           >  2:2362323/1‑100 (MQ=255)
                                                                                    gtgtTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCGATCTTATTCACCGCCTGTTGATATTGGCGGCGTGa          <  2:1809740/100‑1 (MQ=255)
                                                                                    gtgtTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCACCGCCTGTTGATATTGGCGGCGTGa          >  1:3266844/1‑100 (MQ=255)
                                                                                    gtgtTACTGCGCACAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCACCGCCTGTTGATATTGGCGGCGTg           >  2:585636/1‑99 (MQ=255)
                                                                                     tgtTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCACCGCCTGTTGATATTGGCGGCGTGAt         >  1:2740039/1‑100 (MQ=255)
                                                                                      gtTACTGCGCACAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCACCGCCTGTTGATATTGGCGGCGAGAtc        >  1:596571/1‑99 (MQ=255)
                                                                                        tACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCCTCTTATTCACCGCCTGTTGATATTGGCGGCGTGATGGt      >  2:708161/1‑100 (MQ=255)
                                                                                        tACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCACCGCCTGTTGATATTGGCGGCGTGATGGt      <  1:1790683/100‑1 (MQ=255)
                                                                                         aCTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCACCGCCTGTTGATATTGGCGGCGTGATGGTc     >  1:2826111/1‑100 (MQ=255)
                                                                                         aCTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCACCGCCTGTTGATATTGGCGGCGTGATGGTc     >  1:2358481/1‑100 (MQ=255)
                                                                                         aCTGCGCACAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCACCGCCTGTTGATATTGGCGGCGTGATGGTc     >  1:3425237/1‑100 (MQ=255)
                                                                                          cTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCACCGCCTGTTGATATTGGCGGCGTGATGGTcg    >  1:2275371/1‑100 (MQ=255)
                                                                                            gcgcTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCACCGCCTGTTGATATTGGCGGCGTGATGGTcggc  >  1:860545/1‑100 (MQ=255)
                                                                                                |                                                                                               
CTGGCCTCGCGTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCACAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCACCGCCTGTTGATATTGGCGGCGTGATGGTCGGC  >  NC_000913/1519063‑1519254

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