Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,291,923 A→C I138L (ATC→CTC)  nrfE → heme lyase (NrfEFG) for insertion of heme into c552, subunit NrfE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,291,9230AC99.5% 726.0 / ‑4.2 201I138L (ATC→CTC) nrfEheme lyase (NrfEFG) for insertion of heme into c552, subunit NrfE
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base C (92/108);  total (92/109)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.12e-01

Only 100 of 201 total aligned reads displayed.

CCGCTGACGTTAGCCGTTTTATCTCTCATGCTCGCCGCACTGCTACTGTTTGTGGTGCTGTGGTCCGATCCCTTCGTGCGGATATTTCCACCAGCAATCGAAGGCCGCGATCTCAATCCGATGCTGCAACATCCCGGTCTTATCTTTCATCCACCGCTGCTTTACCTTGGCTATGGCGGTTTGATGGTAGC  >  NC_000913/4291827‑4292017
                                                                                                |                                                                                              
ccGCTGACGTTAGCCGTTTTATCTCTCATGCTCGCCGCACTGCTACTGTTTGTGGTGCTGTGGTCCGATCCCTTCGTGCGGATATTTCCACCAGCACTCg                                                                                             <  1:4884372/100‑1 (MQ=255)
 cGCTGACGTTAGCCGTTTTATCTCTCATGCTCGCCGCACTGCTACTGTTTGTGGTGCTGTGGTCCGATCCCTTCGTGCGGATATTTCCACCAGCACTCGa                                                                                            >  1:2905350/1‑100 (MQ=255)
   cTGACGTTAGCCGTTTTATCTCTCATGCTCGCCGCACTGCTACTGTTTGTGGTGCTGTGGTCCGATCCCTTCGTGCGGATATTTCCACCAGCACTCGAAg                                                                                          >  2:863154/1‑100 (MQ=255)
   cTGACGTTAGCCGTTTTATCTCTCATGCTCGCCGCACTGCTACTGTTTGTGGTGCTGTGGTCCGATCCCTTCGTGCGGATATTTCCACCAGCACTCGAAg                                                                                          <  2:5314323/100‑1 (MQ=255)
    tGACGTTAGCCGTTTTATCTCTCATGCTCGCCGCACTGCTACTGTTTGTGGTGCTGTGGTCCGATCCCTTCGTGCGGATATTTCCACCAGCACTCGAAgg                                                                                         <  1:3405487/100‑1 (MQ=255)
     gACGTTAGCCGTTTTATCTCTCATGCTCGCCGCACTGCTACTGTTTGTGGTGCTGTGGTCCGATCCCTTCGTGCGGATATTTCCACCAGCACTCGAAGGc                                                                                        <  1:1227576/100‑1 (MQ=255)
     gACGTTAGCCGTTTTATCTCTCATGCTCGCCGCACTGCTACTGTTTGTGGTGCTGTGGTCCGATCCCTTCGTGCGGATATTTCCACCAGCACTCGAAGGc                                                                                        <  1:1920844/100‑1 (MQ=255)
           aGCCGTTTTATCTCTCATGCTCGCCGCTCTGCTACTGTTTGTGGTGCTGTGGTCCGATCCCTTCGTGCGGATATTTCCACCAGCACTCGAAGGCCGCGAt                                                                                  <  1:543907/100‑1 (MQ=255)
           aGCCGTTTTATCTCTCATGCTCGCCGCACTGCTACTGTTTGTGGTGCTGTGGTCCGATCCCTTCGTGCGGATATTTCCACCAGCACTCGAAGGCCGCGAt                                                                                  <  2:2674183/100‑1 (MQ=255)
            gCCGTTTTATCTCTCATGCTCGCCGCACTGCTACTGTTTGTGGTGCTGTGGTCCGATCCCTTCGTGCGGATATTTCCACCAGCACTCGAAGGCCGCGAtc                                                                                 <  1:1314697/100‑1 (MQ=255)
            gCCGTTTTATCTCTCATGCTCGCCGCACTGCTACTGTTTGTGGTGCTGTGGTCCGATCCCTTCGTGCGGATATTTCCACCAGCACTCGAAGGCCGCGAtc                                                                                 >  1:606797/1‑100 (MQ=255)
             ccGTTTTATCTCTCATGCTCGCCGCACTGCTACTGTTTGTGGTGCTGTGGTCCGATCCCTTCGTGCGGATATTTCCACCAGCACTCGAAGGCCGCGAtct                                                                                >  1:4309510/1‑100 (MQ=255)
               gTTTTATCTCTCATGCTCGCCGCACTGCTACTGTTTGTGGTGCTGTGGTCCGATCCCTTCGTGCGGATATTTCCACCAGCACTCGAAGGCCGCGATCTCa                                                                              <  2:2295352/100‑1 (MQ=255)
                ttttATCTCTCATGCTCGCCGCACTGCTACTGTTTGTGGTGCTGTGGTCCGATCCCTTCGTGCGGATATTTCCACCAGCACTCGAAGGCCGCGATCTCaa                                                                             >  1:3474356/1‑100 (MQ=255)
                  ttATCTCTCATGCTCGCCGCACTGCTACTGTTTGTGGTGCTGTGGTCCGATCCCTTCGTGCGGATATTTCCACCAGCACTCGAAGGCCGCGATCTCAATc                                                                           <  2:582581/100‑1 (MQ=255)
                  ttATCTCTCATGCTCGCCGCACTGCTACTGTTTGTGGTGCTGTGGTCCGATCCCTTCGTGCGGATATTTCCACCAGCACTCGAAGGCCGCGATCTCAATc                                                                           <  2:3240014/100‑1 (MQ=255)
                   tgtctctcATGCTCGCCGCACTGCTACTGTTTGTGGTGCTGTGGTCCGATCGCTTCGTGCGGATATTTCCACCAGCACTCGAAGGCCGCGATCTCAATcc                                                                          <  2:5472859/98‑1 (MQ=255)
                   tATCTCTCATGCTCGCCGCACTGCTACTGTTTGTGGTGCTGTGGTCCGATCCCTTCGTGCGGATATTTCCACCAGCACTCGAAGGCCGCGATCTCAATcc                                                                          <  1:3757174/100‑1 (MQ=255)
                   tATCTCTCATGCTCGCCGCACTGCTACTGTTTGTGGTGCTGTGGTCCGATCCCTTCGTGCGGATATTTCAACCAGCACTCGAAGGCCGCGATCTCAATcc                                                                          >  1:5305189/1‑100 (MQ=255)
                      ctctcATGCTCGCCGCACTGCTACTGTTTGTGGTGCTGTGGTCCGGTCCCTTCGTGCGGATATTTCCACCAGCACTCGAAGGCCGCGATCTCAATCCGAt                                                                       >  2:1181203/1‑100 (MQ=255)
                      ctctcATGCTCGCCGCACTGCTACTGTTTGTGGTGCTGTGGTCCGATCCCTTCGTGCGGATATTTCCACCAGCACTCGAAGGCCGCGATCTCAATCCGAt                                                                       <  2:1494778/100‑1 (MQ=255)
                      ctctcATGCTCGCCGCACTGCTACTGTTTGTGGTGCTGTGGTCCGATCCCCTCGTGCGGATATTTCCACCAGCACTCGAAGGCCGCGATCTCAATCCGAt                                                                       <  1:1146837/100‑1 (MQ=255)
                       tctcATGCTCGCCGCACTGCTACTGTTTGTGGTGCTGTGGTCCGATCCCTTCGTGCGGATATTTCCACCAGCACTCGAAGGCCGCGATCTCAATCCGAtg                                                                      >  2:5263144/1‑100 (MQ=255)
                       tctcATGCTCGCCGCACTGCTACTGTTTGTGGTGCTGTGGTCCGATCCCTTCGTGCGGATATTTCCACCAGCACTCGAAGGCCGCGATCTCAATCCGAtg                                                                      <  1:2142103/100‑1 (MQ=255)
                        ctcATGCTCGCCGCACTGCTACTGTTTGTGGTGCTGTGGTCCGATCCCTTCGTGCGGATATTTCCACCAGCACTCGAAGGCCGCGATCTCAATCCGAtgc                                                                     <  2:3123377/100‑1 (MQ=255)
                        ctcATGCTCGCCGCACTGCTACTGTTTGTGGTGCTGTGGTCCGATCCCTTCGTGCGGATATTTCCACCAGCACTCGAAGGCCGCGATCTCAATCCGAtgc                                                                     >  1:1623099/1‑100 (MQ=255)
                         tcATGCTCGCCGCACTGCTACTGTTTGTGGTGCTGTGGTCCGATCCCTTCGTGCGGATATTTCCACCAGCACTCGAAGGCCGCGATCTCAATCCGAtgct                                                                    >  2:472031/1‑100 (MQ=255)
                         tcATGCTCGCCGCACTGCTACTGTTTGTGGTGCTGTGGTCCGATCCCTTCGTGCGGATATTTCCACCAGCACTCGAAGGCCGCGATCTCAATCCGAtgct                                                                    <  2:2998337/100‑1 (MQ=255)
                         tcATGCTCGCCGCACTGCTACTGTTTGTGGTGCTGTGGTCCGATCCCTTCGTGCGGATATTTCCACCAGCACTCGAAGGCCGCGATCTCAATCCGAtgct                                                                    >  1:1701724/1‑100 (MQ=255)
                           aTGCTCGCCGCACTGCTACTGTTTGTGGTGCTGTGGTCCGATCCCTTCGTGCGGATATTTCCACCAGCACTCGAAGGCCGCGATCTCAATCCGAtgctgc                                                                  >  1:3566478/1‑100 (MQ=255)
                            tGCTCGCCGCACTGCTACTGTTTGTGGTGCTGTGGTCCGATCCCTTCGTGCGGATATTTCCACCAGCACTCGAAGGCCGCGATCTCAATCCGATGCTGCa                                                                 >  1:10932/1‑100 (MQ=255)
                             gCTCGCCGCACTGCTACTGTTTGTGGTGCTGTGGTCCGATCCCTTCGTGCGGATATTTCCACCAGCACTCGAAGGCCGCGATCTCAATCCGATGCTGCaa                                                                >  2:4792986/1‑100 (MQ=255)
                             gCTCGCCGCACTGCTACTGTTTGTGGTGCTGTGGTCCGATCCCTTCGTGCGGATATTTCCACCAGCACTCGAAGGCCGCGATCTCAATCCGATGCTGCaa                                                                >  2:4271933/1‑100 (MQ=255)
                              cTCGCCGCACTGCTACTGTTTGTGGTGCTGTGGTCCGATCCCTTCGTGCGGATATTTCCACCAGCACTCGAAGGCCGCGATCTCAATCCGATGCTGCaan                                                               <  1:1391879/100‑2 (MQ=255)
                              cTCGCCGCACTGCTACTGTTTGTGGTGCTGTGGTCCGATCCCTTCGTGCGGATATTTCCACCAGCACTCGAAGGCCGCGATCTCAATCCGATGCTGCAAc                                                               >  2:1715230/1‑100 (MQ=255)
                              cTCGCCGCACTGCTACTGTTTGTGGTGCTGTGGTCCGATCCCTTCGTGCGGATATTTCCACCAGCACTCGAAGGCCGCGATCTCAATCCGATGCTGCAAc                                                               >  2:4473795/1‑100 (MQ=255)
                                 gccgcACTGCTACTGTTTGTGGTGCTGTGGTCCGATCCCTTCGTGCGGATATTTCCACCAGCACTCGAAGGCCGCGATCTCAATCCGATGCTGCAACATc                                                            >  1:2293801/1‑100 (MQ=255)
                                 gccgcACTGCTACTGTTTGTGGTGCTGTGGTCCGATCCCTTCGTGCGGATATTTCCACCAGCACTCGAAGGCCGCGATCTCAATCCGATGCTGCAACATc                                                            >  1:1688007/1‑100 (MQ=255)
                                  ccgcACTGCTACTGTTTGTGGTGCTGTGGTCCGATCCCTTCGTGCGGATATTTCCACCAGCACTCGAAGGCCGCGATCTCAATCCGATGCTGCAACATcc                                                           >  2:1710765/1‑100 (MQ=255)
                                    gcACTGCTACTGTTTGTGGTGCTGTGGTCCGATCCCTTCGTGCGGATATTTCCACCAGCACTCGAAGGCCGCGATCTCAATCCGATGCTGCAACATCCCg                                                         >  1:4087638/1‑100 (MQ=255)
                                     gcctccTACTGTTGGTGGGGCTGGGGTCGGTTCCCTACCTGCGGGTCTTTCCCCCGGCACTCGAAGGCCGCGATCTCAATCCGATGCTGCAACATCCCgg                                                        <  2:4962117/95‑1 (MQ=255)
                                     cACTGCTACTGTTTGTGGTGCTGTGGTCCGATCCCTTCGTGCGGATATTTCCACCAGCACTCGAAGGCCGCGATCTCAATCCGATGCTGCAACATCCCgg                                                        <  1:964586/100‑1 (MQ=255)
                                     cACTGCTACTGTTTGTGGTGCTGTGGTCCGATCCCTTCGTGCGGATATTTCCACCAGCACTCGAAGGCCGCGATCTCAATCCGATGCTGCAACATCCCgg                                                        <  2:3513511/100‑1 (MQ=255)
                                         gCTACTGTTTGTGGTGCTGTGGTCCGATCCCTTCGTGCGGATATTTCCACCAGCACTCGAAGGCCGCGATCTCAATCCGATGCTGCAACATCCCGGTCtt                                                    >  2:2569373/1‑100 (MQ=255)
                                          cTACTGTTTGTGGTGCTGTGGTCCGATCCCTTCGTGCGGATATTTCCACCAGCACTCGAAGGCCGCGATCTCAATCCGATGCTGCAACATCCCGGTCTTa                                                   <  1:1458265/100‑1 (MQ=255)
                                          cTACTGTTTGTGGTGCTGTGGTCCGATCCCTTCGTGCGGATATTTCCACCAGCACTCGAAGGCCGCGATCTCAATCCGAGGCCGCAACATCCCGGCCTTa                                                   >  2:1924295/1‑100 (MQ=255)
                                            aCTGTTTGTGGTGCTGTGGTCCGATCCCTTCGTGCGGATATTTCCACCAGCACTCGAAGGCCGCGATCTCAATCCGATGCTGCAACATCCCGGTCTTATc                                                 >  2:1785365/1‑100 (MQ=255)
                                             cTGTTTGTGGTGCTGTGGTCCGATCCCTTCGTGCGGATATTTCCACCAGCACTCGAAGGCCGCGATCTCAATCCGATGCTGCAACATCCCGGTCTTATCt                                                >  2:5382/1‑100 (MQ=255)
                                              tGTTTGTGGTGCTGTGGTCCGATCCCTTCGTGCGGATATTTCCACCAGCACTCGAAGGCCGCGATCTCAATCCGATGCTGCAACATCCCGGTCTTATCtt                                               <  1:1548995/100‑1 (MQ=255)
                                              tGTTTGTGGTGCTGTGGTCCGATCCCTTCGTGCGGATATTTCCACCAGCACTCGAAGGCCGCGATCTCAATCCGATGCTGCAACATCCCGGTCTTATCtt                                               <  1:3271693/100‑1 (MQ=255)
                                               gTTTGTGGTGCTGTGGTCCGATCCCTTCGTGCGGATATTTCCACCAGCACTCGAAGGCCGCGATCTCAATCCGATGCTGCAACATCCCGGTCTTATCttt                                              >  1:1516000/1‑100 (MQ=255)
                                               gTTTGTGGTGCTGTGGTCCGATCCCTTCGTGCGGATATTTCCACCAGCACTCGAAGGCCGCGATCTCAATCCGATGCTGCAACATCCCGGTCTTATCttt                                              >  1:4375483/1‑100 (MQ=255)
                                                tttGTGGTGCTGTGGTCCGATCCCTTCGTGCGGATATTTCCACCAGCACTCGAAGGCCGCGATCTCAATCCGATGCTGCAACATCCCGGTCTTATCTTTc                                             <  1:564121/100‑1 (MQ=255)
                                                tttGTGGTGCTGTGGTCCGATCCCTTCGTGCGGATATTTCCACCAGCACTCGAAGGCCGCGATCTCAATCCGATGCTGCAACATCCCGGTCTTATCTTTc                                             <  2:3533508/100‑1 (MQ=255)
                                                  tgtgGTGCTGTGGTCCGATCCCTTCGTGCGGATATTTCCACCAGCACTCGAAGGCCGCGATCTCAATCCGATGCTGCAACATCCCGGTCTTATCTTTCAt                                           <  2:3998129/100‑1 (MQ=255)
                                                   gtggtgCTGTGGTCCGATCCCTTCGTGCGGATATTTCCACCAGCACTCGAAGGCCGCGATCTCAATCCGATGCTGCAACATCCCGGTCTTATCTTTCATc                                          >  2:491829/1‑100 (MQ=255)
                                                       tgCTGTGGTCCGATCCCTTCGTGCGGATATTTCCACCAGCACTCGAAGGCCGCGATCTCAATCCGATGCTGCAACATCCCGGTCTTATCTTTCATCCAcc                                      >  1:4811216/1‑100 (MQ=255)
                                                         cTGTGGTCCGATCCCTTCGTGCGGATATTTCCACCAGCACTCGAAGGCCGCGATCTCAATCCGATGCTGCAACATCCCGGTCTTATCTTTCATCCACCgc                                    <  1:3314329/100‑1 (MQ=255)
                                                         cTGTGGTCCGATCCCTTCGTGCGGATATTTCCACCAGCACTCGAAGGCCGCGATCTCAATCCGATGCTGCAACATCCCGGTCTTATCTTTCATCCACCgc                                    <  1:1302891/100‑1 (MQ=255)
                                                         cTGTGGTCCGATCCCTTCGTGCGGATATTTCCACCAGCACTCGAAGGCCGCGATCTCAATCCGATGCTGCAACATCCCGGTCTTATCTTTCATCCACCgc                                    >  1:843463/1‑100 (MQ=255)
                                                          tgtgGTCCGATCCCTTCGTGCGGATATTTCCACCAGGACTCGAAGGCCGCGATCTCAATCCGATGCTGCAACATCCCGGTCTTATCTTTCATCCACCgct                                   >  2:2037089/1‑100 (MQ=255)
                                                          tgtgGTCCGATCCCTTCGTGCGGATATTTCCACCAGCACTCGAAGGCCGCGATCTCAATCCGATGCTGCAACATCCCGGTCTTATCTTTCATCCACCgct                                   >  1:3861070/1‑100 (MQ=255)
                                                          tgtgGTCCGATCCCTTCGTGCGGATATTTCCACCAGCACTCGAAGGCCGCGATCTCAATCCGATGCTGCAACATCCCGGTCTTATCTTTCATCCACCgct                                   <  1:1485425/100‑1 (MQ=255)
                                                           gtgGTCCGATCCCTTCGTGCGGATATTTCCACCAGCACTCGAAGGCCGCGATCTCAATCCGATGCTGCAACATCCCGGTCTTATCTTTCATCCACCgctg                                  >  2:654799/1‑100 (MQ=255)
                                                             ggTCCGATCCCTTCGTGCGGATATTTCCACCAGCACTCGAAGGCCGCGATCTCAATCCGATGCTGCAACATCCCGGTCTTATCTTTCATCCACCgctgct                                >  2:3684182/1‑100 (MQ=255)
                                                             ggTCCGATCCCTTCGTGCGGATATTTCCACCAGCACTCGAAGGCCGCGATCTCAATCCGATGCTGCAACATCCCGGTCTTATCTTTCATCCACCgctgct                                <  1:5527238/100‑1 (MQ=255)
                                                              gTCCGATCCCTTCGTGCGGATATTTCCACCAGCACTCGAAGGCCGCGATCTCAATCCGATGCTGCAACATCCCGGTCTTATCTTTCATCCACCGCTGCtt                               >  1:655589/1‑100 (MQ=255)
                                                              gTCCGATCCCTTCGTGCGGATATTTCCACCAACACTCGAAGGCCGCGATCTCAATCCGATGCTGCAACATCCCGGTCTTATCTTTCATCCACCGCTGCtt                               >  1:889961/1‑100 (MQ=255)
                                                               tCCGATCCCTTCGTGCGGATATTTCCACCAGCACTCGAAGGCCGCGATCTCAATCCGATGCTGCAACATCCCGGTCTTATCTTTCATCCACCGCTGCttt                              >  1:2240586/1‑100 (MQ=255)
                                                                ccGATCCCTTCGTGCGGATATTTCCACCAGCACTCGAAGGCCGCGATCTCAATCCGATGCTGCAACATCCCGGTCTTATCTTTCATCCACCGCTGCTTTa                             <  1:2088119/100‑1 (MQ=255)
                                                                ccGATCCCTTCGTGCGGATATTTCCACCAGCACTCGAAGGCCGCGATCTCAATCCGATGCTGCAACATCCCGGTCTTATCTTTCATCCACCGCTGCTTTa                             <  2:2830898/100‑1 (MQ=255)
                                                                ccGATCCCTTCGTGCGGATATTTCCACCAGCACTCGAAGGCCGCGATCTCAATCCGATGCTGCAACATCCCGGTCTTATCTTTCATCCACCGCTGCTTTa                             <  1:5415544/100‑1 (MQ=255)
                                                                  gATCCCTTCGTGCGGATATTTCCACCAGCACTCGAAGGCCGCGATCTCAATCCGATGCTGCAACATCCCGGTCTTATCTTTCATCCACCGCTGCTTTAcc                           >  1:1861752/1‑100 (MQ=255)
                                                                   aTCCCTTCGTGCGGATATTTCCACCAGCACTCGAAGGCCGCGATCTCAATCCGATGCTGCAACATCCCGGTCTTATCTTTCATCCACCGCTGCTTTACCt                          <  1:2637763/100‑1 (MQ=255)
                                                                    tCCCTTCGTGCGGATATTTCCACCAGCACTCGAAGGCCGCGATCTCAATCCGATGCTGCAACATCCCGGTCTTATCTTTCATCCACCGCTGCTTTACCtt                         <  2:2781402/100‑1 (MQ=255)
                                                                    tCCCTTCGTGCGGATATTTCCACCAGCACTCGAAGGCCGCGAGCTCAATCCGATGCTGCAACATCCCGGTCTTATCTTTCATCCACCGCTGCTTTACCtt                         >  1:4747753/1‑100 (MQ=255)
                                                                      ccTTCGTGCGGATATTTCCACCAGCACTCGAAGGCCGCGATCTCAATCCGATGCTGCAACATCCCGGTCTTATCTTTCATCCACCGCTGCTTTACCTTgg                       <  1:34420/100‑1 (MQ=255)
                                                                      ccTTCGTGCGGATATTTCCACCAGCACTCGAAGGCCGCGATCTCAATCCGATGCTGCAACATCCCGGTCTTATCTTTCATCCACCGCTGCTTTACCTTgg                       <  1:727964/100‑1 (MQ=255)
                                                                      ccTTCGTGCGGATATTTCCACCAGCACTCGAAGGCCGCGATCTCAATCCGATGCTGCAACATCCCGGTCTTATCTTTCATCCACCGCTGCTTTACCTTgg                       <  2:221831/100‑1 (MQ=255)
                                                                       cTTCGTGCGGATATTTCCACCAGCACTCGAAGGCCGCGATCTCAATCCGATGCTGCAACATCCCGGTCTTATCTTTCATCCACCGCTGCTTTACCTTGGc                      <  2:1372285/100‑1 (MQ=255)
                                                                       cTTCGTGCGGATATTTCCACCAGCACTCGAAGGCCGCGATCTCAATCCGATGCTGCAACATCCCGGTCTTATCTTTCATCCACCGCTGCTTTACCTTGGc                      <  2:2243590/100‑1 (MQ=255)
                                                                          cGTGCGGATATTTCCACCAGCACTCGAAGGCCGCGATCTCAATCCGATGCTGCAACATCCCGGTCTTATCTTTCATCCACCGCTGCTTTACCTTGGCTAt                   <  1:3053999/100‑1 (MQ=255)
                                                                          cGTGCGGATATTTCCACCAGCACTCGAAGGCCGCGATCTCAATCCGATGCTGCAACATCCCGGTCTTATCTTTCATCCACCGCTGCTTTACCTTGGCTAt                   <  1:27006/100‑1 (MQ=255)
                                                                               ggATATTTCCACCAGCACTCGAAGGCCGCGATCTCAATCCGATGCTGCAACATCCCGGTCTTATCTTTCATCCACCGCTGCTTTACCTTGGCTATGGCgg              >  2:769084/1‑100 (MQ=255)
                                                                               ggATATTTCCACCAGCACTCGAAGGCCGCGATCTCAATCCGATGCTGCAACATCCCGGTCTTATCTTTCATCCACCGCTGCTTTACCTTGGCTATGGCgg              <  2:2542480/100‑1 (MQ=255)
                                                                               ggATATTACCACCAGCACTCGAAGGCCGCGATCTCAATCCGATGCTGCAACATGCCGGTCTTATCTTTCATCCACCGCTGCTTTACCTTGGCTATGGCgg              <  2:5115906/100‑1 (MQ=255)
                                                                                  tatTTCCACCAGCACTCGAAGGCCGCGATCTCAATCCGATGCTGCAACATCCCGGTCTTATCTTTCATCCACCGCTGCTTTACCTTGGCTATGGCGGttt           >  1:3054601/1‑100 (MQ=255)
                                                                                    tttCCACCAGCACTCGAAGGCCGCGATCTCAATCCGATGCTGCAACATCCCGGTCTTATCTTTCATCCACCGCTGCTTTACCTTGGCTATGGCGGTTTGa         >  1:5187438/1‑100 (MQ=255)
                                                                                    tttCCACCAGCACTCGAAGGCCGCGATCTCAATCCGATGCTGCAACATCCCGGTCTTATCTTTCATCCACCGCTGCTTTACCTTGGCTATGGCGGTTTGa         <  2:3558244/100‑1 (MQ=255)
                                                                                    tttCCACCAGCACTCGAAGGCCGCGATCTCAATCCGATGCTACAACATCCCGGTCTTATCTTTCATCCACCGCTGCTTTACCTTGGCTATGGCGGTTTGa         >  1:3595148/1‑100 (MQ=255)
                                                                                     ttCCACCAGCACTCGAAGGCCGCGATCTCAATCCGATGCTGCAACATCCCGGTCTTATCTTTCATCCACCGCTGCTTTACCTTGGCTATGGCGGTTTGAt        <  2:1603538/100‑1 (MQ=255)
                                                                                     ttCCACCAGCACTCGAAGGCCGCGATCTCAATCCGATGCTGCAACATCCCGGTCTTATCTTTCATCCACCGCTGCTTTACCTTGGCTATGGCGGTTTGAt        <  2:3359940/100‑1 (MQ=255)
                                                                                      tCCGCCAGCACTCGAAGGCCGCGATCTCAATCCGATGCTGCAACATCCCGGTCTTATCTTTCATCCACCGCTGCTTTACCTTGGCTATGGCGGTTTGATg       <  2:1832965/100‑1 (MQ=255)
                                                                                      tCCACCAGCACTCGAAGGCCGCGATCTCAATCCGATGCTGCAACATCCCGGTCTTATCTTTCATCCACCGCTGCTTTACCTTGGCTATGGCGGTTTGATg       <  1:1355207/100‑1 (MQ=255)
                                                                                      tCCACCAGCACTCGAAGGCCGCGATCTCAATCCGATGCTGCAACATCCCGGTCTTATCTTTCATCCACCGCTGCTTTACCTTGGCTATGGCGGTTTGATg       <  2:3442088/100‑1 (MQ=255)
                                                                                      tCCACCAGCACTCGAAGGCCGCGATCTCAATCCGATGCTGCAACATCCCGGTCTTATCTTTCATCCACCGCTGCTTTACCTTGGCTATGGCGGTTTGATg       <  2:107433/100‑1 (MQ=255)
                                                                                      tCCACCAGCACTCGAAGGCCGCGATCTCAATCCGATGCTGCAACATCCCGGTCTTATCTTTCATCCACCGCTGCTTTACCTTGGCTATGGCGGTTTGATg       <  2:1664544/100‑1 (MQ=255)
                                                                                      tCCACCAGCACTCGAAGGCCGCGATCTCAATCCGATGCTGCAACATCCCGGTCTTATATTTCATCCACCGCTGCTTTACCTTGGCTATGGCGGTTTGATg       <  1:1730805/100‑1 (MQ=255)
                                                                                      tCCACCAGCAATCGAAGGCCGCGATCTCAATCCGATGCTGCAACATCCCGGTCTTATCTTTCATCCACCGCTGCTTTACCTTGGCTATGGCGGTTTGATg       <  1:399114/100‑1 (MQ=255)
                                                                                        caccaGCACTCGAAGGCCGCGATCTCAATCCGATGCTGCAACATCCCGGTCTTATCTTTCATCCACCGCTGCTTTACCTTGGCTATGGCGGTTTGATGGt     <  2:2734326/100‑1 (MQ=255)
                                                                                                |                                                                                              
CCGCTGACGTTAGCCGTTTTATCTCTCATGCTCGCCGCACTGCTACTGTTTGTGGTGCTGTGGTCCGATCCCTTCGTGCGGATATTTCCACCAGCAATCGAAGGCCGCGATCTCAATCCGATGCTGCAACATCCCGGTCTTATCTTTCATCCACCGCTGCTTTACCTTGGCTATGGCGGTTTGATGGTAGC  >  NC_000913/4291827‑4292017

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